Next Article in Journal
Thermal Behavior, Local-Scale Morphology, and Phase Composition of Spherulites in Melt-Crystallized Poly(Vinylidene Fluoride) Films
Previous Article in Journal
Transparent Celadon with Phase-Separated Structure: Study on the Technological Characteristics and Coloring Mechanism of Celadons from the Lieshan Kiln
Previous Article in Special Issue
Highly Sensitive SnS2/rGO-Based Gas Sensor for Detecting Chemical Warfare Agents at Room Temperature: A Theoretical Study Based on First-Principles Calculations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Q Resonances Induced by Toroidal Dipole Bound States in the Continuum in Terahertz Metasurfaces

Electronic Engineering College, Heilongjiang University, Harbin 150080, China
*
Author to whom correspondence should be addressed.
Crystals 2025, 15(1), 96; https://doi.org/10.3390/cryst15010096
Submission received: 3 January 2025 / Revised: 17 January 2025 / Accepted: 17 January 2025 / Published: 20 January 2025
(This article belongs to the Special Issue Organic Photonics: Organic Optical Functional Materials and Devices)

Abstract

:
The radiation mode of the interaction between electromagnetic waves and materials has always been a research hotspot in nanophotonics, and bound states in the continuum (BICs) belong to one of the nonradiative modes. Owing to their high-quality factor characteristics, BICs are extensively employed in nonlinear harmonic generators and sensors. Here, the influence of structural parameters on radiation modes has been systematically analyzed using band theory; the mechanisms of quasi-BIC mode and BIC mode were also analyzed through multipole decomposition of scattered power and near-field distribution. Notably, this study presents the discovery that the toroidal dipole-BIC (TD-BIC) arises from the interference and cancellation of electric and toroidal dipoles. The research results indicate that the structure, which supports symmetry-protected BICs, is sensitive to variations in the concentration of NaCl solution in its surroundings, making it applicable for liquid detection in miniaturized metal sensors. The proposed scheme broadens the applicability of BIC-based sensors and provides a prospective platform for biological and chemical sensing.

1. Introduction

Nonradiative states have been extensively studied in quantum mechanics, including nonradiative motion, extended electron models, radiation processes, and the Aharonov–Bohm phenomenon [1,2,3,4]. Nonradiative states, arising from oscillating charge–current distributions, are electromagnetic states that do not emit energy into the far-field region, instead confining fields to a localized space [5]. Recent advancements in metasurface photonics have led to a growing recognition of the significant potential of nonradiative optical modes for various applications. These include the enhancement of localized electric fields, control of far-field scattering in metasurfaces, and amplification of nonlinear effects [6]. Consequently, research on nonradiative electromagnetic modes has emerged as a prominent topic within the realms of classical electrodynamics and nanophotonics. In Mie-resonant dielectric structures, the coupling and complex interaction of multipole modes can give rise to nonradiative electromagnetic states [7]. These states are characterized by the absence of energy leakage and spatially localized profiles, which are known as BICs [8,9,10]. When a system has reflection or rotational symmetries, bound states and the continuum can belong to separate symmetry classes, which prohibits coupling if the symmetry remains intact [11,12,13,14,15]. When external disturbances affect the structure, such as the breaking of symmetry or material radiation and ohmic losses, symmetry-protected BICs will transition from dark modes to bright modes, becoming quasi-BICs. At this point, the energy confined within the structure is released, leading to Fano resonance with the external environment [16]. The characteristics of BICs enable their application in a wide range of technologies, including ultra-narrowband filters [12,17,18,19,20], laser action [21,22], highly sensitive sensing [23,24,25], enhanced nonlinear effects [26,27,28,29,30], and active devices [31,32,33,34,35].
Regardless of the context of nuclear, atomic, and molecular physics, solid state physics, classical electrodynamics, and photonics, the systems with toroidal topology structures need to consider the existence of toroidal multipole moments [36,37,38]. The existence of toroidal multipole moments was not proposed in Mie theory [38,39,40], but they have been proven to have significant effects in many applications, such as anapole states [41,42,43], lasing spasers [44,45], optical forces [46,47,48], and chirality [49]. Recently, several studies have been presented, including a theoretical analysis of TD-BICs [36], a polarization-insensitive TD-BIC implemented in a permittivity-asymmetric all-dielectric metasurface [50], a symmetry analysis of TD resonance in all-dielectric metasurfaces [51], and a study on high-Q factor-resonant dielectric metasurfaces [52]. Nevertheless, research on toroidal dipole BIC-driven high-Q resonances with a metallic terahertz metasurface is rare, especially when the BIC is mainly composed of interference cancellation between electric dipoles and toroidal dipoles; another mode that supports this special situation is the anapole mode.
In this study, we propose a terahertz metallic metasurface, which supports the BIC mode driven by a TD dipole. The metasurface structure consists of a periodically arranged aluminum double split-rectangular resonator (DSRR) array on a SiO2 substrate. We first demonstrate the BIC characteristics of resonances and choose the combination of an aluminum double split-rectangular as the asymmetric structure. To further explore the radiation modes of the structure, we conduct an in-depth examination of its dipole composition and orientation through multipole decomposition. Finally, perturbation theory is employed to theoretically analyze the sensing performance of the proposed structure. For example, our work proposes a sensor designed to detect the concentration of NaCl solution. Refractive index sensing reveals that the sensitivity of our structure reaches 335 GHz/RIU, higher than that of earlier reported metasurface THz sensors [53,54,55,56,57,58,59]. Our findings facilitate applications in optics and photonics, such as low-threshold lasers and sensing.

2. Theory and Model Design

The simulations are performed using the finite element method with COMSOL Multiphysics 6.2. A unit cell is simulated with periodic boundary conditions along the in-plane x- and y-axes and perfectly matched layers (PML) along the propagation direction (z-axis). Terahertz (THz) plane waves are incident perpendicularly on the metasurface along the z-axis, with the electric field linearly polarized along the x-axis. The electromagnetic and scattering properties of the proposed metasurface are analyzed by decomposing the fields within the metal elements into Cartesian multipole moments.

2.1. Metal Metasurface

The proposed metasurface structure consists of a periodically arranged aluminum double split-rectangular resonator (DSRR) array on a SiO2 substrate, as shown in Figure 1a. Figure 1b presents the top view (x-y plane) of the DSRR: the periodic constants in the x and y directions are Px = Py = 100 μm; the side lengths of the outer and inner rectangular shapes of the DSRR are a1 = 20 μm and a2 = 30 μm, respectively; the widths of the upper and lower gaps of the rectangular body in the DSRR are g2 = 4–10 μm and g1 = 4 μm, respectively; the distance between the two rectangular bodies is t = 6 μm and the height of the rectangular bodies is h = 8 μm; here, the asymmetry parameter α is defined as α = | g 2 g 1 | / ( g 2 + g 1 ) ; the DSRR is built on a SiO2 substrate with a thickness of 100 μm and a refractive index of n = 1.54 [60,61,62]. The DSRR is made of aluminum (Al) and the conductivity of Al is 3.56 × 107  S / m [23].

2.2. Theoretical Formula

To identify the radiation effects between electromagnetic waves and the material, as well as the corresponding dipole composition and orientation, we conducted a multipole expansion analysis of the DSRR based on the following equations. The formulas for multipole decomposition are derived from the current. We first need to calculate the current based on the electric field at the corresponding position [40,63]:
j = i ω [ ε r 1 ] E ( r )
Here, j is the current vector, ω is the frequency of the incident wave, ε is the dielectric constant of the material, and E(r) is the electric field at the corresponding position. Then, the electric dipole moment ED, magnetic dipole moment MD, toroidal dipole moment TD, electric quadrupole moment EQ, and magnetic quadrupole moment MQ can be derived as follows [40,63]:
E D = 1 i ω V j ( r ) d r
M D = 1 2 c V r × j ( r ) d r
T D = 1 10 c V { [ r j ( r ) ] r 2 ( r · r ) j ( r ) } d r
E Q α , β = 1 2 i ω V [ r α j β ( r ) + r β j α ( r ) 2 3 ( r j ( r ) ) δ α , β ] d r
M Q α , β = 1 3 c V { [ r × j ( r ) ] α r β + r α [ r × j ( r ) ] β } d r
Finally, by using the following equation, the corresponding scattering cross-section C s c a s u m can be obtained [64]:
C s c a s u m k 0 4 6 π ε 0 2 E inc 2 E D + i k 0 T D 2 + k 0 4 ε d μ 0 6 π ε 0 E inc 2 M D 2 + k 0 6 ε d 720 π ε 0 2 E inc 2 E Q 2 + k 0 6 ε d 2 μ 0 80 π ε 0 E inc 2 M Q 2
k0 is the wave number in vacuum; E inc is the magnitude of the incident electric field, which is 1 V/m; μ 0 is the magnetic permeability in vacuum; ε d is the dielectric constant of the corresponding material.

3. Results and Discussion

Figure 2a depicts the simulated spectra of the DSRR composed of a perfect electric conductor (PEC) and a realistic metal (aluminum) with different g2; the simulation of the PEC represents the ideal simulation of a metal with no loss in THz waves. As shown in Table 1 and Figure 2a, when g2 is varied from 4 to 10 μm, it can be observed that the Fano resonance gradually strengthens within the structure, and the transmission curve shows a typical Fano resonance curve, and the absorption of the DSRR gradually increases; when g2 = 4 μm, the Fano resonance half-peak full width (FWHM) of the DSRR structure at 2.11–2.15 THz disappears, and the corresponding absorption also tends to a smooth curve, it means a transform in radiation dipole modes in the DSRR. The mode at 2.1 THz is a weakly coupled mode of Fano resonance [8]. Figure 2b shows the two-dimensional transmission spectrum in different g2, and it can be easily observed that Fano resonance FWHM gradually decreases and disappears at 2.112 THz (marked by the green star).
In order to further investigate the properties of BIC mode, a perfect electric conductor (PEC) was used for the DSRR in simulations to calculate characteristic frequency and quality factors. The Q-factor is defined as Q = ω0/2γ, where ω0 represents the resonant characteristic frequency and γ represents the resonance damping [12]. Figure 2c shows the characteristic frequency and radiation Q-factor under the influence of the periodic boundary conditions of the structure, where a represents the period size of the structure. The position of the wave vector is from the M-point to the highly symmetric center point Γ-point and then to the X-point, which mainly indicates whether the disturbance of periodic boundary conditions affects the mirror symmetry and rotational symmetry (C2v symmetry) of the structure. From Figure 2c, the characteristic frequency at the highly symmetrical center point is 2.112 THz, and the radiation Q-factor at this position is at its maximum and close to infinity. When gradually moving away from the highly symmetrical center point, the radiation Q factor sharply decreases. This indicates that the DSRR maintains mirror symmetry and rotational symmetry at the highly symmetrical center Γ-point. Combining Figure 2a,b, it can be observed that when g2 = 4 μm, the structure maintains C2V symmetry. When the symmetry of the structure is disrupted, either by moving away from the highly symmetrical center Γ-point or when the structural parameter g2 exceeds 4 μm, the energy originally confined within the structure is gradually released. It leads to mutual coupling between the structure and the extended wave, resulting in increasingly pronounced Fano resonance. Figure 2d presents a logarithmic plot of the radiation Q-factor as a function of asymmetry α, illustrating the relationship between asymmetry α and the radiation Q-factor as Q α 2 . According to the definition of symmetric protected BIC, two conditions must be met: first, when the structure maintains C2v symmetry, it has an infinitely high Q-factor in an infinitely narrow frequency spectrum [6,8]; second, the relationship between the asymmetric parameter α of the structure and the Q-factor is Q α 2 [11,12]. From Figure 2a–d, it can be inferred that the mode at 2.112 THz is a symmetric protected BIC.
Figure 2 clearly demonstrates that when the structure is in a BIC mode, electromagnetic waves within the structure do not couple with extended waves; the electromagnetic wave absorption within the structure tends toward 0. In contrast, in a quasi-BIC mode, these waves do couple, accompanied by a noticeable Fano resonance, and the electromagnetic wave absorption within the structure gradually increases with the increase of asymmetry. This indicates that in structures with mirror or inversion symmetry, BICs prevent coupling between different modes.
To investigate the reason for the increase in absorption within the DSRR, we performed a multipole decomposition in the Cartesian coordinate system. For specific calculations, please refer to the theoretical discussion mentioned earlier. Figure 3a depicts the total scattering power in different g2. It is worth noting that the minimum total scattering rate is at the BIC mode position (at 2.112 THz), while the maximum total scattering rate is at the quasi-BIC position. Figure 3b illustrates the composition of dipoles at the positions of the extremum, corresponding to Figure 3a. For BICs, it can be observed clearly that the contribution of the ED (at 2.112 THz) is dominant, indicating that BICs are primarily induced by the ED. For quasi-BICs, the scattering is mainly contributed by the superposition of other dipole powers. By decomposing the x and y components of the ED and TD scattering powers in Figure 3a, we find that the x component of the scattering power from the ED and TD predominates and their values are nearly equal. The y and z components of scattering power from the ED and TD approach zero when compared to the x component. This indicates that the ED and TD are aligned along the x direction when the metasurface is illuminated by x-polarized light. For the BIC mode, Figure 3a shows that the total scattering power at 2.112 THz reaches a minimum value, indicating the presence of interference and the cancellation behavior of dipoles. On the contrary, for the quasi-BIC mode, the total scattering power is at a maximum, indicating the presence of dipole superposition behavior. To verify the proposed interference cancellation and superposition mechanisms, we plotted the scattering patterns between EDx and TDx, as well as the scattering power maps of ED-ik0TD at g2 = 4 μm and g2 = 10 μm, respectively, as shown in Figure 3c,d. The blue and green lines represent EDx and −ik0TDx, respectively. When these two lines intersect, it indicates that the angular difference between their directions is π, reflecting opposite orientations and the cancellation of dipole interference. When these two lines do not intersect and differ by an angle of π, it indicates that the directional difference between them is 0, meaning they are aligned in the same direction, which leads to dipole superposition behavior. The red dashed line represents the total radiated power between EDx and TDx. When there is a minimum value of the scattering power, it indicates the interference and cancellation behavior between the EDx and TDx. Conversely, when there is a maximum value of the scattering power, it signifies the superposition behavior of EDx and TDx. This observation is consistent with the solid blue and green lines. For the BIC mode (at g2 = 4 μm), the blue and green lines intersect at the 2.112 THz, while the red dashed line shows a minimum value at 2.112 THz (see Figure 3c). This clearly indicates that the direction of EDx and TDx are opposite, demonstrating that ED and TD undergo interference cancellation behavior. Thus, the generation of the BIC mode at this location arises from the interference cancellation behavior of the ED and TD. For the quasi-BIC mode (at g2 = 10 μm), the difference in value between the blue and green lines at 2.242 THz is π, while the red dashed line exhibits a maximum value at 2.242 THz (see Figure 3d). This clearly indicates that the direction of EDx and TDx are aligned, demonstrating that the ED and TD undergo superposition behavior. Thus, the generation of the quasi-BIC mode at this location arises from the superposition behavior of the ED and TD.
To further verify the directions of the ED and TD in the BIC and quasi-BIC modes, we plotted the electromagnetic field diagrams for g2 = 4 μm and g2 = 10 μm, as shown in Figure 4a–d. The direction of the ED is defined as the path from the negative charge to the positive charge, indicating that the opposite direction of the electric field vector corresponds to the direction of the ED. The generation of the TD arises from the MD, which in turn is derived from the surface current. The directions of the TD and MD can be determined using the right-hand screw rule. According to Figure 4a,c, the electric field vector of the DSRR is oriented from left to right in the central gap, while the direction of the ED is from right to left. In the BIC mode, the dipole direction at the upper and lower open gaps is from right to left; in contrast, for the quasi-BIC mode, the dipole direction at these gaps is from left to right. As shown in Figure 4b, the surface current in the upper cube flows in a left-handed circular direction, indicating that the direction of the MD is oriented upward. In the lower cube, the surface current flows in a right-handed circular direction, resulting in a downward orientation of the MD. Therefore, the direction of the TD is to the right. In contrast, as shown in Figure 4c, the direction of the TD is oriented to the left. Figure 4a–d confirm the observations presented in Figure 3c,d, indicating that the BIC mode induced by the DSRR structure results from destructive interference between the ED and TD, while the quasi-BIC mode arises from the constructive interference between dipoles.
The increase in the scattering power of the TD significantly enhances the absorption of the structure, which is also reflected in the transmission dips observed in Figure 2a,b. In the BIC mode, due to the destructive interference between the ED and TD, the interaction between the THz electromagnetic waves and the structure does not exhibit a pronounced Fano resonance, with the primary contribution to radiation arising from the ED mode. On the other hand, in the quasi-BIC mode, the ED and TD are in a state of constructive interference, and the relatively strong scattering power of the TD significantly influences the interaction between the THz electromagnetic waves and the structure. This results in a pronounced Fano resonance, with the primary contributions to radiation arising from the ED, TD, and MQ modes.
The Fano resonance characteristics of metasurfaces are influenced by the surrounding background medium, making them promising for liquid sensing applications [65]. We designed a metasurface liquid sensor featuring Al cuboid dimer clusters, a SiO2 substrate, and an integrated microfluidic channel, as shown in Figure 5a. The sensitivity (S) is defined as the resonant frequency shift per unit refractive index change, using the formula S = Δ f / Δ n [66], where Δ f represents the resonance peak frequency shift and Δ n is the refractive index change of the analyte liquid. The figure of merit (FOM) is defined as FOM = S/FWHM, where FWHM is the full width at half maximum. As shown in Figure 5b,c, increasing the refractive index of liquid from 1.3333 to 1.351 results in a redshift of the resonant frequency. When g2 is 12, 10, and 8 μm, the maximum absorption rates of the DSRR in quasi-BIC mode are 47%, 44%, and 19%, respectively. The overall structure of the DSRR is immersed in water at room temperature and pressure, and the thermal radiation generated by THz electromagnetic waves can be ignored [67]. The refractive index corresponding to NaCl solutions at concentrations of 0%, 1%, 5%, and 10% are 1.3333, 1.3352, 1.3424, and 1.351, respectively [67]. When g2 = 12 μm, the Fano resonance dips across different background media occur between 2.06 and 2.07 THz, with an average FWHM of 64.53 RIU−1. When g2 = 10 μm, the Fano resonance dips occur between 2.054 and 2.064 THz, with an average FWHM of 99.04 RIU−1. When g2 = 8 μm, the Fano resonance dips occur between 2.047 and 2.057 THz, with an average FWHM of 159.52 RIU−1. As shown in Figure 5c,d, when g2 = 12, 10, 8 μm, the average values of refractive index sensitivity are S = 332 ± 3, 272 ± 2, and 230 ± 2 GHz/RIU−1, respectively. Notably, when g2 = 12 μm, the refractive index sensitivity better linearity compared to other g2. Additionally, we compared our work with others, showing better performance (as shown in Table 2). This enhanced linearity ensures more accurate and consistent detection, making the metasurface sensor highly suitable for liquid sensing applications, where precise measurement of refractive index changes is crucial.

4. Conclusions

In summary, a THz metal liquid sensor capable of generating toroidal dipole BIC-driven has been proposed, and the influence of structural parameters on radiation modes has been systematically analyzed using band theory through the relationship between structural characteristic frequency and radiation Q-factor. The mechanisms of the quasi-BIC mode and BIC mode were also analyzed through multipole decomposition of scattered power and near-field distribution. With sensitivities of 335 GHz/RIU and an FOM of 64.53 RIU⁻¹, the quasi-BIC resonant sensor outperforms previously reported designs, demonstrating promising liquid sensing capabilities for precision sensing applications.

Author Contributions

Conceptualization, L.G. and Y.G.; methodology, L.G. and Y.G.; software, L.G. and Y.G.; validation, L.G.; formal analysis, L.G.; investigation, L.G.; writing—original draft preparation, L.G.; writing—review and editing, L.G. and Y.G.; funding acquisition, Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in this study are included in the article. Further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bohm, D.; Weinstein, M. The Self-Oscillations of a Charged Particle. Phys. Rev. 1948, 74, 1789–1798. [Google Scholar] [CrossRef]
  2. Davidson, M.P. Quantum Wave Equations and Non-Radiating Electromagnetic Sources. Ann. Phys. 2007, 322, 2195–2210. [Google Scholar] [CrossRef]
  3. Goedecke, G.H. Classically Radiationless Motions and Possible Implications for Quantum Theory. Phys. Rev. 1964, 135, B281–B288. [Google Scholar] [CrossRef]
  4. Nemkov, N.A.; Basharin, A.A.; Fedotov, V.A. Nonradiating Sources, Dynamic Anapole, and Aharonov-Bohm Effect. Phys. Rev. B 2017, 95, 165134. [Google Scholar] [CrossRef]
  5. Koshelev, K.; Favraud, G.; Bogdanov, A.; Kivshar, Y.; Fratalocchi, A. Nonradiating Photonics with Resonant Dielectric Nanostructures. Nanophotonics 2019, 8, 725–745. [Google Scholar] [CrossRef]
  6. Baryshnikova, K.V.; Smirnova, D.A.; Luk’yanchuk, B.S.; Kivshar, Y.S. Optical Anapoles: Concepts and Applications. Adv. Opt. Mater. 2019, 7, 1801350. [Google Scholar] [CrossRef]
  7. Azzam, S.I.; Kildishev, A.V. Photonic Bound States in the Continuum: From Basics to Applications. Adv. Opt. Mater. 2021, 9, 2001469. [Google Scholar] [CrossRef]
  8. Limonov, M.F. Fano Resonances in Photonics. Nat. Photonics 2017, 11, 543–554. [Google Scholar] [CrossRef]
  9. Hsu, C.W.; Zhen, B.; Stone, A.D.; Joannopoulos, J.D.; Soljačić, M. Bound States in the Continuum. Nat. Rev. Mater. 2016, 1, 16048. [Google Scholar] [CrossRef]
  10. Friedrich, H.; Wintgen, D. Interfering Resonances and Bound States in the Continuum. Phys. Rev. A 1985, 32, 3231–3242. [Google Scholar] [CrossRef]
  11. Li, S.; Zhou, C.; Liu, T.; Xiao, S. Symmetry-Protected Bound States in the Continuum Supported by All-Dielectric Metasurfaces. Phys. Rev. A 2019, 100, 063803. [Google Scholar] [CrossRef]
  12. Koshelev, K.; Lepeshov, S.; Liu, M.; Bogdanov, A.; Kivshar, Y. Asymmetric Metasurfaces with High-Q Resonances Governed by Bound States in the Continuum. Phys. Rev. Lett. 2018, 121, 193903. [Google Scholar] [CrossRef] [PubMed]
  13. Overvig, A.C.; Malek, S.C.; Carter, M.J.; Shrestha, S.; Yu, N. Selection Rules for Quasibound States in the Continuum. Phys. Rev. B 2020, 102, 035434. [Google Scholar] [CrossRef]
  14. Sadrieva, Z.; Frizyuk, K.; Petrov, M.; Kivshar, Y.; Bogdanov, A. Multipolar Origin of Bound States in the Continuum. Phys. Rev. B 2019, 100, 115303. [Google Scholar] [CrossRef]
  15. Sadrieva, Z.F.; Belyakov, M.A.; Balezin, M.A.; Kapitanova, P.V.; Nenasheva, E.A.; Sadreev, A.F.; Bogdanov, A.A. Experimental Observation of a Symmetry-Protected Bound State in the Continuum in a Chain of Dielectric Disks. Phys. Rev. A 2019, 99, 53804. [Google Scholar] [CrossRef]
  16. Wang, Y.; Han, Z.; Du, Y.; Qin, J. Ultrasensitive Terahertz Sensing with High-Q Toroidal Dipole Resonance Governed by Bound States in the Continuum in All-Dielectric Metasurface. Nanophotonics 2021, 10, 1295–1307. [Google Scholar] [CrossRef]
  17. Kawanishi, K.; Shimatani, A.; Lee, K.J.; Inoue, J.; Ura, S.; Magnusson, R. Cross-Stacking of Guided-Mode Resonance Gratings for Polarization-Independent Flat-Top Filtering. In Proceedings of the 2019 24th Microoptics Conference (MOC), Toyama, Japan, 17–20 November 2019. [Google Scholar] [CrossRef]
  18. Liu, W.; Li, Y.; Jiang, H.; Lai, Z.; Chen, H. Controlling the Spectral Width in Compound Waveguide Grating Structures. Opt. Lett. 2013, 38, 163. [Google Scholar] [CrossRef]
  19. Wang, X.; Li, S.; Zhou, C. Polarization-Independent Toroidal Dipole Resonances Driven by Symmetry-Protected BIC in Ultraviolet Region. Opt. Express 2020, 28, 11983. [Google Scholar] [CrossRef]
  20. Ramesh, R.; Qing, N.; Alshehri, H.; Azeredo, B.; Liping, W. Design of Selective Metasurface Filter for Thermophotovoltaic Energy Conversion. ES Energy Environ. 2023, 22, 999. [Google Scholar] [CrossRef]
  21. Wu, M.; Ha, S.T.; Shendre, S.; Durmusoglu, E.G.; Koh, W.-K.; Abujetas, D.R.; Sánchez-Gil, J.A.; Paniagua-Domínguez, R.; Demir, H.V.; Kuznetsov, A.I. Room-Temperature Lasing in Colloidal Nanoplatelets via Mie-Resonant Bound States in the Continuum. Nano Lett. 2020, 20, 6005–6011. [Google Scholar] [CrossRef]
  22. Zhang, H.; Wang, T.; Tian, J.; Sun, J.; Li, S.; De Leon, I.; Zaccaria, R.P.; Peng, L.; Gao, F.; Lin, X.; et al. Quasi-BIC Laser Enabled by High-Contrast Grating Resonator for Gas Detection. Nanophotonics 2022, 11, 297–304. [Google Scholar] [CrossRef] [PubMed]
  23. Yang, F.; Li, J.; Wu, L.; Zheng, C.; Yue, Z.; Li, H.; Song, C.; Ding, X.; Zhang, Y.; Yao, J. Bending Sensing Based on Quasi Bound States in the Continuum in Flexible Terahertz Metasurface. Adv. Opt. Mater. 2023, 11, 2300909. [Google Scholar] [CrossRef]
  24. Wu, F.; Wu, J.; Guo, Z.; Jiang, H.; Sun, Y.; Li, Y.; Ren, J.; Chen, H. Giant Enhancement of the Goos-Hänchen Shift Assisted by Quasibound States in the Continuum. Phys. Rev. Appl. 2019, 12, 14028. [Google Scholar] [CrossRef]
  25. Srivastava, Y.K.; Ako, R.T.; Gupta, M.; Bhaskaran, M.; Sriram, S.; Singh, R. Terahertz Sensing of 7 Nm Dielectric Film with Bound States in the Continuum Metasurfaces. Appl. Phys. Lett. 2019, 115, 151105. [Google Scholar] [CrossRef]
  26. Koshelev, K.; Tang, Y.; Li, K.; Choi, D.-Y.; Li, G.; Kivshar, Y. Nonlinear Metasurfaces Governed by Bound States in the Continuum. ACS Photonics 2019, 6, 1639–1644. [Google Scholar] [CrossRef]
  27. Liu, Z.; Xu, Y.; Lin, Y.; Xiang, J.; Feng, T.; Cao, Q.; Li, J.; Lan, S.; Liu, J. High- Q Quasibound States in the Continuum for Nonlinear Metasurfaces. Phys. Rev. Lett. 2019, 123, 253901. [Google Scholar] [CrossRef]
  28. Li, N.; Xu, Z.; Dong, Y.; Hu, T.; Zhong, Q.; Fu, Y.H.; Zhu, S.; Singh, N. Large-Area Metasurface on CMOS-Compatible Fabrication Platform: Driving Flat Optics from Lab to Fab. Nanophotonics 2020, 9, 3071–3087. [Google Scholar] [CrossRef]
  29. Mikheeva, E.; Koshelev, K.; Choi, D.-Y.; Kruk, S.; Lumeau, J.; Abdeddaim, R.; Voznyuk, I.; Enoch, S.; Kivshar, Y. Photosensitive Chalcogenide Metasurfaces Supporting Bound States in the Continuum. Opt. Express 2019, 27, 33847. [Google Scholar] [CrossRef]
  30. Chen, W.; Yongzheng, W.; Jingbo, S.; Ji, Z. Recent Progress on Optical Frequency Conversion in Nonlinear Metasurfaces and Nanophotonics. ES Mater. Manuf. 2022, 17, 1–13. [Google Scholar] [CrossRef]
  31. Han, S.; Cong, L.; Srivastava, Y.K.; Qiang, B.; Rybin, M.V.; Kumar, A.; Jain, R.; Lim, W.X.; Achanta, V.G.; Prabhu, S.S.; et al. All-Dielectric Active Terahertz Photonics Driven by Bound States in the Continuum. Adv. Mater. 2019, 31, 1901921. [Google Scholar] [CrossRef]
  32. Karl, N.; Vabishchevich, P.P.; Liu, S.; Sinclair, M.B.; Keeler, G.A.; Peake, G.M.; Brener, I. All-Optical Tuning of Symmetry Protected Quasi Bound States in the Continuum. Appl. Phys. Lett. 2019, 115, 141103. [Google Scholar] [CrossRef]
  33. Li, J.; Li, J.; Zheng, C.; Yue, Z.; Wang, S.; Li, M.; Zhao, H.; Zhang, Y.; Yao, J. Free Switch between Bound States in the Continuum (BIC) and Quasi-BIC Supported by Graphene-Metal Terahertz Metasurfaces. Carbon 2021, 182, 506–515. [Google Scholar] [CrossRef]
  34. Tan, T.C.; Srivastava, Y.K.; Ako, R.T.; Wang, W.; Bhaskaran, M.; Sriram, S.; Al-Naib, I.; Plum, E.; Singh, R. Active Control of Nanodielectric-Induced THz Quasi-BIC in Flexible Metasurfaces: A Platform for Modulation and Sensing. Adv. Mater. 2021, 33, 2100836. [Google Scholar] [CrossRef]
  35. Ma, W.; Zhou, C.; Chen, D.; You, S.; Wang, X.; Wang, L.; Jin, L.; Huang, L.; Wang, D.; Miroshnichenko, A.E. Active Quasi-BIC Metasurfaces Assisted by Epsilon-near-Zero Materials. Opt. Express 2023, 31, 13125. [Google Scholar] [CrossRef]
  36. He, Y.; Guo, G.; Feng, T.; Xu, Y.; Miroshnichenko, A.E. Toroidal Dipole Bound States in the Continuum. Phys. Rev. B 2018, 98, 161112. [Google Scholar] [CrossRef]
  37. Papasimakis, N.; Fedotov, V.A.; Savinov, V.; Raybould, T.A.; Zheludev, N.I. Electromagnetic Toroidal Excitations in Matter and Free Space. Nat. Mater. 2016, 15, 263–271. [Google Scholar] [CrossRef]
  38. Radescu, E.E.; Vaman, G. Exact Calculation of the Angular Momentum Loss, Recoil Force, and Radiation Intensity for an Arbitrary Source in Terms of Electric, Magnetic, and Toroid Multipoles. Phys. Rev. E 2002, 65, 046609. [Google Scholar] [CrossRef]
  39. Evlyukhin, A.B.; Chichkov, B.N. Multipole Decompositions for Directional Light Scattering. Phys. Rev. B 2019, 100, 125415. [Google Scholar] [CrossRef]
  40. Evlyukhin, A.B.; Fischer, T.; Reinhardt, C.; Chichkov, B.N. Optical Theorem and Multipole Scattering of Light by Arbitrarily Shaped Nanoparticles. Phys. Rev. B 2016, 94, 205434. [Google Scholar] [CrossRef]
  41. Savinov, V.; Papasimakis, N.; Tsai, D.P.; Zheludev, N.I. Optical Anapoles. Commun. Phys. 2019, 2, 69. [Google Scholar] [CrossRef]
  42. Díaz-Escobar, E.; Bauer, T.; Pinilla-Cienfuegos, E.; Barreda, Á.I.; Griol, A.; Kuipers, L.; Martínez, A. Radiationless Anapole States in On-Chip Photonics. Light Sci. Appl. 2021, 10, 204. [Google Scholar] [CrossRef] [PubMed]
  43. Miroshnichenko, A.E.; Evlyukhin, A.B.; Yu, Y.F.; Bakker, R.M.; Chipouline, A.; Kuznetsov, A.I.; Luk’yanchuk, B.; Chichkov, B.N.; Kivshar, Y.S. Nonradiating Anapole Modes in Dielectric Nanoparticles. Nat. Commun. 2015, 6, 8069. [Google Scholar] [CrossRef] [PubMed]
  44. Huang, Y.-W.; Chen, W.T.; Wu, P.C.; Fedotov, V.A.; Zheludev, N.I.; Tsai, D.P. Toroidal Lasing Spaser. Sci. Rep. 2013, 3, 1237. [Google Scholar] [CrossRef] [PubMed]
  45. Zdagkas, A.; McDonnell, C.; Deng, J.; Shen, Y.; Li, G.; Ellenbogen, T.; Papasimakis, N.; Zheludev, N.I. Observation of Toroidal Pulses of Light. Nat. Photon. 2022, 16, 523–528. [Google Scholar] [CrossRef]
  46. Hong, I.; Hong, C.; Tutanov, O.S.; Massick, C.; Castleberry, M.; Zhang, Q.; Jeppesen, D.K.; Higginbotham, J.N.; Franklin, J.L.; Vickers, K.; et al. Anapole-Assisted Low-Power Optical Trapping of Nanoscale Extracellular Vesicles and Particles. Nano Lett. 2023, 23, 7500–7507. [Google Scholar] [CrossRef]
  47. Hong, I.; Anyika, T.; Hong, C.; Yang, S.; Ndukaife, J.C. Hybrid Optical and Diffusiophoretic Nanomanipulation Using All-Dielectric Anapole-Enhanced Thermonanophotonics. ACS Photonics 2023, 10, 4038–4044. [Google Scholar] [CrossRef]
  48. Conteduca, D.; Brunetti, G.; Barth, I.; Quinn, S.D.; Ciminelli, C.; Krauss, T.F. Multiplexed Near-Field Optical Trapping Exploiting Anapole States. ACS Nano 2023, 17, 16695–16702. [Google Scholar] [CrossRef]
  49. Raybould, T.A.; Fedotov, V.A.; Papasimakis, N.; Kuprov, I.; Youngs, I.J.; Chen, W.T.; Tsai, D.P.; Zheludev, N.I. Toroidal Circular Dichroism. Phys. Rev. B 2016, 94, 35119. [Google Scholar] [CrossRef]
  50. Yu, S.; Wang, Y.; Gao, Z.; Li, H.; Song, S.; Yu, J.; Zhao, T. Dual-Band Polarization-Insensitive Toroidal Dipole Quasi-Bound States in the Continuum in a Permittivity-Asymmetric All-Dielectric Meta-Surface. Opt. Express 2022, 30, 4084. [Google Scholar] [CrossRef]
  51. Dmitriev, V.; Kupriianov, A.S.; Silva Santos, S.D.; Tuz, V.R. Symmetry Analysis of Trimer-Based All-Dielectric Metasurfaces with Toroidal Dipole Modes. J. Phys. D Appl. Phys. 2021, 54, 115107. [Google Scholar] [CrossRef]
  52. Algorri, J.F.; Zografopoulos, D.C.; Ferraro, A.; García-Cámara, B.; Beccherelli, R.; Sánchez-Pena, J.M. Ultrahigh-Quality Factor Resonant Dielectric Metasurfaces Based on Hollow Nanocuboids. Opt. Express 2019, 27, 6320. [Google Scholar] [CrossRef] [PubMed]
  53. Ma, T.; Huang, Q.; He, H.; Zhao, Y.; Lin, X.; Lu, Y. All-Dielectric Metamaterial Analogue of Electromagnetically Induced Transparency and Its Sensing Application in Terahertz Range. Opt. Express 2019, 27, 16624. [Google Scholar] [CrossRef] [PubMed]
  54. Xie, Q.; Dong, G.-X.; Wang, B.-X.; Huang, W.-Q. High-Q Fano Resonance in Terahertz Frequency Based on an Asymmetric Metamaterial Resonator. Nanoscale Res. Lett. 2018, 13, 294. [Google Scholar] [CrossRef]
  55. Yahiaoui, R.; Tan, S.; Cong, L.; Singh, R.; Yan, F.; Zhang, W. Multispectral Terahertz Sensing with Highly Flexible Ultrathin Metamaterial Absorber. J. Appl. Phys. 2015, 118, 83103. [Google Scholar] [CrossRef]
  56. Gupta, M.; Srivastava, Y.K.; Manjappa, M.; Singh, R. Sensing with Toroidal Metamaterial. Appl. Phys. Lett. 2017, 110, 121108. [Google Scholar] [CrossRef]
  57. Xu, J.; Liao, D.; Gupta, M.; Zhu, Y.; Zhuang, S.; Singh, R.; Chen, L. Terahertz Microfluidic Sensing with Dual-torus Toroidal Metasurfaces. Adv. Opt. Mater. 2021, 9, 2100024. [Google Scholar] [CrossRef]
  58. He, X.; Zhang, Q.; Lu, G.; Ying, G.; Wu, F.; Jiang, J. Tunable Ultrasensitive Terahertz Sensor Based on Complementary Graphene Metamaterials. RSC Adv. 2016, 6, 52212–52218. [Google Scholar] [CrossRef]
  59. Zhong, Y.; Du, L.; Liu, Q.; Zhu, L.; Meng, K.; Zou, Y.; Zhang, B. Ultrasensitive Specific Sensor Based on All-Dielectric Metasurfaces in the Terahertz Range. RSC Adv. 2020, 10, 33018–33025. [Google Scholar] [CrossRef]
  60. Beruete, M.; Jáuregui-López, I. Terahertz Sensing Based on Metasurfaces. Adv. Opt. Mater. 2020, 8, 1900721. [Google Scholar] [CrossRef]
  61. Kim, J.S.; Kwon, Y.S.; Ahn, S.; Choi, B.K.; Cho, S.Y.; Jeon, M.Y. Measurement of Refractive Indices and Absorption Coefficients for Glass Materials and Nematic Liquid Crystals in THz Frequency Band. J. Korean Phys. Soc. 2024, 84, 750–757. [Google Scholar] [CrossRef]
  62. Palik, E.D.; Ghosh, G. (Eds.) Handbook of Optical Constants of Solids; Academic Press: San Diego, CA, USA, 1998; ISBN 978-0-12-544420-0. [Google Scholar]
  63. Savinov, V.; Fedotov, V.A.; Zheludev, N.I. Toroidal Dipolar Excitation and Macroscopic Electromagnetic Properties of Metamaterials. Phys. Rev. B 2014, 89, 205112. [Google Scholar] [CrossRef]
  64. Evlyukhin, A.B.; Reinhardt, C.; Evlyukhin, E.; Chichkov, B.N. Multipole Analysis of Light Scattering by Arbitrary-Shaped Nanoparticles on a Plane Surface. J. Opt. Soc. Am. B 2013, 30, 2589. [Google Scholar] [CrossRef]
  65. Algorri, J.F.; Zografopoulos, D.C.; Ferraro, A.; García-Cámara, B.; Vergaz, R.; Beccherelli, R.; Sánchez-Pena, J.M. Anapole Modes in Hollow Nanocuboid Dielectric Metasurfaces for Refractometric Sensing. Nanomaterials 2018, 9, 30. [Google Scholar] [CrossRef]
  66. Li, H.; Yu, S.; Yang, L.; Zhao, T. High Q-Factor Multi-Fano Resonances in All-Dielectric Double Square Hollow Metamaterials. Opt. Laser Technol. 2021, 140, 107072. [Google Scholar] [CrossRef]
  67. Li, S.; Jiang, H.; Zhu, X.; Shi, Y.; Han, Z. A High-Sensitivity Refractive Index Sensor with Period-Doubling Plasmonic Metasurfaces to Engineer the Radiation Losses. ACS Appl. Opt. Mater. 2023, 1, 736–744. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic of the DSRR under the THz wave illumination. (b) Detailed diagram of the unit structure of metamaterials.
Figure 1. (a) Schematic of the DSRR under the THz wave illumination. (b) Detailed diagram of the unit structure of metamaterials.
Crystals 15 00096 g001
Figure 2. (a) Simulated spectra of the ideal DSRR composed of PEC and realistic metallic (Al) with varying structural parameters g2. (b) Simulated transmittance spectra map of the DSRR for various g2. (c) The band structure and radiative Q-factor of the DSRR structure at gap sizes of 4 μm. (d) The curve shows the change in the radiative Q-factor with variations in the asymmetry parameter α at the Γ-point, depicted in logarithmic plots.
Figure 2. (a) Simulated spectra of the ideal DSRR composed of PEC and realistic metallic (Al) with varying structural parameters g2. (b) Simulated transmittance spectra map of the DSRR for various g2. (c) The band structure and radiative Q-factor of the DSRR structure at gap sizes of 4 μm. (d) The curve shows the change in the radiative Q-factor with variations in the asymmetry parameter α at the Γ-point, depicted in logarithmic plots.
Crystals 15 00096 g002
Figure 3. (a) The total scattering power. (b) The scattering powers of different multipoles including ED, MD, TD, EQ, and MQ in the Cartesian coordinate system. In the cases of (c) g2 = 4 μm and (d) g2 = 10 μm, scattering direction diagrams between ED and TD and radiation power diagrams between EDx−ik0TDx were plotted, respectively.
Figure 3. (a) The total scattering power. (b) The scattering powers of different multipoles including ED, MD, TD, EQ, and MQ in the Cartesian coordinate system. In the cases of (c) g2 = 4 μm and (d) g2 = 10 μm, scattering direction diagrams between ED and TD and radiation power diagrams between EDx−ik0TDx were plotted, respectively.
Crystals 15 00096 g003
Figure 4. Electromagnetic field distribution of BIC (a,b) and quasi-BIC (c,d). The green, red, blue, and black arrows represent the directions of ED, MD, TD, and surface current, respectively.
Figure 4. Electromagnetic field distribution of BIC (a,b) and quasi-BIC (c,d). The green, red, blue, and black arrows represent the directions of ED, MD, TD, and surface current, respectively.
Crystals 15 00096 g004
Figure 5. (a) Schematic diagram of NaCl concentration solution sensor. (b,c) Simulated spectra versus the refractive index of liquid. (d) Refractive index sensitivity curve of frequency versus surrounding liquid concentration. (e) Figure of merit as a function of the liquid’s refractive index.
Figure 5. (a) Schematic diagram of NaCl concentration solution sensor. (b,c) Simulated spectra versus the refractive index of liquid. (d) Refractive index sensitivity curve of frequency versus surrounding liquid concentration. (e) Figure of merit as a function of the liquid’s refractive index.
Crystals 15 00096 g005
Table 1. The Fano resonance FWHM of BICs and quasi-BICs under different g2 values.
Table 1. The Fano resonance FWHM of BICs and quasi-BICs under different g2 values.
Material
Simulation
FWHM (THz)
The Value of g2 (μm)
46810
Aluminum00.011620.024110.02996
PEC00.003590.011320.01966
Table 2. Comparing THz sensing performance achieved in different works.
Table 2. Comparing THz sensing performance achieved in different works.
Types of
Resonance
MaterialsSensitivity (GHz/RIU)FOM (RIU−1)References
EITSilicon23164.7[53]
BICGold1057.501[54]
FanoAluminum139.2-[55]
TDAluminum186-[56]
TDAluminum124.3-[57]
EITGraphene177.759.3[58]
BICSilicon7711.1[59]
TD-BICAluminum33564.53This work
EIT, electromagnetically induced transparency.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, L.; Gao, Y. High-Q Resonances Induced by Toroidal Dipole Bound States in the Continuum in Terahertz Metasurfaces. Crystals 2025, 15, 96. https://doi.org/10.3390/cryst15010096

AMA Style

Guo L, Gao Y. High-Q Resonances Induced by Toroidal Dipole Bound States in the Continuum in Terahertz Metasurfaces. Crystals. 2025; 15(1):96. https://doi.org/10.3390/cryst15010096

Chicago/Turabian Style

Guo, Lincheng, and Yachen Gao. 2025. "High-Q Resonances Induced by Toroidal Dipole Bound States in the Continuum in Terahertz Metasurfaces" Crystals 15, no. 1: 96. https://doi.org/10.3390/cryst15010096

APA Style

Guo, L., & Gao, Y. (2025). High-Q Resonances Induced by Toroidal Dipole Bound States in the Continuum in Terahertz Metasurfaces. Crystals, 15(1), 96. https://doi.org/10.3390/cryst15010096

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop