Hybrid bound states in the continuum in terahertz metasurfaces
Resonant cavities can effectively confine light and enhance light-matter interactions, which are of great importance to technologies and applications including lasers1-4, filters5, 6, harmonic generation7-11, and sensors12, 13. Quality factor (Q) is one of the most important parameters to characterize the strength of light-matter interactions. Different schemes have been proposed to improve Q in microcavities14, such as microdisks15, 16, Bragg reflector microcavities17, and photonic crystals18. In the same context, a generalized concept of bound states in the continuum (BIC) was raised19-21, which is in fact a topological defect in the momentum space and localizes in the continuous spectrum of extended states but is unable to couple to free space. Therefore, an ideal BIC will theoretically exhibit an infinite lifetime. Since the first demonstration, BICs have been applied to improve the performance in a plethora of optical applications by opening a coupling channel to free space via symmetry breaking of unit cells. Although any value of Q could be theoretically obtained, the measured values are commonly much lower than their theoretical predictions due to unavoidable scattering losses by fabrication defects, finite size of resonator array, and Ohmic losses in practice22-24. One solution to reduce the susceptibility of BIC to defect-induced scattering losses is by merging multiple BICs in the vicinity of Г point in the Brillouin zone (BZ) that in fact decreases the slope of radiative Q versus wavevector (k) from an inverse quadratic relationship (
Here, we introduce a generalized scheme to access robust and high-Q BICs by decreasing the radiation density in a hybrid metasurface lattice. For a common lattice supporting symmetry-protected BICs, the leaky channel is opened by uniformly breaking the symmetry of resonators in the unit cell, while a hybrid lattice indicates that half or a quarter (or less) of resonators in a supercell are symmetry breaking so that the radiation density reduces in the lattice. In the course of radiation suppression, the radiative Q versus k will theoretically have a 16-time larger coefficient in a 1/4 hybrid lattice compared to the common BICs that would enable a robust high-Q resonance. The idea was numerically and experimentally demonstrated by an array of classical double-gap split ring resonators (DSRR) with C2 symmetry in terahertz regime that supports symmetry-protected BICs28. Band folding analysis in BZ was performed to visualize the evolution of modes and quality factors between uniform and hybrid lattices. In the hybrid lattice, the high-Q portion of band was retained from the uniform lattice, while the low-Q portion was discarded in the band folding process resulting in a relatively high Q (more than 14.6 times higher than a uniform lattice in simulations) and robust BIC. Accompanying with the band folding, certain inaccessible modes in a uniform lattice become accessible at Г point, and multiple Fano resonances were captured in the hybrid lattice which could enable an alternative approach for broadband molecular fingerprint sensing12, 13.
2 Materials and Methods
Numerical simulations were carried out using commercially available software (COMSOL Multiphysics) with RF module of the finite-element frequency-domain solver. Periodic boundary conditions were employed for the unit cell, and perfectly matching layer (PML) was applied at the input and output ports. A nondispersive refractive index n = 3.45 was set for silicon without loss as substrate. For the calculations of eigenvalues and radiative quality factors, perfect electric conductor (PEC) was used for DSRR. In simulation of transmission spectra, DSRR was set as aluminum with a conductivity of 3.56×107 S·m−1.
Prior to fabrication, a 500 μm thick high-resistivity silicon (resistivity > 10000 Ω·cm) wafer was cleaned in an ultrasonic bath with acetone for 10 min and rinsed with isopropanol followed by baking on a hot plate at 120 °C for 180 s. Afterwards, a 2 μm RZJ 304.50 photoresist was spin-coated on the silicon at a speed of 5000 r/s for 30 s. The substrate with photoresist was then baked on a hot plate (100 °C, 180 s). Conventional UV photolithography (SUSS-MA6) was used to transfer DSRR pattern on photoresist, and then the sample was developed with a RZX3038 developer for ~ 30 s. The patterned sample was then baked on a hot plate (120 °C, 90 s). Finally, electron beam evaporation (TF500) was used to deposit 200 nm thick aluminum, and liftoff of the remaining photoresist was done in bath with acetone (60 °C, 30 min).
Transmission spectra were measured with a commercially available terahertz time-domain spectroscopic system29. After Fourier transform, we obtain transmission spectra of samples (ts) and references (bare silicon substrate, tr), and get normalized transmission spectra (t=ts/tr). With measured transmission spectra, the total Qtot
3 Results and discussion
We demonstrate the idea with array of classical DSRRs which could be generalized to other symmetry-protected BIC scenarios without prerequisite of multiple BICs in the band. In the symmetric scenario of DSRR, symmetry-protected BIC is uncoupled to free space exhibiting an infinite radiative Q as illustrated in
Hybrid BIC lattices. ( a– c) Schematic diagram of a symmetry-protected BIC lattice without radiation channel (a), a uniform quasi-BIC lattice with radiation channel open by breaking symmetry of all the resonators (b), and a hybrid quasi-BIC lattice with C2 symmetry preserved in the neighboring resonators along x-axis in a supercell (c). ( d) A double gap split ring resonator as the unit cell of the metasurface. ( e, f) Simulated (e) and experimental (f) transmission amplitude spectra for the three-type lattices. The same asymmetry degree (α = 4.95%) was applied for U-qBIC and Hx-BIC metasurfaces.
In the scenario of U-qBIC, all the symmetry-breaking resonators contribute to the far-field radiation that follows the inverse quadratic law in
The intuitive evolution of BIC from U-qBIC to Hx-BIC is unambiguously interpreted from the eigenvalue analysis in the reciprocal space. As shown in
Interpretation of hybrid BIC from reciprocal space. ( a, b) Brillouin zones of monoatomic and diatomic supercells at α = 0 when the periods were chosen with a and 2a in the x direction. ( c) Illustration of Brillouin zones for monoatomic and diatomic supercells showing the band folding operation where X and M points in the BZ of a monoatomic supercell are folded to X′ and M′ points in the BZ of a diatomic supercell, and Г(Г′) point is fixed. ( d) Band diagrams of monoatomic (black circles) and diatomic (orange lines) supercells showing the folding behavior where all the modes of a monoatomic supercell in unshaded region are reflected into the shaded region representing modes of a diatomic supercell. ( e) Radiative quality factors of monoatomic and diatomic supercells. The same folding behavior of Q is inherited from the eigenmodes. ( f) Comparison of radiative Q versus k between monoatomic and diatomic supercells. Circles are simulated values and solid lines are fitting curves with Eq. (2) whose coefficient is 4-times larger in a diatomic supercell than that of a monoatomic supercell. Here, px is period of supercell along x direction. Perfect electric conductor (PEC) was used for DSRRs in simulations to calculate eigenvalues and quality factors.
The mode evolution in the band folding process was verified with eigenmode analysis in simulations with DSRRs at
The above discussions could be numerically and experimentally demonstrated via far-field measurements by breaking the symmetry of the resonators (i.e.,
Experimental demonstration of the high-Q hybrid BIC. ( a, b) Microscopic images of U-qBIC and Hx-BIC metasurfaces. Supercells of U-qBIC and Hx-BIC metasurfaces are shown in the inset. Scale bar, 20 μm. ( c) Simulated transmission amplitude spectra of U-qBIC (left) and Hx-BIC (right) metasurfaces at an asymmetry degree of 2.97% with excitation electric field polarized along y-axis. Band diagrams of U-qBIC and Hx-BIC supercells are shown in the middle. ( d) Simulated radiative Q (circles) versus asymmetry degree (α) with inverse quadratic fitting curves (solid lines) of U-qBIC (black) and Hx-BIC (orange) supercells. Here, an additional constant of β is necessary to account for the nonuniform asymmetry in the hybrid lattice. ( e) Experimental transmission amplitude spectra of U-qBIC and Hx-BIC metasurfaces at an asymmetry degree of 7.42% with excitation electric field polarized along y-axis. The linewidth of Fano resonances is larger than that of simulations due to Ohmic loss in metallic resonators (aluminum) and finite number of supercells.
We fabricated the samples with conventional photolithography, and the microscopic images of samples are shown in
With the interpretation of band folding in the hybrid lattice, we could further expand the scheme to high order in the BZ by introducing a diagonal nonradiative resonators in the 2×2 supercells. In this scenario, two configurations, where one (
Generalized high-order hybrid BICs. ( a, b) Microscopic images of Ht-BIC and Hq-BIC metasurfaces with three and one asymmetric resonators out of four in a 2×2 supercell, respectively, and the period is 2a along both x and y axes. Scale bar, 20 μm. ( c) Schematic diagram of band folding from U-qBIC lattice (black) to Ht-BIC/Hq-BIC (red) in the Brillouin zone. ( d) Simulated transmission amplitude spectra of the Ht-BIC (left) and Hq-BIC (right) metasurfaces at an asymmetry degree of 2.97%. The band structure of Ht-BIC/Hq-BIC is shown in the middle, and the modes at the Г point marked with different colored circles are folded from X (red), Y (orange), and M (blue) points in the Brillouin zone of U-qBIC lattice, respectively. The highlighted resonances show the original modes inherited from U-qBIC lattice. ( e) Experimental (orange) and simulated (black) transmission amplitude spectra of Ht-BIC and Hq-BIC metasurfaces at an asymmetry degree of 7.42%. The overall linewidth of Fano resonances is larger than that of simulations due to Ohmic loss in metallic resonators (aluminum) and finite number of supercells.
Symmetry-protected BICs have found a plethora of important applications in lasing, nonlinear optics, terahertz generation, and biosensors. Significant improvement of light-matter interactions could be obtained with these BICs that have enabled lower lasing threshold, higher-efficiency generation of harmonics and terahertz radiations, and hyperspectral sensing. Traditionally, most common logic to access a quasi-BIC is to break the symmetry in the level of resonator itself. In this work, we have introduced a different approach to tailor the radiative losses from the level of the entire lattice. Via selectively preserving the C2 symmetry of resonators in the supercell in a hybrid BIC lattice, we can reduce the radiation density and thus effectively improve the overall quality factors. We numerically and experimentally investigated four-type typical hybrid lattices via investigating 0/4 (U-qBIC), 1/4 (Ht-BIC), 2/4 (Hx/y/d-BIC), and 3/4 (Hq-BIC) supercells, and found a progressive increase of the overall quality factors (
Significant Q improvement in hybrid BIC supercells and robustness against fabrication imperfections. ( a) Evolution of radiative Q versus AD for U-qBIC, Ht-BIC, Hx-BIC, and Hq-BIC supercells. The overall quality factors are improved in hybrid supercells with a lower radiation density. ( b) Influences of fabrication imperfection on quality factors in the four scenarios. Imperfection is introduced by adjusting the sharp right-angle of square in the resonators to rounded angles indicated by radius r. Q and Q' indicate radiative quality factors for lattices with right-angle and rounded-angle resonators, respectively.
A common fabrication imperfection in photolithography and E-beam lithography is the rounded angles in the square or rectangular resonators instead of the designed right angles. The radius (r) of rounded angles is usually impossible to be accurately predicted which will thus deteriorate the very sensitive high-Q resonances as shown in
In summary, we report a generalized approach to improve quality factors of symmetry-protected BICs with hybrid BIC supercells. A progressive improvement of overall quality factors is observed in the four-type hybrid supercells, and high-quality factors insensitive to asymmetry degree were realized in the Hq-BIC supercell which guarantees robust resonances with high-quality factors to fabrication defects and disorders. The underlying physics is interpreted from reciprocal space, which uncovers the origin of the stable and improved quality factors of hybrid BIC supercells from band folding process in the Brillouin zone. The consequent multiple Fano resonances at Г point accompanied with band folding provide an excellent solution for hyperspectral sensing and would be useful for high quality applications in optoelectronic devices. The hybrid BIC supercells exhibit great robustness to fabrication imperfections and disorders and would release the rigid requirements of fabrication accuracy, especially for applications requiring extremely precise quality factors.
This work was supported by the National Natural Science Foundation of China (Award No. 62175099), Guangdong Basic and Applied Basic Research Foundation (Award No. 2023A1515011085), Stable Support Program for Higher Education Institutions from Shenzhen Science, Technology & Innovation Commission (Award No. 20220815151149004), Global recruitment program of young experts of China, and startup funding of Southern University of Science and Technology. The authors acknowledge the assistance of SUSTech Core Research Facilities and thank Yao Wang for helpful discussions on fabrication.
L. Q. Cong initiated the idea and supervised the project. J. X. Fan fabricated samples, performed experiments and simulations, analyzed data, and initiated the manuscript. Z. L. Li performed partial experiments. All authors read and commented on the manuscript. L. Q. Cong, J. Q. Gu, and J. G. Han supervise the project.
The authors declare no competing financial interests.
 ZhaoAK, JiangN, PengJF, LiuSQ, ZhangYQet alParallel generation of low-correlation wideband complex chaotic signals using CW laser and external-cavity laser with self-phase-modulated injectionOpto-Electron Adv2022520002610.29026/oea.2022.200026
 LiuZJ, WangJY, ChenB, WeiYM, LiuWJ, et alGiant enhancement of continuous wave second harmonic generation from few-layer GaSe coupled to high-Q quasi bound states in the continuum Nano Lett2021217405741010.1021/acs.nanolett.1c01975
 ShiT, DengZL, GengGZ, ZengXZ, ZengYX, et alPlanar chiral metasurfaces with maximal and tunable chiroptical response driven by bound states in the continuumNat Commun202213411110.1038/s41467-022-31877-1
 MoiseevEI, KryzhanovskayaN, PolubavkinaYS, MaximovMV, KulaginaMM, et alLight outcoupling from quantum dot-based microdisk laser via plasmonic nanoantennaACS Photonics2017427528110.1021/acsphotonics.6b00552
Junxing Fan, Zuolong Li, Zhanqiang Xue, Hongyang Xing, Dan Lu, Guizhen Xu, Jianqiang Gu, Jiaguang Han, Longqing Cong. Hybrid bound states in the continuum in terahertz metasurfaces[J]. Opto-Electronic Science, 2023, 2(4): 230006.