Abstract
1. INTRODUCTION
Trapping electromagnetic waves in optics systems is critical for enabling strong light–matter interactions that facilitate high-performance optical devices. To achieve high quality () factors so that electromagnetic waves can be trapped for a longer period before decaying, a variety of resonant structures have been proposed [1–7]. Nevertheless, the factors of these localized states are fundamentally limited by inevitable radiative losses and defects in the fabrication process [8]. Thus, there has been a persistent desire for achieving a perfect bound state where electromagnetic waves can be localized for an infinite time.
Recently, a novel class of localized states, so-called bound states in the continuum (BICs), has greatly attracted researchers’ attention [8–10]. A surprising feature of BICs is that they theoretically process infinite factors even lying inside continuous spectra. Empowered by their unique features to confine light, BICs have been widely observed for various waves, including quantum waves [11–13], water waves [14,15], electromagnetic waves [16–21], and acoustic waves [22–26]. To date, researches on BICs have emerged in various optical platforms such as metamaterials [27–30] and photonic crystals (PhCs) [31–44], relying on their flexible block construction and size-dependent spectral scalability. BICs are generally divided into two categories: symmetry-protected and accidental. Symmetry-protected BICs are usually fixed at high-symmetry points in the Brillouin zone [27,28], i.e., the point, where the states are completely decoupled from continuous spectra due to incompatible symmetry with outgoing waves in the radiation channels. Nevertheless, accidental BICs originate from the destructive interference between radiation channels of the system [45–48]. Therefore, accidental BICs, which usually exist at off-high symmetry points, are no longer dependent on structural symmetry but are sensitive to their spatial and geometric parameters.
In practical applications, the factors of BICs are intensively limited due to inevitable fabrication imperfections. An effective approach to overcome this issue is to exploit the topological properties of BICs [49–51]. BICs have been verified to be singular points of polarizations in momentum space, known as eigenmodes of an optical system whose far-field polarization direction is undefined. As a proven momentum-space vortex center for linearly polarized far fields, BICs have an integer topological charge (). The topological nature ensures the robust existence of BICs as long as and symmetries of the system are preserved [51]. Topological charge conservation indicates that BICs can be moved in reciprocal space with the symmetry of the system unchanged [52]. Thus, multiple isolated BICs can be tuned to the same wave vector to construct merging BICs by carefully varying structural parameters. For example, Jin
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
On the other hand, though and off- merging BICs can significantly enhance the factors of quasi-BICs in a wide range of wave vectors, the resonance wavelength will shift obviously with the increase of incident angle due to the dispersive band. Therefore, in practical applications, tiny dispersion is meaningful to stabilize the operating wavelength of quasi-BICs, especially for those evolved from off- merging BICs. As far as we know, the current research on flatband BICs [56] provides an excellent mechanism for combining a moiré flatband with a quasi-BIC by controlling the twist angle and the distance of bilayer PhC slabs. Such single moiré BICs offer a new avenue to overcome the problem of angular sensitivity. However, the performance always suffers from scattering loss caused by fabrication imperfections. Thus, utilizing flatband merging BICs to realize angle and disturbance robustness is significant and useful.
In this paper, we propose a new scheme to construct merging BICs at almost arbitrary wave vectors on a quasi-flatband without breaking the structural symmetry. A symmetry-protected BIC can be designed at point by using a PhC with symmetry. At off- points, a few accidental BICs can be formed by adjusting the structural parameters. Since the symmetry of the structure is preserved when only the periodic constant a is changed, the symmetry-protected BIC is fixed at point, while accidental BICs will move their positions in momentum space owing to their sensitivity to structural parameters. Thus, as we first increase the periodic constant , the accidental BICs near will move towards the center of the Brillouin zone () along and directions, and merge with the symmetry-protected BICs pinned at the point. Since these accidental BICs have the same topological charge, they cannot annihilate each other. Owing to the restriction of the square lattice symmetry, they are also unable to further form higher-order BICs, and will therefore deflect in direction. A perfect accidental BIC does not exhibit the Fano feature at a certain oblique incidence angle since it is decoupled from far-field radiation. However, incident light will cause sharp Fano resonances when we slightly change the structural parameters, such as periodic constant . The variation of structural parameters will affect the wave vector and amplitude of the resonance, resulting in the moving of accidental BICs. Next, we show that one accidental BIC and another accidental BIC gradually approach each other along direction by decreasing and finally merge at off- points. Since the topological charges of the two accidental BICs are opposite, they will eventually annihilate into a quasi-BIC, as governed by charge conservation. Furthermore, we demonstrate that the merging at off- points can occur at any arbitrary wave vector by adjusting structural parameters without breaking symmetry. Compared with an isolated BIC, merging multiple BICs at point or off- points can improve the scaling property from to or , respectively. Especially, such and off- merging BICs occur on an isolated band with flat dispersion, which means that they have a more stable frequency under varying oblique incidence angles. Finally, as an example of application, efficient angle-insensitive second-harmonic generation (SHG) assisted by merging BICs is demonstrated. Such steerable merging BICs on a quasi-flatband provide an excellent mechanism to design a variety of optical devices that require high factors, including low-threshold lasers [34,57,58], ultrasensitive sensing [59,60], nonlinear optics [29,61–63], and narrowband filters [64].
2. DESIGN AND MECHANISM
To demonstrate that merging BICs can be manipulated without reducing structural symmetry, we consider a PhC slab that consists of dielectric pillars (relative permittivity ) arranged in a square lattice [see Fig. 1(a)], where , , and . The energy bands of structures can be classified into TE-like and TM-like bands, corresponding to and , respectively at the mirror plane in the direction. The calculated TE-like bands using COMSOL Multiphysics are shown in Fig. 1(b), where we mainly focus on flatband TE2 in the middle, marked in red. The other two adjacent bands, TE1 and TE3, correspond to blue and black solid lines, respectively. The mode profiles () at point for the three bands are shown in the insets of Fig. 1(b). These bands are inside the light cone where there is a continuum of radiative modes in the surrounding material. Since they are below the diffraction limit, they can radiate into free space only through zeroth-order diffraction. As shown by the red circle in Fig. 1(c), TE2 shows five BICs with infinite factors along high-symmetry direction , one of which is a symmetry-protected BIC fixed at point and the others are accidental BICs formed by destructive interference. In contrast, for band TE3, there is only one symmetry-protected BIC at point [blue circle in Fig. 1(c)]. For more clarity, we give the 3D band structure of the system as shown in Fig. 1(d). The colors represent the factors of the eigenmodes, which approach infinite at these BICs. As long as the symmetry of the structure is maintained, the symmetry-protected BIC will be pinned at point, while the position of the accidental BIC is sensitive to structural parameters and can be moved easily in momentum space. Thus, the simultaneous existence of multiple different types of BICs on the flatband makes it possible to construct merging BICs at arbitrary wave vectors with nearly constant frequency. It should be noted that since the isolated BICs are all located in and , the formation of merging BICs will also occur in and . Therefore, the discussion on flatbands is considered along the direction of high symmetry.
Figure 1.(a) Schematic of the proposed PhC slab with lattice constant
A. Photonic Crystal Exhibiting a
First, we show that four accidental BICs near the center of the Brillouin zone can be tuned to merge with the symmetry-protected BIC at point. As the vortex polarization singularities ( points), the topological properties of BICs can be revealed by calculating the polarization distribution of far-field radiation, which is the projection of the Bloch wave function on a plane basis. The topological charge () of the polarization vortex can be defined as the winding number of the polarization vectors around a BIC in reciprocal space. That is,
Here, is a closed simple path in space along the counterclockwise direction that surrounds the vortex center, and is the angle of polarization where and denote the electric fields along and directions, respectively. Thus, the polarization field distribution of the vortex (antivortex) corresponds to a positive (negative) topological charge. Figure 2(a) depicts the evolution of several BICs close to the center of the Brillouin zone on band TE2 as the periodic constant increases gradually. When , five divergent points of the factor can be observed, representing one symmetry-protected BIC and four accidental BICs. Polarization vortices emerge around each BIC, and topological charges of the symmetry-protected BIC and surrounding accidental BICs are and , respectively. The symmetry-protected BIC is locked at point, while the topologically protected accidental BIC will move with the change of structural parameters. As the periodic constant increases to 1101 nm, accidental BICs are close to the center along the direction and merge with symmetry-protected BICs. The merging BICs formed at point significantly improve the factors of surrounding modes in a wide wave vector range. By further increasing , BICs with the same topological charge of cannot annihilate and will deflect in the direction as shown in the lower panel of Fig. 2(a).
Figure 2.(a) Simulated polarization vectors (white arrows) around the BICs near
To confirm the above BIC merging phenomenon, the factors along high-symmetry directions and at different periodic constants are provided in Fig. 2(b). The red, blue, and magenta lines represent factors on the TE2 band for , 1101, and 1261 nm, respectively. Before charges merge at the center of the Brillouin zone (, red), the factors decay roughly as away from . Herein, and are wave vectors of two accidental BICs in the direction. The contribution of accidental BICs that are relatively far from point to factors in the area discussed above cannot be ignored, so it is taken into account when fitting curves. By increasing , accidental BICs closest to the center gradually approach and finally meet symmetry-protected BIC at (red). This corresponds to the merging BICs case that factors decay as and are sufficiently large over a broader wave vector range around the diverging point compared with those of isolated BICs before merging. If we further increase the period to 1261 nm, the accidental BICs will be kicked away from point along the direction, leaving only one symmetry-protected BIC at point. In that case, the factor is almost not affected by accidental BICs in a small wave vector range and exhibits the rule as [see the magenta line in Fig. 2(b)].
So far, we have investigated the merging of BICs at the point by tuning accidental BICs to merge with symmetry-protected BICs. It is relatively easy to achieve merging BICs at high-symmetry points in reciprocal space, such as point. In contrast, merging multiple BICs at an off- point is rare, especially if structural symmetry is preserved. In previous work, the implementation of merging BICs at an off- point often needs to reduce the symmetry of the structure [54,55], which will undoubtedly bring inevitable difficulties in the fabrication process. To date, there is not yet a feasible mechanism for implementing tunable merging BICs at and off- points while maintaining structural symmetry. In the following, we demonstrate that merging BICs at an off- point can be realized by merging two accidental BICs without breaking structural symmetry.
B. Merging BICs at Off-
We still consider the PhC slab and focus on the Brillouin zone that deviates from . As shown in Fig. 3(a), in addition to the above-mentioned accidental BIC that participates in the merging BICs of point, there is another accidental BIC farther from located at on the TE2 band. The far-field polarization of these two accidental BICs exhibits characteristics of the vortex (charge ) and antivortex (charge ), so they have opposite topological charges. By keeping other parameters unchanged and only gradually decreasing the periodic constant , these two accidental BICs are tuned to approach each other along direction. When decreases to 1083 nm, they meet at an off- point and form merging BICs. Owing to their opposite topological charges, the two accidental BICs will annihilate each other if we further decrease , and will evolve into a quasi-BIC as shown by the dashed circle in the last panel of Fig. 3(a), where the factors are not infinite but remain large and can be regarded as a supercavity. If the structure possesses higher symmetry, it is possible to merge the basic BICs into a higher-order BIC at an off- point. In the above evolutionary process, the symmetry-protected BIC has been fixed at the point since the structural symmetry has not been changed. We also give the factors distribution in the direction as shown in Fig. 3(b), indicating that the factors near merging BICs (red) have been significantly enhanced over a broad wave vector compared with either isolated symmetry-protected BICs or isolated accidental BICs (blue). The factors decay as at , while the factors of merging BICs constructed from two accidental BICs at follow the scaling rule as , which greatly improves the performance of potential applications requiring momentum selection. The factors are dominated by the intrinsic loss of materials at merging BICs, and thus high-refractive-index dielectric material with ultralow intrinsic loss is an excellent choice to create a high- cavity (see Appendix C). It should be noted that accidental BICs originate from destructive interference between radiation channels, which is common in various optical systems. Therefore, this mechanism of using multiple accidental BICs to realize the merging phenomenon is also applicable to other systems, such as typical air-hole PhCs.
Figure 3.(a) Simulated polarization vectors (white arrow) around the accidental BICs with
C. Merging BICs at an Arbitrary Location
Next, we demonstrate that merging BICs can be achieved at an arbitrary point in reciprocal space by changing periodic constant and height simultaneously but keeping symmetry. Whether decreasing (increasing) the periodic constant or increasing (decreasing) the height , the above two accidental BICs located at off- point will approach (away from) each other along direction. For the change of periodic constant , the moving speed of the two accidental BICs in momentum space is almost the same. However, the accidental BIC with topological charge is much more sensitive to the height than the one with topological charge , which means that by changing the height , the different accidental BICs have different moving scales (see Appendix A). This phenomenon of asymmetric shift may be explained by Fabry–Perot (FP) mode with different orders participating in the formation of different accidental BICs, which generally exists in the PhC slab depending on the height . Owing to the ultralow factor of FP resonance, they are often regarded as the background spectrum, while the resonance properties are ignored [40]. Because the types and number of BICs have not changed during the formation of merging BICs, the phases (, , ) of BICs defined in Ref. [40] always belong to areas of (2, 1, 0), where , , and denote the number of accidental BICs, symmetry-protected BICs, and Friedrich–Wintgen BICs (FW-BICs), respectively. It is worth noting that this is different from previous work [55] based on FW-BIC with (1, 0, 1) phase. When the height changes, the frequency shift of lower-order FP mode is slower, so that the accidental BICs resulting from the coupling of the magnetic dipole (MD) mode and different orders of FP mode have asymmetric moving speeds. Thus, merging BICs can be designed to appear at an arbitrary point in direction as long as we choose an appropriate periodic constant and height simultaneously. For example, as shown in Fig. 4(a), when () is equal to (1048.7 nm, 845 nm), (1083.2 nm, 850 nm), and (1115 nm, 857 nm), we tune merging BICs to , , and , respectively. The factors satisfy the scaling rule shown in Fig. 4(b), which further confirms the existence of merging BICs.
Figure 4.(a) Simulated polarization vectors (white arrows) around accidental BICs at different wave vectors with
D. Quasi-Flatband
Above, we have proved that merging BICs can be constructed at arbitrary wave vectors while maintaining structural symmetry. The unique topological property of BICs ensures that factors can be greatly improved in a larger wave vector range by merging several BICs that are robust against fabrication imperfections (see Appendix B) and improve the performance of optoelectronic devices. However, merging BICs at either point or off- point usually occur in a band with large dispersion [54,55,65], which limits its performance in many practical applications based on angle selectivity and stable operating frequency. Next, by comparing the resonance modes near point on TE1 and TE3 energy bands, we prove that the quasi-BIC evolved from a merging BIC on flatband TE2 is insensitive to the incident angle, not only in terms of factors but also the resonance wavelength. At , a merging BIC is formed at point on the TE2 band [red line in Fig. 2(b)], while there is always an isolated symmetry-protected BIC on the TE3 band [blue circle in Fig. 2(b)] due to the maintenance of structural symmetry. We consider that the PhC slab is illuminated by a TE obliquely incident plane wave polarized along the axis. As the incident angle increases from 0°, the two BICs transform into quasi-BICs. We call the quasi-BIC that evolves from symmetry-protected BIC as quasi-BIC1, and that from merging BICs as quasi-BIC2. The simulated transmittance spectra of two quasi-BICs under different incident angles are shown in Figs. 5(a) and 5(b). At normal incidence (), one symmetry-protected BIC and one merging BIC exist in TE3 and TE2, respectively, having vanishing linewidths, and the electric dipole (ED) resonance mode located on TE1 can be observed. The isolated BIC on the TE3 band can be tailored, and the dispersion can be flattened to reduce the sensitivity to the angle of incidence; see Figs. 1(b) and 5(a). Nevertheless, the factor of the flattened quasi-BIC1 drops rapidly as the incident angle increases, as shown in Fig. 5(a). Interestingly, by precise design, we can achieve merging BICs in a flattened band that possesses an ultrahigh resonance with a certain degree of incident-angle robustness. One can see from Fig. 5(b) that the resonance wavelength and the factors of quasi-BIC2 are nearly unchanged even at a large incident angle. As increases from 0° to 12°, the relative shift of ED mode is equal to 1%, where is the resonance peak wavelength shift, and is the wavelength when , while that of quasi-BIC2 is only 0.1%.
Figure 5.(a) Transmission spectra of the quasi-BIC1 originated from isolated BIC at TE3 under different incident angles
To gain deep insight into the spectral response of the two distinct quasi-BICs and the ED mode, we first calculated the scattering powers of different multipoles of three resonance modes in the Cartesian coordinate system based on multipole scattering theory [66], as visualized in Figs. 5(d)–5(f). It can be observed from Fig. 5(e) that the contribution of the electric quadrupole (EQ) moment (orange line) is dominant, which indicates that the strong Fano resonance (quasi-BIC1) is mainly induced by the EQ moment. Different from quasi-BIC1, quasi-BIC2 stems from MD resonance since the contribution of the MD moment (red line) is dominated, as shown in Fig. 5(f). To further understand that the merging BIC occurs in the flatband, we choose quasi-BIC2 at 1551 nm as the operating wavelength and give the magnetic field distribution of the structure under different incident angles, as shown in Fig. 5(c). One can see that the magnetic fields can be strongly enhanced as the incident angle changes because the factors and resonance wavelength of quasi-BIC2 are almost unchanged as the incident angle increases. Namely, quasi-BIC2 is located in a flatband. For comparison, we also choose a conventional ED mode at 1567 nm as the operating wavelength and also give the magnetic field distribution of the structure under different incident angles, as shown in Fig. 5(d). It can be seen that, compared with the flat BIC2, the enhancement of magnetic fields assisted by conventional ED mode at 1567 nm is rapidly changed under different incident angles. Since the formation of resonance-trapped BICs in metamaterial originates from local resonance in structural elements that can be excited by plane waves in different directions, they are insensitive to the incident angle and therefore correspond to almost flat bands [67,68]. Note that merging BICs within a flatband would be very useful in some applications as shown below.
E. Efficient SHG Enhanced by Flat Merging BICs
Finally, we investigated the superior performance in the enhancement of SHG by flatband merging BICs. AlGaAs is chosen as the material of a dielectric cylinder, which has been widely used for SHG owing to its large second-order nonlinear coefficient and excellent broadband transparency. For the dispersion of the refractive index of AlGaAs, we used the analytical model proposed in Ref. [69]. We first consider merging BICs at point. At , , and , the numerically simulated transmission spectra of quasi-BIC1 and quasi-BIC2 by sweeping incident angle are given in Figs. 6(a) and 6(b), respectively. It is clear that quasi-BIC1 is almost obscured by the surrounding background when reaches 15°, while quasi-BIC2 still maintains an extremely narrow linewidth, and its relative shift is only approximately 0.12%, which exhibits an excellent ability to enhance SHG. The SHG efficiency is defined as
Figure 6.(a), (b) Simulated transmission spectra of different quasi-BICs by sweeping
Generally, using a quasi-BIC to enhance SHG is based on symmetry-protected BICs [29,30,71–73]. This means that the conversion efficiency can be significantly improved only at a small incident angle, as for the quasi-BIC1 mentioned above. BICs merged at off- provide an excellent way to overcome the limitation of incident angle. When () are (1083.2 nm, 850 nm), the linewidth of the transmission spectrum vanishes at a certain incident angle as shown in Fig. 7(a). The simulated SHG with different incident angles in Figs. 7(b)–7(d) clearly shows the unique advantages of merging BICs at off-; even at a large oblique incidence angle, it can still maintain an efficient and stable SHG in a certain angle range. It is worth mentioning that no matter whether SHG is enhanced by merging BICs at point or off- point, the wavelength with maximum values of basically does not shift. As the angle increases by 8°, the relative shifts are all less than 0.3%. Hence, we can deduce that quasi-BICs originated from flatband merging BICs exhibit highly efficient performances in various applications requiring field enhancement.
Figure 7.(a) Simulated transmission spectra of merging BICs at off-
3. CONCLUSIONS
In summary, without breaking structural symmetry, we have proposed a feasible avenue to construct merging BICs at almost arbitrary points in momentum space. The proposed structure satisfies symmetry and supports the concurrence of symmetry-protected BICs and a few accidental BICs. The topological property of BICs ensures that accidental BICs can be moved in reciprocal space. By increasing the period of the PhC, the BIC near the point will converge to the center of the Brillouin zone, and finally merge with the symmetry-protected BIC at point. If the period decreases, these two accidental BICs carrying opposite topological charges approach each other and then form a merging BIC at off- point. We further showed that merging BICs can be constructed at almost arbitrary wave vectors by selecting appropriate heights and periods. Moreover, the band on which BICs are manipulated to realize merging BICs possesses quasi-flat dispersion. Thus, these flatband merging BICs eliminate the performance limitation of quasi-BICs caused by dispersive effects under wide-angle illuminations. Finally, compared with an isolated symmetry-protected BIC, the efficiency and robustness of SHG assisted by the designed flatband merging BICs can be greatly enhanced. When merging BICs appear at off- point, stable and efficient SHG at a larger incident angle is also proved. Our finding may have many potential applications based on local field enhancement, such as nonlinear, sensing, and quantum effects.
APPENDIX A: EVOLUTION OF BICs WITH THE VARIATION OF STRUCTURAL PARAMETERS
The evolution of merging BICs at with the variation of the periodic constant and height is shown in Fig.
Figure 8.Evolution of merging BICs at
Figure 9.Evolution of merging BICs at off-
APPENDIX B: ROBUSTNESS OF MERGING BICS AGAINST FABRICATION IMPERFECTIONS
The influence of three typical types of fabrication imperfections is considered in the factor. In the first type, the diameter of the cylinder has a deviation from the design value, as shown in Fig.
Figure 10.Influence of fabrication defect on
Figure 11.Influence of fabrication defect on
APPENDIX C: INTRINSIC LOSS OF MATERIALS
The loss of the system is mainly composed of intrinsic loss caused by materials and radiation loss, following . For merging BICs, the radiation loss is eliminated due to the destructive interference between radiation channels, which means that factors are dominated by the intrinsic loss of materials, as shown in Fig.
Figure 12.Influence of intrinsic loss on merging BICs. (a) Simulation for the evolution of
References
[4] S. Jahani, Z. Jacob. All-dielectric metamaterials. Nat. Nanotechnol., 11, 23-36(2016).
[9] A. F. Sadreev. Interference traps waves in an open system: bound states in the continuum. Rep. Prog. Phys., 84, 055901(2021).
Set citation alerts for the article
Please enter your email address