Metasurface engineering through bound states in the continuum

Metasurface engineering through bound states in the continuum

Anton S. Kupriianov    Yi Xu    Andrey Sayanskiy    Victor Dmitriev    Yuri S. Kivshar    Vladimir R. Tuz; College of Physics, Jilin University, Changchun 130012, China Department of Electronic Engineering, College of Information Science and Technology, Jinan University, Guangzhou 510632, China ITMO University, St. Petersburg 197101, Russia Electrical Engineering Department, Federal University of Para, PO Box 8619, Agencia UFPA, CEP 66075-900 Belem, Para, Brazil Nonlinear Physics Center, Australian National University, Canberra ACT 2601, Australia International Center of Future Science, State Key Laboratory of Integrated Optoelectronics, College of Electronic Science and Engineering, Jilin University, Changchun 130012, China Institute of Radio Astronomy of the National Academy of Science, Kharkiv 61002, Ukraine

Metasurfaces have attracted a lot of attention in recent years due to novel ways they provide for the efficient wavefront control and engineering of the resonant transmission. We introduce a new concept for effective engineering of the sharp Fano resonances in metasurfaces associated with a new type of bound states in the continuum. We demonstrate that by breaking the symmetry transversely, in the direction perpendicular to a metasurface with a complex unit cell, we can control the number, frequency, and type of high- resonances originating from bound states in the continuum. As example, we demonstrate experimentally the metasurfaces with magnetic dipole and toroidal dipole responses governed by the physics of multipolar bound states.

Bound states in the continuum (BICs) are waves that remain localized in dynamic systems even though they coexist with a continuous spectrum of radiating waves. Although BICs were first introduced in quantum mechanics Von Neumann and Wigner (1929); Stillinger and Herrick (1975), it is a general wave phenomenon that has been found in many physical systems, such as electromagnetic, acoustic, water and elastic waves. Specifically, these states have also been identified in a wide range of optical systems, including dielectric gratings Marinica et al. (2008), optical waveguides Plotnik et al. (2011); Weimann et al. (2013), photonic crystals Lee et al. (2012); Hsu et al. (2013), graphene quantum-dot structures González et al. (2010), and hybrid plasmonic-photonic systems Azzam et al. (2018) (see a comprehensive review on BICs in Ref. Hsu et al., 2016 and references therein).

One interesting type of BICs relies on the symmetry properties of both discrete mode and coexisted radiative continuum. When a discrete mode is characterized by a particular symmetry class, radiative continuum belonging to a different symmetry classes completely decouples from such a discrete state Lee et al. (2012); He et al. (2018). The coupling between the bound states and continuum band is forbidden as long as their symmetry properties are preserved, and they are classified as symmetry-protected BICs Lee et al. (2012).

Figure 1: (a) A new concept to manipulate the symmetry breaking in all-dielectric metasurfaces composed of cluster-based unit cells. Two typical examples are (b) dimer and (c) trimer clusters of dissimilar subwavelength dielectric disks.

A particular signature of symmetry-protected BICs in the field of optical metasurfaces is the observation of exceptionally narrow resonances (sometimes called trapped modes) in the optical spectra of metasurfaces composed of in-plane asymmetric structural elements Prosvirnin and Zouhdi (2003); Fedotov et al. (2007). The appearance of ultrahigh- resonances are attributed to the excitation of asymmetric modes whose electromagnetic field distributions are slightly deviated from the, otherwise inaccessible, symmetric modes. The corresponding quality factors of leaky resonances depend on the degree of the introduced structural asymmetry Koshelev et al. (2018). The concept of manipulating the transition between the symmetry-protected BICs and leaky resonances was first proposed for plasmonic metasurfaces composed of metallic split-ring resonators Khardikov et al. (2010), and then it was developed for all-dielectric metasurfaces Khardikov et al. (2012); Zhang et al. (2013); Jain et al. (2015); Khardikov et al. (2016); Kruk and Kivshar (2017); Tuz et al. (2018a, b); Sayanskiy et al. (2019); Yu et al. (2019).

In general, BIC modes transform into leaky resonances when the symmetry of either eigenmodes of the structure or radiative continuum is broken. As a result, one needs to break the symmetry of either discrete mode or radiative continuum in order to access such symmetry protected multipolar BICs. To date, all the mechanisms reported so far for manipulating the transition between BICs and radiative continuum can be classified as in-plane symmetry breaking. No effort has been paid to manipulate the transition between BICs and leaky resonances utilizing the out-of-plane symmetry, which will provide more flexible approach to engineer and control the resonant properties of metasurfaces.

In this Letter, we introduce a new concept for effectively controlling the high- Fano resonances, which are originated from the BICs of optical metasurfaces, by utilizing out-of-plane symmetry breaking in the unit cell of an all-dielectric metasurface. We demonstrate that our idea is quite general, which has been validated in metasurfaces composed of lattices of particle clusters corresponding to different symmetries. By introducing this unique degree of freedom beyond the in-plane symmetry, one can efficiently control the symmetry of the discrete mode and thus manipulate its interaction with the continuum. Employing this concept, we realize accesses to various kinds of multipolar BIC associated leaky modes existed in the cluster-based all-dielectric metasurfaces. We further perform proof-of-principle microwave experiments by demonstrating leaky resonances associated with the toroidal dipole BICs in both far-field and near-field measurements.

Figure 2: (a) A set of the BIC-type eigenstates and (b) the corresponding electromagnetic near-field distribution for metasurfaces with the unit cell composed of a pair of disks (dimer). (c) Access to leaky resonance originated from the BIC Mode D2 of the dimer metasurface via out-of-plane symmetry breaking. (d) Dependence of the quality factors of leaky resonances originated from BIC modes D1-D4 on the out-of-plane asymmetry parameter . Parameters of the metasurface are:  mm,  mm,  mm,  mm, , .

Without loss of generality, we consider an all-dielectric metasurface whose unit cell consists of different types of particle clusters [Fig. 1(a)]. The simplest configuration of the unit cell is a disk dimer arranged in a square lattice, as shown in [Fig. 1(b)]. Therefore, we start from the demonstration of our concept by evaluating how the out-of-plane symmetry breaking will affect on the transition between various BIC modes and the corresponding leaky resonances in a dimer-based metasurface. Then it is generalized to a more complex trimer-based metasurface composed of equilateral triangular unit cells [Fig. 1(c)].

All notations and parameters related to the geometry of the metasurface under study are summarized in Fig. 1. The disks are made of a non-magnetic material with relative permittivity and the metasurface is placed on a dielectric substrate with relative permittivity . We consider that the metasurface is illuminated by a normally incident () linearly polarized plane wave with the electric field vector directed either along the axis (, -polarized wave) or the axis (, -polarized wave).

Figure 3: (a) A set of BIC-type eigenstates and (b) their electromagnetic near-field distributions for metasurfaces with the twin-trimer unit cell. (c) Schemes of access to Mode T1 in the trimer-based metasurface. (d) Evolution of the transmission spectra vs the asymmetry parameter . (e) Electromagnetic near-field distribution of the leaky mode at the resonant frequency corresponding to the manifestation of Mode T1.  mm, all other parameters are the same as in Fig. 2.

For the metasurface under study we can recognize the BIC modes by performing an eigenvalue analysis using the RF module of commercial COMSOL Multiphysics® finite-element electromagnetic solver com (). Such BIC modes can be distinguished by their real-valued eigenfrequencies. In accordance with our experimental possibilities, we consider the BICs mode in the microwave part of the spectrum ( GHz). From the known set of modes we selected several BIC modes which appear in the chosen frequency band for the given parameters of the structure. We denoted the modes of interest by symbols D1-D4 [Fig. 2(a)].

These modes manifest themselves by distinct displacement current geometries, which can be regarded as specific combinations of magnetic dipole (MD) and electric dipole (ED) moments inherent to an individual disk [Fig. 2(b)]. In particular, for Mode D2 the magnetic field vector is anti-parallel located in the - plane, whereas displacement current flows inside the particles in the - plane and has a solenoid-like behavior enclosing the magnetic field. It resembles the typical out-of-plane toroidal dipole (TD) BIC mode which was theoretically studied recently He et al. (2018). Rather than enabling the coupling between TD BIC mode and the obliquely incident plane wave, we show that this TD BIC mode can be turned to a TD leaky resonance by breaking the out-of-plane symmetry, as shown in Fig. 2(c). Specifically, such a coupling to Mode D2 can be realized for the -polarized incident wave as soon as a difference in the thicknesses of disks forming a dimer is introduced [we denote this difference as the asymmetry parameter , see Fig. 1(b)]. Notice that the group of symmetry in the - plane of the non-perturbed dimer with is , and of the perturbed dimer with is .

Figure 4: Experimental results. (a) Transmission coefficient of the trimer-based metasurface. The inset presents a fragment of the metasurface prototype and several sets of particles. In the simulation actual material losses () in ceramic particles are taken into account, while the substrate is modeled as a loss-less dielectric. The asymmetry factor is  mm. All other material and geometrical parameters of the metasurface are the same as in Fig. 3. (b) Simulated and (c) measured near-field distribution of the real part of the normal component of the electric field () at the resonant frequency . The measurement plane is 1 mm away from the metasurface. The field is normalized to the maximum value of .

Indeed, while the symmetry of dimer is preserved, there is no peculiarities in the transmitted spectra of the metasurface under study at the frequencies of the BIC modes [the eigenfrequency of Mode D2 is denoted as , see Fig. 2(c)]. As soon as the asymmetry is introduced, there is a leaky resonance (the resonant frequency of the leaky mode originated from Mode D2 is denoted as ) in the transmission spectrum and such resonance acquires a peak-and-trough (Fano) profile which is typical for the metasurfaces composed of in-plane asymmetric structural elements Fedotov et al. (2007). By controlling the asymmetry factor, we can manipulate the -factors of the leaky resonances. The less out-of-plane asymmetry is, the higher the quality factor of the leaky mode resonance is [Fig. 2(d)]. At the same time, the corresponding frequencies and related to Mode D2 are found to be different but they are rather closely spaced to each other. The difference between them also depends on the degree of the out-of-plane symmetry breaking. Therefore, breaking the out-of-plane symmetry of the unit cell in the metasurface presents an efficient way to manipulate the quality factor of resonances related to the BIC modes.

It should be pointed out that breaking the out-of-plane symmetry of the cluster composed of disks has fundamentally different effects on the quality factor of the BIC modes compared with the case of breaking the in-plane symmetry. One can anticipate that introduced asymmetry in the thicknesses of disks forming the dimer has a stronger effect on its out-of-plane ED moment and in-plane MD moment, since the polarization currents are parallel to such plane. In contrast, breaking the in-plane symmetry by changing the radius of one disk has more pronounced effect on the in-plane ED moment and out-of-plane MD moment. For instance, among the selected BIC modes of dimer, Mode D2 is considered as a TD BIC composed of two anti-parallel in-plane MD moments. Therefore, symmetry breaking in the out-of-plane direction turns such a TD-BIC into a leaky resonance much more efficiently compared with other types of BIC modes. At the same time, the decrease of -factor for Mode D4 consisted in the in-plane MD moments is also faster than that for Modes D1 and D3, as shown in Fig. 2(d). As a result, the out-of-plane symmetry breaking can be regarded as a new degree of freedom and an efficient way to tailor the resonant properties of the symmetry protected BIC modes, especially for the out-of-plane TD BIC and in-plane magnetic BIC modes.

We further present an example utilizing the out-of-plane symmetry breaking induced transition between the symmetry protected BIC modes and the corresponding leaky resonances in an all-dielectric metasurface based on trimers [Fig. 1(c)]. The existence of the TD mode in the metasurface whose lattice is formed by trimers arranged in square unit cells was examined in Xu et al. (2019). Here we introduce some modifications to the structure of unit cell studied earlier to provide further degree of freedom for resonance manipulation in both near- and far-fields. Firstly, we consider that the unit cell has the form of an equilateral triangle. Secondly, two triangle unit cells are paired together in such a way that a rhomboid cluster appears, i.e., generally, the metasurface under study is constructed from the twin-trimer super-cells.

After the eigenvalue analysis of this metasurface, a set of BIC modes existed in the frequency band of interest are selected. We denote these modes by symbols T1-T3 [Fig. 3(a)]. All these modes present fruitful electromagnetic near-field distributions manifested as electric and magnetic hot-spots in the near-field [Fig. 3(b)]. Specifically, these BIC modes are composed of various MD moments in each disk, which are arranged in the - plane. When all disks are identical, the symmetry of the twin-trimer super-cell is and it belongs to the rhomboid system (see Fig. 3(c), upper part). As can be seen from this figure, neither the -polarized nor the -polarized normally incident plane wave can access to these BIC modes. As the in-plane MD moment is more sensitive to the out-of-plane asymmetry for the all-dielectric metasurface under consideration, these BIC modes can be effectively turned into leaky resonance by increasing the thickness of certain disks. For example, perturbation of the unit cell by two disks with different thicknesses reduces the symmetry to of the monoclinic system (Fig. 3(c), middle part). This perturbation allows one to excite the MD leaky resonance originated from BIC Mode T1. It should be pointed out that the metasurface has a polarization-dependent response, where the transmissions of the metasurface are different under the excitation by the -polarized and -polarized plane waves [Fig. 3(c)].

At the same time, the out-of-plane symmetry breaking can also be used to tailor the frequency and -factor of the leaky Fano resonance. By changing the asymmetry factor from to  mm, one can observe the shift of the resonant frequency of the leaky mode compared with the eigenfrequency of Mode T1, as shown by Figs. 3(d). It can be observed that the -factor of this resonance can also be manipulated accordingly by using different degrees of out-of-plane symmetry breaking. More importantly, the electromagnetic near-field can also be manipulated by using the asymmetry parameter . The electromagnetic field distributions at the resonant frequency under the excitation of the -polarized plane wave are shown in Figs. 3(e) for two typical asymmetry parameters and , respectively. As can be seen from this figure, they are distinct by their electromagnetic near-field distributions though they are originated from the same BIC Mode T1, validating the concept of out-of-plane symmetry breaking for the manipulation of the vector near field. It should be pointed out that our concept is quite general and it can be easily generalized to resonant metasurfaces operating at the optical spectrum. (see two examples in Supplemental Material Sup ()). The ability to control the electromagnetic near-field by utilizing out-of-plane symmetry breaking provides possibility to precisely tailor both the emission of the ED and MD types of dipole emitters Feng et al. (2018).

To experimentally verify the out-of-plane symmetry breaking effect as an effective mean to trigger the transition between BIC modes and high- Fano resonance, an all-dielectric metasurface is fabricated and studied in the microwave range. All details regarding to the samples preparation and the measurement setup can be found in Refs. Sayanskiy et al. (2019) and Xu et al. (2019) (see also Supplemental Material Sup ()). The manifestation of the corresponding leaky resonance is detected in the far-field measurement of the transmission coefficient of the metasurface [Fig. 4(a)]. In order to provide a reference for experiment, we present the simulated electric near-field distribution [Real] of the metasurface at the frequency . As can be seen from Fig. 4(b), the electric field is mostly concentrated in the center of trimers, being out of phase in the pair of trimers forming the super-cell. The electromagnetic near field distribution are then experimentally measured for the all-dielectric metasurface at the same resonant frequency . Indeed, the simulated [Fig. 4(b)] and experimental [Fig. 4(c)] results show excellent agreement with each other, validating the excitation of MD leaky resonance originated from BIC Mode T1.

In summary, we have introduced a new mechanism for a highly efficient engineering and control of resonant metasurfaces via symmetry-breaking in the out-of-plane direction, which essentially trigger the transition between symmetry protected BIC modes and leaky resonances. We show that both the resonance frequency and -factor of leaky resonance can be manipulated by the asymmetry factor . At the same time, the electromagnetic near field can also be tailored simultaneously to address the electric and magnetic hot-spots locally. As an example, we have presented a microwave experiment to validate the concept of out-of-plane symmetry breaking for metasurface engineering.

A.S. and V.R.T. acknowledge the Jilin University for hospitality and financial support; V.D. thanks the Brazilian agency CNPq for a financial support; Y.S.K. acknowledges a support from the Strategic Fund of the Australian National University.


  • Von Neumann and Wigner (1929) J. Von Neumann and E. Wigner, “Über merkwürdige diskrete eigenwerte (about strange discrete eigenvalues),” Phys. Zschr. 30, 465–467 (1929).
  • Stillinger and Herrick (1975) F. H. Stillinger and D. R. Herrick, “Bound states in the continuum,” Phys. Rev. A 11, 446–454 (1975).
  • Marinica et al. (2008) D. C. Marinica, A. G. Borisov,  and S. V. Shabanov, “Bound states in the continuum in photonics,” Phys. Rev. Lett. 100, 183902 (2008).
  • Plotnik et al. (2011) Y. Plotnik, O. Peleg, F. Dreisow, M. Heinrich, S. Nolte, A. Szameit,  and M. Segev, “Experimental observation of optical bound states in the continuum,” Phys. Rev. Lett. 107, 183901 (2011).
  • Weimann et al. (2013) S. Weimann, Y. Xu, R. Keil, A. E. Miroshnichenko, A. Tünnermann, S. Nolte, A. A. Sukhorukov, A. Szameit,  and Y. S. Kivshar, “Compact surface Fano states embedded in the continuum of waveguide arrays,” Phys. Rev. Lett. 111, 240403 (2013).
  • Lee et al. (2012) J. Lee, B. Zhen, S.-L. Chua, W. J. Qiu, J. D. Joannopoulos, M. Soljačić,  and O. Shapira, “Observation and differentiation of unique high- optical resonances near zero wave vector in macroscopic photonic crystal slabs,” Phys. Rev. Lett. 109, 067401 (2012).
  • Hsu et al. (2013) C. W. Hsu, B. Zhen, J. Lee, S.-L. Chua, S. G Johnson, J. D Joannopoulos,  and M. Soljačić, “Observation of trapped light within the radiation continuum,” Nature 499, 188 (2013).
  • González et al. (2010) J. W. González, M. Pacheco, L. Rosales,  and P. A. Orellana, “Bound states in the continuum in graphene quantum dot structures,” EPL 91, 66001 (2010).
  • Azzam et al. (2018) S. I. Azzam, V. M. Shalaev, A. Boltasseva,  and A. V. Kildishev, “Formation of bound states in the continuum in hybrid plasmonic-photonic systems,” Phys. Rev. Lett. 121, 253901 (2018).
  • Hsu et al. (2016) C. W. Hsu, B. Zhen, A. D. Stone, J. D. Joannopoulos,  and M. Soljačić, “Bound states in the continuum,” Nat. Rev. Mater. 1, 16048 (2016).
  • He et al. (2018) Y. He, G. Guo, T. Feng, Y. Xu,  and A. E. Miroshnichenko, “Toroidal dipole bound states in the continuum,” Phys. Rev. B 98, 161112 (2018).
  • Prosvirnin and Zouhdi (2003) S. Prosvirnin and S. Zouhdi, “Resonances of closed modes in thin arrays of complex particles,” in Advances in Electromagnetics of Complex Media and Metamaterials, edited by S. Zouhdi and M. Arsalane (Kluwer Academic Publishers, Dordrecht, 2003) pp. 281–290.
  • Fedotov et al. (2007) V. A. Fedotov, M. Rose, S. L. Prosvirnin, N. Papasimakis,  and N. I. Zheludev, “Sharp trapped-mode resonances in planar metamaterials with a broken structural symmetry,” Phys. Rev. Lett. 99, 147401 (2007).
  • Koshelev et al. (2018) K. Koshelev, S. Lepeshov, M. Liu, A. Bogdanov,  and Y. Kivshar, “Asymmetric metasurfaces with high- resonances governed by bound states in the continuum,” Phys. Rev. Lett. 121, 193903 (2018).
  • Khardikov et al. (2010) V. V. Khardikov, E. O. Iarko,  and S. L. Prosvirnin, “Trapping of light by metal arrays,” J. Opt. 12, 045102 (2010).
  • Khardikov et al. (2012) V. V. Khardikov, E. O. Iarko,  and S. L. Prosvirnin, “A giant red shift and enhancement of the light confinement in a planar array of dielectric bars,” J. Opt. 14, 035103 (2012).
  • Zhang et al. (2013) J. Zhang, K. F. MacDonald,  and N. I. Zheludev, “Near-infrared trapped mode magnetic resonance in an all-dielectric metamaterial,” Opt. Express 21, 26721–26728 (2013).
  • Jain et al. (2015) A. Jain, P. Moitra, T. Koschny, J. Valentine,  and C. M. Soukoulis, “Electric and magnetic response in dielectric dark states for low loss subwavelength optical meta atoms,” Adv. Opt. Mater. 3, 1431–1438 (2015).
  • Khardikov et al. (2016) V. Khardikov, P. Mladyonov, S. Prosvirnin,  and V. Tuz, “Electromagnetic wave diffraction by periodic planar metamaterials with nonlinear constituents,” in Contemporary Optoelectronics: Materials, Metamaterials and Device Applications, edited by O. Shulika and I. Sukhoivanov (Springer Netherlands, Dordrecht, 2016) Chap. 5, pp. 81–98.
  • Kruk and Kivshar (2017) S. Kruk and Y. Kivshar, “Functional meta-optics and nanophotonics governed by Mie resonances,” ACS Photonics 4, 2638–2649 (2017).
  • Tuz et al. (2018a) V. R. Tuz, V. V. Khardikov,  and Y. S. Kivshar, “All-dielectric resonant metasurfaces with a strong toroidal response,” ACS Photonics 5, 1871–1876 (2018a).
  • Tuz et al. (2018b) V. R. Tuz, V. V. Khardikov, A. S. Kupriianov, K. L. Domina, S. Xu, H. Wang,  and H.-B. Sun, “High-quality trapped modes in all-dielectric metamaterials,” Opt. Express 26, 2905–2916 (2018b).
  • Sayanskiy et al. (2019) A. Sayanskiy, A. S. Kupriianov, S. Xu, P. Kapitanova, V. Dmitriev, V. V. Khardikov,  and V. R. Tuz, “Controlling high- trapped modes in polarization-insensitive all-dielectric metasurfaces,” Phys. Rev. B 99, 085306 (2019).
  • Yu et al. (2019) P. Yu, A. Kupriianov, V. Dmitriev,  and V. Tuz, “All-dielectric metasurfaces with trapped modes: Group-theoretical description,” J. Appl. Phys. 125, 143101 (2019).
  • (25) Comsol Application Gallery ID: 15711.
  • Xu et al. (2019) S. Xu, A. Sayanskiy, A. S. Kupriianov, V. R. Tuz, P. Kapitanova, H.-B. Sun, W. Han,  and Y. S. Kivshar, “Experimental observation of toroidal dipole modes in all-dielectric metasurfaces,” Adv. Opt. Mater. 7, 1801166 (2019).
  • (27) See Supplemental Material at (url) for information on samples preparation and characterization as well as description of metasurfaces targeted on optical range (wavelength of nm).
  • Feng et al. (2018) T. Feng, W. Zhang, Z. Liang,  and Y. Xu, “Unidirectional emission in an all-dielectric nanoantenna,” J. Phys. Condens. Matter 30, 124002 (2018).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description