• Photonics Research
  • Vol. 7, Issue 4, 457 (2019)
Ke Bi1, Daquan Yang1, Jia Chen1, Qingmin Wang1, Hongya Wu2, Chuwen Lan1、*, and Yuping Yang3
Author Affiliations
  • 1State Key Laboratory of Information Photonics and Optical Communications, School of Science, Beijing University of Posts and Telecommunications, Beijing 100876, China
  • 2School of Materials Science and Engineering, Shijiazhuang Tiedao University, Shijiazhuang 050043, China
  • 3School of Science, Minzu University of China, Beijing 100081, China
  • show less
    DOI: 10.1364/PRJ.7.000457 Cite this Article Set citation alerts
    Ke Bi, Daquan Yang, Jia Chen, Qingmin Wang, Hongya Wu, Chuwen Lan, Yuping Yang, "Experimental demonstration of ultra-large-scale terahertz all-dielectric metamaterials," Photonics Res. 7, 457 (2019) Copy Citation Text show less

    Abstract

    All-dielectric metamaterials have emerged as a promising platform for low-loss and highly efficient terahertz devices. However, existing fabrication methods have difficulty in achieving a good balance between precision and cost. Here, inspired by the nano-template-assisted self-assembly method, we develop a micro-template-assisted self-assembly (MTAS) method to prepare large-scale, high-precision, and flexible ceramic microsphere all-dielectric metamaterials with an area exceeding 900 cm×900 cm. Free from organic solvents, vacuum, and complex equipment, the MTAS method ensures low-cost and environmentally friendly fabrication. The ceramic microsphere resonators can be readily assembled into nearly arbitrary arrangements and complex aggregates, such as dimers, trimers, quadrumers, and chains. Finally, using the heat-shrinkable substrate and dipole coupling effect, a broadband reflector with a bandwidth of 0.15 THz and a reflection of up to 95% is demonstrated. This work provides a versatile and powerful platform for terahertz all-dielectric metamaterials, with potential to be applied in a wide variety of high-efficiency terahertz devices.

    1. INTRODUCTION

    Manipulating light at a subwavelength level is the cornerstone for many optical applications, such as compact photonic devices and all-optical chips. This demand has driven the rapid development of metamaterials comprising arrayed artificial atoms. Although having exotic properties and wide applications, most of the metamaterials still remain at the laboratory stage. One of the main limitations is the significant ohmic loss associated with plasma resonances in conventional metamaterials, limiting their efficiency for a wide range of applications. Thus, reducing this loss is essential for highly efficient photonic devices. Recently, it was demonstrated that limitations could be addressed by using high-index dielectric resonators, because dielectric resonances are induced by displacement current rather than conduction current. Over the past several years, all-dielectric metamaterials have been investigated to exhibit various novel properties, including zero index [1], optical magnetism [2], local field enhancement [35], and directional scattering [6], which highlight a wide variety of applications in wavefront control [712], sensing [13,14], surface-enhanced spectroscopy [15], reflectors [1618], and absorbers [19].

    Whereas all-dielectric metamaterials have become a hotbed of research, ranging from the infrared to optical regions, they are seldom reported in the terahertz range, which is an important spectral region that has wide applications in wireless communication, sensing, security, and chemical identification. This is different from those made of metallic structures, which have been widely reported in terahertz ranges [2031]. One of the main limitations is the high cost of fabrication techniques. Over the past few years, attempts have been made to obtain terahertz all-dielectric metamaterials. There are three methods: direct laser patterning, deep etching, and a chemical method. The direct laser-patterning method has proven to be a high-precision and material-free fabrication approach to all-dielectric metamaterials [32,33]. However, this method suffers from expensive fabrication, inevitably blurred boundaries, and complex thickness reduction processes. Using deep-etching methods, silicon-based all-dielectric metamaterials in the terahertz range have been demonstrated in several studies [3437]. Whereas these metamaterials show high-precision patterns and high-quality Mie resonances, reducing the high fabrication costs remains an immense challenge. Recently, a cost-effective method was reported, based on high-throughput titanium dioxide spherical resonators using a chemical method. However, the resonators are difficult to assemble, and applications were impeded in some cases [38]. Clearly, the above fabrication methods have difficulty in achieving a good balance between precision and cost.

    In this work, inspired by the template-assisted self-assembly method, we address this limitation by developing a novel approach, called the micro-template-assisted self-assembly (MTAS) method. We demonstrate that this method is capable of fabricating ultra-large-scale, high-precision, and flexible ceramic microsphere all-dielectric metamaterials with an area exceeding 900  cm×900  cm. Additionally, by using proper templates, the ceramic microsphere resonators can be readily assembled into nearly any arrangement, including complex aggregates. Next, the all-dielectric metamaterial is prepared on an active material, with which the microsphere resonators can be assembled into closely packed metasurfaces under heat treatment (HT), leading to a broadband reflector with a bandwidth of 0.15 THz and a reflection up to 95%. This work provides a versatile and powerful platform for terahertz all-dielectric metamaterials, which should find a wide variety of applications in low-loss and highly efficient terahertz devices.

    2. METHOD

    Before introducing the MTAS method, the nano-template-assisted self-assembly (NTSA) method is illustrated in Fig. 1(a) for comparison [39]. The fundamental process for this technology includes several steps. First, a template with holes or other patterns is fabricated and mounted on a motorized stage. A droplet of particles, usually a concentrated monodispersed solution, is sandwiched between the template and a fixed glass slide. Then, the template is slowly moved relatively to the slide, drawing the droplet meniscus across the template. When the meniscus moves over the holes, nanoparticles at the template–meniscus interface are pushed into the holes and are assembled by using capillary and gravity forces. Based on this process, by choosing appropriate templates and nanoparticles, complex aggregates of nanoparticles with well-defined shapes, sizes, and structures can be readily obtained.

    Illustrations of the (a) NTSA [39] and (b) MTAS methods.

    Figure 1.Illustrations of the (a) NTSA [39] and (b) MTAS methods.

    Inspired by NTSA technology, we developed the MTAS method, illustrated in Fig. 1(b). This method requires only a few steps. First, a template with the desired geometry is attached to the target substrate. Then, the microspheres are placed on the template. A soft slicker is employed to blow the microspheres from left to right. The microspheres are pushed into the holes and assembled with external forces (Fe) and gravity (Fg). Here, Fe functions as the capillary force in the NTSA method. Note that when using the NTSA method, the nanoparticles can be in physical contact because of the attractive capillary force caused by solvent evaporation. To fix the microspheres, a thin adhesive layer is employed between the template and target substrate. Finally, the template is removed and a well-controlled microsphere-based all-dielectric metamaterial is obtained.

    In the fabrication, the stainless steel membrane template is fabricated by a solid-state laser (Enpon-Nano-H532, a nanosecond laser), which has a work wavelength, output power, and pulse duration of 532 nm, 7 W, and 10 ns, respectively. First, the stainless steel membrane with a thickness of 150 μm is placed on the workbench. Holes with various geometries are etched with proper AutoCAD files. The measured terahertz transmission spectra are obtained by using the photoconductive switch-based terahertz time-domain spectroscopy. The reflection of the fabricated sample is achieved by using terahertz spectroscopic/imaging system TAS7500SL, ADVANTEST.

    3. RESULTS AND DISCUSSION

    To demonstrate its capability and feasibility, we fabricated a single-layer ZrO2-based all-dielectric metamaterial on a flexible substrate. The micro-ZrO2 ceramic spheres were prepared according to our previous work [40]. Here, smooth, uniform, and spherical ZrO2 ceramic spheres with a diameter of d=150  μm were employed. In the experiment, a 50-μm-thick polyimide (PI) film was used as a substrate, because of its low permittivity (ε=3.0) and low loss (loss tangent=0.001) in the terahertz range. To obtain an adhesive layer, 2-μm-thick glue was cast on the PI film. A 150-μm-thick stainless steel membrane with micro through-holes (diameter=160  μm) was chosen as the template, and attached to the adhesive layer. The stainless steel membrane template was fabricated by laser etching, which ensures high precision and low cost. To assemble the microspheres on the substrate, they were placed on one side of the template. A slicker, used in screen printing, was employed to softly blow the microspheres from the left to the right. Using an appropriate strength, single-layer spheres were pushed into the holes and firmly attached to the substrate. Then, the template was carefully peeled off, and a single-layer all-dielectric metamaterial was obtained.

    In this method, no organic solvents, vacuum, or complex equipment are required, ensuring low-cost and environmentally friendly fabrication. Additionally, this method is fast and effective. The microscope image of the fabricated sample is shown in Fig. 2(a), where micro ZrO2 ceramic spheres are highly arrayed on the substrate. Thus, a single-layer all-dielectric metamaterial, based on micro ZrO2 ceramic spheres is obtained. Because the arrangement of the ceramic spheres is highly dependent on the stainless steel membrane template, one can prepare metamaterials with the desired structures using proper templates. Furthermore, this method is a generic way that can be extended to other terahertz all-dielectric metamaterials based on microspheres. To confirm these two advantages, an Al2O3 microsphere-based all-dielectric metamaterial was prepared with this method. Al2O3 microspheres with a diameter of 80 μm were employed, whereas the template had holes of diameter 85 μm and a lattice constant (P) of 200 μm. The fabricated metamaterial sample is shown in Fig. 2(b), where uniform Al2O3 microspheres are arrayed on the substrate with the desired arrangement. In this method, the fabricated area is also dependent on the area of the employed template. By using a large-scale template, one can obtain an ultra-large-scale metamaterial, which is highly desirable in practical applications. In the experiment, a 90-mesh sieve with a diameter of 40 cm was utilized. By using the MTAS method, we fabricated an ultra-large-scale and rolled metamaterial with an area exceeding 900  cm×900  cm [Figs. 2(c) and 2(d)]. This sample, to the best of our knowledge, is the largest all-dielectric metamaterial operating in the terahertz range. It should be noted that large-scale metallic metamaterials in the terahertz range have been reported [41,42]. However, large-scale all-dielectric metamaterials are rarely found.

    Microscope images of the fabricated (a) ZrO2 and (b) Al2O3 all-dielectric metamaterials. Photographs of (c) the fabrication process and (d) the fabricated ultra-large-scale flexible all-dielectric metamaterial using the MTAS method. (e) Simulated and measured transmissions for ZrO2 and Al2O3 all-dielectric metamaterials. The insets are simulated magnetic field intensity distributions at the corresponding resonance dips in the H–k plane.

    Figure 2.Microscope images of the fabricated (a) ZrO2 and (b) Al2O3 all-dielectric metamaterials. Photographs of (c) the fabrication process and (d) the fabricated ultra-large-scale flexible all-dielectric metamaterial using the MTAS method. (e) Simulated and measured transmissions for ZrO2 and Al2O3 all-dielectric metamaterials. The insets are simulated magnetic field intensity distributions at the corresponding resonance dips in the Hk plane.

    The resonances can be predicted by using Mie theory. Here, a nonmagnetic sphere with relative permittivity εs=ns2, relative permeability μh=1, and radius rs, in air or vacuum (nh=1), is applied. The first magnetic resonance (namely, the magnetic dipole) is considered. According to Wheeler et al. [43], the magnetic dipole portion of the scattered field is proportional to where rs=nsk0rs, rh=nhk0rs, ψ(x), and ξ(x) are the Riccati–Bessel functions. By equating the denominator of this equation to zero and employing the long-wavelength condition, one can obtain the first magnetic resonance (f), approximatively where c is the velocity of light in vacuum. The permittivity and loss of the ZrO2 ceramic spheres were 35.0 and 0.01, whereas those for Al2O3 were 10.0 and 0.01, respectively. By applying this equation to the ZrO2 and Al2O3 metamaterials, one can obtain the first magnetic resonances at frequencies of 0.337 and 1.181 THz, respectively.

    Using a terahertz time-domain spectroscopy system, the corresponding transmission spectra for the above samples were measured [Fig. 2(e)]. Regarding the ZrO2 ceramic sphere all-dielectric metamaterial, two remarkable resonance dips were observed at 0.366 and 0.509 THz. In the Al2O3 sphere all-dielectric metamaterial, the first two resonance dips were observed at 1.228 and 1.333 THz. Compared to the ZrO2 all-dielectric metamaterial, the resonances of the Al2O3 all-dielectric metamaterial shift to higher frequencies due to the lower permittivity and smaller size of Al2O3 spheres. It can be found that the theoretical predictions for the magnetic resonances are in good coincidence with experimental results. The deviation can be explained in two ways. First, in the theory prediction, the substrate has not been considered. Second, the permittivity in the experiment is frequency dependent. Simulations based on CST Microwave Studio (commercial software) were also conducted to obtain the transmission spectrum. In the simulation, unit cell boundary conditions are employed. Good agreement can be observed between the simulation and experimental results, confirming the high quality of the fabricated samples. It should be noted that the measured resonances have a lower Q-factor value, which is associated with the deviation of the sphere diameter and placing locations.

    To further demonstrate the robustness of this method, we assembled these microspheres into a variety of complex meta-aggregates, such as touching dimers, trimers, quadrumers, and chains, which are difficult or impossible to obtain with other approaches. These meta-aggregates are of great interest, because they exhibit many unique properties, which would find plenty of applications. For example, the all-dielectric dimer generates a remarkable electric-field enhancement between two resonators; this can be applied to high surface-enhanced fluorescence and surface-enhanced Raman scattering [44]. Meanwhile, the touching all-dielectric trimer has been demonstrated to support electromagnetically induced transparency [45], whereas the all-dielectric trimer can produce a new type of Fano resonance [46]. However, this meta-aggregate seems difficult to obtain by using conventional methods. In fact, this limitation can be addressed by adjusting the shape and dimensions of the templates. In the experiment, the touching dimer, trimer, quadrumer, and chain were obtained by using proper templates [Figs. 3(a)3(d)]. Here, the holes should be just larger than the meta-aggregates. This ensures that the meta-aggregates can pass through the holes, and thus prevent redundant spheres from passing through. The fabricated touching dimer, trimer, quadrumer, and chain are shown in Figs. 3(e)3(h), where one can see that the meta-aggregates are well-assembled and arrayed. We have also provided the simulated transmission of the dimer all-dielectric metamaterial (see Fig. 4). Here, the electric field direction of the wave is depicted by the inset. When the spacing between these two spheres is 20 μm, the first electric mode can be observed at 0.461 THz. In addition, at this resonance, remarkable electric-field enhancement between these two resonators is obtained (about 1000 times, as depicted by the inset). When the spacing between these two spheres is decreased to 0 μm, the first electric mode is split into two modes (w+ and w). For mode w+, great electric-field enhancement between these two resonators is obtained (about 3000 times). The structures with electric-field enhancement can be used in sensor and chemical identification.

    Illustrations of preparation for a touching (a) dimer, (b) trimer, (c) quadrumer, and (d) chain. The upper insets are templates and those below are the resultant metamaterials. (e)–(h) Corresponding fabricated samples.

    Figure 3.Illustrations of preparation for a touching (a) dimer, (b) trimer, (c) quadrumer, and (d) chain. The upper insets are templates and those below are the resultant metamaterials. (e)–(h) Corresponding fabricated samples.

    Simulated transmissions of dimer all-dielectric metamaterial with different spacings. The insets are the simulated electric field intensities.

    Figure 4.Simulated transmissions of dimer all-dielectric metamaterial with different spacings. The insets are the simulated electric field intensities.

    The introduction of an adhesion layer is of crucial importance for the successful fabrication of high-quality all-dielectric metamaterials. The adhesion layer ensures that the microspheres can be readily fabricated on other substrates. For example, all-dielectric metamaterials are fabricated on solid substrates, such as glass and silicon, by using a poly(methylmethacrylate) (PMMA) adhesion layer. In the experiment, PMMA was heated to 120°C, its glass-transition temperature. After the metamaterials were fabricated, PMMA and the substrates were cooled, and the metamaterials were firmly attached to the substrate. The stability of these metamaterials was tested by tape-peeling and ultrasound, both of which confirmed their high stability. More importantly, by using this method, the all-dielectric metamaterials could be easily transferred to the same active substrates, which provides more freedom for new applications. For example, the all-dielectric metamaterial can be prepared on an active material (heat-shrinkable) to obtain a broadband all-dielectric metamaterial reflector, which is difficult to achieve by using conventional methods. Here, we demonstrate that by using active material and the coupling effect, one can achieve a broadband all-dielectric metamaterial reflector.

    All-dielectric metamaterial reflectors are of interest because of their considerable potential applications (e.g., low-loss mirrors, surface protection, bio-imaging, and nano-sensing [1618]). The conventional metallic reflector is shown in Fig. 5(a), where a continuous metallic film is used to reflect an incident wave [it should be noted that the incident angle of a wave throughout this paper is 0° and Figs. 5(a) and 5(b) are just schematic illustrations]. The all-dielectric metamaterial reflector is made of separated dielectric resonators, as shown in Fig. 5(b). To obtain a perfect reflection, we first derive the corresponding conditions. For a semi-infinite medium, the reflectance R at normal incidence from vacuum can be expressed as [18] where z and z are the real and imaginary components of complex impedance (z=z+iz). Here, z=με, where ε and μ are the complex permittivity and the complex permeability, respectively. To obtain a perfect reflection (R=1), we need to create proper ε and μ. This equation can be solved when z=0, which leads to the following conditions: The effective permittivity ε and permeability μ should have opposite signs to meet these two equations. In such a reflector, the effective permittivity and permeability with opposite signs can be produced near the Mie resonances, leading to perfect reflection. Whereas almost all all-dielectric metamaterial perfect reflections have been investigated in the infrared and optical ranges [1618], they have rarely been reported in the terahertz range. In addition, whereas all-dielectric metamaterials can generate perfect reflection, they usually suffer from limited bandwidth. From the viewpoint of practical applications, broadening the working bandwidth becomes very important.

    Schematic illustrations of (a) a metallic reflector and (b) an all-dielectric metamaterial reflector. (c) Simulated reflection of the all-dielectric metamaterial with different lattice constants. Calculated effective permittivity and permeability with different lattice constants: (d) P=160 μm, (e) P=120 μm, and (f) P=90 μm. The regions where the effective permittivity and permeability have opposite signs are marked in blue.

    Figure 5.Schematic illustrations of (a) a metallic reflector and (b) an all-dielectric metamaterial reflector. (c) Simulated reflection of the all-dielectric metamaterial with different lattice constants. Calculated effective permittivity and permeability with different lattice constants: (d) P=160  μm, (e) P=120  μm, and (f) P=90  μm. The regions where the effective permittivity and permeability have opposite signs are marked in blue.

    To accomplish this goal, one can use the coupling effect between resonators. This has been investigated to broaden the bandwidth of near-zero permeability [47]. Simulations were carried out to investigate the influence of P on the reflection. Originally, the metamaterial was arranged in a hexagonal array and the geometrical parameters were P=160  μm and d=90  μm, where P is the lattice constant and d is the diameter of dielectric resonators. The thickness and permittivity of the substrate were 50 μm and 3.0, respectively. The simulated reflection is shown in Fig. 5(c), where two reflection peaks reaching nearly 100% can be observed at 0.67 and 0.89 THz. Simulated electric fields in the insets suggest that the first magnetic resonance and the first electric resonance were induced near those peaks. When P decreased from 160 to 90 μm, the magnetic resonance moved to a higher frequency, and the electric resonance moved to a slightly lower frequency. Additionally, the reflection valley between the two peaks vanished slowly with the reduction of P. Finally, when P=90  μm, about 0.17 THz (from 0.70 to 0.87 THz) bandwidth near-perfect reflection (>97%) could be obtained. Thus, a broadband THz all-dielectric metamaterial perfect reflector can be achieved when the single-layer sphere resonators are closely packed. To provide deeper insight into this phenomenon, S-parameter retrieval technology was employed to obtain the effective parameters based on the simulated complex transmission and reflection coefficients [Figs. 5(d)5(f)]. The Mie resonances produced regions where the effective permittivity and permeability had opposite signs (marked in blue). This property would result in a near-perfect reflection. It is also shown that these regions were broadened greatly, providing a broad bandwidth near-perfect reflection when P decreases. This phenomenon can be attributed to the coupling effect between resonances, which can make the magnetic plasma frequency shift to higher frequencies [47].

    As shown above, the key to a broadband perfect reflector is the preparation of single-layer close-packed sphere resonators. However, high-quality single-layer closely packed sphere resonators with a size of 100  μm cannot be easily fabricated with conventional methods. In our study, we demonstrate that this goal can be reached with a combination of MTAS and shrinkable film, as depicted in Fig. 6(a). The ZrO2 ceramic spheres were arrayed on a shrinkable film using the MTAS method (hexagonal lattice). Then, a shrinkable film with the spheres was heated, resulting in a metamaterial with closely packed spheres. The microscope images of the prepared sample before and after HT are shown in Figs. 6(b) and 6(c), respectively. Clearly, HT led to the reduction of P, and single-layer closely packed sphere resonators were obtained. Using the terahertz spectroscopic imaging system TAS7500SL, ADVANTEST, we measured the reflection of the prepared sample before and after HT [Fig. 6(d)]. As shown in this figure, before HT, two reflection peaks, reaching nearly 75%, are observed at 0.69 and 0.93 THz. However, the bandwidth of these reflection peaks is limited. After HT, a broadband (0.15  THz, from 0.74 to 0.89 THz) near-perfect reflection (>95%) can be observed. The experiment and simulation results showed good agreement, confirming the robustness of the proposed method. It should be noted that this method can also extend to other active composites [48,49], which would provide much more freedom for metamaterials.

    (a) Fabrication process of the broadband terahertz all-dielectric metamaterial reflector based on shrinkable film. Microscope images of the all-dielectric metamaterial (b) before and (c) after HT. (d) Measured reflections of the fabricated terahertz all-dielectric metamaterial before and after HT.

    Figure 6.(a) Fabrication process of the broadband terahertz all-dielectric metamaterial reflector based on shrinkable film. Microscope images of the all-dielectric metamaterial (b) before and (c) after HT. (d) Measured reflections of the fabricated terahertz all-dielectric metamaterial before and after HT.

    4. CONCLUSION

    A novel and generic approach called template-assisted fabrication, used to obtain ultra-large-scale, high-precision, and flexible ceramic microsphere all-dielectric metamaterials, is developed. Using this method, ceramic microsphere resonators (ZrO2 and Al2O3 resonators) could be readily assembled into nearly any arrangement, including complex aggregates, such as dimers, trimers, quadrumers, and chains, providing more freedom for manipulating terahertz radiation. The resultant all-dielectric metamaterials were facilely prepared on other active media. A broadband reflector, with a bandwidth of 0.15 THz and a reflection up to 95%, was obtained. This work provides a versatile and powerful platform for terahertz all-dielectric metamaterials, which should find a wide variety of applications in low-loss and high-efficiency terahertz devices.

    References

    [1] P. Moitra, Y. Yang, Z. Anderson, I. I. Kravchenko, D. P. Briggs, J. Valentine. Realization of an all-dielectric zero-index optical metamaterial. Nat. Photonics, 7, 791-795(2013).

    [2] J. C. Ginn, I. Brener, D. W. Peters, J. R. Wendt, J. O. Stevens, P. F. Hines, L. I. Basilio, L. K. Warne, J. F. Ihlefeld, P. G. Clem. Realizing optical magnetism from dielectric metamaterials. Phys. Rev. Lett., 108, 097402(2012).

    [3] A. S. Shorokhov, E. V. Melikgaykazyan, D. A. Smirnova, B. Hopkins, K. E. Chong, D. Y. Choi, M. R. Shcherbakov, A. E. Miroshnichenko, D. N. Neshev, A. A. Fedyanin. Multifold enhancement of third-harmonic generation in dielectric nanoparticles driven by magnetic Fano resonances. Nano Lett., 16, 4857-4861(2016).

    [4] P. Albella, R. A. de la Osa, F. Moreno, S. A. Maier. Electric and magnetic field enhancement with ultralow heat radiation dielectric nanoantennas: considerations for surface-enhanced spectroscopies. ACS Photon., 1, 524-529(2014).

    [5] M. R. Shcherbakov, D. N. Neshev, B. Hopkins, A. S. Shorokhov, I. Staude, E. V. Melikgaykazyan, M. Decker, A. A. Ezhov, A. E. Miroshnichenko, I. Brener. Enhanced third-harmonic generation in silicon nanoparticles driven by magnetic response. Nano Lett., 14, 6488-6492(2014).

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

    [7] Y. F. Yu, A. Y. Zhu, R. Paniagua-Domínguez, Y. H. Fu, B. Luk’Yanchuk, A. I. Kuznetsov. High-transmission dielectric metasurface with 2π phase control at visible wavelengths. Laser Photon. Rev., 9, 412-418(2015).

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

    [9] F. Aieta, M. A. Kats, P. Genevet, F. Capasso. Multiwavelength achromatic metasurfaces by dispersive phase compensation. Science, 347, 1342-1345(2015).

    [10] R. C. Devlin, M. Khorasaninejad, 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).

    [11] F. Capasso, F. Aieta, M. Khorasaninejad, P. Genevet, R. Devlin. Recent advances in planar optics: from plasmonic to dielectric metasurfaces. Optica, 4, 139-152(2017).

    [12] K. E. Chong, L. Wang, I. Staude, A. R. James, J. Dominguez, S. Liu, G. S. Subramania, M. Decker, D. N. Neshev, I. Brener. Efficient polarization-insensitive complex wavefront control using Huygens’ metasurfaces based on dielectric resonant meta-atoms. ACS Photon., 3, 514-519(2016).

    [13] Y. Yang, I. I. Kravchenko, D. Briggs, J. Valentine. Dielectric metasurface analogue of electromagnetically induced transparency. Nat. Commun., 5, 5753(2014).

    [14] E. Semouchkina, R. Duan, G. Semouchkin, R. Pandey. Sensing based on Fano-type resonance response of all-dielectric metamaterials. Sensors, 15, 9344-9359(2015).

    [15] Z. Huang, J. Wang, Z. Liu, G. Xu, Y. Fan, H. Zhong, B. Cao, C. Wang, K. Xu. Strong-field-enhanced spectroscopy in silicon nanoparticle electric and magnetic dipole resonance near a metal surface. J. Phys. Chem. C, 20, 47-50(2015).

    [16] P. Moitra, B. A. Slovick, W. Li, I. I. Kravchencko, D. P. Briggs, S. Krishnamurthy, J. Valentine. Large-scale all-dielectric metamaterial perfect reflectors. ACS Photon., 2, 692-698(2015).

    [17] P. Moitra, B. A. Slovick, Z. Gang Yu, S. Krishnamurthy. Experimental demonstration of a broadband all-dielectric metamaterial perfect reflector. Appl. Phys. Lett., 104, 171102(2014).

    [18] S. Krishnamurthy, B. Slovick, Z. G. Yu, M. Berding. Perfect dielectric-metamaterial reflector. Phys. Rev. B, 88, 5514-5518(2013).

    [19] C. Y. Yang, J. H. Yang, Z. Y. Yang, Z. X. Zhou, M. G. Sun, V. E. Babicheva, K. P. Chen. Nonradiating silicon nanoantenna metasurfaces as narrow-band absorbers. ACS Photon., 5, 2596-2601(2018).

    [20] H. T. Chen, W. J. Padilla, M. J. Cich, A. K. Azad, R. D. Averitt, A. J. Taylor. A metamaterial solid-state terahertz phase modulator. Nat. Photonics, 3, 148-151(2009).

    [21] M. Choi, S. H. Lee, Y. Kim, S. B. SKang, J. Shim, M. H. Kwak, K. Y. Kang, Y. H. Lee, N. Park, B. Min. A terahertz metamaterial with unnaturally high refractive index. Nature, 470, 369-373(2011).

    [22] H. T. Chen, W. J. Padilla, J. M. O. Zide, A. C. Gossard, A. J. Taylor, R. D. Averitt. Active terahertz metamaterial devices. Nature, 444, 597-600(2006).

    [23] S. Y. Chiam, R. Singh, C. Rockstuhl, F. Lederer, W. Zhang, A. A. Bettiol. Analogue of electromagnetically induced transparency in a terahertz metamaterial. Phys. Rev. B, 80, 153103(2009).

    [24] R. Singh, I. A. I. Al-Naib, M. Koch, W. Zhang. Sharp Fano resonances in THz metamaterials. Opt. Express, 19, 6312-6319(2011).

    [25] T. Hu, C. M. Bingham, A. C. Strikwerda, D. Pilon, D. Shrekenhamer, N. I. Landy, K. Fan, X. Zhang, W. J. Padilla, R. D. Averitt. Highly-flexible wide angle of incidence terahertz metamaterial absorber. Phys. Rev. B, 78, 1879-1882(2008).

    [26] L. H. Gao, Q. Cheng, J. Yang, S. J. Ma, J. Zhao, S. Liu, H. B. Chen, Q. He, W. X. Jiang, H. F. Ma. Broadband diffusion of terahertz waves by multi-bit coding metasurfaces. Light Sci. Appl., 4, e324(2015).

    [27] L. Liu, X. Zhang, K. Mitchell, X. Su, N. Xu, C. Ouyang, Y. Shi, J. Han, W. Zhang, S. Zhang. Broadband metasurfaces with simultaneous control of phase and amplitude. Adv. Mater., 26, 5031-5036(2014).

    [28] X. Zhang, Z. Tian, W. Yu, J. Gu, S. Zhang, J. Han, W. Zhang. Broadband terahertz wave deflection based on C-shape complex metamaterials with phase discontinuities. Adv. Mater., 25, 4567-4572(2013).

    [29] H. Xu, K. Bi, Y. Hao, J. Zhang, J. Xu, J. Dai, K. Xu, J. Zhou. Switchable complementary diamond-ring-shaped metasurface for radome application. IEEE Antennas Wireless Propag. Lett., 17, 2494-2497(2018).

    [30] X. Wang, Y. Cui, T. Li, M. Lei, J. Li, Z. Wei. Recent advances in the functional 2D photonic and optoelectronic devices. Adv. Opt. Mater., 7, 1801274(2018).

    [31] Q. Wang, X. Li, L. Wu, P. Lu, Z. Di. Electronic and interface properties in graphene oxide/hydrogen-passivated Ge heterostructure. Phys. Status Solidi (RRL), 13, 1800461(2019).

    [32] Y. Yang, B. Cui, Z. Geng, S. Feng. Terahertz magnetic and electric Mie resonances of an all-dielectric one-dimensional grating. Appl. Phys. Lett., 106, 111106(2015).

    [33] H. Němec, P. Kužel, F. Kadlec, C. Kadlec, R. Yahiaoui, P. Mounaix. Tunable terahertz metamaterials with negative permeability. Phys. Rev. B, 79, 241108(2009).

    [34] I. V. Shadrivov, K. Fan, W. J. Padilla, X. Liu. Experimental realization of a terahertz all-dielectric metasurface absorber. Opt. Express, 25, 191-201(2017).

    [35] C. Ouyang, C. Hu, H. Zhang, J. Han, J. Gu, M. Wei, Q. Wang, Q. Xu, W. Zhang, X. Zhang. Polarization-independent all-silicon dielectric metasurfaces in the terahertz regime. Photon. Res., 6, 24-29(2017).

    [36] Z. Ma, S. M. Hanham, P. Albella, B. Ng, H. T. Lu, Y. Gong, S. A. Maier, M. Hong. Terahertz all-dielectric magnetic mirror metasurfaces. ACS Photon., 3, 1010-1018(2016).

    [37] D. Headland, S. Nirantar, W. Withayachumnankul, P. Gutruf, D. Abbott, M. Bhaskaran, C. Fumeaux, S. Sriram. Terahertz magnetic mirror realized with dielectric resonator antennas. Adv. Mater., 27, 7137-7144(2015).

    [38] H. Němec, C. Kadlec, F. Kadlec, P. Kužel, R. Yahiaoui, U. C. Chung, C. Elissalde, M. Maglione, P. Mounaix. Resonant magnetic response of TiO2 microspheres at terahertz frequencies. Appl. Phys. Lett., 100, 061117(2012).

    [39] J. A. Fan, K. Bao, L. Sun, J. Bao, V. N. Manoharan, P. Nordlander, F. Capasso. Plasmonic mode engineering with templated self-assembled nanoclusters. Nano Lett., 12, 5318-5324(2012).

    [40] J. Gao, C. Lan, Q. Zhao, B. Li, J. Zhou. Experimental realization of Mie-resonance terahertz absorber by self-assembly method. Opt. Express, 26, 13001-13011(2018).

    [41] M. Walther, A. Ortner, H. Meier, U. Loffelmann, P. J. Smith, J. G. Korvink. Terahertz metamaterials fabricated by inkjet printing. Appl. Phys. Lett., 95, 251107(2009).

    [42] K. Takano, T. Kawabata, C. F. Hsieh, K. Akiyama, F. Miyamaru, Y. Abe, Y. Tokuda, R. P. Pan, C. L. Pan, M. Hangyo. Fabrication of terahertz planar metamaterials using a super-fine ink-jet printer. Appl. Phys. Express, 3, 016701(2010).

    [43] M. S. Wheeler, J. S. Aitchison, M. Mojahedi. Coupled magnetic dipole resonances in sub-wavelength dielectric particle clusters. J. Opt. Soc. Am. B, 27, 1083-1091(2010).

    [44] M. Caldarola, P. Albella, E. Cortés, M. Rahmani, T. Roschuk, G. Grinblat, R. F. Oulton, A. V. Bragas, S. A. Maier. Non-plasmonic nanoantennas for surface enhanced spectroscopies with ultra-low heat conversion. Nat. Commun., 6, 7915(2015).

    [45] J. H. Yan, P. Liu, Z. Y. Lin, H. Wang, H. J. Chen, C. X. Wang, G. W. Yang. Magnetically induced forward scattering at visible wavelengths in silicon nanosphere oligomers. Nat. Commun., 6, 7042(2015).

    [46] B. Hopkins, D. S. Filonov, A. E. Miroshnichenko, F. Monticone, A. Alu, Y. S. Kivshar. Interplay of magnetic responses in all-dielectric oligomers to realize magnetic Fano resonances. ACS Photon., 2, 724-729(2015).

    [47] J. Sun, L. Kang, R. Wang, L. Liu, L. Sun, J. Zhou. Low loss negative refraction metamaterial using a close arrangement of split-ring resonator arrays. New J. Phys., 12, 083020(2010).

    [48] G. Zhang, Q. Liao, Z. Zhang, Q. Liang, Y. Zhao, X. Zheng, Y. Zhang. Novel piezoelectric paper-based flexible nanogenerators composed of BaTiO3 nanoparticles and bacterial cellulose. Adv. Sci., 3, 1500257(2016).

    [49] K. Bi, M. Bi, Y. Hao, W. Luo, W. Cai, X. Wang, Y. Huang. Ultrafine core-shell BaTiO3@SiO2 structures for nanocomposite capacitors with high energy density. Nano Energy, 51, 513-523(2018).

    Ke Bi, Daquan Yang, Jia Chen, Qingmin Wang, Hongya Wu, Chuwen Lan, Yuping Yang, "Experimental demonstration of ultra-large-scale terahertz all-dielectric metamaterials," Photonics Res. 7, 457 (2019)
    Download Citation