• Photonics Research
  • Vol. 7, Issue 12, 1501 (2019)
Chang Liu1、2、†, Lei Chen1、3、†, Tiesheng Wu2, Yumin Liu1、*, Jing Li1, Yu Wang1, Zhongyuan Yu1, Han Ye1、4, and Li Yu1
Author Affiliations
  • 1State Key Laboratory of Information Photonics and Optical Communications, Beijing University of Posts and Telecommunications, Beijing 100876, China
  • 2Guangxi Key Laboratory of Wireless Wideband Communication and Signal Processing, Guilin 541004, China
  • 3College of Optical Sciences, The University of Arizona, Tucson, Arizona 85721, USA
  • 4e-mail: dabombyh@aliyun.com
  • show less
    DOI: 10.1364/PRJ.7.001501 Cite this Article Set citation alerts
    Chang Liu, Lei Chen, Tiesheng Wu, Yumin Liu, Jing Li, Yu Wang, Zhongyuan Yu, Han Ye, Li Yu, "All-dielectric three-element transmissive Huygens’ metasurface performing anomalous refraction," Photonics Res. 7, 1501 (2019) Copy Citation Text show less

    Abstract

    Metasurfaces have pioneered a new avenue for advanced wave-front engineering. Among the various types of metasurfaces, Huygens’ metasurfaces are thought to be a novel paradigm for flat optical devices. Enabled by spectrally overlapped electric resonance and magnetic resonance, Huygens’ metasurfaces are imparted with high transmission and full phase coverage of 2π, which makes them capable of realizing high-efficiency wave-front control. However, a defect of Huygens’ metasurfaces is that their phase profiles and transmissive responses are often sensitive to the interaction of neighboring Huygens’ elements. Consequently, the original assigned phase distribution can be distorted. In this work, we present our design strategy of transmissive Huygens’ metasurfaces performing anomalous refraction. We illustrate the investigation of Huygens’ elements, realizing the overlapping between an electric dipole and magnetic dipole resonance based on cross-shaped structures. We find that the traditional discrete equidistant-phase design method is not enough to realize a transmissive Huygens’ surface due to the interaction between neighboring Huygens’ elements. Therefore, we introduce an extra optimization process on the element spacing to palliate the phase distortion resulting from the element interaction. Based on this method, we successfully design unequally spaced three-element transmissive metasurfaces exhibiting anomalous refraction effect. The anomalous refractive angle of the designed Huygens’ metasurface is 30°, which exceeds the angles of most present transmissive Huygens’ metasurfaces. A transmissive efficiency of 83.5% is numerically derived at the operating wavelength. The far-field electric distribution shows that about 93% of transmissive light is directed along the 30° refractive direction. The deflection angle can be tuned by adjusting the number of Huygens’ elements in one metasurface unit cell. The design strategies used in this paper can be inspiring for other functional Huygens’ metasurface schemes.

    1 INTRODUCTION

    The shaping and engineering of electromagnetic wavefronts at subwavelength scale have attracted much attention due to their extensive applications in recent years. The control of propagating wavefront enables functional optical devices such as lens, anomalous reflectors or refractors, vortex beam exciters, cloaking, and polarization convertors to name a few [18]. These applications are designed to be imparted with particular phase profiles to realize abundant functions. In optical scale, phase engineering can be achieved by the imposition of discrete phase changes on a wavefront using arranged micro- or nano-antenna arrays. Each element inside these antenna arrays acts as an individual phase convertor and changes an incident phase response into a required reflected or transmitted phase response, utilizing their intrinsic electromagnetic resonances or optical modes.

    According to the literature, metasurfaces can be roughly divided into two groups, namely plasmonic and dielectric metasurfaces. Plasmonic metasurfaces are plagued by huge intrinsic ohmic loss [911], which leads to a decrease in their working efficiency and prohibits their application in high-efficiency systems. Dielectric metasurfaces, however, are characterized by very low ohmic loss and are preferred in high-efficiency systems. Generally speaking, there are three groups of dielectric metasurfaces. The first group employs Fabry–Perot elements [1215]. These elements can be regarded as resonating Fabry–Perot cavities, where light reflects back and forth at both ends of these structures. For the Fabry–Perot type metasurfaces, their elements are often required to have enough height to support Fabry–Perot resonances. As a result, they often possess large aspect ratios which make them relatively difficult to produce at nano-scale. The second group of metasurfaces often utilizes circularly polarized incident-light illuminated resonant elements to realize Pancharatnam–Berry phase shift [1619]. These types of metasurfaces often have high efficiency [20], while the incident light is often limited only to circularly polarized light. The last group of metasurfaces, and the one discussed in this paper, is Huygens’ metasurfaces. Huygens’ metasurfaces adopt the spectral and spatial overlapping between an electric dipole (ED) and a magnetic dipole (MD) resonance to achieve phase change [2126,27]. The accumulation of π phase shift resulting from ED and MD resonance together leads to a 2π phase shift at the ED–MD overlapping wavelength. According to Kerker effect, the overlapping between the ED and MD resonance results in forward-scattering performance [28,29]. Consequently, the transmission of the Huygens’ elements is enhanced to near unity. Compared with the other two types of metasurfaces, Huygens’ metasurfaces often have low aspect ratios and can be designed to work in any incident conditions, which can compensate for the shortages of the other two types of metasurfaces.

    However, Huygens’ metasurfaces are plagued by other problems. For example, the introduction of substrates to the Huygens’ elements leads to a variation in MD and ED resonant radiative patterns [30]. As a result, the Kerker condition is no longer satisfied, and the spectral phase response and transmittance of Huygens’ elements will change. Therefore, some publications have designed substrate-free Huygens’ metasurfaces [31,32], which, however, are hard to apply to real optical systems. The interaction between Huygens’ elements can be another crucial obstacle to the realization of functional metasurfaces. It is well known that Mie resonances arise from induced polarized charges and circular currents. The interaction between neighboring Huygens’ elements will no doubt change the distributions of the induced charges and currents. Consequently, the phase response and transmission of Huygens’ elements inside a designed metasurface can be different from those of isolated Huygens’ elements. Therefore, some publications only introduce how to realize ED–MD spectral overlapping or the design of isolated Huygens’ elements, but they do not realize the design of functional Huygens’ metasurfaces [21,22,26]. Some other publications adopt large element spacing to decrease the influence of element interaction or use a large number of Huygens’ elements to maintain phase continuity. However, their functionalities or efficiency is hence limited due to these arrangements. For example, the anomalous refractive angles of present Huygens’ metasurfaces are often less than 20° [2325], which greatly limits their real applications.

    In this paper, we aim at both the design of Huygens’ elements endowed with spectrally overlapped ED and MD resonance, and the realization of efficient transmissive Huygens’ metasurfaces performing anomalous refraction. We will start with the analysis of the ED and MD resonance inside dielectric cross-shaped particles to derive Huygens’ elements performing full 2π phase response. We will discuss the spectral reflection response and the phase response of these elements, and identify the corresponding elements’ parameters required for full phase coverage. Then we take a step further and discuss the design of anomalous refracting Huygens’ metasurfaces. Considering that the inevitable interaction between these elements can have an undesired impact on the phase and refractive response of the metasurface, we will make an improvement on the traditional Huygens’ metasurface design process and realize functional Huygens’ metasurface design. In the end, we will present that our method is also efficient in designing metasurfaces with more elements, through which the anomalous refractive angles of the metasurfaces can be controlled.

    2 DESIGN STRATEGY

    A. ED and MD inside an Isolated Cross-Shaped Particle

    For a single element of Huygens’ metasurfaces, one of the major requisites is at least one structural variable closely related with the ED and MD resonances inside the structure. Specifically speaking, the ED and MD resonant wavelengths should spectrally overlap with the variation of a certain structural parameter. Some regular-shaped structures, such as spheres and cubes, cannot be applied to the Huygens’ elements because the ED and MD resonances inside these structures are both directly proportional to the size parameters and cannot overlap without any other manipulating processes. The basic element structure chosen for metasurface design in this paper is depicted in Fig. 1(a). Instead of commonly used cuboid- or cylinder-shaped structures, the cross-shaped structure is applied as the fundamental element of our metasurface. Compared with other structures, the cross-shaped structure has certain particular characteristics. At normal TM incidence, its two orthogonal rectangle arms direct towards the incident electric and magnetic fields, respectively. The induced electric dipole moment is along the x direction while the induced magnetic dipole moment is along the y direction, as shown in Fig. 1(c). The orthogonality of the two arms ensures relative independent variation for the ED and MD moments. The variation of one arm of the cross-shaped structure will lead to a stronger influence for one dipole moment and a weaker influence for the other. The inconsistent influence is very convenient for designing ED–MD overlapping elements.

    (a) Sketch of an isolated cross-shaped silicon structure. Fixed parameters include ly=700 nm, wy=150 nm, h=400 nm, and wx=140 nm. (b) Simulation unit cell of a periodical cross-shaped structure array, consisting of a cross-shaped particle, substrate, and coating layer. The thicknesses of the substrate and coating layer are 200 and 500 nm, respectively, and Px=Py=1 μm. (c) Sketch of the dipolar moments induced in the two arms of the cross-shaped particle. The gray-colored area depicts the top view of the particle. (d) Sketch of the proposed refracting metasurface. The three cross-shaped elements in the structure are marked as C1, C2, and C3, respectively. Fixed parameters are the same as in (a). The refractive index of the substrate is set as 1.45, while the refractive index of the coating layer is set as 1.4. The refractive index of the silicon is set as 3.5. The ordinate origin is set as the center of the cross particle for (a) and (b), while it is set as the center of the cross particle C2 for the structure in (d).

    Figure 1.(a) Sketch of an isolated cross-shaped silicon structure. Fixed parameters include ly=700  nm, wy=150  nm, h=400  nm, and wx=140  nm. (b) Simulation unit cell of a periodical cross-shaped structure array, consisting of a cross-shaped particle, substrate, and coating layer. The thicknesses of the substrate and coating layer are 200 and 500 nm, respectively, and Px=Py=1  μm. (c) Sketch of the dipolar moments induced in the two arms of the cross-shaped particle. The gray-colored area depicts the top view of the particle. (d) Sketch of the proposed refracting metasurface. The three cross-shaped elements in the structure are marked as C1, C2, and C3, respectively. Fixed parameters are the same as in (a). The refractive index of the substrate is set as 1.45, while the refractive index of the coating layer is set as 1.4. The refractive index of the silicon is set as 3.5. The ordinate origin is set as the center of the cross particle for (a) and (b), while it is set as the center of the cross particle C2 for the structure in (d).

    For an isolated particle, the ED and MD scattering cross sections (SCSs) can be used to present the ED and MD resonant wavelengths and strengths. For the cross-shaped structures, the commonly known Mie theory is no longer efficient in acquiring the desired scattering cross sections. Here, we use the electromagnetic multipole expansion (EME) method to calculate the corresponding dipolar coefficients [33]. Electric and magnetic multipolar coefficients inside an arbitrary particle can be expressed as aE(l,m)=(i)(l1)k2ηOlmE0[π(2l+1)]1/2exp(imϕ){[Ψl(kr)+Ψl(kr)]·Plm(cosθ)r^·JS,j(r)+Ψl(kr)kr[τlm(θ)θ^·JS,j(r)iπlm(θ)ϕ^·JS,j(r)]}d3raM(l,m)=(i)(l1)k2ηOlmE0[π(2l+1)]1/2exp(imϕ)jl(k,r)[τlm(θ)ϕ^·JS,j(r)+iπlm(θ)θ^·JS,j(r)]}d3r,where l represents the multipole order; m is an integer; Ψl(kr) and Ψl(kr) are the first and second derivatives of the Riccati–Bessel functions Ψl(kr)=krjl(kr) with respect to the argument kr; Olm, πlm, and τlm are parameters [33,34]; and Plm are the associated Legendre polynomials. JS,j(r) is the effective current density that creates the scattered field of the j-th particle in the self-consistent solution of Maxwell equations. The integrals in Eq. (1) are operated only inside the particle. Out of the particle, integrands are equal to zero. The current density distributions and the integrals in Eq. (1) are calculated by the numerical simulation software, COMSOL Multiphysics. As for simulation, the particles are placed in the air, out of which it is a thick perfectly matched layer (PML). The boundary between air and PML is set as scattering boundary condition (SBC). The mesh is manually set with a minimum mesh size of 5 nm. The incident field is set as a normal plane field. The operating condition is set as a built-in scattering-field environment, where the incident electric field can be omitted from the total field. The scattering cross section for a particle can be derived through [33] Cs=πk2l=1m=ll(2l+1)[|aE(l,m)|2+|aM(l,m)|2].The terms of the series in Eq. (2) allow one to determine the contribution of each multipole excitation to the overall SCS. For example, scattering contribution of ED and MD can be calculated through σED=3πk2m=11|aE(1,m)|2,σMD=3πk2m=11|aM(1,m)|2,namely when l=1 for Eq. (2).

    In Fig. 2(a), we plot the calculated scattering cross sections of ED and MD components inside isolated cross-shaped structures as a function of the x arm length lx. The intensities of the ED and MD resonances both increase with the increment of lx, while their resonant wavelengths both redshift with the increment of lx. Such effect is known as the size effect [35]. The intensity of the ED resonance increases more prominently than the intensity of MD resonance when lx increases. In Fig. 2(b), we plot the ED and MD resonant wavelengths as a function of lx. Apparently, the ED resonant wavelength redshifts relatively quicker than the MD resonant wavelength, resulting in spectral approaching between the ED and MD resonance with the increment of lx, whereby the spectral spacing between ED and MD reduces from 200 to about 100 nm with the variation of the parameter lx. This result proves that ED resonance is more active to the change of the x arm length lx, compared with MD resonance. Such effect can be further applied to realize ED–MD overlapping.

    (a) Scattering cross sections of ED and MD components varying with x arm length lx. (b) ED and MD resonant wavelengths as a function of lx for an isolated cross-shaped particle.

    Figure 2.(a) Scattering cross sections of ED and MD components varying with x arm length lx. (b) ED and MD resonant wavelengths as a function of lx for an isolated cross-shaped particle.

    B. ED and MD inside Periodically Arranged Cross-Shaped Structures

    In periodically arranged array, ED and MD resonances exist in the form of coupled dipole resonances called the array effect. When the particle size is far smaller than the incident wavelength, the particle array can be regarded as a dipolar moment array with efficient polarizabilities. The efficient polarizability can be calculated as [36] ε0/αeffED=ε0/αEDks2Gxx0/2,1/αeffMD=1/αMDks2Gyy0/2,where ks is the wave vector in the surrounding environment; Gxx0 and Gyy0 are the diagonal elements of a Green dyadic G0 [37,38]; αED and αMD are the electric and magnetic polarizabilities, respectively, of an isolated particle extracted from the periodical array. The terms ks2Gxx0/2 and ks2Gyy0/2 depict the influence of the ED–ED and MD–MD couplings inside an array, and are mainly determined by the lattice constants Px and Py. Equation (4) makes a connection between the dipolar resonances of an isolated particle and the dipolar resonances inside its periodical array. The array effect can bring about spectral shifting for ED and MD resonances [36]. Ignoring the influence of the higher-order multipoles, the spectral reflective and transmissive coefficients for a lossless particle array can be expressed by the efficient ED and MD polarizabilities as [36] r=iks2DL(αeffEDαeffMD),t=1iks2DL(αeffEDαeffMD),where DL is the area of one unit cell. When the ED and MD resonances spectrally overlap with the same intensity, namely αeffED=αeffMD, the reflection is zero while the transmission is unity. This condition corresponds to the Kerker effect for a single particle, where ε=μ is required. The particle inside the array can be then used to build a Huygens’ metasurface. Equation (5) also indicates that a single dipole resonance in the periodical element array brings about a reflection peak, while the overlapping between an ED and MD resonance leads to a transmission peak. However, elements inside a real metasurface cannot all be designed to be in ED–MD overlapping state, considering that various phase responses are required. Therefore, some elements are working at the quasi-ED–MD overlapping state or single dipolar resonant state. As a result, the efficiency of a transmissive Huygens’ metasurface is hard to reach unity.

    For applicable metasurfaces, substrates are often required to support periodical elements and keep them flush with surrounding elements. However, the introduction of substrates will result in the distortion of dipolar resonant radiative patterns [30]. The original ED or MD polarizability αED or αMD will change, and the overlapping condition αeffED=αeffMD will no longer be satisfied. Consequently, the directional or quasi-directional scattering performances vanish, and the transmissive efficiency decreases. To eliminate the influence of the substrate as much as possible, the particles above the substrate can be coated with another dielectric layer, with a close refractive index to the substrate. The sketch of the particle attached with a substrate and coating layer is plotted in Fig. 1(b). The refractive index of the substrate is set as 1.45, while the refractive index of the coating layer is set as 1.4. With this process, the particle can be regarded as being embedded inside a homogeneous layer. Such arrangement only changes the background environment of a particle, namely only the wave vector ks changes in Eq. (4). Consequently, the ED and MD resonances inside the array both spectrally shift while their radiative pattern changes slightly, and directional scattering can still be realized. Other publications [21,24] also adopt this arrangement, but the implied physics mechanism is not demonstrated.

    In Figs. 3(a) and 3(b), the spectral reflection and phase response of the cross-shaped structure imposed with the substrate and coating layer in Fig. 1(b) are plotted as a function of the x arm length lx. As mentioned before, the ED and MD resonances induce reflection peaks in periodical particle arrays [36,39]. Therefore, we use spectral reflection responses here to address the ED and MD resonant spectral wavelengths inside the periodical structure. As for the simulation of the transmission and phase responses of periodical arrays, we use the software Lumerical FDTD here to derive the required data. Periodical boundary conditions (PBCs) are adopted at both the x and y boundaries of the simulated structures. Perfectly matched layers are imposed on the top and bottom of the simulated unit cell. A minimum mesh size of 5 nm is applied, ensuring that the calculation is independent of the mesh size.

    (a) Spectral reflection and (b) spectral phase as a function of the structure parameter lx in Fig. 1(b).

    Figure 3.(a) Spectral reflection and (b) spectral phase as a function of the structure parameter lx in Fig. 1(b).

    In Fig. 3(a), there are two distinguishable resonant peaks, which are spectrally close to the resonant peaks of ED and MD in Fig. 2(b). To prove these resonant peaks to be Mie resonances, we plot the corresponding electric field distributions in Fig. 4. The electric field vector distributions in Figs. 4(c) and 4(d) confirm the resonant modes to be ED and MD resonances. Referring back to Fig. 3(a), the ED and MD resonances both spectrally redshift with the increment of the parameter lx. However, the ED resonance shows a more active spectral shifting effect than the MD counterpart. Consequently, the two resonant peaks merge at around 1.49 μm. According to the Kerker effect, the spectral overlapping between ED and MD resonance leads to directional forward scattering. Therefore, the two reflection peaks change into a transmission peak at the ED–MD overlapping wavelength in Fig. 3(a). In Fig. 3(b), the phase of the transmitted light is plotted as a function of the structural parameter lx. According to the numerical results, there is a π phase-shifting at both ED and MD resonances. When the two resonances spectrally overlap, the total phase shifting of the transmitting light changes to 2π. Applying the low reflection and the swift phase change at around 1.5 μm in Figs. 3(a) and 3(b), we can further design the desired metasurfaces.

    Electric field magnitude distributions |E|2 at the (a) first reflection peak and (b) second reflection peak in Fig. 3(a). Electric field vector distribution at the (c) first reflection peak and (b) second reflection peak in Fig. 3(a). The incident electric field is along the negative x direction. The variable lx is set as 0.27 μm. (a)–(d) are the electric field distributions at the y=0 plane.

    Figure 4.Electric field magnitude distributions |E|2 at the (a) first reflection peak and (b) second reflection peak in Fig. 3(a). Electric field vector distribution at the (c) first reflection peak and (b) second reflection peak in Fig. 3(a). The incident electric field is along the negative x direction. The variable lx is set as 0.27 μm. (a)–(d) are the electric field distributions at the y=0 plane.

    C. Metasurface Design

    In general, anomalous refracting angle decreases with increment of the elements inside a periodical unit cell of a metasurface. To derive a large reflecting angle, we set the number of elements inside one unit cell as 3. Traditional Huygens’ metasurface design requires discrete equidistant-phase Huygens’ elements to realize certain phase profiles. Therefore, three cross-shaped structures with an equidistant phase discretization are selected to form the metasurface. The variables of the three cross-shaped structures and the spectral responses are listed in Table 1. The three elements are designed to together realize a 2π phase shifting while maintaining a relatively high transmission (>87%). The total periodicity of the metasurface is set as P3x=3  μm, which is three times the periodicity of a single element (1 μm). The phase and transmission are measured at a selected wavelength λ=1.507  μm. With these arrangements, the metasurface is supposed to realize a linear phase profile at around 1.507 μm. In Fig. 1(d), we plot the sketch of the designed metasurface. The parametric configuration here is marked as configuration 1.

    Elementlx (μm)TransmissionPhase (rad)Ideal Phase (rad)Real Phase conf. 1Real Phase conf. 2Real Phase conf. 3
    E10.240.872.102.09 (2π/3)3.42 (rad)2.17 (rad)2.12 (rad)
    E20.390.974.214.19 (4π/3)0.575 (rad)5.10 (rad)4.31 (rad)
    E30.430.976.236.28 (2π)0.81 (rad)6.20 (rad)6.15 (rad)

    Table 1. First Configuration of the Phase Gradient Metasurfacea

    In Fig. 5(e), we plot the phase distribution of the designed metasurface at 1.493 μm, where the phase profile is closest to the predesigned linear phase profile. Apparently, the phase of the transmitted light distributes irregularly along the x direction. The real phase response of each element inside the unit cell for configuration 1 is calculated and listed in Table 1. The real phase responses of these elements in the metasurface unit cell turn out to be quite different from the phase responses of these elements in their corresponding periodical arrays. This result can be attributed to the interaction between these elements. In Figs. 5(a)5(d), we respectively plot the electric field distributions of the three periodically arranged elements, together with the electric field distribution of the metasurface unit cell. Apparently, the interaction between the three particles changes the original field distribution. Consequently, the phase response and transmission of the metasurface will both get changed. Instead of a linear phase profile, the phase profile in Fig. 5(e) shows an irregular change. The transmissions of the original three elements are 87%, 97%, and 97%, respectively, while the transmission of the metasurface plotted in Fig. 5(f) (conf. 1) shows only 78% transmission at 1.493 μm. Such variation indicates that undesired changes occur to both the phase response and transmission of the metasurface due to element–element interaction. A prominent transmission dip, resulting from the Rayleigh anomaly [40], occurs at around 1.5 μm, where the periodicity of the metasurface and the incident wavelength fulfill the relationship P3x=2λ. Figures 5(a)5(e) prove that the traditional design method is no longer efficient in our design scheme. Other methods need to be utilized to realize the desired metasurfaces.

    (a)–(c) Electric field distributions of the periodically arranged cross-shaped structures E1, E2, and E3, respectively, at 1.507 μm. (d) Electric field distribution of one unit cell of the metasurface configuration 1, consisting of three elements, namely E1, E2, and E3. (a)–(d) depict the electric field distributions at the y=0 plane. (e) Phase distribution of the metasurface configuration 1 at 1.493 μm. Three unit cells (UCs), marked as UC1, UC2, and UC3, are plotted to have a comfortable aspect ratio for the figure. (f) Spectral transmission for the proposed three configurations.

    Figure 5.(a)–(c) Electric field distributions of the periodically arranged cross-shaped structures E1, E2, and E3, respectively, at 1.507 μm. (d) Electric field distribution of one unit cell of the metasurface configuration 1, consisting of three elements, namely E1, E2, and E3. (a)–(d) depict the electric field distributions at the y=0 plane. (e) Phase distribution of the metasurface configuration 1 at 1.493 μm. Three unit cells (UCs), marked as UC1, UC2, and UC3, are plotted to have a comfortable aspect ratio for the figure. (f) Spectral transmission for the proposed three configurations.

    In periodical structures, the interaction is mainly determined by the distance between neighboring elements, the geometries of the elements, and the periodicity of the unit cell. In most metasurface design schemes [1125], a stereotype is that metasurface elements have to be linearly arranged with equal spacing in the phase–gradient shifting plane. When the spacing between elements is unequal, the phase profile is distorted. However, the phase profile for the Huygens’ metasurface designed here is already distorted due to the intrinsic interaction between elements. Instead, the variation of the spacing between elements may have an unexpected influence on the phase profile. With this hypothesis, we then consider how to realize specific phase profiles through element spacing control. However, the interaction between the three non-simple elements is complex, and there is no analytical theory illustrating how to adjust these parameters on purpose. Therefore, we use a global optimization method here to investigate our assumption.

    To reduce the complexity of calculation, some main parameters lx1, lx2, lx3, d1, and d2 in Fig. 1(d) are chosen to be the optimizing parameters. The periodicity is fixed as P3x=3 and Py=1  μm. Of particular importance is that the y-arm parameters wy and ly are fixed during optimization, since these two parameters seriously affect the MD resonant response of the structure. Considering the extensive possible combinations of these parameters, we take the step-by-step optimization method to orderly optimize the surface. In the first optimization process, the element spaces d1 and d2 are optimized first. We consider both the equal spacing condition, d1=d2, and the unequal spacing condition, namely d1d2. After numerical calculation, the best result of linear phase profile occurs when d1=d2=0.86  μm. The derived parameter configuration is represented as configuration 2. The sketches of configurations 1 and 2 are plotted in Figs. 6(a) and 6(b), respectively, for intuitive understanding. The phase map of configuration 2 is plotted in Fig. 6(d). Compared with the phase profile of configuration 1 in Fig. 5(e), Fig. 6(d) shows a better phase gradient changing profile. A roughly linear phase-changing profile can be observed. In Fig. 5(f), we plot the transmission of configuration 2. At 1.495 μm, 82% transmission is achieved for configuration 2.

    (a) Top view of the metasurface configuration 1. (b) Top view of the metasurface configuration 2. The parameters of each element stay unchanged, while only the spaces d1 and d2 between three elements change from 1 to 0.86 μm. (c) Top view of the metasurface configuration 3. Element parameters lx are optimized. The optimized parameters lx of the three elements are 0.3, 0.4, and 0.56 μm, respectively. (d) Calculated phase distribution of the second configuration of the metasurface; the operating wavelength is 1.495 μm. (e) Calculated phase distribution of the third configuration of the metasurface; the incident wavelength is 1.521 μm.

    Figure 6.(a) Top view of the metasurface configuration 1. (b) Top view of the metasurface configuration 2. The parameters of each element stay unchanged, while only the spaces d1 and d2 between three elements change from 1 to 0.86 μm. (c) Top view of the metasurface configuration 3. Element parameters lx are optimized. The optimized parameters lx of the three elements are 0.3, 0.4, and 0.56 μm, respectively. (d) Calculated phase distribution of the second configuration of the metasurface; the operating wavelength is 1.495 μm. (e) Calculated phase distribution of the third configuration of the metasurface; the incident wavelength is 1.521 μm.

    Based on configuration 2, the parameter lx of the three elements is further optimized. A sketch in Fig. 6(c) summarizes the difference between configurations 2 and 3. The optimized phase profile is plotted in Fig. 6(e). Compared with the two former results, this latest result has a better linear phase-changing effect. In Fig. 5(f), the transmission of configuration 3 is 83.5% at 1.521 μm, where the best linear phase profile is achieved. These results indicate that the optimization processes are very effective and necessary in our designing scheme. According to our simulation, the two-step optimization processing can be taken inversely. The same result as Fig. 6(e) can also be realized by inversing the two optimization processes. The real phase response of the elements under configurations 2 and 3 are listed in Table 1. Clearly, the real phase response approaches the ideal phase response when the optimization processes are imposed on the metasurface from configuration 1 to configurations 2 and 3.

    In our two-step optimization process, the element–parameter optimization corrects the originally distorted phase response of the three elements inside the metasurface unit cell. The optimized three elements present near-ideal phase responses. The element-spacing optimization can improve the linearity of metasurface phase profile. To better show the effect of the element-spacing optimization, we plot the transmissive response and phase profile of the surface varying with the element spacing d in Figs. 7(a) and 7(b). For simplicity, we only show the result at the condition d=d1=d2. According to Fig. 7(b), the structure can realize a 2π phase shift at any element spacing d, but the corresponding linearity of the 2π phase shift is varied. This indicates that the element-spacing variation has limited change on the phase profile of the designed metasurface. However, in Fig. 7(a), we can see that the transmission of the metasurface changes from 40% to near 90% with the variation of the element spacing d. This result implies that the efficiency of the metasurface can be efficiently optimized by adjusting the element spacing. The adjusting of the element spacing leads to constructive or destructive interference to the elements. According to Figs. 7(a) and 7(b), the interference seriously affects the transmissive response for the metasurface.

    (a) Transmissive response of the metasurface varying with the element spacing d. Black dotted line marks the applied working wavelength. (b) Phase response along the x direction inside a metasurface unit cell varying with the element spacing d at working wavelength. The parameters lx of the three elements are 0.3, 0.4, and 0.56 μm, respectively.

    Figure 7.(a) Transmissive response of the metasurface varying with the element spacing d. Black dotted line marks the applied working wavelength. (b) Phase response along the x direction inside a metasurface unit cell varying with the element spacing d at working wavelength. The parameters lx of the three elements are 0.3, 0.4, and 0.56 μm, respectively.

    According to numerical calculation, the anomalous deflecting angle of the proposed metasurface is about 30°, and the transmission is 83.5%. In Fig. 8, far-field transmitting intensities are plotted under the three configurations. Apparently, the light deflected into the ±2, 0, and 1 diffraction orders is well suppressed with the optimizing processing from configuration 1 to 3. The final result, namely the third configuration in Fig. 8(c), contains a major +1 ordered diffracted light and minor 1 and 0 ordered diffracted light. In Figs. 8(a)8(c), the transmitting intensity is normalized to the total transmitting intensity. Our calculation shows that the sum of the transmitting intensity at around +1 ordered diffraction possesses 93% of the total transmitting intensity, while the total transmitting intensity of 1 and 0 ordered diffraction together possesses only 7% of the total transmitting intensity.

    Far-field transmitting intensity (T.I.) for the three configurations at different diffraction angles. The transmitting intensity is normalized to the total transmission intensity.

    Figure 8.Far-field transmitting intensity (T.I.) for the three configurations at different diffraction angles. The transmitting intensity is normalized to the total transmission intensity.

    According to the generalized Snell’s law [1], the angle of anomalous refraction can be theoretically predicted by θt=arcsin{1nt[λ0P+(sinθi)ni]},where nt and ni are the refractive index of the transmitting and incident environments, respectively; P is the periodicity of the array in the phase-gradient direction; and λ0 is the operating wavelength. In our situation, the incident and transmitting environments are set as air, namely nt=ni=1. The influence of the substrate and coating layer on the diffraction angle can be neglected, as their thicknesses are less than half of the incident wavelength. At the normal incidence condition, the incident angle θi is 0°, the operating wavelength is 1.521 μm, and the periodicity is 3 μm in the x direction. The theoretically calculated anomalous deflection angle is 30.5°, which is in good agreement with the simulation result (30°). Compared with many other Huygens’ phase-gradient metasurfaces, the proposed metasurface has a larger deflecting angle while the phase-change linearity is well. The excellent performance is likely owing to two main reasons. The first is that the number of the surface unit cell elements is small; only three elements are applied to build the phase gradient surface. Another reason is that the element-spacing optimization process is used to better control the phase profile of the metasurface.

    According to Eq. (6), the anomalous refraction angle can be tuned by changing the periodicity of the metasurface. However, the periodicity of Huygens’ metasurface influences the interaction between Huygens’ elements. A large adjustment of periodicity can easily bring about distortion for the phase distribution of the metasurfaces. Although minor adjustment of periodicity maintains the functionality of the metasurface, the anomalous refraction angle can only be tuned in a very limited region (less than 1°). An efficient way to realize different anomalous refraction angles is to add elements to expand the periodicity of the metasurface. Therefore, we further design metasurfaces with four, five, and six Huygens’ elements. The optimized results are plotted in Fig. 9, showing that near-linear phase profiles can be realized in four, five, and six-element Huygens’ metasurfaces. The corresponding parameters and detailed numerical results are listed in Table 2. The refraction angles range from 30° to 15.8°. In our optimizing process, the element geometric parameter lx, periodicity Px, and space d between neighboring elements are the main optimizing parameters. A shortcoming of the proposed optimization method is that the optimizing parameters increase with the number of elements. Therefore, the optimizing process is too complex when all the optimizing parameters are considered. These results prove that our design strategy is effective for more than one specific case.

    Phase distribution of the Huygens’ metasurfaces consisting of different element numbers. (a) Four-element Huygens’ metasurface. (b) Five-element Huygens’ metasurface. (c) Six-element Huygens’ metasurface. Each figure consists of two unit cells to better show the anomalous deflection effect. The spacing between neighboring elements in these three configurations is set as 0.86 μm for simplicity. Other parameters are listed in Table 2.

    Figure 9.Phase distribution of the Huygens’ metasurfaces consisting of different element numbers. (a) Four-element Huygens’ metasurface. (b) Five-element Huygens’ metasurface. (c) Six-element Huygens’ metasurface. Each figure consists of two unit cells to better show the anomalous deflection effect. The spacing between neighboring elements in these three configurations is set as 0.86 μm for simplicity. Other parameters are listed in Table 2.

    E. N.lx1lx2lx3lx4lx5lx6Pλ0θt
    30.30.40.5631.52130°
    40.270.380.420.563.81.50623.4°
    50.230.350.380.40.551.49517.4°
    60.230.330.380.420.460.565.61.52415.8°

    Table 2. Crucial Parameters of the Huygens’ Metasurfaces Consisting of Different Numbers of Huygens’ Elementsa

    3 CONCLUSION

    In this paper, we propose an all-dielectric three-element Huygens’ metasurface realizing anomalous refraction based on cross-shaped structures. Utilizing the multipole decomposition method, we derive the ED and MD resonant wavelengths and intensities inside cross-shaped structures. We find that the ED resonance shows a more active response to the change of the x-directed arm of the cross particle on resonant intensity and wavelengths than the MD resonance. Utilizing this performance, we realize the spectral overlapping between the ED and MD resonances. We prove that the traditionally fixed element spacing inside a metasurface can be an efficient phase-profile-adjusting parameter. Through a two-step optimization process, the traditionally designed Huygens’ metasurface can realize desired anomalous refraction. The proposed three-element metasurface can refract a normal TM incident beam into a 30° oblique transmissive beam with a transmissive efficiency of 83.5%.

    We also design metasurfaces with different element numbers according to our design strategy. By increasing the element number in the metasurface, the anomalous refractive angle can be tuned in the range from 15.8° to 30°. Our results prove that functional Huygens’ metasurfaces with small element numbers can be efficiently realized through element–element spacing control. The methods and solutions used in this paper can help to design other high-efficiency metasurfaces.

    References

    [1] N. Yu, P. Genevet, M. A. Kats, F. Aieta, J. P. Tetienne, F. Capasso, Z. Gaburro. Light propagation with phase discontinuities: generalized laws of reflection and refraction. Science, 334, 333-337(2011).

    [2] L. Philippe, C. Pierre. Metalenses at visible wavelengths: past, present, perspectives. Laser Photon. Rev., 11, 1600295(2017).

    [3] A. E. Minovich, A. E. Miroshnichenko, A. Y. Bykov, T. V. Murzina, D. N. Neshev, Y. S. Kivshar. Functional and nonlinear optical metasurfaces. Laser Photon. Rev., 9, 195-213(2015).

    [4] S. Sun, K. Y. Yang, C. M. Wang, T. K. Juan, W. T. Chen, C. Y. Liao, Q. He, S. Xiao, W. T. Kung, G. Y. Guo. High-efficiency broadband anomalous reflection by gradient meta-surfaces. Nano Lett., 12, 6223-6229(2012).

    [5] A. Amir, H. Yu, J. B. Alexander, B. Mahmood, F. Andrei. Subwavelength-thick lenses with high numerical apertures and large efficiency based on high-contrast transmitarrays. Nat. Commun., 6, 7069(2015).

    [6] X. Ni, Z. J. Wong, M. Mrejen, Y. Wang, X. Zhang. An ultrathin invisibility skin cloak for visible light. Science, 349, 1310-1314(2015).

    [7] Y. Huang, Q. Zhao, S. K. Kalyoncu, R. Torun, Y. Lu, F. Capolino, O. Boyraz. Phase-gradient gap-plasmon metasurface based blazed grating for real time dispersive imaging. Appl. Phys. Lett., 104, 161106(2014).

    [8] Z. Wang, X. Ding, K. Zhang, B. Ratni, S. N. Burokur, X. Gu, Q. Wu. Huygens metasurface holograms with the modulation of focal energy distribution. Adv. Opt. Mater., 6, 1800121(2018).

    [9] N. Yu, F. Capasso. Flat optics with designer metasurfaces. Nat. Mater., 13, 139-150(2014).

    [10] Y. Zhao, X.-X. Liu, A. Alù. Recent advances on optical metasurfaces. J. Opt., 16, 035403(2014).

    [11] V. N. Gururaj, M. S. Vladimir, B. Alexandra. Alternative plasmonic materials: beyond gold and silver. Adv. Mater., 25, 3264-3294(2013).

    [12] A. Amir, H. Yu, B. Mahmood, F. Andrei. A dielectric metasurfaces for complete control of phase and polarization with subwavelength spatial resolution and high transmission. Nat. Nanotechnol., 10, 937-944(2015).

    [13] M. Khorasaninejad, A. Y. Zhu, C. Roques-Carmes, W. T. Chen, J. Oh, I. Mishra, R. C. Devlin, F. Capasso. Polarization-insensitive metalenses at visible wavelength. Nano Lett., 16, 7229-7234(2016).

    [14] E. Arbabi, A. Arbabi, S. M. Kamali, Y. Horie, A. Faraon. High efficiency double wavelength dielectric metasurface lenses with dichroic birefringent meta-atoms. Opt. Express, 24, 18468-18477(2016).

    [15] S. S. Kruk, B. Hopkins, I. I. Kravchenko, A. Miroshnichenko, D. N. Neshev, Y. S. Kivshar. Broadband highly efficient dielectric metadevices for polarization control. APL Photon., 1, 030801(2016).

    [16] Z. Guoxing, M. Holger, K. Mitchell, L. Guixin, Z. Thomas, Z. Shuang. Metasurface holograms reaching 80% efficiency. Nat. Nanotechnol., 10, 308-312(2015).

    [17] L. Dianmin, F. Pengyu, H. Erez, L. B. Mark. Dielectric gradient metasurface optical elements. Science, 345, 298-302(2016).

    [18] R. C. Devlin, K. Mohammadreza, W. T. Chen, J. Oh, F. Capasso. Broadband high efficiency dielectric metasurfaces for the visible spectrum. Proc. Natl. Acad. Sci. USA, 113, 10473-10478(2016).

    [19] E. Maguid, I. Yulevich, M. Yannai, V. Kleiner, M. L. Brongersma, E. Hasman. Multifunctional interleaved geometric-phase dielectric metasurfaces. Light Sci. Appl., 6, e17027(2017).

    [20] M. Khorasaninejad, W. T. Chen, R. C. Devlin, J. Oh, A. Y. Zhu, F. Capasso. Metalenses at visible wavelengths: diffraction-limited focusing and subwavelength resolution imaging. Science, 352, 1190-1194(2016).

    [21] M. Decker, I. Staude, M. Falkner, J. Dominguez, D. N. Neshev, I. Brener, P. Thomas, Y. S. Kivshar. High-efficiency dielectric Huygens’ surfaces. Adv. Opt. Mater., 3, 813-820(2015).

    [22] D. Arslan, K. E. Chong, A. Miroshnichenko, D. Y. Choi, D. N. Neshev, T. Pertsch, Y. S. Kivshar, I. Staude. Angle-selective all-dielectric Huygens’ metasurfaces. J. Phys. D, 50, 434002(2017).

    [23] L. Zhang, J. Ding, H. Y. Zheng, S. S. An, H. T. Lin, B. W. Zheng, D. Y. Du, G. F. Yin, J. Michon, Y. F. Zhang. Ultra-thin high-efficiency mid-infrared transmissive Huygens meta-optics. Nat. Commun., 9, 1481(2018).

    [24] A. J. Ollanik, J. A. Smith, M. J. Belue, M. D. Escarra. High efficiency all-dielectric Huygens metasurfaces from the ultraviolet to the infrared. ACS Photon., 5, 1351-1358(2018).

    [25] M. I. Shalaev, J. B. Sun, A. Tsukernik, A. Pandey, K. Nikolskly, N. M. Litchinister. High-efficiency all-dielectric metasurfaces for ultracompact beam manipulation in transmission mode. Nano Lett., 15, 6261-6266(2015).

    [26] S. Liu, A. Vaskin, S. Campione, O. Wolf, M. B. Sinclair, J. L. Reno, G. A. Keeler, I. Staude, I. Brener. Huygens’ metasurfaces enabled by magnetic dipole resonance tuning in split dielectric nanoresonators. Nano Lett., 17, 4297-4303(2017).

    [27] R. Zhao, Z. Zhu, G. Dong, T. Lv, Y. Li, C. Guan, J. Shi, H. Zhang. High-efficiency Huygens’ metasurface for terahertz wave manipulation. Opt. Lett., 44, 3482-3485(2019).

    [28] M. Kerker, D. S. Wang, C. L. Giles. Electromagnetic scattering by magnetic spheres. J. Opt. Soc. Am., 73, 765-767(1983).

    [29] I. Staude, A. E. Miroshnichenko, M. Decker, N. T. Fofang, S. Liu, E. Gonzales, J. Domingguez, T. S. Luk, D. N. Neshev, I. Brener, Y. Kivshar. Tailoring directional scattering through magnetic and electric resonances in subwavelength silicon nanodisks. ACS Nano, 7, 7824-7832(2013).

    [30] J. S. Eismann, M. Neugebauer, P. Banzer. Exciting a chiral dipole moment in an achiral nanostructure. Optica, 5, 954-959(2018).

    [31] W. Y. Zhao, H. Jiang, B. Y. Liu, J. Song, Y. Y. Jiang. High-efficiency beam manipulation combining geometric phase with anisotropic Huygens surface. Appl. Phys. Lett., 108, 181102(2016).

    [32] Y. Z. Chen, J. Mei. All-dielectric two-dimensional metasurfaces based on electric and magnetic dipolar Mie resonances. Europhys. Lett., 122, 54002(2018).

    [33] P. Grahn, A. Shevchenko, M. Kaivola. Electromagnetic multipole theory for optical nanomaterials. New J. Phys., 14, 093033(2012).

    [34] J. D. Jackson. Classical Electrodynamics(1999).

    [35] E. Ringe, M. R. Langille, K. Sohn, J. Zhang, J. Huang, C. A. Mirkin, R. P. Van Duyne, L. D. Marks. Plasmon length: a universal parameter to describe size effects in gold nanoparticles. J. Phys. Chem. Lett., 3, 1479-1483(2012).

    [36] A. B. Evlyukhin, C. Reinhardt, A. Seidel, B. S. Luk’yanchuk, B. N. Chichkov. Optical response features of Si-nanoparticle arrays. Phys. Rev. B, 82, 45404(2010).

    [37] V. E. Babicheva, A. B. Evlyukhin. Resonant lattice Kerker effect in metasurfaces with electric and magnetic optical responses. Laser Photon. Rev., 11, 1700132(2017).

    [38] S. D. Swiecicki, J. E. Sipe. Surface-lattice resonances in two-dimensional arrays of spheres: multipolar interactions and a mode analysis. Phys. Rev. B, 95, 195406(2017).

    [39] C. Liu, T. Wu, Y. Liu, J. Li, Y. Wang, Z. Yu, H. Ye, L. Yu. Realization of perfect selective absorber based on multipole modes in all-dielectric moth-eye structure. Opt. Express, 27, 5703-5718(2019).

    [40] L. Rayleigh. Note on the remarkable case of diffraction spectra described by Prof. Wood. Philos. Mag., 14, 60-65(1907).

    Chang Liu, Lei Chen, Tiesheng Wu, Yumin Liu, Jing Li, Yu Wang, Zhongyuan Yu, Han Ye, Li Yu, "All-dielectric three-element transmissive Huygens’ metasurface performing anomalous refraction," Photonics Res. 7, 1501 (2019)
    Download Citation