Abstract
Gradient metasurfaces have been extensively applied in the unprecedented control of light beams over thin optical components. However, these metasurfaces suffer from low efficiency when bending light through large angles and high fabrication demand when it requires fine discretion. In this work, we investigate alldielectric metagratings based on the generalized Kerker effect induced by interference between Mietype resonances. It allows extraordinary optical diffraction for beam steering through ultralarge angles. The coupling inside and between the lattices in the metagrating can be used to tune the excited states of the electric and magnetic resonances, including both the fundamental dipoles and highorder multipoles, leading to an ideal asymmetrical scattering pattern that redistributes the energy between the diffraction channels as required. The quadrupole and hexadecapole not only significantly enhance the working efficiency but also enable distinctive possibilities for wave manipulation that cannot be reached by dipoles. Utilizing a thin array of silicon rods, largeangle negative refraction and reflection are realized with almost unity efficiency under both transverse magnetic and transverse electric polarization. Compared with conventional metasurfaces, such an alldielectric metagrating has the merits of high flexibility, high efficiency, and low fabrication requirements. The coupling and interactions among the multipoles may serve as a foundation for various forms of onchip optical wave control.
1 Introduction
The rapid development of exotic metamaterials has enabled the unprecedented manipulation of electromagnetic waves which cannot be attained by natural materials [1, 2]. Recently, dielectric nanoparticles with low dissipative losses and high permittivity, which support Mietype resonance, have been developed as favorable candidates for various nanophotonics applications [3], [4], [5]. Compared with their metallic counterparts, dielectric nanoparticles can excite a rich range of electric and magnetic multipoles without needing geometries such as splitring resonators or Ushaped particles [6]. The low losses in combination with the localized multipolar nature can support highquality resonances with the field increased by several orders of magnitude [7, 8]. In addition, the resonances are sustained inside the dielectric particles [9], unlike resonances for plasmonic structures, which mainly rely on surface modes. These properties produce to a strong light–matter interaction and are promising for nonlinear phenomena [10, 11], Fano resonances [12, 13], sensing [14], guiding [15], spectroscopy [16] and magnetic mirrors [17].
Most importantly, by tuning the shape or geometric parameters of dielectric particles, the interplay among the electric and magnetic multipoles can be used to realize unconventional scattering [3], [4], [5]. For instance, by overlapping the electric dipole (ED) and magnetic dipole (MD) with identical values in the spectra, nobackward or noforward scattering can be achieved [18]. Unidirectional scattering has significant potential in various optical applications [19, 20]. As one of the most promising concepts, Huygens metasurfaces based on the Kerker effect can freely modulate the electromagnetic characteristics with an ultrahigh efficiency [17, 21], [22], [23], [24], [25], [26], [27], [28], [29], [30]. By accumulating the gradient changes of periodic elements, the phase difference between ED and MD is able to achieve a full 2π phase shift [25]. The abrupt phase changes induced by Huygens metasurfaces enable intriguing wavefront control, leading to various applications, such as abnormal refraction [22], [23], [24], [25], [26], [27], [28], focusing [17, 22], [23], [24, 26], absorption [29], and artificial colors [30]. Compared with conventional optical devices, a Huygens metasurface has a high flexibility, a low profile, and a high transmission. However, Huygens metasurfaces also have some limitations. Achieving sufficient spatial resolution for the wavefront control requires fine spatial discretization and elaborate fabrication. The corresponding manufactured resolution has a considerable influence on the working performance of the optical devices. When routing beams through large angles, owing to the intrinsic drawbacks of gradient metasurfaces [31, 32], complex geometric configurations are needed to avoid the dramatic decrease of efficiency [33]. The anomalous refractive angles achieved with Huygens metasurfaces are often smaller than 30° [17, 24], [25], [26], [27], [28], which is insufficient for many applications.
Due to diffraction and resonance, gratings can transfer the incident energy into a discrete set of diffraction channels according to the Floquet theory. Under these circumstances, the light is bent through a large angle by the uniform elements, which are not so small as to require an extremely high fabrication resolution. A metasurface adopting the diffraction effect is also known as the metagrating which has advantages on achieving extreme wave manipulation [34], [35], [36]. The essential principle for the metagrating is completely suppressing all the unselected diffraction channels using asymmetrical scattering elements. It is still difficult to realize this suppression since each channel is not independently controllable [37]. Recently, it has been shown that the interference between Mie resonances in periodic structures with defined multipolar characteristics is a novel way to form scattering patterns, which is known as the generalized Kerker effect [38]. This designed scattering pattern can be used to choose the diffraction order and eliminate all the unwanted channels. Utilizing simple dielectric grating configurations, some pioneering researchers have achieved anomalous refraction with angles of up to 90° [39], [40], [41], [42], [43], [44, 47]. However, there remain several important issues plaguing the development of highindex dielectric metagratings based on the generalized Kerker effect.
First, the influence of the lattice coupling between the periodic nanoparticles is comparable to the effect of the resonance inside an individual element [45, 46]. Therefore, the scattering by periodic structures can be totally different from an isolated lattice (IL) [47]. Some studies also utilize the lattice coupling to overlap the ED and MD spectra, but relevant designs for abnormal wavefront control have not yet been fully investigated [45]. Many other investigations just neglect the influence of lattice coupling [39, 48] or do not discuss it in sufficient detail [40], [41], [42], [43]. Second, as well as the fundamental ED and MD, the highorder multipoles offer new possibilities in the generalized Kerker effect for shaping complex scattering patterns [49]. Although some works have recognized the importance of these modes, such as in enhancing the scattering directionality or building magnetic mirrors [47, 50, 51], there is still no method for applying or optimizing highorder multipoles in diffraction selection. In addition, periodic doublet configurations have been employed for improving working performances [41], [42], [43]. The interparticle coupling inside doublets is also important for the generalized Kerker effect, but there are few analytical works on this.
In this work, we analytically investigate the efficient manipulation of waves through large angles, as induced by the generalized Kerker effect inside periodic highindex nanoparticles. Based on the multiple scattering theory, both the simple grating (SG), formed by a single dielectric rod in one period, and the compound grating (CG), formed by a dimer in one period, are investigated in detail. The coupling inside and between the periodic dielectric lattices allows not only EDs and MDs but also highorder multipoles to manipulate waves. With appropriate designs, the scattering by the nanoparticles can be elaborately tailored by tuning the interference between the multipoles, allowing the light energy to be harnessed in both the transmission and reflection diffraction channels. Thus, the highorder multipoles not only help to improve the working performance of the dipoles but also play a dominant role in wave control. They have different features from dipoles, which can increase flexibility and working efficiency. Efficiencies close to 100% for the largeangle negative reflection and refraction of incident waves are specifically demonstrated for both transverse magnetic (TM) and transverse electric (TE) polarization. Thus, these lowloss alldielectric structures may stimulate the design of compact optical components with a low profile and high efficiency.
2 Theory for the interaction of light with dielectric cylinders
Let us start by considering a onedimensional (1D) array of rods with a refraction index of 3.5 (approaching that of silicon). As shown in Figure 1a and b, the lattices in the array have one or two dielectric cylinders, denoted as SG or CG, respectively. The period of the lattice is L and cylinders are uniform in the zdirection. For the SG, all the rods have the same radius r. A CG lattice is composed of two cylinders with radiuses r _{1} and r _{2} and a centertocenter distance d. In the local coordinate system for a specific lattice, the line connecting the two rods has an angle ϕ with respect to the xaxis.
Figure 1:
To simplify the derivation, we first consider twodimensional (2D) scattering cases when a plane wave with incident angle θ impinges on the grating. Realizable threedimensional (3D) cases are calculated and discussed in the last section. The polarization of the incident waves can be TM or TE when the electric or magnetic field is polarized along the zaxis. Under a polar coordinate (ρ, φ), the scattering fields of an isolated cylinder are expanded into a series of cylindrical harmonics:
where P = [i ^{ m } H _{ m } ^{(1)}(k _{ 0 } ρ)e ^{ imφ }] and σ = [σ_{ m }] are column vectors. F is the E _{ z } or H _{ z } component for TE or TM polarization. H _{ m } ^{(1)} is the mthorder Hankel function of the first kind and σ _{ m } is the scattering coefficient of the harmonic of the mth cylinder normalized by 2λ/π, where λ is the working wavelength [52]:
in which J _{ m } is the mthorder Bessel function, J′ _{ m } and H′ _{ m } ^{(1)} are the derivatives of the original functions, N is the refractive index of the rod. Then, the total scattering crosssection σ _{s} of a single rod is the sum of all scattering orders:
For TMpolarized waves, m = 0, 1, 2, and 3 correspond to the MD, ED, electric quadrupole (EQ), and electric hexadecapole (EH). For TE polarization, m = 0, 1, 2, and 3 are the ED, MD, magnetic quadrupole (MQ), and magnetic hexadecapole (MH). The scattering spectra of an isolated silicon rod under the two types of polarization are plotted in Figure S1a and b (Section 1 of the supplementary material). The field distributions and scattering properties of these multipoles are also calculated and illustrated in the supplementary material. It is clear that there are electric and magnetic multipoles under both types of polarization, which enriches the scattering phenomenon of the dielectric rod. The interplay among these multipoles can lead to the generalized Kerker effect. Therefore, it is possible to tailor the excited states of the multipoles through proper geometric design to realize unique scattering properties.
Note that, owing to the rotational symmetry, opposite orders always have the same scattering coefficient and phase difference, i.e., σ _{ m } = σ _{−m }. When the cylinders form a 1D SG, the detailed scattering properties can be analytically investigated by the multiple scattering theory. Thus, the scattering field of an SG can be expressed as [52]:
Here,
where
For the CG, each periodic lattice consists of two dielectric rods. The fields scattered by one rod become the incident waves for the other. Thus, an extra interlattice coupling must be considered. To see the influence of this interlattice coupling, we first analyze the scattering of an isolated doublet. By considering the multiple couplings and the boundary conditions, the coefficients for the multipoles inside each rod (
where T is a square matrix describing the interlattice effects, including electric–electric coupling, magnetic–magnetic coupling, and electric–magnetic coupling. For example,
where ϕ ^{
j
} = ϕ or –ϕ for j = 1 or 2, representing the angle of the rod center with respect to the xaxis in its local coordinates. After setting the values of the T matrices (see Section 2 of the supplementary material),
When the doublet is periodically arranged as a CG, both the interlattice coupling and lattice coupling affect the generalized Kerker effect. To characterize the scattering properties quantitatively, each unit cell of the CG is considered to be a single composite structure, such that its scattering features are also described by the multipoles with coefficient
where
The total transmission and reflection coefficients of the mthorder diffraction harmonics T _{ m } and R _{ m } in Cartesian coordinates are
where
The term
3 Wave steering under TM polarization
As mentioned above, to realize abnormal wave control with extra parameters, the dielectric grating must work in the metagrating regime
In the metagrating regime, at least four possible diffraction channels are available, i.e., T
_{0} with θ
_{T0} = θ, T
_{−1} with θ
_{
T−1} = arc cos(−λ/L + cos(θ)), R
_{0} with θ
_{
R0} = −θ, and R
_{−1} with θ
_{
R−1} = −arc cos(−λ/L + cos(θ)) − π. The negative refraction under TM polarization when preserving the T
_{−1} channel is considered first. We start by considering an SG with L = 4.3r illuminated by a TMpolarized wave at λ =
Figure 2:
As discussed, the cancellation of T _{0} is caused by destructive interference between the incident waves and the waves scattered by the cylinders. The angular scattering pattern of each rod is compared with two approximations, the IR approximation and the dipole approximation (DA), for ED + MD, to illustrate the contributions of highorder multipoles and lattice coupling. The importance of lattice coupling for the generalized Kerker effect is obvious as the unique scattering behavior disappears under the IR approximation. However, the DA can reproduce the scattering pattern well, indicating that highorder multipoles do not participate in this steering process. To enhance the working performance, it is necessary to optimize the angular scattering shape by tuning the states of the dipoles or taking highorder multipoles into consideration, but the limited tuning freedom of SG constrains the further improvement.
The extra coupling inside the lattices in a CG provides more possibilities than an SG for the generalized Kerker effect. The first CG (CG1) is also taken from reference [41]. Its parameters are r _{1} = r, r _{2} = 0.4r, d = r _{1} + r _{2} + 0.05r, L = 4.3r, and ϕ = 90°. The results, shown in Figure 2(b1)–(b4), indicate a working efficiency of 92%. Unlike the analysis in reference [41], this doublet dimer mainly contributes to enhancing the coefficient of MD rather than the ED. Compared with the SG, the R _{−1} channel is further suppressed due to the enhanced interference between the dipoles. This improvement is induced by the interaction inside the doublet, as shown in Figure 2(b3). However, there are still no highorder multipoles. As above, the farfield scattering pattern of the DA is slightly different from that for CG1. In addition, the scattering pattern of an IL is like that for an IR. Therefore, the negative refraction of CG1 is also mainly caused by lattice coupling.
To take the highorder multipoles into the generalized Kerker effect, one can increase the effective size of the composite unit, which is achieved by increasing the radius of each cylinder or enlarging the separation between them (see Section 2 of the supplementary material). After optimization, a new CG (CG2) is chosen with parameters r _{1} = 1.06r, r _{2} = 0.64r, ϕ = 90°, and d = r _{1} + r _{2} + 0.2r. As illustrated in Figure 2(c1), more than 99% of the incident energy is transferred into the T _{−1} channel. Note that this highefficiency negative refraction with a bending angle of 90° is hard to realize with common gradient metasurfaces, verifying the advantages of a CG for largeangle steering. The nearfield distribution in Figure 2(c3) shows that strong ED and MD are excited in the bigger and smaller rods, unlike Figure 2(b3), where most of the energy is confined in the bigger rods. More importantly, as shown in Figure 2(c2), CG2 takes EQs into the interference process. According to the generalized Kerker effect, the participation of EQs produces a highly directional scattering pattern, which cannot be approximated by the DA, as indicated in Figure 3(c3). Moreover, the R _{0} and R _{−1} channels are eliminated well owing to the zero backward scattering. The totally symmetrical forward scattering cancels the T _{0} channel and generates negative refraction with nearunity efficiency. Also, notice that the scattering pattern for the IL is very similar to that for a CG, but different from IR1 and IR2 (Figure S3b in the supplementary material). Therefore, the desired scattering pattern is mainly caused by the interlattice coupling, unlike CG1. The lattice coupling helps to optimize the scattering shape and further improves the working efficiency. Thus, we have demonstrated that the extra coupling inside the lattice increases the potential of the metagrating and excites highorder multipoles for the generalized Kerker effect.
Figure 3:
To further show the promising ability of highorder multipoles, we present an extra example (CG3) with r _{1} = 1.3r, r _{2} = 0.71r, ϕ = 90°, and d = r _{1} + r _{2} + 0.08r. Its working efficiency also exceeds 99%. The corresponding results are illustrated in Figure 2(d1)–(d4). In this case, the intriguing negative refraction is caused by the combination of the first four scattering modes. In contrast to the former designs, the highorder multipolar modes play a dominant role, as shown by the excited coefficients and farfield angular scattering patterns. From the nearfield map of the lattice, we also see features of highorder multipoles. There are several zero points in the angular scattering pattern. Although CG3 has nozero backward scattering, the directions of R _{0} and R _{−1} coincide with the zero points. This is clear evidence that the highorder multipoles provide extra flexibility for the generalized Kerker effect, leading to the enhancement of working efficiency. Figure 2(d4) clearly indicates that the interlattice coupling is very important for forming the required scattering pattern. The lattice coupling is also necessary for the total suppression of unwanted diffraction channels. Therefore, by utilizing the coupling effect inside and among lattices, the highorder multipoles which are usually hard to control, can be excited efficiently and applied in wavefront engineering with impressive performances.
Next, we investigate the dependence of the inplane negative refraction on the wavelength and the period of the grating, the results of which are plotted in Figure 3a and b, respectively. The theoretical calculations are accompanied by the simulations performed by COMSOL Multiphysics. Periodic ports and boundaries are utilized in COMSOL for mimicking the incident plane wave and infinite metagratings. The theoretical results match the simulations well, verifying the validity of the theory. In Figure 3a, near the investigated wavelengths, Fanoshaped profiles can be recognized in the spectra of CG2 and CG3. The peaks and dips are associated with constructive and destructive interference between the dipoles and the highorder multipoles. This further indicates that the unity working efficiency is attributed to the participation of highorder multipoles. SG and CG1 have smooth wavelength responses because only dipoles are excited in both structures. In addition, the high efficiency of CG2 and CG3 is sustained over a large range of periods because the generalized Kerker effect for the two structures is mainly caused by the interlattice coupling. The variation of the period only has influences on the lattice coupling. As a comparison, SG and CG1 which are mainly based on the lattice coupling, have low tolerance with the change of period.
According to the above analysis for a metagrating, each value of the period is associated with a bending angle. Therefore, we can carefully modify the period to change the wavefront, thus achieving a complex wave control. As a specific realization, we designed a focusing lens. Its working schematic is illustrated in Figure 3c. According to the phase profile for focusing light [54], the period of the metagrating should be set as
where
By setting a focal distance f = 12λ, we can get the distribution of period L(x). After considering the value of local L(x) and the global space of this metagrating, the position of each metaatom is obtained. The parameters of this focusing lens are put in Section 5 of the supplementary material, including the geometric size, phase response, and transmission coefficient of the metaatoms. The numerical simulation of a Gaussian beam impinging on the grating is presented in Figure 3d. The light is bent and focused to a specific position with an efficiency of over 85%. The full width at half maximum (FWHM) of the focusing point is about 0.98λ, as shown in Figure 3d. The focusing can be further optimized as Eqs. (12) and (13) do not consider the influence of the gradient period. The achieved result demonstrates the potential of these periodic nanoparticles. Of course, more applications based on this wavefrontshaping method are also possible.
The inplane negative reflection means that the incident and reflected waves are on the same side of the normal. According to the above analysis, the waves scattered from the multipolar modes must undergo destructive interference with the incident waves in the 0th transmission channel and constructive interference in the −1st reflection channel. Moreover, the scattering due to specular reflection and −1st transmission should be suppressed. According to the generalized Kerker effect, the required angular scattering pattern cannot be achieved just by using the interaction between the MD and the ED.
An optimized SG for negative reflection is set as L = 3.1r and λ =
Figure 4:
The negative reflection of the chosen SG is very sensitive to the wavelength. The corresponding working bandwidth is of about 5%, as plotted in Figure 5a. These attributes to the properties of highorder multipoles are inherited in the negative refraction as the contribution of the ED is dominant. Due to the importance of lattice coupling and the high sensitivity of the ED, the working state strongly depends on the period of this grating, as shown in Figure 5b. These factors may create great challenges for practical fabrication. However, the high sensitivity of the working wavelength and geometry also offers significant potential for many applications [55], such as biomedical sensing, wavelength filtering, and hyperspectral imaging. As a specific example for use in a narrowband filter or multiplexer, the field distribution of this SG is put in Figure 5c whose calculated frequency is set as λ = 4r. More than 90% of the incident energy directly passes through the structure with little diffraction.
Figure 5:
The performances of CG4 are similar to structures whose dipoles are dominant, as the highorder multipoles help to form the scattering pattern. Therefore, CG4 is tolerant to changes in the wavelength. The interlattice coupling also increases the stability with changes to the grating period. The slight deviation between theory and simulation when L < 2.9r is caused by that the size of an IL (r _{1} + r _{2} + d) is much larger than the period of grating L. In the deduction of the T matrix for CGs, the lattice is assumed as a composite structure. This assumption is inaccurate when the period is much smaller than the lattice. However, the trend of the theoretical curve is still well matched with the simulation in Figure 5b. By carefully tailoring the geometry, CG4 can be applied for complex wave steering. Figure 5d depicts an example by setting L = 3.5r. The angle of reflection is −53° for an efficiency of 85%. Therefore, the interlattice coupling is able to provide further potential for realizing the generalized Kerker effect and achieving highefficient wave steering. In addition, the highorder multipoles can not only improve the working performance formed by dipoles but also dominantly achieve unique scattering patterns which are distinctive from dipoles.
4 Wave steering under TE polarization
In the following, we will show that the highefficiency negative refraction and reflection realized for TMpolarized waves are also achievable for TE polarization. As above, we first consider a CG with L = 4.3r illuminated by 45° TEpolarized waves with λ =
Figure 6:
Inspired by this result, we further explore the potential of the higher scattering mode MH using an SG with r = 1.4r _{1} (r _{1} is the same as for CG5) and L = 3.06r. The corresponding results are depicted in Figure 6(c1)–(c4), respectively. Negative refraction with efficiency of over 89% is obtained and the contribution of hexadecapole appears in the nearfield snapshot. Because of the multipole scattering and coupling, there are also significant excitation of ED, MD, and MQ, resulting in negative refraction through interference. The directions of the R _{0} and R _{−1} channels coincide with the zero points in the angular scattering pattern. As above, the generalized Kerker condition cannot be fulfilled without the MH, fully verifying the necessity of the highorder multipoles. These designs clearly indicate that intriguing wave steering can be induced by the higherorder scattering multipoles. The working efficiency is plotted against the working wavelength in Figure 6d. Here, we emphasize that, although the quadrupole participates in the negative refraction for both TM and TE polarization, the situations are quite different. For the results in Figure 2c, the quadrupole is applied to enhance the working performance induced by the dipoles. For TE polarization, the quadrupole dominates the wave control. A similar conclusion is also valid for the MH in Figure 6c. Therefore, both structures have a high sensitivity to the working parameters. The relative working bandwidths are 7 and 6% for the wire gratings in Figure 6b and c.
Negative reflection is produced under TE polarization in the same way as for TM polarization, i.e., the highorder multipoles must be excited to realize the generalized Kerker effect. Based on this, two optimized CGs with r
_{1} = 0.22r, r
_{2} = r, ϕ = 90°, d = r
_{1} + r
_{2} + 0.05r, and L = 4.3r (CG6) and r
_{1} = 0.42r, r
_{2} = 0.86r, ϕ = 0°, d = r
_{1} + r
_{2} + 0.44r, and L = 4.3r (CG7) are chosen. Their efficiencies are respectively 92 and 94%. Interestingly, the parameters of CG6 are the same as those in Figure 6b, when the incident angle is set as 135°. The calculated results for this CG are shown in Figure 7(a1)–(a4) with θ = 45° and λ =
Figure 7:
5 Discussion and conclusion
Note that, although the promising wave steering is realized by using the metagrating profiles, the ability of dielectric rods is obviously not limited to the several functions demonstrated. Depending on the requirements, a rod grating can work in transmission or reflection mode and the working performance can be sensitive or robust to the parameters. The multipolar interference and collective lattice coupling provide flexible and complex wavefront control for TM and TE polarization. The working efficiency and fabrication demand of the dielectric rods can exceed the conventional gradient metasurfaces especially bending light with large angles. Note that, by using an asymmetrical CG, the abnormal control of waves with normal incidence is also realizable. We also present CG8, which has r _{1} = 0.9r, r _{2} = 0.58r, ϕ = 67°, d = r _{1} + r _{2} + 0.05r, and L = 7r and CG9, which has r _{1} = 0.81r, r _{2} = 0.54r, ϕ = 67°, d = r _{1} + r _{2} + r, and L = 7r. These can achieve abnormal refraction of 50° for TM and TE polarization with θ = 90°. The calculated wavelength is L sin(50°), leading to six diffraction channels: R _{0}, R _{−1}, R _{1}, T _{0}, T _{1}, and T _{−1}. The efficiencies achieved for the two structures are 92 and 90%, respectively, as shown in Figure 8a and b. For this design, the scattering patterns must have high directionality to suppress the four unwanted diffraction channels and cancel out the 0th transmission, which is more complex than previous investigations. From the farfield angular scattering patterns plotted in the insets, we can conclude that the backward scattering is suppressed well and the forward scattering is asymmetrical and highly directional, consistent with our predictions. The multiple couplings inside and among the lattice help to efficiently produce the first three multipoles, which lead to the required scattering properties.
Figure 8:
Considering practical fabrication, the metagratings must be truncated in the zdirection. Therefore, we performed serval fullwave 3D simulations to study the working performances of the metagratings. The height of the rods is set as 4 μm, located on silica (refractive index is 1.5) substrate with 300 nm thickness. The metagratings under these parameters are totally achievable through microfabrication, such as the Ebeam lithography and the reactive ion etch [41, 43]. The configuration of the simulations is presented in Figure 8c, and more details are put in Section 6 of the supplement. Here, an individual lattice is chosen for calculation and periodic conditions are applied in the xdirection to mimic infinite gratings. Simulated results for CG2, CG4, CG5, and CG7 are illustrated in Figure 8d–g, respectively, corresponding to the negative refraction and reflection under both types of polarization. Further, other mentioned structures in Figures 2 –8 are also calculated and put in the supplementary material. Clearly, the field maps in Figure 8 are fully consistent with the corresponding 2D results. Therefore, although the working properties of multipoles are influenced by the height of the rods [4], the highefficiency wave steering can be achieved under realizable situations. The grating parameters can be further optimized to compensate for the finite thickness. In addition, the influence of loss in the infrared spectra is also investigated in detail (Section 7 of supplementary material). Although the working efficiency is slightly reduced, the unique scattering patterns due to interference among the multipoles are conserved when the loss is taken into account. In addition, compared with dipoles, the highorder multipoles are more sensitive to losses or the height of rods. For these designs (CG5 and CG6) which are mainly based on the highorder multipoles, the decrease of working efficiency is inevitable. There are still many challenges in efficiently applying the highorder multipoles.
In conclusion, in this work, we show that highefficiency wave steering is achievable by the generalized Kerker effect which is based on interactions and coupling between multipole resonant modes of alldielectric nanoparticles. Based on multipole scattering theory, both an SG formed by one rod per period and a CG formed by two rods per period are analyzed. The coupling inside and among the lattices can efficiently tailor the multipole states, leading to extensive scattering properties. By using the simple lattice configuration at the metagrating regime, highefficiency negative refraction and reflection through large bending angles are demonstrated for both types of polarization. These effects are realized by controlling the diffraction channels through multipole interference. In addition to the fundamental dipoles, highorder multipoles with distinctive properties are studied and extensively applied for wave manipulation. These provide extra possibilities for nanoparticlebased beam control. This work may provide support for onchip optical components, with functions such as bending, sensing, filtering, and focusing.
Supplement
The supplementary material is available for some supporting information. Section 1 investigates the scattering properties of cylindrical multipoles; In Section 2, we present the scattering properties of an isolated doublet; In Section 3, the T matrix for a CG is derived and the coupling inside the lattice is analyzed; Section 4 shows the scattering properties IRs in CG; Section 5 gives the parameters of the focusing lens; Section 6 describes the 3D simulations and some simulated results for the rods with a finite thickness; Section 7 investigates the influence of loss on the working performance.
Funding source: National Key Research and Development Program of China10.13039/501100012166
Award Identifier / Grant number: 2019YFA0210203
Funding source: National Natural Science Foundation of China10.13039/501100001809
Award Identifier / Grant number: 61971013

Author contribution: All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.

Research funding: The study was supported by the National Key Research and Development Program under Grant No. 2019YFA0210203 and the National Natural Science Foundation of China under Grant No. 61971013.

Conflict of interest: The authors declare no conflict of interest.
References
[1] T. J. Cui, D. R. Smith, and R. Liu, Metamaterials, Spring Street, NY, Springer, 2010. Search in Google Scholar
[2] J. M. Hao, M. Qiu, and L. Zhou, “Manipulate light polarizations with metamaterials: from microwave to visible,” Front. Phys. China, vol. 5, pp. 291–307, 2010. https://doi.org/10.1007/s114670100005y. Search in Google Scholar
[3] N. Yu, P. Genevet, M. A. Kats, et al., “Optically resonant dielectric nanostructures,” Science, vol. 354, p. 6314, 2016. Search in Google Scholar
[4] Y. H. Ko and R. Magnusson, “Wideband dielectric metamaterial reflectors: Mie scattering or leaky Bloch mode resonance?,” Optica, vol. 5, pp. 289–294, 2018. https://doi.org/10.1364/optica.5.000289. Search in Google Scholar
[5] I. Staude and J. Schilling, “Metamaterialinspired silicon nanophotonics,” Nat. Photonics, vol. 11, pp. 274–284, 2017. https://doi.org/10.1038/nphoton.2017.39. Search in Google Scholar
[6] L. Zhou, X. Huang, Y. Zhang, et al.., “Resonance properties of metallic ring systems,” Mater. Today, vol. 12, pp. 52–59, 2009. https://doi.org/10.1016/s13697021(09)703177. Search in Google Scholar
[7] E. Calandrini, A. Cerea, F. De Angelis, et al.., “Magnetic hotspot generation at optical frequencies: from plasmonic metamolecules to alldielectric nanoclusters,” Nanophotonics, vol. 8, pp. 45–62, 2018. https://doi.org/10.1515/nanoph20180138. Search in Google Scholar
[8] Y. Yang, V. A. Zenin, and S. I. Bozhevolnyi, “Anapoleassisted strong field enhancement in individual alldielectric nanostructures,” ACS Photonics, vol. 5, pp. 1960–1966, 2018. https://doi.org/10.1021/acsphotonics.7b01440. Search in Google Scholar
[9] K. Koshelev, G. Favraud, A. Bogdanov, et al.., “Nonradiating photonics with resonant dielectric nanostructures,” Nanophotonics, vol. 8, pp. 725–745, 2019. https://doi.org/10.1515/nanoph20190024. Search in Google Scholar
[10] L. Carletti, A. Locatelli, D. Neshev, et al.., “Shaping the radiation pattern of secondharmonic generation from AlGaAs dielectric nanoantennas,” ACS Photonics, vol. 3, pp. 1500–1507, 2016. https://doi.org/10.1021/acsphotonics.6b00050. Search in Google Scholar
[11] D. A. Smirnova, A. B. Khanikaev, L. A. Smirnov, et al.., “Multipolar thirdharmonic generation driven by optically induced magnetic resonances,” ACS Photonics, vol. 3, pp. 1468–1476, 2016. https://doi.org/10.1021/acsphotonics.6b00036. Search in Google Scholar
[12] B. Hopkins, D. S. Filonov, A. E. Miroshnichenko, et al.., “Interplay of magnetic responses in alldielectric oligomers to realize magnetic Fano resonances,” ACS Photonics, vol. 2, pp. 724–729, 2015. https://doi.org/10.1021/acsphotonics.5b00082. Search in Google Scholar
[13] J. Yan, P. Liu, Z. Lin, et al.., “Directional F‵ano resonance in a silicon nanosphere dimer,” ACS Nano, vol. 9, pp. 2968–2980, 2015. https://doi.org/10.1021/nn507148z. Search in Google Scholar
[14] O. Yavas, M. Svedendahl, and R. Quidant, “Unravelling the role of electric and magnetic dipoles in biosensing with Si nanoresonators,” ACS Nano, vol. 13, pp. 4582–4588, 2019. https://doi.org/10.1021/acsnano.9b00572. Search in Google Scholar
[15] R. M. Bakker, Y. F. Yu, R. PaniaguaDomínguez, et al.., “Resonant light guiding along a chain of silicon nanoparticles,” Nano Lett., vol. 17, pp. 3458–3464, 2017. https://doi.org/10.1021/acs.nanolett.7b00381. Search in Google Scholar
[16] A. Krasnok, M. Caldarola, N. Bonod, et al.., “Spectroscopy and biosensing with optically resonant dielectric nanostructures,” Adv. Opt. Mater., vol. 6, 2018, Art no. 1701094. https://doi.org/10.1002/adom.201701094. Search in Google Scholar
[17] Z. Ma, S. M. Hanham, P. Albella, et al.., “Terahertz alldielectric magnetic mirror metasurfaces,” ACS Photonics, vol. 3, pp. 1010–1018, 2016. https://doi.org/10.1021/acsphotonics.6b00096. Search in Google Scholar
[18] Y. H. Fu, A. I. Kuznetsov, A. E. Miroshnichenko, et al.., “Directional visible light scattering by silicon nanoparticles,” Nat. Commun., vol. 4, pp. 1–6, 2013. https://doi.org/10.1038/ncomms2538. Search in Google Scholar
[19] F. Todisco, R. Malureanu, C. Wolff, et al.., “Magnetic and electric Mieexciton polaritons in silicon nanodisks,” Nanophotonics, vol. 9, pp. 803–814, 2020. https://doi.org/10.1515/nanoph20190444. Search in Google Scholar
[20] A. Furasova, E. Calabró, E. Lamanna, et al.., “Resonant silicon nanoparticles for enhanced light harvesting in halide perovskite solar cells,” Adv. Opt. Mater., vol. 6, 2018, Art no. 1800576. https://doi.org/10.1002/adom.201800576. Search in Google Scholar
[21] S. Xiao, J. Wang, F. Liu, et al.., “Spindependent optics with metasurfaces,” Nanophotonics, vol. 6, pp. 215–234, 2017. https://doi.org/10.1515/nanoph20160121. Search in Google Scholar
[22] O. Tsilipakos, A. C. Tasolamprou, T. Koschny, et al.., “Pairing toroidal and magnetic dipole resonances in elliptic dielectric rod metasurfaces for reconfigurable wavefront manipulation in reflection,” Adv. Opt. Mater., vol. 6, 2018, Art no. 1800633. https://doi.org/10.1002/adom.201800633. Search in Google Scholar
[23] A. Forouzmand and H. Mosallaei, “Alldielectric Cshaped nanoantennas for light manipulation: tailoring both magnetic and electric resonances to the desire,” Adv. Opt. Mater., vol. 5, 2017, Art no. 1700147. https://doi.org/10.1002/adom.201700147. Search in Google Scholar
[24] J. Engelberg, C. Zhou, N. Mazurski, et al.., “NearIR widefieldofview Huygens metalens for outdoor imaging applications,” Nanophotonics, vol. 9, pp. 361–370, 2020. https://doi.org/10.1515/nanoph20190177. Search in Google Scholar
[25] C. Liu, L. Chen, T. Wu, et al.., “Alldielectric threeelement transmissive Huygens’ metasurface performing anomalous refraction,” Photon. Res., vol. 7, pp. 1501–1510, 2019. https://doi.org/10.1364/prj.7.001501. Search in Google Scholar
[26] L. Zhang, J. Ding, H. Zheng, et al.., “Ultrathin highefficiency midinfrared transmissive Huygens metaoptics,” Nat. Commun., vol. 9, p. 1481, 2018. https://doi.org/10.1038/s41467018038317. Search in Google Scholar
[27] M. I. Shalaev, J. Sun, A. Tsukernik, et al.., “Highefficiency alldielectric metasurfaces for ultracompact beam manipulation in transmission mode,” Nano Lett., vol. 15, pp. 6261–6266, 2015. https://doi.org/10.1021/acs.nanolett.5b02926. Search in Google Scholar
[28] S. M. Kamali, E. Arbabi, A. Arbabi, et al.., “A review of dielectric optical metasurfaces for wavefront control,” Nanophotonics, vol. 7, pp. 1041–1068, 2018. https://doi.org/10.1515/nanoph20170129. Search in Google Scholar
[29] J. Tian, Q. Li, P. A. Belov, et al.., “HighQ alldielectric metasurface: super and suppressed optical absorption,” ACS Photonics, vol. 7, pp. 1436–1443, 2020. https://doi.org/10.1021/acsphotonics.0c00003. Search in Google Scholar
[30] S. Q. Li, W. Song, M. Ye, et al.., “Generalized method of images and reflective color generation from ultrathin multipole resonators,” ACS Photonics, vol. 5, pp. 2374–2383, 2018. https://doi.org/10.1021/acsphotonics.8b00161. Search in Google Scholar
[31] N. M. Estakhri and A. Alu, “Wavefront transformation with gradient metasurfaces,” Phys. Rev. X, vol. 6, 2016, Art no. 041008. https://doi.org/10.1103/physrevx.6.041008. Search in Google Scholar
[32] A. DíazRubio, V. S. Asadchy, A. Elsakka, et al.., “From the generalized reflection law to the realization of perfect anomalous reflectors,” Sci. Adv., vol. 3, 2017, Art no. e1602714. https://doi.org/10.1126/sciadv.1602714. Search in Google Scholar
[33] A. Epstein and G. V. Eleftheriades, “Synthesis of passive lossless metasurfaces using auxiliary fields for reflectionless beam splitting and perfect reflection,” Phys. Rev. Lett., vol. 117, 2016, Art no. 256103. https://doi.org/10.1103/physrevlett.117.256103. Search in Google Scholar
[34] Y. Ra’di, D. L. Sounas, and A. Alù, “Metagratings: beyond the limits of graded metasurfaces for wave front control,” Phys. Rev. Lett., vol. 119, 2017, Art no. 067404. Search in Google Scholar
[35] E. Khaidarov, H. Hao, R. PaniaguaDomínguez, et al.., “Asymmetric nanoantennas for ultrahigh angle broadband visible light bending,” Nano Lett., vol. 17, pp. 6267–6272, 2017. https://doi.org/10.1021/acs.nanolett.7b02952. Search in Google Scholar
[36] L. Z. Yin, T. J. Huang, F. Y. Han, et al.., “Terahertz multichannel metasurfaces with sparse unit cells,” Opt. Lett., vol. 44, pp. 1556–1559, 2019. https://doi.org/10.1364/ol.44.001556. Search in Google Scholar
[37] X. Song, L. Huang, C. Tang, et al.., “Selective diffraction with complex amplitude modulation by dielectric metasurfaces,” Adv. Opt. Mater., vol. 6, 2018, Art no. 1701181. https://doi.org/10.1002/adom.201701181. Search in Google Scholar
[38] W. Liu and Y. S. Kivshar, “Generalized Kerker effects in nanophotonics and metaoptics,” Opt. Express, vol. 26, pp. 13085–13105, 2018. https://doi.org/10.1364/oe.26.013085. Search in Google Scholar
[39] R. PaniaguaDominguez, Y. F. Yu, E. Khaidarov, et al.., “A metalens with a nearunity numerical aperture,” Nano Lett., vol. 18, pp. 2124–2132, 2018. https://doi.org/10.1021/acs.nanolett.8b00368. Search in Google Scholar
[40] D. Lin, M. Melli, E. Poliakov, et al.., “Optical metasurfaces for high angle steering at visible wavelengths,” Sci. Rep., vol. 7, pp. 1–8, 2017. https://doi.org/10.1038/s41598017021674. Search in Google Scholar
[41] J. Du, Z. Lin, S. T. Chui, et al.., “Optical beam steering based on the symmetry of resonant modes of nanoparticles,” Phys. Rev. Lett., vol. 106, 2011, Art no. 203903. https://doi.org/10.1103/physrevlett.106.203903. Search in Google Scholar
[42] H. Li, C. Qiu, J. Du, et al., “Broadband inplane light bending with a doublet silicon nanopost array,” IEEE Photon. J., vol. 8, pp. 1–5, 2015. Search in Google Scholar
[43] A. Wu, H. Li, J. Du, et al.., “Experimental demonstration of inplane negativeangle refraction with an array of silicon nanoposts,” Nano Lett., vol. 15, pp. 2055–2060, 2015. https://doi.org/10.1021/nl5049516. Search in Google Scholar
[44] H. Chalabi, Y. Ra’Di, D. L. Sounas, et al.., “Efficient anomalous reflection through nearfield interactions in metasurfaces,” Phys. Rev. B, vol. 96, 2017, Art no. 075432. https://doi.org/10.1103/physrevb.96.075432. Search in Google Scholar
[45] V. E. Babicheva and A. B. Evlyukhin, “Resonant lattice Kerker effect in metasurfaces with electric and magnetic optical responses,” Laser Photonics Rev, vol. 11, 2017, Art no. 1700132. https://doi.org/10.1002/lpor.201700132. Search in Google Scholar
[46] V. E. Babicheva and J. V. Moloney, “Lattice effect influence on the electric and magnetic dipole resonance overlap in a disk array,” Nanophotonics, vol. 7, pp. 1663–1668, 2018. https://doi.org/10.1515/nanoph20180107. Search in Google Scholar
[47] W. Liu and A. E. Miroshnichenko, “Beam steering with dielectric metalattices,” ACS Photonics, vol. 5, pp. 1733–1741, 2017. https://doi.org/10.1021/acsphotonics.7b01217. Search in Google Scholar
[48] H. Song, L. Sun, and G. P. Wang, “Tunable perfect magnetic mirrors and retroreflectors in terahertz band,” Opt Express, vol. 28, pp. 753–759, 2020. https://doi.org/10.1364/oe.382965. Search in Google Scholar
[49] I. M. Hancu, A. G. Curto, M. CastroLópez, et al.., “Multipolar interference for directed light emission,” Nano Lett., vol. 14, pp. 166–171, 2014. https://doi.org/10.1021/nl403681g. Search in Google Scholar
[50] W. Liu, “Generalized magnetic mirrors,” Phys. Rev. Lett., vol. 119, 2017, Art no. 123902. https://doi.org/10.1103/physrevlett.119.123902. Search in Google Scholar
[51] P. D. Terekhov, A. B. Evlyukhin, D. Redka, et al.., “Magnetic octupole response of dielectric quadrumers,” Laser Photonics Rev, vol. 14, 2020, Art no. 1900331. https://doi.org/10.1002/lpor.201900331. Search in Google Scholar
[52] C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles, Germany, John Wiley & Sons, 2008. Search in Google Scholar
[53] K. Yasumoto and K. Yoshitomi, “Efficient calculation of lattice sums for freespace periodic Green’s function,” IEEE Trans. Antenn. Propag., vol. 47, pp. 1050–1055, 1999. https://doi.org/10.1109/8.777130. Search in Google Scholar
[54] X. Yin, H. Zhu, H. Guo, et al.., “Hyperbolic metamaterial devices for wavefront manipulation,” Laser Photonics Rev, vol. 13, 2019, Art no. 1800081. https://doi.org/10.1002/lpor.201800081. Search in Google Scholar
[55] C. Yan, K. Y. Yang, and O. J. Martin, “Fanoresonanceassisted metasurface for color routing,” Light Sci. Appl., vol. 6, p. e17017, 2017. https://doi.org/10.1038/lsa.2017.17. Search in Google Scholar
Supplementary Material
The online version of this article offers supplementary material (https://doi.org/10.1515/nanoph20210158).
© 2021 TieJun Huang et al., published by De Gruyter, Berlin/Boston
This work is licensed under the Creative Commons Attribution 4.0 International License.