Photonics Research, Volume. 13, Issue 9, 2497(2025)

Toroidal dipole Fabry–Perot bound states in the continuum metasurfaces for ultrasensitive chiral detection

Chengfeng Li1,2,3,4,5、†, Tao He1,2,3,4、†, Xiaofeng Rao1,2,3,4,5, Chao Feng1,2,3,4,5, Jingyuan Zhu1,2,3,4,5, Siyu Dong1,2,3,4,5, Zeyong Wei1,2,3,4,5, Hongfei Jiao1,2,3,4,5, Yuzhi Shi1,2,3,4,5,7、*, Zhanshan Wang1,2,3,4,5, and Xinbin Cheng1,2,3,4,5,8、*
Author Affiliations
  • 1Institute of Precision Optical Engineering, School of Physics Science and Engineering, Tongji University, Shanghai 200092, China
  • 2MOE Key Laboratory of Advanced Micro-Structured Materials, Shanghai 200092, China
  • 3Shanghai Frontiers Science Center of Digital Optics, Shanghai 200092, China
  • 4Shanghai Professional Technical Service Platform for Full-Spectrum and High-Performance Optical Thin Film Devices and Applications, Shanghai 200092, China
  • 5Shanghai Institute of Intelligent Science and Technology, Tongji University, Shanghai 200092, China
  • 7e-mail: yzshi@tongji.edu.cn
  • 8e-mail: chengxb@tongji.edu.cn
  • show less

    Circular dichroism (CD) spectroscopy, widely used for chiral sensing, has been limited by the detection sensitivity. Enhancing optical chirality in the light fields interacting with chiral molecules is crucial for achieving ultrasensitive chiral detection. Here, we present a new paradigm for ultrasensitive chiral detection by creating accessible chiral hotspots using a toroidal dipole Fabry–Perot bound state in the continuum (TD FP-BIC) metasurface. BIC resonance is achieved by controlling the coupling between the TD resonance and its multilayer reflector-induced perfect mirror image. This method enables unprecedented local maximum and average optical chirality enhancements of up to 6×104-fold and 2×103-fold, respectively, within non-structured regions, resulting in an 866-fold increase in CD signals compared to chiral molecules alone without nanostructures. Our results pave the way for enhanced light–matter interactions and ultrasensitive enantiomeric operation.

    1. INTRODUCTION

    Chirality, or the inability of objects to align with their mirror counterparts, is prevalent at all scales, from the quantum [1] to the macroscopic [2]. In biology, many essential biomolecules, such as DNA and amino acids, exist exclusively as right-handed or left-handed enantiomers. Alternations in the chirality of these biomolecules can lead to serious diseases. Consequently, the study of chirality has garnered significant attention for its applications in various fields, including the food industry [3], drug development [4,5], medical diagnostics [6,7], materials science [8], and lateral optical forces [9,10]. Chiral molecules often exhibit optical activity, such as circular dichroism (CD), which reflects the differential absorption of right-handed and left-handed circularly polarized (RCP and LCP) light. Based on this principle, CD spectroscopy is considered to be a valuable tool for chiroptical analysis. However, due to the scale mismatch between chiral molecules and wavelength, CD signals are typically significantly weaker (103106) than the absorbance signals [11], limiting the sensitivity of chiral detection.

    Pioneering works [12,13] have indicated that CD signals are proportional to the optical chirality C, defined as C=ω2c2Im(E*·H), where ω and c denote the angular frequency and speed of light, and E and H are the complex electric and magnetic field vectors, respectively. Therefore, improving chiral detection sensitivity, or the magnitude of the CD signal, lies in generating enhanced optical chirality [1416], particularly in the light field interacting with chiral molecules. Recently, metasurfaces [1722] have been widely used to amplify CD signals by manipulating the electromagnetic field to generate enhanced optical chirality. Plasmonic metasurfaces [2325] were initially proposed to improve chiral sensing through the localized surface plasmon resonance of metals. However, concerns arose in the photothermal heating effects potentially damaging biological molecules. Subsequently, dielectric Mie metasurfaces [2630] were proposed to create strong chiral hotspots, offering a promising alternative to plasmonic metasurfaces. Nonetheless, the enhancements of optical chirality achieved through these fundamental multipoles are typically limited to an order of magnitude and primarily occur within high-index nanostructures, hindering the full interaction between chiral molecules and chiral fields, thereby reducing detection sensitivity.

    Here, we introduce a new paradigm for ultrasensitive chiral detection by proposing a toroidal dipole Fabry–Perot bound state in the continuum (TD FP-BIC) metasurface, aimed at accessing significant optical chirality enhancements outside the nanostructures, as illustrated in Fig. 1(a). The polarization-insensitive TD resonance is induced through dual symmetry breaking, creating strong chiral hotspots outside the clustered metasurface. By controlling the coupling between the TD resonance and its multilayer reflector-induced perfect mirror image, we achieve a metasurface operating at FP-BIC [31,32]. Unlike metallic mirrors that introduce significant absorption losses, the dielectric reflector provides high reflectivity with minimal dissipation, enabling low-loss interference. When the spacer thickness satisfies the FP-BIC phase condition, destructive interference between the TD mode and its mirror image suppresses radiation leakage, giving rise to a high-Q resonance with strongly enhanced electromagnetic fields and optical chirality. Under these conditions, the local maximum enhancement and average optical chirality enhancement (Cmax and Cavg), normalized to circularly polarized light in air, reach approximately 6×104-fold and 2×103-fold, respectively, within the near-infrared (NIR) spectral regime. Notably, the choice of the NIR regime is primarily motivated by its superior compatibility with practical sensing requirements, including higher detector quantum efficiency, lower material absorption losses, and more mature nanofabrication techniques. This marks significant progress compared to state-of-the-art dielectric metasurfaces [11,2628,30,3337], as shown in Fig. 1(b). As a proof of concept, the TD FP-BIC metasurface enables an 866-fold enhancement of CD signals compared to the chiral molecules alone without nanostructures. More detailed comparisons between this work and previously reported studies on chiral sensing are provided in Appendix B. Our results present a novel strategy for manipulating the chiral electromagnetic field, revealing new possibilities in ultrasensitive sensing and chiral photonics.

    (a) Schematic diagram illustrating ultrasensitive chiral sensing based on a TD FP-BIC metasurface, which consists of a dielectric metasurface and a multilayer reflector separated by a dielectric spacer. The emergence of TD FP-BIC results from destructive interference between TD resonance and the multilayer reflector-induced perfect mirror image. (b) Comparison among this work and previously reported metasurfaces on Cmax and Cavg for chiral detection.

    Figure 1.(a) Schematic diagram illustrating ultrasensitive chiral sensing based on a TD FP-BIC metasurface, which consists of a dielectric metasurface and a multilayer reflector separated by a dielectric spacer. The emergence of TD FP-BIC results from destructive interference between TD resonance and the multilayer reflector-induced perfect mirror image. (b) Comparison among this work and previously reported metasurfaces on Cmax and Cavg for chiral detection.

    2. RESULT

    A. Design of the Polarization-Independent TD Metasurface

    TD [3841], the third family of electromagnetic dipoles alongside the electric dipole (ED) and magnetic dipole (MD), has garnered considerable attention due to its significant field enhancements and unique exposed field characteristic. A straightforward method to induce TD resonance [42] involves creating the symmetry breaking within the dielectric cluster, as shown in Figs. 2(a) and 2(b). The unit cell of the metasurface is immersed in water and composed of four Si nanofins on a SiO2 substrate. The metasurface has a period (p) of 600 nm in both x and y directions, with an intra-cluster distance (d) of 220 nm and height (h) of 500 nm. The square nanofins have lengths of a=180  nm and b=120  nm. For a symmetric metasurface with four identical nanofins [Fig. 2(a)], only one resonance appears in the transmission spectra, but not the TD resonance. Analysis of multipole decomposition and electromagnetic field distributions (see Fig. 5 in Appendix C) suggests this resonance arises from a collective response of four longitudinal MDs. Due to the confinement of field enhancements within nanostructures, symmetric metasurfaces with strong MD responses are not suitable for chiral sensing applications. Introducing symmetry breaking results in the splitting of the previous single resonance in transmission spectra [Fig. 2(b)], revealing additional resonances. Multipole decomposition and electromagnetic field distributions (see Fig. 6 in Appendix C) show a pronounced TD response at approximately 1090 nm exclusively under the TM excitation. However, polarization-dependent TD responses in the metasurface hinder the generation of enhanced optical chirality necessary for high-sensitivity chiral detection.

    Design of the polarization-independent TD metasurface. The transmission spectra of symmetric metasurfaces (a), metasurfaces with single symmetry breaking (b), and asymmetric metasurfaces with dual symmetry breaking (c). Here, TE refers to the electric field polarized along the y direction, and TM refers to the electric field polarized along the x direction. (d) The multipole decomposition of the TD metasurface. More details about multipole decomposition can be found in Appendix A. (e) The Cavg of the TD metasurface under RCP and LCP incidence, respectively. (e) The calculated Cavg of the TD metasurface under RCP and LCP incidence, respectively. The light orange area over the metasurface shown in the inset corresponds to the actual volume used for calculating Cavg. The electric field (f), magnetic field (g), and optical chirality enhancements (h) at z=h/2 under RCP excitation. The electric field (i), magnetic field (j), and optical chirality enhancements (k) at z=h/2 under LCP excitation. The white arrows in (f)–(k) represent the electric displacement fields.

    Figure 2.Design of the polarization-independent TD metasurface. The transmission spectra of symmetric metasurfaces (a), metasurfaces with single symmetry breaking (b), and asymmetric metasurfaces with dual symmetry breaking (c). Here, TE refers to the electric field polarized along the y direction, and TM refers to the electric field polarized along the x direction. (d) The multipole decomposition of the TD metasurface. More details about multipole decomposition can be found in Appendix A. (e) The Cavg of the TD metasurface under RCP and LCP incidence, respectively. (e) The calculated Cavg of the TD metasurface under RCP and LCP incidence, respectively. The light orange area over the metasurface shown in the inset corresponds to the actual volume used for calculating Cavg. The electric field (f), magnetic field (g), and optical chirality enhancements (h) at z=h/2 under RCP excitation. The electric field (i), magnetic field (j), and optical chirality enhancements (k) at z=h/2 under LCP excitation. The white arrows in (f)–(k) represent the electric displacement fields.

    In response to this challenge, we propose the concept of dual symmetry breaking to obtain a polarization-independent TD response. Although the metasurface exhibits asymmetries in both x and y directions, as shown in Fig. 2(c), identical additional resonances are observed in the transmission spectra for both TM and TE polarization excitations. Multipole decomposition and electromagnetic field distributions for these excitations are detailed in Fig. 7 (Appendix C), confirming the polarization-independent TD responses. This design evolution, as detailed in Appendix C, enables the stepwise activation of high-Q TD resonances under orthogonal polarizations. As a result, the final dual-symmetry-breaking configuration achieves polarization-insensitive high-Q responses, extending previous single-axis designs toward more versatile bidirectional control [42]. Multipole decomposition under RCP excitation, illustrated in Fig. 2(d), highlights the predominance of the TD resonance. To assess chiral detection enhancement, we calculate the Cavg under RCP and LCP excitations, as illustrated in Fig. 2(e). Here, the simulation domain for calculating Cavg is complementary to the nanostructures, as indicated by the light orange region in the inset of Fig. 2(e). It corresponds to the non-structured portion of a 600  nm×600  nm×500  nm volume along the x, y, and z directions, respectively, and is normalized to the incident circularly polarized light. Remarkably, the Cavg reaches nearly a 100-fold enhancement. This achievement is attributed to the significant electromagnetic field enhancements and unique field configuration of TD resonance.

    For an intuitive understanding of the significant enhancements in optical chirality, Figs. 2(f)–2(h) depict the electromagnetic field distributions at z=h/2 under the RCP illumination. The z=h/2 plane corresponds to the midpoint of the nanofin’s height. Within the unit cell, torus-like electromagnetic field lines and ring-like electric fields closely resemble the theoretical patterns of TD responses [4244], affirming their presence and significant contribution. Therefore, although both MD and TD modes are present in the multipole decomposition, the dominant toroidal features support the identification of this structure as a TD metasurface. The electric and magnetic fields exhibit maximum enhancements of approximately 35 times and 45 times, respectively, while Cmax reaches up to 650 times outside the nanostructures, providing a favorable condition for detection. Similarly, Figs. 2(i)–2(k) illustrate the electromagnetic field distributions at z=h/2 under the LCP illumination. It is noteworthy that the TD metasurface exhibits identical electromagnetic field distributions but opposite enhancements in optical chirality under RCP and LCP excitation, respectively. These results confirm the intrinsically achiral nature of the metasurface and ensure that the subsequently observed CD signals originate solely from the chiral molecules, allowing the metasurface to function as a high-efficiency and background-free CD signal amplifier.

    B. Demonstration of the TD FP-BIC Metasurface

    To further enhance optical chirality, we create multilayer reflector-induced TD resonance to couple with the TD resonance metasurface. This controlled coupling between the TD resonance and its multilayer reflector-induced counterpart gives rise to TD FP-BIC phenomena. Figure 3(a) illustrates that the TD FP-BIC metasurface consists of the previously designed TD metasurface and an all-dielectric multilayer reflector separated by a SiO2 spacer, with additional details provided in Appendix E. To determine the BIC condition, the dielectric spacer is meticulously designed to control the coupling between the TD resonance and its mirror image. By varying the thickness of the spacer (t), we can observe clear transitions in the transmission spectra from quasi-TD FP-BIC to TD FP-BIC and back to quasi-TD FP-BIC in Fig. 3(a). Moreover, the Cavg exhibits a similar variation pattern to the transmission spectra, as shown in Fig. 3(b), suggesting that the significant enhancement in optical chirality stems from the BIC phenomena. Opting for a spacer thickness of 843 nm, as indicated by the green dashed rectangular frame, allows us to excite quasi-TD FP-BIC and achieve significant optical chirality enhancement. We speculate that the low transmission of this quasi-TD FP-BIC resonance arises from the reflectance mismatch between the upper metasurface and the underlying multilayer reflector, with reflectance of 95% and 100%, respectively. Further analysis and evidence can be found in Appendix D.

    Design and characterization of the TD FP-BIC metasurface. (a) The transmission spectra of the TD FP-BIC metasurface with respect to the thickness of spacer t. The illustrations depict the schematic diagram of the metasurface. (b) The Cavg of the TD FP-BIC metasurface versus the thickness of spacer t. (c) The multipole decomposition of the TD FP-BIC metasurface. (d) The Cavg of the TD FP-BIC metasurface under RCP excitation. The optical chirality enhancements at z=h/2 for wavelength 1152.44 nm (e), 1152.45 nm (f), and 1152.46 nm (g) under RCP incidence. The white arrows in (e)–(g) represent the magnetic fields.

    Figure 3.Design and characterization of the TD FP-BIC metasurface. (a) The transmission spectra of the TD FP-BIC metasurface with respect to the thickness of spacer t. The illustrations depict the schematic diagram of the metasurface. (b) The Cavg of the TD FP-BIC metasurface versus the thickness of spacer t. (c) The multipole decomposition of the TD FP-BIC metasurface. (d) The Cavg of the TD FP-BIC metasurface under RCP excitation. The optical chirality enhancements at z=h/2 for wavelength 1152.44 nm (e), 1152.45 nm (f), and 1152.46 nm (g) under RCP incidence. The white arrows in (e)–(g) represent the magnetic fields.

    Figure 3(c) illustrates the multipole decomposition of the TD FP-BIC metasurface, indicating that the TD responses remain dominant. The difference is that, in this case, the MTD dominates instead of the ETD. Under these conditions, the Cavg in Fig. 3(d) shows a clear sign reversal from negative to positive around the FP-BIC resonance, achieving a remarkable enhancement of 2200-fold. To gain a better insight into the enhancement of optical chirality, we demonstrate the optical chirality enhancement distributions at z=h/2 for three different wavelengths, as shown in Figs. 3(e)–3(g). Here, the blue, green, and yellow pentagons are consistent with those in Fig. 3(d), representing wavelengths of 1152.44 nm, 1152.45 nm, and 1152.46 nm, respectively. Strikingly, the local maximum enhancement of optical chirality reaches about 6×104-fold for a wavelength of 1152.45 nm. To the best of our knowledge, this is the highest value within non-structured regions reported to date. Strikingly, the local maximum enhancement of optical chirality reaches approximately 6×104-fold at a wavelength of 1152.45 nm. To the best of our knowledge, this represents one of the highest values reported to date, even when accounting for the slight material absorption due to fabrication imperfections. Further details on the impact of fabrication-induced variations and material absorption on the metasurface performance are provided in Fig. 9 of Appendix F, highlighting the robustness of our design strategy.

    Additionally, cut planes demonstrating optical chirality enhancements at various z heights under RCP illumination are illustrated in Appendix G. Remarkably, these enhancements maintain consistently high levels across different heights, ranging from z=0 to z=h. These results underscore the significant potential of the proposed TD FP-BIC metasurface to enhance CD signals for chiral molecules across a wide volume range, emphasizing its utility in chiral sensing.

    C. Ultrasensitive Chiral Detection Based on the TD FP-BIC Metasurface

    To explore the practical applications of the TD FP-BIC metasurface in ultrasensitive chiral detection, we consider the scenario where the metasurface is uniformly immersed in a homogenous chiral molecule layer with a thickness of 500 nm. The chiral molecules exhibit a refractive index of n=1.33+105i and a fixed Pasteur parameter of κ=105i, consistent with the parameters used in previous works [45] and the realistic value for chiral biomolecules [46,47], such as aqueous solutions of monosaccharides [48,49].

    Chiral detection based on the TD FP-BIC metasurface. The absorption spectra of chiral molecules alone (a), chiral molecules enhanced by TD metasurface (b), and chiral molecules enhanced by the TD FP-BIC metasurface (c) under RCP and LCP incidence, respectively. The insets show the simulation schematics. (d) CD signals for the three scenarios mentioned above [shown in (a)–(c)]. (e) CD signal enhancements with the assistance of the TD metasurface and TD FP-BIC metasurface.

    Figure 4.Chiral detection based on the TD FP-BIC metasurface. The absorption spectra of chiral molecules alone (a), chiral molecules enhanced by TD metasurface (b), and chiral molecules enhanced by the TD FP-BIC metasurface (c) under RCP and LCP incidence, respectively. The insets show the simulation schematics. (d) CD signals for the three scenarios mentioned above [shown in (a)–(c)]. (e) CD signal enhancements with the assistance of the TD metasurface and TD FP-BIC metasurface.

    3. CONCLUSION

    In conclusion, we have introduced an ultrasensitive chiral detection platform utilizing a TD FP-BIC metasurface. The TD FP-BIC metasurface leverages controlled coupling between TD resonance and its mirrored counterpart, achieving record-high Cmax and Cavg enhancements of 6×104-fold and 2×103-fold, respectively, in the non-structured region. With the TD FP-BIC metasurface, we realize a 2-order-of-magnitude enhancement in the CD signal. To our knowledge, this marks a significant advancement over prior studies. Additionally, our proposed paradigm is compatible with resonance-gradient metasurfaces [53,55,56] for broad spectral range chiral detection. Compared with the actively tunable-material-based metasurfaces [5759], this approach avoids additional fabrication complexities, reduces costs, and eliminates response time delays, while maintaining high-Q resonances and strong optical chirality. This strategy offers a practical, scalable solution for real-world chiroptical sensing, as shown in Fig. 12 of the Appendix I. Our results pave a new avenue toward high-performance enantiomeric detection and separation, providing a novel design strategy to manipulate the near-field electromagnetic fields, which may inspire more innovative photonic devices for enhanced light–matter interaction in the future.

    APPENDIX A: NUMERICAL SIMULATIONS AND MULTIPOLE DECOMPOSITION

    The numerical simulations of the spectra, multipole decomposition, electromagnetic field distributions, and CD signals presented in this paper are performed by commercially available software (COMSOL Multiphysics 6.0) using the frequency domain solver in the Wave Optics Module. Periodic boundary conditions are used in x and y directions, while PML boundary conditions are used in the z direction. It is worth mentioning that although our simulations are based on infinite periodicity, previous studies have shown that quasi-BIC resonances can persist in finite arrays as small as 10×10 or 27×27 unit cells [60,61], depending on the measurement technique. This suggests that the strong chiral field enhancement demonstrated in our design remains practically achievable in realistic device implementations. The maximum size of the mesh was set to λ/10, with λ being the incident wavelength. Linearly and circularly polarized incident plane waves are defined at the incident port. The refractive indices of Si, SiO2, and water in the manuscript are 3.5, 1.5, and 1.33, respectively. The chiral molecules can be modeled as a homogeneous layer whose constitutive relations can be expressed as D=εmE+(iκ/c)H,B=(iκ/c)E+μmH,where εm and μm are the permittivity and permeability, and κ is the Pasteur parameter.

    The scattering intensity of each electromagnetic multipole can be expressed as Isca=2ω43c3(|p|2+|m|2)+2ω6εm3c5(ε0μ0εm|T(e)|2+|T(m)|2)+ω65c5|Qαβ(e)|2+ω640c5|Qαβ(m)|2,where α,β=x,y,z, εm is the permittivity of the medium, and p, m, T(e), T(m), Qαβ(e), and Qαβ(m) denote the electric dipole (ED), magnetic dipole (MD), electric toroidal dipole (ETD), magnetic toroidal dipole (MTD), electric quadrupole (EQ), and magnetic quadrupole (MQ) in Cartesian coordinate, respectively.

    APPENDIX B: COMPARISON AMONG THIS WORK AND PREVIOUSLY REPORTED STUDIES

    In this section, we provide detailed comparison among our work and previously reported studies on chiral sensing in Table 1.

    APPENDIX C: EVOLUTION OF THE POLARIZATION-INDEPENDENT TD METASURFACE

    In this section, we calculate the multipole decomposition and field distribution of the symmetric metasurface (Fig. 5), the metasurface with single symmetry breaking (Fig. 6), and the metasurface with double symmetry breaking (Fig. 7) to illustrate the detailed design process of the polarization-independent TD metasurface.

    Analysis of the symmetric metasurface. (a) The multipole decomposition of the symmetric metasurface under TM polarization. The electric field (b) and magnetic field (c) of the metasurface at z=h/2 under TM polarization. (d) The multipole decomposition of the symmetric metasurface under TE polarization. The electric field (e) and magnetic field (f) of the metasurface at z=h/2 under TE polarization. The white arrows in figures represent the electric displacement fields.

    Figure 5.Analysis of the symmetric metasurface. (a) The multipole decomposition of the symmetric metasurface under TM polarization. The electric field (b) and magnetic field (c) of the metasurface at z=h/2 under TM polarization. (d) The multipole decomposition of the symmetric metasurface under TE polarization. The electric field (e) and magnetic field (f) of the metasurface at z=h/2 under TE polarization. The white arrows in figures represent the electric displacement fields.

    Analysis of the metasurface with single symmetry breaking. (a) The multipole decomposition of the metasurface with single symmetry breaking under TM polarization. The electric field (b) and magnetic field (c) of the metasurface at z=h/2 under TM polarization. (d) The multipole decomposition of the metasurface with single symmetry breaking under TE polarization. The electric field (e) and magnetic field (f) of the metasurface at z=h/2 under TE polarization. The white arrows in figures represent the electric displacement fields.

    Figure 6.Analysis of the metasurface with single symmetry breaking. (a) The multipole decomposition of the metasurface with single symmetry breaking under TM polarization. The electric field (b) and magnetic field (c) of the metasurface at z=h/2 under TM polarization. (d) The multipole decomposition of the metasurface with single symmetry breaking under TE polarization. The electric field (e) and magnetic field (f) of the metasurface at z=h/2 under TE polarization. The white arrows in figures represent the electric displacement fields.

    Analysis of the metasurface with dual symmetry breaking. (a) The multipole decomposition of the metasurface with dual symmetry breaking under TM polarization. The electric field (b) and magnetic field (c) of the metasurface at z=h/2 under TM polarization. (d) The multipole decomposition of the metasurface with dual symmetry breaking under TE polarization. The electric field (e) and magnetic field (f) of the metasurface at z=h/2 under TE polarization. The white arrows represent the electric displacement fields.

    Figure 7.Analysis of the metasurface with dual symmetry breaking. (a) The multipole decomposition of the metasurface with dual symmetry breaking under TM polarization. The electric field (b) and magnetic field (c) of the metasurface at z=h/2 under TM polarization. (d) The multipole decomposition of the metasurface with dual symmetry breaking under TE polarization. The electric field (e) and magnetic field (f) of the metasurface at z=h/2 under TE polarization. The white arrows represent the electric displacement fields.

    APPENDIX D: REFLECTANCE MISMATCH BETWEEN THE METASURFACE AND MULTILAYER REFLECTOR

    To illustrate that the low transmittance of the TD FP-BIC metasurface arises from the reflectance mismatch between the upper metasurface and underlying multilayer reflectors (Fig. 3), we used a multilayer film with a lower reflectance to achieve the BIC resonance. As illustrated in Fig. 8, due to the decrease in the reflectance of the multilayer film, the reflectance of the multilayer film (97%) and metasurface (95%) match better, resulting in a higher transmittance of the TD FP-BIC, but a lower optical chirality enhancement.

    Analysis of the reflectance mismatch between the metasurface and multilayer reflector. The reflectance properties of the multilayers with varying numbers of Si and SiO2 pairs: eight pairs used in the main text with reflectance of 100% (a), and three pairs with reflectance of 97% (b). (c) The transmission spectra of the TD FP-BIC metasurface with respect to the thickness of spacer t. The illustrations depict the schematic diagram of the metasurface, which consists of multilayer reflector with lower reflectance. (d) The Cavg of the TD FP-BIC metasurface versus the thickness of spacer t.

    Figure 8.Analysis of the reflectance mismatch between the metasurface and multilayer reflector. The reflectance properties of the multilayers with varying numbers of Si and SiO2 pairs: eight pairs used in the main text with reflectance of 100% (a), and three pairs with reflectance of 97% (b). (c) The transmission spectra of the TD FP-BIC metasurface with respect to the thickness of spacer t. The illustrations depict the schematic diagram of the metasurface, which consists of multilayer reflector with lower reflectance. (d) The Cavg of the TD FP-BIC metasurface versus the thickness of spacer t.

    APPENDIX E: PARAMETERS OF MULTILAYER FILMS

    In this section, we provide detailed information about Si/SiO2 multilayer films in Table 2.

    The Thicknesses of the Si(H) and SiO2(L) Multilayer Films (in nm)

    Layerd1/Hd2/Ld3/Hd4/Ld5/Hd6/Ld7/Hd8/Ld9/H
    Thickness7517575175751757517575
    Layerd10/Ld11/Hd12/Ld13/Hd14/Ld15/Hd16/LSpacer/L
    Thickness175751757517575175843

    APPENDIX F: FABRICATION TOLERANCE ANALYSIS OF THE TD FP-BIC METASURFACE

    Figure 9 shows the fabrication tolerance analysis of the TD FP-BIC metasurface. The analysis includes the following cases. System.Xml.XmlElementSystem.Xml.XmlElement

    (a) Cavg for the TD FP-BIC metasurface with a −2% deviation in both width and length of the nanostructures. (b) Cavg for the TD FP-BIC metasurface with a +2% deviation in both width and length of the nanostructures. (c) Cavg for the TD FP-BIC metasurface with a −2% variation in the multilayer reflector thickness, spacer thickness, and the height of the nanostructures. (d) Cavg for the TD FP-BIC metasurface with a +2% variation in the multilayer reflector thickness, spacer thickness, and the height of the nanostructures.

    Figure 9.(a) Cavg for the TD FP-BIC metasurface with a 2% deviation in both width and length of the nanostructures. (b) Cavg for the TD FP-BIC metasurface with a +2% deviation in both width and length of the nanostructures. (c) Cavg for the TD FP-BIC metasurface with a 2% variation in the multilayer reflector thickness, spacer thickness, and the height of the nanostructures. (d) Cavg for the TD FP-BIC metasurface with a +2% variation in the multilayer reflector thickness, spacer thickness, and the height of the nanostructures.

    APPENDIX G: DETAILED OPTICAL CHIRALITY ENHANCEMENTS

    Figure 10 shows the detailed optical chirality enhancements of the TD FP-BIC metasurface at various heights for different wavelengths. All the enhancements maintain consistently high levels across different heights ranging from z=0 to z=h, emphasizing its utility in chiral sensing.

    Detailed optical chirality enhancements of the TD FP-BIC metasurface at various heights for different wavelengths. The optical chirality enhancements at different z positions: at z=0 for wavelengths 1152.44 nm (a), 1152.45 nm (b), and 1152.46 nm (c); at z=h/2 for wavelengths 1152.44 nm (d), 1152.45 nm (e), and 1152.46 nm (f); at z=h for wavelengths 1152.44 nm (g), 1152.45 nm (h), and 1152.46 nm (i), respectively. All conditions are under RCP excitation.

    Figure 10.Detailed optical chirality enhancements of the TD FP-BIC metasurface at various heights for different wavelengths. The optical chirality enhancements at different z positions: at z=0 for wavelengths 1152.44 nm (a), 1152.45 nm (b), and 1152.46 nm (c); at z=h/2 for wavelengths 1152.44 nm (d), 1152.45 nm (e), and 1152.46 nm (f); at z=h for wavelengths 1152.44 nm (g), 1152.45 nm (h), and 1152.46 nm (i), respectively. All conditions are under RCP excitation.

    APPENDIX H: THE IMPACT OF CHIRAL MOLECULES ON CD ENHANCEMENT

    Figure 11 illustrates chiral detection based on the TD FP-BIC metasurface for chiral molecules with n=1.33+104i. The degeneration in magnitude between the CD enhancement and Cavg can be attributed to the perturbation caused by chiral molecules on the resonance of the BIC, which can be further alleviated by considering the coupling between the chiral molecules and metasurfaces.

    Chiral detection based on the TD FP-BIC metasurface for chiral molecules with n=1.33+10−4i. The absorption spectra of chiral molecules alone (a), chiral molecules enhanced by the TD metasurface (b), and chiral molecules enhanced by the TD FP-BIC metasurface (c) under RCP and LCP incidence, respectively. The insets show the simulation schematics. (d) CD signals for the three scenarios mentioned above [shown in (a)–(c)]. (e) CD signal enhancements with the assistance of the TD metasurface and TD FP-BIC metasurface.

    Figure 11.Chiral detection based on the TD FP-BIC metasurface for chiral molecules with n=1.33+104i. The absorption spectra of chiral molecules alone (a), chiral molecules enhanced by the TD metasurface (b), and chiral molecules enhanced by the TD FP-BIC metasurface (c) under RCP and LCP incidence, respectively. The insets show the simulation schematics. (d) CD signals for the three scenarios mentioned above [shown in (a)–(c)]. (e) CD signal enhancements with the assistance of the TD metasurface and TD FP-BIC metasurface.

    APPENDIX I: BROADBAND DETECTION BASED ON THE TD FP-BIC METASURFACE

    Figure 12 presents the concept of a resonance gradient metasurface tailored for broadband spectral detection based on our proposed strategy, demonstrating its strong potential for practical chiral sensing.

    Concept of the gradient metasurface for broad spectral range detection. (a) Schematic representation of a series of gradient metasurfaces in which the in-plane dimensions of the resonators increase gradually from left to right. The structural parameters, length a and length b, are identical to those in Fig. 2, and the factor s denotes the applied scaling factor. (b) Numerical results for Cavg of the TD FP-BIC metasurfaces with scaling factors ranging from 0.98 to 1.02. The peak resonance of Cavg shifts from 1130 nm to 1170 nm, while maintaining high enhancement values, showcasing the potential for border spectral range detection.

    Figure 12.Concept of the gradient metasurface for broad spectral range detection. (a) Schematic representation of a series of gradient metasurfaces in which the in-plane dimensions of the resonators increase gradually from left to right. The structural parameters, length a and length b, are identical to those in Fig. 2, and the factor s denotes the applied scaling factor. (b) Numerical results for Cavg of the TD FP-BIC metasurfaces with scaling factors ranging from 0.98 to 1.02. The peak resonance of Cavg shifts from 1130 nm to 1170 nm, while maintaining high enhancement values, showcasing the potential for border spectral range detection.

    [28] J. Garcia-Guirado, M. Svedendahl, J. Puigdollers. Enhanced chiral sensing with dielectric nanoresonators. Nano Lett., 20, 585-591(2020).

    Tools

    Get Citation

    Copy Citation Text

    Chengfeng Li, Tao He, Xiaofeng Rao, Chao Feng, Jingyuan Zhu, Siyu Dong, Zeyong Wei, Hongfei Jiao, Yuzhi Shi, Zhanshan Wang, Xinbin Cheng, "Toroidal dipole Fabry–Perot bound states in the continuum metasurfaces for ultrasensitive chiral detection," Photonics Res. 13, 2497 (2025)

    Download Citation

    EndNote(RIS)BibTexPlain Text
    Save article for my favorites
    Paper Information

    Category: Nanophotonics and Photonic Crystals

    Received: May. 1, 2025

    Accepted: Jun. 28, 2025

    Published Online: Aug. 25, 2025

    The Author Email: Yuzhi Shi (yzshi@tongji.edu.cn), Xinbin Cheng (chengxb@tongji.edu.cn)

    DOI:10.1364/PRJ.559587

    CSTR:32188.14.PRJ.559587

    Topics