Multiple symmetry protected BIC lines in two dimensional synthetic parameter space

Fengyuan Zhang; Qiongqiong Chu; Qiang Wang; Shining Zhu; Hui Liu

doi:10.1515/nanoph-2022-0781

Article Open Access

Multiple symmetry protected BIC lines in two dimensional synthetic parameter space

Fengyuan Zhang , Qiongqiong Chu , Qiang Wang , Shining Zhu and Hui Liu

Published/Copyright: February 6, 2023

Published by

Become an author with De Gruyter Brill

Submit Manuscript Author Information

From the journal Nanophotonics Volume 12 Issue 13

Abstract

Bound states in the continuum (BICs) have attracted significant interest in recent years due to their unique optical properties, such as infinite quality factor and wave localization. In order to improve the optical performance of BICs based devices, more degrees of freedom are required to tune BICs in high-dimension parameter space for practical applications. To effectively tune more BICs, we form a 2D synthetic parameter space based on a nanohole metasurface array. Multiple symmetry protected BIC modes with high Q factors can be achieved at high-order symmetry point. Through manipulating asymmetry parameters, BIC lines formed by a series of BIC modes can be found in the 2D synthetic parameter space. Moreover, the electric field distributions are investigated to demonstrate the generation and evolution of BICs. By measuring the absorption spectra, the tuning of multiple BICs with synthetic asymmetry parameters is experimentally explored, which agrees well with theoretical results. Therefore, our design can provide new insight for a variety of on-chip applications, such as nonlinear devices, integrated nanolasing array, and high-resolution sensors for infrared molecular detection.

Keywords: BIC lines; bound states in the continuum; nanohole metasurface; synthetic parameter space

1 Introduction

BICs are eigenmodes embedded in the continuum radiation spectrum, generally characterized by infinite Q factor and perfect optical mode confinement [1–3]. This intriguing phenomenon was first proposed in quantum mechanics by von Neumann and Wigner [1], and has been found in many physical systems, such as electromagnetic, acoustic, and water waves [4–6]. When small perturbations (such as structural parameters or incident angle) are introduced into the resonant system, symmetry protected BICs will transformation into quasi-BICs, exhibiting as high-Q resonances with tight optical confinement. The corresponding Q factors of leaky quasi-BIC modes depend on the asymmetry degrees caused by introduced perturbations [7].

In recent years, BICs have been extensively studied in various optical designs such as gratings [8–12], photonic crystals [13–16], waveguides [17–21], dielectric metasurfaces [22–25], and plasmonic systems [26–30], enabling numerous applications in nonlinear effect enhancement [31, 32], lasers [33–35], sensors [36, 37], and filters [9]. System with multiple BICs can provide more perspectives to flexibly manipulate the optical responses. New method for effectively tuning BICs in high-dimension parameter space is strongly needed for the development of various BIC based optical devices. Currently, the studies of optical mode modulation in high-dimension synthetic parameter space have become a hot topic due to the greatly increased tuning degrees of freedom. Based on 2D or 3D synthetic parameter space, our research group has achieved BIC mode modulation [38], charge-2 Dirac point in topological superlattices [39] and rotated Weyl physics probing [40, 41].

On the basis of above researches, by introducing the concept of synthetic parameter space, which is formed by two independent asymmetry parameters, we demonstrate a system with multiple BICs in 2D synthetic parameter space through designed nanohole metasurface array. At high-order symmetry point, there are two BIC modes controlled by the x-direction asymmetry parameter and one BIC mode controlled by the y-direction asymmetry parameter. These BIC modes can achieve high Q factors which satisfies the common properties of symmetry-protected BICs. Away from this point, corresponding BIC mode will transform into quasi-BIC mode due to symmetry broken, manifesting as absorption peaks. Furthermore, BIC lines formed by a series of BIC modes can be achieved through adjusting the asymmetry parameters in proposed 2D synthetic parameter space. In experiments, measured absorption spectra under varied asymmetry parameters show the same trend as theoretical evolution of BICs. Benefiting from the continuously tunable optical responses of BICs, this proposed nanohole metasurface design can further enhance the light–matter interaction, opening a new avenue for plenty of physics fields in laser, sensing, and nonlinear optics.

2 Results and discussion

2.1 Multiple BICs in 1D parameter space

As shown in Figure 1a, our proposed multiple BICs system is composed of nanohole metasurface array, a gold mirror, and a dielectric Si layer sandwiched between them. In particular, each unit cell of artificially designed metasurface contains two identical nanoholes with fixed length L and width W as 2.2 μm and 0.2 μm, as shown in Figure 1c. For better understanding, we consider two nanoholes as two individual parts with the period along x (y) axis Px(Py) fixed as Px = 1.5 μm (Py = 3.6 μm). If we move one nanohole while the other keeps still, the symmetry of nanoholes will be broken. According to the moving distances Δx and Δy, we define two asymmetry parameters p = Δx/(Px/2) and q = Δy/(Py/2) to form the 2D synthetic parameter space (p − q Space), as shown in Figure 1d. The SEM picture of one kind of asymmetry nanohole metasurfaces (p = 0.4, q = 0) is displayed in Figure 1b. In p − q space, the point of p = 0, q = 0 is the high-order symmetry point where the nanoholes simultaneously possess the σx and σy symmetry, which means that the structure is invariant under the mirror reflection σx and σy, as illustrated in Figure 1c. With the increase of the parameters p and q, the electric field distribution of nanoholes changes accordingly to the increased asymmetry. When two nanoholes possess two out-of-phase electric dipole field distributions, they can be regarded as a nonradiative electric quadrupole. Specifically, when the electric field of two nanoholes manifests as two out-of-phase x-polarized (y-polarized) electric dipoles, we call it the x-polarized (y-polarized) quadrupole mode.

Figure 1:

Schematic diagram of multiple BICs system based on nanohole metasurface array. (a) Schematic of designed nanohole metasurface. (b) The SEM picture of fabricated one kind of asymmetric nanohole metasurfaces (p = 0.4, q = 0) the unit cell of nanohole metasurface at high-order symmetry point where nanoholes simultaneously possess σx and σy symmetry, (c) and the asymmetry points where the symmetry of nanoholes is broken (d).

The mode properties of proposed multiple BICs system in 1D parameter space are firstly analyzed. We calculated the eigenfrequency variations of nanohole metasurfaces with varied asymmetry parameters p and q, respectively, by COMSOL Multiphysics (COMSOL Inc.), as depicted in Figure 2a and b. In the simulation, the gold material is described by Drude model with plasma frequency set as 1.37 × 10¹⁶ rad/s and the collision frequency Г set as 0 to realize lossless simulation. Here, we first focus on four modes at high-order symmetry point, two of which manifest as x-polarized electric quadrupole mode, that is, the BIC modes protected by σx symmetry, denoted as X1, X2. One mode manifests as y-polarized electric quadrupole, that is, the BIC mode protected by σy symmetry, denoted as Y1. The last ordinary mode is denoted as N1. Figure 2c shows the Ex distribution of four modes. The above three BIC modes and N1 mode are indicated by red and green circles, respectively.

Figure 2:

Multiple BIC modes investigation in 1D parameter space. (a) Simulated eigenfrequency variations of four modes X1, X2, Y1, and N1 and corresponding absorption spectra variations with varied asymmetry parameter p while q fixed as 0. The position of the BICs (non-BIC) is marked by red (green) circles. (b) Simulated eigenfrequency variations of four modes and corresponding absorption spectra variations with varied asymmetry parameter q while p fixed as 0. (c) The Ex distribution of four modes at high-order symmetry point (p = 0, q = 0). Variations of Q factors with varied asymmetry parameter p (q = 0), (d) and asymmetry parameter q (p = 0), (e), respectively, for four modes. (f) Linear fit of the Q factor of modes X1 and X2 to the p ⁻². (g) Linear fit of the Q factor of mode Y1 to the q ⁻².

To verify the existence of above three BIC modes, we further calculate the absorptivity of designed metasurfaces as a function of the asymmetry parameters p and q, respectively, as shown in Figure 2a and b. Here, the collision frequency Г of gold is set to 4.05 × 10¹⁰ rad/s to obtain the absorption results, since the absorption of lossless system is nearly 0. It should be noted that this introduced metal loss is only one thousandth of the actual metal loss (Г_{Normal Au} = 4.05 × 10¹³ rad/s), which is basically negligible. In Figure 2a, modes X1 and X2 have no corresponding absorption peaks at high-order symmetry point (indicated by the red circle), performing as nonradiative modes. However, when the parameter p changes continuously with q remains at 0, the absorption peaks corresponding to these two modes appear performing as radiative modes. As parameter p moves away from the high-order symmetry point, the linewidth of absorption peaks gradually broaden, meaning that these two radiative modes become leakier. These results show that X1 and X2 modes are two nonradiative BIC modes at high-order symmetry point. Away from this point, X1 and X2 modes gradually evolve into quasi-BIC modes, manifesting as leaky modes due to broken σx symmetry. When p is fixed at 0 and q changes continuously, X1 and X2 modes have no corresponding absorption peaks and remain as BIC modes, as shown in Figure 2b. Clearly, these two BIC modes at high-order symmetry point are protected by σx symmetry (p = 0). The slight discrepancies between the eigenfrequency and absorption peaks of modes X1 and X2 in Figure 2a are caused by the different Г set in the simulations. As shown in Figure 2c, the Ex of modes X1 and X2 at high-order symmetry point exhibit an odd symmetry distribution (x-polarized quadrupole), satisfying the central inverse symmetry (x, y) → (−x, −y). At this point, two nanoholes act as two radiation components that can cancel each other and then the radiation to the outside space is consequently forbidden. When σx symmetry is broken (p ≠ 0), the symmetry of Ex distributions change from odd to even and the radiation channel to the outside space opens up.

Different from X1, X2 modes, the absorption peak of Y1 mode disappears at the high-order symmetry point and it will appear when q changes continuously while the parameter p is fixed as 0, as shown in Figure 2b. Therefore, this mode is a BIC mode at the high-order symmetry point and will evolve into quasi-BIC mode after moving away from that point. Similarly, when p changes continuously and q remains at 0, Y1 mode has no corresponding absorption peak and remains as BIC mode, indicating that this BIC mode is protected by σy symmetry (q = 0), as shown in Figure 2a. In contrast, mode N1 always exhibits varying absorption peaks when p or q change continuously, which proves that this mode is a non-BIC mode (green circle), as shown in Figure 2a and b. The Ex of mode Y1 at high-order symmetry point similarly exhibits an odd symmetry distribution (y-polarized quadrupole) with no radiation to the outside space, as shown in Figure 2c. When σy symmetry is broken (q ≠ 0), the symmetry of Ex distribution changes from odd to even. However, the Ex of mode N1 always exhibits an even symmetry distribution, representing a radiative mode.

In order to reveal the underlying physical mechanism of these three BIC modes mentioned above, we calculated the variation of Q factors of them as the asymmetry parameters p or q change, as shown in Figure 2d and e. We can see that when q is fixed as 0 and p gradually varies, the Q factors of X1 and X2 modes reach the maximum values at the high-order symmetry point. After moving away from this point, these two modes transform into quasi-BIC modes and the Q factors of them drop rapidly. Moreover, the Q factors of two modes are stabilized near the maximum value when p is fixed as 0 and q gradually varies. In contrast, when p is fixed as 0 and q gradually varies, the Q factor of the Y1 mode reaches its maximum value at the high-order symmetry point and decreases rapidly after moving away from the point, indicating the evolution from BIC mode to quasi-BIC mode. Similarly, the Q factor of Y1 mode is stabilized near its maximum value when q is fixed as 0 and p gradually varies. For comparison, the N1 mode shows relatively low Q factor during the variation of parameters p or q. Then, through linear fitting, we find that the Q factors of X1, X2, and Y1mdoes satisfy the relationship Q ∝ p ⁻²or Q ∝ q ⁻², further verifying that these three modes are symmetry-protected BIC modes at high-order symmetry point, as shown in Figure 2f and g. The Q factor of N1 mode does not satisfy the linear relationship with the negative quadratic of the asymmetry parameter p or q. The corresponding results for these modes of nanohole metasurfaces under considering metal loss are given in Figure S1.

2.2 BIC lines in 2D synthetic parameter space

Based on the above results in 1D parameter space, we continue to investigate the optical properties of nanohole metasurfaces in 2D parameter p − q space. Here, we call the symmetry points merely protected by the σx (σy) symmetry where p = 0, q ≠ 0 (q = 0, p ≠ 0) as x-direction (y-direction) symmetry points. While the points where p ≠ 0, q ≠ 0, we call them as asymmetry points. Figure 3 shows the calculated Q factor variations of four modes X1, X2, Y1, and N1 in the p − q space, respectively.

Figure 3:

Multiple BIC modes investigation in synthetic p − q space. Distribution of Q factors in p − q space for (a) X1, (b) X2, (c) Y1, and (d) N1 modes. BIC lines are marked by red lines. Four picked points for field analysis are marked by colored dots.

In 2D parameter p − q space, it is clear that the maximum values of Q factors of X1 and X2 modes appear on the parameter line (p = 0) connected by x-direction symmetry points. Away from this line, these two BIC modes transform into quasi-BIC modes, and the Q factors of them decrease continuously with the increase of |p|. Clearly, a BIC line is formed by a series of BIC modes which exist on the parameter line (p = 0). Similarly, the maximum value of Q factor of Y1 mode appears on the parameter line connected by y-direction symmetry points. Away from this line, the Q factor of Y1 mode decreases continuously as |q| increases due to mode transformation. A new BIC line is formed by a series of BIC modes which exist on the parameter line (q = 0). These BIC lines are indicated by red lines in Figure 3.

The Q factor distribution in p − q space demonstrates that these BIC lines are protected by σx or σy symmetry. In addition, the Q factors of N1 mode present low values at all symmetry points, which is inconsistent with the BIC mode properties. In addition, we investigate the Q factors and the real part of the eigenfrequency variations of these four modes in p − q space with considering metal loss (as shown in Figure S2 and Figure S3). Along the BIC lines, BIC mode with highest Q factor can be rapidly found in p − q space, enabling the high-dimension tuning of optical performance. Benefiting from the increased tuning degrees of freedom provided by BIC lines, proposed nanohole metasurfaces can facilitate the developing of BIC based optical devices.

To better understand the physical mechanism of multiple BICs system, we select four points on the Q factor surface in p − q space for each of the four modes, namely, the high-order symmetry point (p = 0, q = 0), the x-direction symmetry point (p = 0, q = −0.2), the y-direction symmetry point (p = −0.2, q = 0) and the asymmetry point (p = −0.2, q = −0.2), for in-plane electric field distribution analysis, as shown in Figure 4. At high-order symmetry point, the Ex distribution of X1 and X2 modes in two nanoholes can be regarded as two out-of-phase electric dipoles, performing as x-polarized electric quadrupole. Similarly, the Ex distribution of Y1 mode performs as y-polarized electric quadrupoles. Those nonradiative quadrupole distributions demonstrate the simultaneous existence of three BIC modes. At x-direction symmetry point, only the Ex distributions of X1, X2 modes are maintained as electric quadrupoles while the Ex distribution of Y1 mode in two nanoholes is in-phase, indicating the transformation from BIC mode to quasi-BIC mode. On the contrary, at y-direction symmetry point, only Y1 mode maintains as electric quadrupole while the X1 and X2 modes are no longer electric quadrupole distribution. However, at asymmetry point where the symmetry is broken in both directions, the Ex distributions of X1, X2, and Y1 modes exhibit same in-phase properties and all transform into radiative quasi-BIC mode. Differently, the Ex distribution of N1 mode at four picked points exhibits the same in-phase distribution, performing as radiative non-BIC mode. To clearly show the radiation mechanism of BIC modes, we also provide the out-of-plane (x − z plane) electric field Ex and |E| distributions for each of the four modes at four points in p − q space, shown as Figure S4 and S5.

Figure 4:

Electric field distribution of modes in p-q space. The in-plane Ex distribution for X1, X2, Y1, and N1 modes at four picked points where p = 0, q = 0 (a), p = 0, q = −0.2 (b), p = −0.2, q = 0 (c), and p = −0.2, q = −0.2 (d). The red dashed box represents the BIC modes.

2.3 Experiment results and analysis

In the following, we introduced metal loss into nanohole metasurfaces and simulated the absorption spectra of metasurfaces under varied asymmetry parameters to verify the existence of BIC lines in lossy system. The absorption spectra under varied asymmetry parameter p (q = 0) were firstly simulated, as shown in Figure 5a. It can be seen that there is no resonant absorption peak for the X1 and X2 modes at the high-order symmetry point. After moving away from this point, the absorption peaks of X1, X2 modes appear and show redshift and blue shift, respectively, indicating the evolution from BICs to quasi-BICs. The absorption peaks of two modes are marked by the red dots in the figure. In contrast, the absorption peaks of N1 mode always present and its resonance wavelength remains almost unchanged. Above results are consistent with that in lossless system in Figure 2a. Note that, compared to lossless system, the absorption peaks of lossy system are broadened due to the presence of metal loss.

Figure 5:

Experimentally demonstrated multiple BICs system in p − q space. Simulated absorption spectra in lossy system (a) and experimentally measured results (b) of nanohole metasurfaces under varied asymmetry parameter p (q = 0). The red line indicates the nonradiative BIC modes (X1 X2) and the red dots indicate the absorption peaks corresponding to the quasi-BIC modes. Comparison of Q factors obtained from simulation and experiment for X1 (c) and X2 (d) modes.

Specifically, we fabricated designed nanohole metasurface array with varied asymmetry parameter p and measured the absorption spectra, as shown in Figure 5b. The variation tendency of absorption peaks of X1, X2 modes show good agreement with the simulation results, verifying the existence of two BIC modes in experiments. At the same time, the Y1 mode shows no corresponding absorption peak and thus maintain as BIC mode, which means that there is a BIC line formed by a series of BICs under varied p. Moreover, the Q factors of X1, X2 modes in experimental and simulated results are compared, as shown in Figure 5c and d. Both of them show the same trend and are consistent with the properties of the σx symmetry protected BIC modes. The slight discrepancies between the Q factors of simulation and experimental results are caused by fabrication errors.

Then, we investigate the absorption spectra variations under varied asymmetry parameter q (p = 0), as shown in Figure 6a and b. It can be seen that the experiment results agree well with the simulation results. There is no absorption peak for the Y1 mode at high-order symmetry point, but away from this point, corresponding absorption peaks appear and gradually shift, indicating the appearance of quasi-BIC. The variation tendency is same as that of lossless system in Figure 2b. Meanwhile, the X1, X2 modes show no corresponding absorption peak and thus both maintain as BIC modes, which means that there are BIC lines formed by a series of BIC modes under varied q. In Figure 6c, the Q factors of experiment results show the same tendency as that of simulation results and are consistent with the properties of the σy symmetry protected BIC modes.

Figure 6:

Above results clearly show the generation and evolution of BIC lines in p − q space at the condition of considering metal loss. Therefore, our proposed multiple BICs system provides more degrees of freedom for BICs tuning, offering a new platform for multifunctional integrated optical devices. For example, by introducing this design method into semiconductor nanolasing regime, low-threshold integrated nanolasing array are expected to be realized. In addition, most studied BIC-based sensors have focused on the wavelength range from visible to near-infrared. However, our design can basically cover the long-wave infrared range of 8–14 μm by asymmetry parameter manipulation, which can be expected to be applied for integrated high-resolution sensors for infrared molecular detection.

3 Conclusions

In summary, based on nanohole metasurface array, we have realized multiple BIC lines formed by a series of BICs in 2D synthetic parameter space. Desired BIC mode with highest Q factor can be found along the BIC lines. Through adjusting the asymmetry parameters in p − q space, the continuously tuning of BICs are simultaneously realized in simulation and experiments. Specifically, the physical mechanisms underlying the generation and evolution of BICs have been carefully investigated. This proposed system with multiple BICs can find many potential optical devices in nonlinear optics, nanolasing and infrared sensing. Meanwhile, our proposed new method for effectively tuning BICs in high-dimension synthetic parameter space can be applied to various BIC systems, paving the way for multidimensional manipulation.

4 Methods

4.1 Nanohole metasurface array fabrication

Firstly, three layers of Au (70 nm)/Si(300 nm)/Au(100 nm) are sequentially deposited on the Si substrate by electron beam evaporation (AdNaNotek EBS-150U). Then nanohole metasurface is etched from the top gold film by a focused ion beam (FIB dual-beam FEI Helios 600i, 30 keV, 100 pA). By adjusting the etching parameters, the asymmetric parameters p and q of nanohole metasurface array can be varied from −1 to 1 while the etching depth is set as 70 nm. Each sample of the metasurface array in p–q space is 100 μm × 100 μm sized.

4.2 Optical characterization

The absorption spectra (A) are derived based on the measured reflection spectra (R) by A = 1 − R. Specifically, each reflection spectrum (R) of the nanohole metasurface array was measured by a Fourier transform infrared (FTIR) spectrometer. The reflection signals in the spectral range of 4–16 μm were collected using a Hyperion 2000 IR microscope with a liquid-nitrogen-cooled HgCdTe (MCT) detector. Measured reflection spectra were normalized with respect to a gold mirror.

4.3 Numerical simulations and analysis

To find eigenmodes of multiple BICs system, numerical simulations are performed to calculate the eigenfrequency of nanohole metasurface array by the eigenfrequency solver in COMSOL Multiphysics. The Q factor is obtained from the real and imaginary parts of the eigenfrequency:

Q = r e a l E i g e n f r e q u e n c y 2 × i m a g E i g e n f r e q u e n c y

The absorption of each individual metasurface can also be simulated. During the process, periodic boundary conditions are applied in both x and y directions for mimicking the periodic nanohole array.

Corresponding author: Hui Liu, National Laboratory of Solid State Microstructures, School of Physics, Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing, Jiangsu 210093, China, E-mail: liuhui@nju.edu.cn

Fengyuan Zhang and Qiongqiong Chu are contributed equally to this work.

Funding source: National Natural Science Foundation of China

Award Identifier / Grant number: 92150302

Award Identifier / Grant number: 92163216

Author contributions: All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.
Research funding: This work was financially supported by the National Natural Science Foundation of China (Nos. 92163216 and 92150302).
Conflict of interest statement: The authors declare no conflicts of interest regarding this article.

References

[1] J. von Neumann and E. P. Wigner, “Über merkwürdige diskrete Eigenwerte,” in Phys. Z., A. S. Wightman, Ed., Berlin, Heidelberg, Springer, 1993, pp. 291–293.10.1007/978-3-662-02781-3_19Search in Google Scholar

[2] H. Friedrich and D. Wintgen, “Interfering resonances and bound-states in the continuum,” Phys. Rev. A, vol. 32, no. 6, pp. 3231–3242, 1985. https://doi.org/10.1103/physreva.32.3231.Search in Google Scholar PubMed

[3] C. W. Hsu, B. Zhen, A. D. Stone, J. D. Joannopoulos, and M. Soljacic, “Bound states in the continuum,” Nat. Rev. Mater., vol. 1, no. 9, p. 16048, 2016. https://doi.org/10.1038/natrevmats.2016.48.Search in Google Scholar

[4] F. Ursell, “Trapping modes in the theory of surface waves,” Math. Proc. Cambridge Philos. Soc., vol. 47, no. 2, pp. 347–358, 1951. https://doi.org/10.1017/s0305004100026700.Search in Google Scholar

[5] M. Robnik, “A simple separable Hamiltonian having bound-states in the continuum,” J. Phys. A: Math. Gen., vol. 19, no. 18, pp. 3845–3848, 1986. https://doi.org/10.1088/0305-4470/19/18/029.Search in Google Scholar

[6] F. Ursell, “Trapped modes in a circular cylindrical acoustic waveguide,” Proc. R. Soc. London, Ser. A, vol. 435, no. 1895, pp. 575–589, 1991.10.1098/rspa.1991.0162Search in Google Scholar

[7] K. Koshelev, S. Lepeshov, M. K. Liu, A. Bogdanov, and Y. Kivshar, “Asymmetric metasurfaces with high-Q resonances governed by bound states in the continuum,” Phys. Rev. Lett., vol. 121, no. 19, p. 193903, 2018. https://doi.org/10.1103/physrevlett.121.193903.Search in Google Scholar

[8] D. C. Marinica, A. G. Borisov, and S. V. Shabanov, “Bound states in the continuum in photonics,” Phys. Rev. Lett., vol. 100, no. 18, p. 183902, 2008. https://doi.org/10.1103/physrevlett.100.183902.Search in Google Scholar

[9] J. M. Foley, S. M. Young, and J. D. Phillips, “Symmetry-protected mode coupling near normal incidence for narrow-band transmission filtering in a dielectric grating,” Phys. Rev. B, vol. 89, no. 16, p. 165111, 2014. https://doi.org/10.1103/physrevb.89.165111.Search in Google Scholar

[10] H. Hemmati and R. Magnusson, “Resonant dual-grating metamembranes supporting spectrally narrow bound states in the continuum,” Adv. Opt. Mater., vol. 7, no. 20, p. 1900754, 2019. https://doi.org/10.1002/adom.201900754.Search in Google Scholar

[11] E. A. Bezus, D. A. Bykov, and L. L. Doskolovich, “Integrated diffraction gratings on the Bloch surface wave platform supporting bound states in the continuum,” Nanophotonics, vol. 10, no. 17, pp. 4331–4340, 2021. https://doi.org/10.1515/nanoph-2021-0352.Search in Google Scholar

[12] Z. L. Deng, F. J. Li, H. A. Li, X. P. Li, and A. Alu, “Extreme diffraction control in metagratings leveraging bound states in the continuum and exceptional points,” Laser Photonics Rev., vol. 16, no. 6, p. 2100617, 2022. https://doi.org/10.1002/lpor.202100617.Search in Google Scholar

[13] M. Minkov, D. Gerace, and S. H. Fan, “Doubly resonant χ2 nonlinear photonic crystal cavity based on a bound state in the continuum,” Optica, vol. 6, no. 8, pp. 1039–1045, 2019. https://doi.org/10.1364/optica.6.001039.Search in Google Scholar

[14] S. Vaidya, W. A. Benalcazar, A. Cerjan, and M. C. Rechtsman, “Point defect localized bound states in the continuum in photonic crystals and structured fibers,” Phys. Rev. Lett., vol. 127, no. 2, p. 023605, 2021. https://doi.org/10.1103/physrevlett.127.023605.Search in Google Scholar PubMed

[15] S. Neale and E. A. Muljarov, “Accidental and symmetry-protected bound states in the continuum in a photonic-crystal slab: a resonant-state expansion study,” Phys. Rev. B, vol. 103, no. 15, p. 155112, 2021. https://doi.org/10.1103/physrevb.103.155112.Search in Google Scholar

[16] M. Kang, S. P. Zhang, M. Xiao, and H. X. Xu, “Merging bound states in the continuum at off-high symmetry points,” Phys. Rev. Lett., vol. 126, no. 11, p. 117402, 2021. https://doi.org/10.1103/physrevlett.126.117402.Search in Google Scholar

[17] C. L. Zou, J. M. Cui, F. W. Sun, et al.., “Guiding light through optical bound states in the continuum for ultrahigh-Q microresonators,” Laser Photonics Rev., vol. 9, no. 1, pp. 114–119, 2015. https://doi.org/10.1002/lpor.201400178.Search in Google Scholar

[18] X. W. Gao, B. Zhen, M. Soljacic, H. S. Chen, and C. W. Hsu, “Bound states in the continuum in fiber bragg gratings,” ACS Photonics, vol. 6, no. 11, pp. 2996–3002, 2019. https://doi.org/10.1021/acsphotonics.9b01202.Search in Google Scholar

[19] A. F. Sadreev, “Interference traps waves in an open system: bound states in the continuum,” Rep. Prog. Phys., vol. 84, no. 5, p. 055901, 2021. https://doi.org/10.1088/1361-6633/abefb9.Search in Google Scholar PubMed

[20] V. Ardizzone, F. Riminucci, S. Zanotti, et al.., “Polariton Bose-Einstein condensate from a bound state in the continuum,” Nature, vol. 605, no. 7910, pp. 447–452, 2022. https://doi.org/10.1038/s41586-022-04583-7.Search in Google Scholar PubMed

[21] F. Ye, Y. Yu, X. Xi, and X. K. Sun, “Second-Harmonic generation in etchless lithium niobate nanophotonic waveguides with bound states in the continuum,” Laser Photonics Rev., vol. 16, no. 3, p. 2100429, 2022. https://doi.org/10.1002/lpor.202100429.Search in Google Scholar

[22] F. Yesilkoy, E. R. Arvelo, Y. Jahani, et al.., “Ultrasensitive hyperspectral imaging and biodetection enabled by dielectric metasurfaces,” Nat. Photonics, vol. 13, no. 6, pp. 390–396, 2019. https://doi.org/10.1038/s41566-019-0394-6.Search in Google Scholar

[23] A. S. Kupriianov, Y. Xu, A. Sayanskiy, V. Dmitriev, Y. S. Kivshar, and V. R. Tuz, “Metasurface engineering through bound states in the continuum,” Phys. Rev. Appl., vol. 12, no. 1, p. 014024, 2019. https://doi.org/10.1103/physrevapplied.12.014024.Search in Google Scholar

[24] D. R. Abujetas, A. Barreda, F. Moreno, A. Litman, J. Geffrin, and J. A. Sanchez-Gil, “High-Q transparency band in all-dielectric metasurfaces induced by a quasi bound state in the continuum,” Laser Photonics Rev., vol. 15, no. 1, p. 2000263, 2021. https://doi.org/10.1002/lpor.202000263.Search in Google Scholar

[25] J. F. Algorri, F. Dell’Olio, P. Roldan-Varona, et al.., “Strongly resonant silicon slot metasurfaces with symmetry-protected bound states in the continuum,” Opt. Express, vol. 29, no. 7, pp. 10374–10385, 2021. https://doi.org/10.1364/oe.415377.Search in Google Scholar

[26] W. Zhang, A. Charous, M. Nagai, D. M. Mittleman, and R. Mendis, “Extraordinary optical reflection resonances and bound states in the continuum from a periodic array of thin metal plates,” Opt. Express, vol. 26, no. 10, pp. 13195–13204, 2018. https://doi.org/10.1364/oe.26.013195.Search in Google Scholar PubMed

[27] L. Q. Cong and R. Singh, “Symmetry-protected dual bound states in the continuum in metamaterials,” Adv. Opt. Mater., vol. 7, no. 13, p. 1900383, 2019. https://doi.org/10.1002/adom.201900383.Search in Google Scholar

[28] X. G. Zhao, C. X. Chen, K. Kaj, et al.., “Terahertz investigation of bound states in the continuum of metallic metasurfaces,” Optica, vol. 7, no. 11, pp. 1548–1554, 2020. https://doi.org/10.1364/optica.404754.Search in Google Scholar

[29] S. Y. Sun, Y. F. Sun, H. Z. Li, et al.., “Tunable plasmonic bound states in the continuum in the visible range,” Phys. Rev. B, vol. 103, no. 4, p. 045416, 2021. https://doi.org/10.1103/physrevb.103.045416.Search in Google Scholar

[30] S. Q. Li, C. H. Du, L. Z. Yin, Z. C. Gao, and F. H. Li, “Quasi-Bound states in the continuum of localized spoof surface plasmons,” ACS Photonics, vol. 9, no. 10, pp. 3333–3340, 2022. https://doi.org/10.1021/acsphotonics.2c00673.Search in Google Scholar

[31] A. P. Anthur, H. Z. Zhang, R. Paniagua-Dominguez, et al.., “Continuous wave second harmonic generation enabled by quasi- bound-states in the continuum on gallium phosphide metasurfaces,” Nano Lett., vol. 20, no. 12, pp. 8745–8751, 2020. https://doi.org/10.1021/acs.nanolett.0c03601.Search in Google Scholar PubMed

[32] Z. J. Huang, M. J. Wang, Y. Li, et al.., “Highly efficient second harmonic generation of thin film lithium niobate nanograting near bound states in the continuum,” Nanotechnology, vol. 32, no. 32, p. 325207, 2021. https://doi.org/10.1088/1361-6528/abfe23.Search in Google Scholar PubMed

[33] A. Kodigala, T. Lepetit, Q. Gu, B. Bahari, Y. Fainman, and B. Kante, “Lasing action from photonic bound states in continuum,” Nature, vol. 541, no. 7636, pp. 196–199, 2017. https://doi.org/10.1038/nature20799.Search in Google Scholar PubMed

[34] S. T. Ha, Y. H. Fu, N. K. Emani, et al.., “Directional lasing in resonant semiconductor nanoantenna arrays,” Nat. Nanotechnol., vol. 13, no. 11, pp. 1042–1047, 2018. https://doi.org/10.1038/s41565-018-0245-5.Search in Google Scholar PubMed

[35] R. M. Ma and R. F. Oulton, “Applications of nanolasers,” Nat. Nanotechnol., vol. 14, no. 1, pp. 12–22, 2019. https://doi.org/10.1038/s41565-018-0320-y.Search in Google Scholar PubMed

[36] Y. H. Liu, W. D. Zhou, and Y. Z. Sun, “Optical refractive index sensing based on high-Q bound states in the continuum in free-space coupled photonic crystal slabs,” Sensors, vol. 17, no. 8, p. 1861, 2017. https://doi.org/10.3390/s17081861.Search in Google Scholar PubMed PubMed Central

[37] Y. L. Wang, Z. H. Han, Y. Du, and J. Y. Qin, “Ultrasensitive terahertz sensing with high-Q toroidal dipole resonance governed by bound states in the continuum in all-dielectric metasurface,” Nanophotonics, vol. 10, no. 4, pp. 1295–1307, 2021. https://doi.org/10.1515/nanoph-2020-0582.Search in Google Scholar

[38] F. Zhong, K. Ding, Y. Zhang, S. Zhu, C. Chan, and H. Liu, “Angle-resolved thermal emission spectroscopy characterization of non-hermitian metacrystals,” Phys. Rev. Appl., vol. 13, no. 1, p. 014071, 2020. https://doi.org/10.1103/physrevapplied.13.014071.Search in Google Scholar

[39] M. Y. Hu, K. Ding, T. Qiao, et al.., “Realization of photonic charge-2 Dirac point by engineering super-modes in topological superlattices,” Commun. Phys., vol. 3, no. 1, p. 130, 2020. https://doi.org/10.1038/s42005-020-00395-1.Search in Google Scholar

[40] Q. Wang, M. Xiao, H. Liu, S. N. Zhu, and C. T. Chan, “Optical interface states protected by synthetic Weyl points,” Phys. Rev. X, vol. 7, no. 3, p. 031032, 2017. https://doi.org/10.1103/physrevx.7.031032.Search in Google Scholar

[41] Z. W. Yan, Q. Wang, M. Xiao, Y. L. Zhao, S. N. Zhu, and H. Liu, “Probing rotated Weyl physics on nonlinear lithium niobate-on-insulator chips,” Phys. Rev. Lett., vol. 127, no. 1, p. 013901, 2021. https://doi.org/10.1103/physrevlett.127.013901.Search in Google Scholar