Excitation of tunable dual quasi-bound states in the continuum in graphene metasurface and terahertz sensing application

Excitation of tunable dual quasi-bound states in the continuum in graphene metasurface and terahertz sensing application | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Excitation of tunable dual quasi-bound states in the continuum in graphene metasurface and terahertz sensing application Zhong Huang, Qing Chen, Xiaoxiao Cao, Ji Ding, Chaojun Tang, Zhendong Yan, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-3087509/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 14 Jul, 2023 Read the published version in Plasmonics → Version 1 posted 7 You are reading this latest preprint version Abstract This study proposes an all-graphene metasurface supporting dual band symmetric bound state in the continuum (BIC) in the terahertz (THz) range for the first time. The structure consists of a unit cell containing a double-gap split ring resonator periodically on a dielectric substrate. By introducing symmetry breaking, two plasmonic quasi-BICs (Q-BICs) transmission dips are observed with finite Q-factors and high modulation depth. Simulation and analysis results simultaneously exhibit that the Q-factors of Q-BICs follow an inverse square dependency with the asymmetry degree. Via changing the graphene’s chemical potential, the Q-BICs’ operating frequency range can be actively expanded. The dual Q-BICs are immune to the variation of incident angle and has a significant slow light effect with time delay up to 23.8 ps. In addition, the sensing performance in THz region is investigated. A maximum sensitivity of 267.5 GHz/RIU is obtained with a FOM of 24.08 RIU − 1 . Our work shows an alternative way to design the class of tunable Q-BIC metasurface, which will provide a valuable reference for future dynamic sensor and other fields. Graphene metasurface Bound state in the continuum (BIC) Plasmonic quasi-BIC Terahertz sensing Tunable Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 1. Introduction In 1929, the idea of bound states in the continuum (BIC) was firstly reported by J. V. Neumann and E. Wigner in quantum mechanics during the process of unique Schrödinger equation calculation [ 1 ]. The BIC, however, is often brittle, necessitates uniquely tuned potentials, and is hence challenging to observe in the quantum experiments. Because the Maxwell equation and the Schrodinger equation are formal equals, the emergence of these exotic non-radiating states in optics, is strongly influenced by the quantum mechanical system [ 2 ]. Similar to quantum systems, the concept of BIC in optics is the result of two or more waves interfering destructively to reduce radiative loss, where radiation resonances can generally carry power from or to infinity. At the operating frequency of BIC, the embedded eigenvalues vanish with the leakage power and cause eigenstates having an infinitely high-quality factor (Q-factor) [ 3 ]. Generally, there are two ways to excite BIC resonance. One is tuning the parameters of optic systems to make outgoing waves interfering destructively, and this type of BIC is named as accidental BIC [ 4 ]. The other way is breaking the in-plane symmetry to generate mismatch between the symmetry of outgoing wave and the resonance mode, and this type is symmetry protected BIC [ 5 ]. Due to its novel ability to have the high Q-factor, BIC has been widely studied in various optical systems of quantum dots lasers [ 6 ], filters [ 7 ], sensors [ 8 ] and nonlinear optics [ 9 ]. Recently, BIC has been introduced into the metasurfaces [ 10 ] owing to the advantage of electromagnetic wave confinement in sub-wavelength, precise wavefront control, and the simplicity of fabrication in two-dimensional planes [ 11 ]. The periodic arrangements of artificial sub-wavelength unit cells in metasurfaces result in some diffraction open channels. The symmetry protected BIC is easier to be existed in the absence of coupling to these open channels by engineering the periodic sub-wavelength structures [ 12 ]. Interestingly, true optical BIC (ideal BIC) is a mathematically defined concept with the infinite value of Q-factor and its resonance linewidth infinitely tends towards 0. This is only possible in an ideal lossless infinite optical system with extreme structure parameters. In practice, suffering from material absorption, finite lateral size of structures and other leakage, the observable BIC often refers to quasi-BIC (Q-BIC) with high, yet finite Q-factor [ 13 ]. Although Q-BIC has received extensive attention in metasurfaces, most of the devices reported in the literatures are lack of dynamic capability. Once the manufacturing parameters are fixed, the operating frequency of Q-BIC is synchronous fixed. So the ability to expand the Q-BIC’s operating frequency range without changing the geometrical parameters of the structure is very important and attractive. Currently, graphene has exhibited powerful tuning ability for active devices [ 14 , 15 ]. By controlling the chemical potential E f of the graphene, the target operating frequency region can be dynamically shifted [ 16 ]. However, the active tuning of Q-BIC reported in the literatures are mainly focused on the switching characteristics controlled by graphene absorption. For example, the authors in ref. [ 17 ] proposed dual Q-BICs in a silicon stair-like metasurface, where the transmittance can be controlled with respect to a change in the chemical potential of graphene layer. Similarly, the authors in ref. [ 18 ] have reported graphene-metal hybrid terahertz metasurfaces with a similar kind of control. So far, to the best knowledge of us, only two groups have reported dynamic tuning the operating frequency of Q-BIC via the changing of the chemical potential E f of the graphene. Wang et al. have investigated the symmetry-protected Q-BIC by using the structure of two contacting graphene nanoribbons [ 19 ], and Roy et al. have proposed rectangular dimer-shaped graphene metasurfaces that supported asymmetric Q-BIC [ 20 ]. Therefore, the tunable aspects of graphene in Q-BIC are still lack of research. Besides that, tunable Q-BIC with only one operating band also faces time delay introduced by modulation, especially in sensing applications in detecting a variety of gas and chemical solvents. In this article, the exciting of symmetry-protected dual Q-BICs in an all graphene-based metasurface is theoretically demonstrated for the first time. The metasurface is made of a square lattice of which unit cell consists of a double-gap split ring resonator (DSRR). By displacing the two gaps from the center, the in-plane C 2 rotational symmetry is broken, as shown in Fig. 1 . This symmetry breaking provides an efficient way to transform BICs into Q-BICs, manifested as two new sharp dips in the transmission spectrum. The Q factor of dual Q-BICs can be well evaluated by the theoretical normalized Fano line shape equation, and follow the inverse square dependence on the asymmetry degree. The dual Q-BICs are further analyzed by the near-field electric field distributions and surface current flows, indicating the leakage of an ideal bound state can be caused by a slight perturbation on the effective electric dipole moments. We show the operating frequency of the dual Q-BICs can be actively tuned by varying the chemical potential E f of the graphene. Our Q-BIC metasurface is immune to the variation of incident angle, and has a significant slow light effect with time delay up to 23.8 ps. Moreover, our Q-BIC metasurface can operate as a refractive index sensor in the THz region. A maximum sensitivity of 267.5 GHz/RIU is obtained with a FOM of 24.08 RIU − 1 . Our work shows an alternative way to design the class of tunable Q-BIC metasurface, which will provide a valuable reference for future dynamic sensor and other fields. 2. Design and simulation of dual Q-BICs metasurface The schematic of the proposed Q-BIC metasurface is shown in Fig. 1 (a). On a thick SiO 2 substrate, a monolayer graphene is patterned into periodic DSRRs along x and y directions. The relative permittivity of the THz transparent SiO 2 substrate is taken as 3.9. Figure 1 (b) shows the unit cell structure. Here, the period along x and y directions are P x and P y , and P x = P y = 50 µm. The arm length of the DSSR is l = 31 µm, and the arm width is w = 4 µm. The gap length g = 2 µm. In order to break the in plane C 2 rotational symmetry of the DSSR, we displace both two gaps from the center by a distance “ δ ”. In this work, we use open souce software “East FDTD” [ 21 , 22 ] based on the finite-difference time-domain method to design our structure and calculate the transmission spectrum numerically. Along the x - and y -axis, we choose periodic boundary conditions to simulate a system of infinite size. Along the z -axis direction, the matesurface is excited by a normally incident broadband terahertz plane wave polarizing along x -direction, and two perfectly matched layers (PMLs) behind Floquet periodic ports are used to absorb the scattered light. For the numerical calculations, graphene is treated as a surface conductive layer with the thickness of 1 nm, and its surface conductivity is defined by analytical Kubo’s formula [ 21 , 22 ]: $${\sigma _g}=\frac{{2{e^2}{k_B}T}}{{\pi {\hbar ^2}}}Ln\left[ {2\cosh (\frac{{{E_f}}}{{2{k_B}T}})} \right]\frac{i}{{\omega +i{\tau ^{ - 1}}}}$$ 1 , where ω represents the frequency of incident terahertz wave, T represents the environmental temperature, e represents an electron charge, k B represents the Boltzmann's constant, ħ represents the reduced Planck constant, and E f represents the chemical potential of graphene DSSRs. In our numerical calculations, room temperature is taken as T = 300 K, and the relaxation time of graphene DSSRs is taken as 𝜏 = 2 ps. 3. Results and discussions We start with the investigation of transmission spectra from Q-BIC metasurface with in plane C 2 rotational symmetry, which two split gaps of the graphene DSRR are located at the center of the left and right arms and corresponding to asymmetric parameter δ = 0. One can see from the black line in Fig. 2 (a) that the graphene DSRRs show an obvious plasmonic resonance at 0.59 THz when excited by the incident THz wave along x -direction. The transmittance at this frequency point is 0.1. In this simulation, the chemical potential of the graphene is E f = 0.4 eV. The form of the electric field distribution at the frequency of the minimum value of transmittance demonstrates the trapped mode nature of the resonance, as shown in Fig. 3 (c). The electric field of plasmonic resonance mode is mostly concentrated in the left and right branches of the graphene DSRR. There are the in-phased surface current flows parallel to the x- and y -axis of the DSRR elements (white arrows in Fig. 3 (c)). The amplitude of the surface current flow parallel to the x -axis is greater than that parallel to the y -axis. This amplitudes distribution of currents provides a strong coupling between the trapped mode and the incident THz wave [ 23 ]. We then plot the transmission spectrum of the periodic graphene DSRRs with distinct asymmetric parameters δ in Fig. 2 (a). Besides the plasmonic resonance at 0.59 THz, there are two new transmission dips appear with the introduction of asymmetric parameter δ = 0.5 µm. The exhibited two new narrow resonances are strongly altered with increasing δ from 0 to 4 µm. The two new transmission dips become prominent with the increase of δ , and the linewidth gradually increase. Meantime, the two new transmission dips are slightly red shifted. As the non-zero value of asymmetric parameters δ breaks the C 2 rotational symmetry of the graphene DSRR, the leaky energy channel between the resonances and the free space continuum can be built. At δ = 0, the two new transmission dips vanish in the transmission spectrum, which implies that the energy leakage channel is closed by the in plane C 2 rotational symmetry. This is the ideal BIC point as marked by red dashed circles in Fig. 2 (a). The two new resonances have finite linewidth induced by the leakage of ideal BICs are two quasi-BICs, and they are named as Q-BIC-I in the low frequency and Q-BIC-II in the high frequency, respectively. To further confirm, the relationship between Q-factors and asymmetry degree α are discussed in Fig. 2 (b). The dimensionless asymmetry degree is defined to be, α = 2 δ / ( l 1 + l 2 ) × 100%. Here, l 1 and l 2 are the total length of upper and lower branches of the graphene DSRR, respectively. When δ = 0.5 µm, and l 2 = l 1 − 2 δ = 14 µm, α equals to 3.45%. The value of Q-factors of Q-BIC-I and Q-BIC-II are 349.75 and 1456.46, respectively. With the decrease in the asymmetry degree α , the increasing Q-factors of modes Q-BIC-I and Q-BIC-II are observed. The fitting solid lines in Fig. 2 (b) indicate the Q-factors follows the asymmetry degree as Q∝ α −2 , and diverge at the ideal BIC points. All the Q-factor values of modes Q-BIC-I and Q-BIC-II discussed above are analyzed by using the following normalized Fano’s equation [ 24 ]: $${T_{Fano}}=a+\frac{{{{(b+c)}^2}}}{{({b^2}+1)(1+{c^2})}}$$ 2 , where a donates the real Fano fitting number. b = 2( f - f 0 )/ Γ tot , where f donates the broadband frequency of incident THz wave, f 0 donates the central plasmonic resonant frequency, and Γ tot is the total loss. c is a non-zero Fano parameter for Q-BICs. Since the transmission dips are due to plasmonic Q-BICs, the Q-factor is given by [ 10 ]: $$Q=\frac{{{f_0}}}{{{\Gamma _{tot}}}}$$ 3 . The match of fitting results via Eq. 2 with the numerical spectra is shown in Fig. 3 (b) by choosing the asymmetry parameter δ = 4 µm. The fitting spectral line (the red solid line) is consistent widely with the simulation curve (the black solid curve). The full transmission spectra of Q-BICs ( δ = 4 µm) in comparison with ideal BIC ( δ = 0) is shown in Fig. 4 (a). The spectral contrast ratio is defined as: $$\frac{{{T_{\hbox{max} }} - {T_{\hbox{min} }}}}{{{T_{\hbox{max} }}+{T_{\hbox{min} }}}} \times 100\%$$ 4 , where T max and T min donate the maximum and minimum values of transmittance. Mode Q-BIC-I at 0.44 THz has a spectral contrast ratio larger than 70%, and that of Q-BIC-II at 1.19 THz is close to 50%. In order to clarify the transition from ideal BICs to Q-BICs more clearly, the electric field distributions and surface current flow distributions at xoy sections at each Q-BIC resonance frequency are presented at Fig. 3 (d) and (e). Under linearly polarized THz plane wave incidence, for DSRRs with C 2 rotational symmetry, only collective modes that have net effective electric dipole moments can be excited due to the symmetry incompatibility [ 25 ]. Anti-phased and equal magnitude currents parallel to the x - and y -axis will result in the zero effective electric dipole moments, which cannot cause the coupling of ideal BIC to the incident THz wave. When the gap of the graphene DSRR is displaced ( δ ≠ 0), the symmetry of surface current flow is broken in terms of magnitude, as shown in Fig. 3 (d). The currents at 0.44THz parallel to the x - and y -axis are anti-phased while the magnitude of current in lower branch is obviously higher than that of upper branch. Therefore, the symmetry of electric dipole moments will be broken not only in terms of magnitude, but also in terms of the orientation, and leads to non-zero net effective electric dipole moments which can cause the coupling of ideal BIC to the incident THz wave. As a result, the ideal BIC transforms into Q-BIC which reveals the characteristic of symmetry protection. Similar trend is observed for mode Q-BIC-II at 1.19 THz. Compared with mode Q-BIC-I which has two branches of current flows, there are four asymmetric branches of current flows are observed in Fig. 3 (e). Here, more nodes of asymmetric electric field distribution and more branches of current distributions also explain why the mode Q-BIC-II can be excited at the higher frequency. The most attractive property of graphene plasmonic resonance is its tunable chemical potential without structure reconfiguration via chemical doping, electrical gating, and ion gel as top gates [ 26 ]. To further expand the Q-BIC’s operating frequency range, the transmission spectra of our Q-BIC metasurface with different chemical potentials are calculated in Fig. 4 . For the graphene DSRRs with asymmetric parameter δ = 4 µm, we change the chemical potential from 0.4 eV to 0.55 eV. With the increase of chemical potential, there is a blueshift of both two Q-BICs resonances. The spectral contrast ratio of mode Q-BIC-I and Q-BIC-II are also increased. At higher chemical potential, graphene becomes more absorbing because of the band filling. Therefore, the total incident wave energy loss due to transmission is less. To proceed, we perform an angle-dependent calculation for our proposed graphene metasurface. The asymmetric parameter is δ = 4 µm and the graphene’s chemical potential is E f = 0.4 eV. The 2D mesh plot in Fig. 5 shows that the mode Q-BIC-I is constantly maintained at 0.44 THz for the incident angle from 0° to 70°, with a spectral contrast ratio larger than 70%. The mode Q-BIC-II is constantly maintained at 0.19 THz as well when changing the incident angle. At higher angles > 60°, its spectral contrast ratio begins to decrease, and the Q-BIC point is slight shift in the spectrum. The operating frequencies of Q-BICs are loosely dependent on incident angle because the geometric dimensions of the graphene DSRR are much smaller than the incident wavelength, thus the phase difference at neighboring graphene DSRR is rather independent of incident angle. Figure 5 shows our dual band Q-BIC metasurface are immune to the variation of incident angle. We also perform the calculation of slow light effect for our proposed Q-BIC metasurface. In this simulation, the asymmetric parameter is δ = 4 µm and the graphene’s chemical potential is E f = 0.4 eV. The slow light effect can be evaluated by group time delay τ g , which is given by: $${\tau _g}= - \frac{{d\phi }}{{d\omega }}$$ 5 , where ф donates the phase change of the transmission spectrum. In Fig. 6 (a), we show the phase oscillates rapidly from − 1.39 rad to -0.42 rad and from − 2.41 rad to -1.98 rad as close to the Q-BIC-I and Q-BIC-II region, respectively. The phase ф changes monotonously like usual plane wave while away from the Q-BIC region. Figure 6 (b) shows the largest positive time delay of 23.8 ps and 14.5 ps, which represent the remarkable slow light effect, can be achieved by the mode Q-BIC-I and Q-BIC-II, respectively. The very slow group velocity when light passing through the metasurface will extremely enhance the light-matter interaction for various applications. Thereafter, we investigate the sensing characteristics of our proposed Q-BIC metasurface. The asymmetric parameter and the graphene’s chemical potential are fixed to δ = 4 µm and E f = 0.4 eV, respectively. We calculate the change of transmission spectrum by varying the refractive index ( n ) of external liquid from 1.33 to 1.37. As shown in Fig. 7 (a) and (c), both two Q-BICs’ resonant frequency has a redshift when increasing the refractive index. As we know that, the important performance index of the refractive index sensor is its sensitivity, which is defined as [ 27 ]: $$S=\frac{{df}}{{dn}}$$ 6 . Figure 7 (b) and (d) show the frequency shifts of the mode Q-BIC-I and Q-BIC-II corresponding to different refractive indexes, respectively. For mode Q-BIC-I in Fig. 7 (b), the curve shows a relatively linear dependence, and the slope is 96 GHz/RIU. Similar trend is observed for mode Q-BIC-II in Fig. 7 (d), while the slope is relative larger and can reach up to 267.5 GHz/RIU. The measurable frequency offset caused by a slight change in the refractive index of external liquid indicates that our proposed Q-BIC metasurface can operate as a THz liquid sensor. The sensing characteristics can also be quantified by calculating the figure of merit (FOM). The FOM is defined as [ 27 ]: $$FOM=\frac{S}{{{\Gamma _{tot}}}}$$ 7 . The maximum FOM is calculated to be 5.87 RIU − 1 and 24.08 RIU − 1 for the mode Q-BIC-I and Q-BIC-II, respectively. In addition, the sensing performance of the proposed graphene metasurface are compared with the other types of recently reported THz sensors as shown in Table 1 . The sensitivities and FOMs of several kinds of sensors are listed. Obviously, our work having high sensitivity is superior to the previous works [ 27 – 32 ] in sensing applications. It is worth noting that the sensitivities of Ref. [ 30 ] and [ 32 ] are higher than that of our proposed Q-BIC metasurface, however, their FOMs are much lower. When the analyte concentration of external liquid is very low, a sensing platform with high FOM is typically desired. Table 1 shows that our proposed Q-BIC metasurface exhibits improved sensing performance over most of the recently reported THz sensors. The simplicity of our tunable design will provide a valuable reference for future dynamic sensor. Table 1 Comparison of sensing performance of the proposed Q-BIC metasurface. Reference Sensitivity (GHz/RIU) FOM (RIU − 1 ) [ 27 ] 170.58 - [ 28 ] 187.1 285.8 [ 29 ] 105 7.501 [ 30 ] 533.7 6.03 [ 31 ] 266 - [ 32 ] 597 9.48 This work 96 267.5 5.87 24.08 4. Conclusions In summary, an all-graphene metasurface capable of simultaneously exciting dual Q-BICs is presented for the first time. The dual ideal BIC points are established by choosing periodic graphene DSRRs on a terahertz transparent SiO 2 substrate. This happens at the asymmetric parameter δ = 0. Dual Q-BICs are observed by varying the asymmetric parameter δ to break the in-plane C 2 rotational symmetry, which manifested as two new resonance dips in the transmission spectrum. We calculate the finite values of two Q-BICs’ Q-factors, which follow the asymmetry degree as Q∝ α −2 and diverge at the ideal BIC points. Furthermore, the Q-BICs’ operating frequency range can be actively expanded via tuning graphene’s chemical potential, and the dual Q-BICs are immune to the variation of incident angle. Meanwhile, a significant slow light effect is reported with time delay up to 23.8 ps. Our proposed Q-BIC metasurface can operate as a THz sensor. A maximum sensitivity of 267.5 GHz/RIU is obtained with a FOM of 24.08 RIU − 1 . We envision that the proposed graphene metasurface empowered by Q-BIC resonances will exhibit applications in actively tunable slow light and sensing devices, and will promote the development of the class of graphene Q-BIC metasurfaces. Declarations Competing interests All authors have no competing interests Authors' contributions Zhong Huang processed the data and wrote the manuscript. Qing chen and Ji Ding give the support of numerical calculation in software. Chaojun Tang, Zhendong Yan and Jing Chen helped with analysis and revised the manuscript. Xiaoxiao cao revised the manuscript. All authors reviewed the paper and approved the final version to be published. Funding This work is financially supported by National Natural Science Foundation of China (11904139), and The Science and Technology Project of Jiangsu Education Department (ZD2020410). Availability of data and materials All the data and materials are in the manuscript. References Von Neumann J, Wigner E P (1929) Über merkwürdige diskrete Eigenwerte. Phys. Z. 30:465-467. Hsu, C W, Zhen, B, Stone, A D, et al (2016) Bound states in the continuum. Nat. Rev. Mater. 1:1-13. https://doi.org/10.1038/natrevmats.2016.48 Fan, K., Shadrivov, I V, Padilla, W J (2019) Dynamic bound states in the continuum. Optica 6:169-173. https://doi.org/10.1364/OPTICA.6.000169 Abujetas D R, Olmos-Trigo J, Sánchez-Gil J A (2022) Sánchez-Gil, Tailoring Accidental Double Bound States in the Continuum in All-Dielectric Metasurfaces. Adv. Opt. Mater. 10:2200301. https://doi.org/10.1002/adom.202200301 Li S, Zhou C, Liu T, et al (2019) Symmetry-protected bound states in the continuum supported by all-dielectric Metasurfaces. Phys. Rev. A 100:063803. https://doi.org/10.1103/PhysRevA.100.063803 Zhong H, Yu Y, Zheng Z, et al (2023) Ultra-low threshold continuous-wave quantum dot mini-BIC lasers. Light: Science & Applications 12:100. https://doi.org/10.1038/s41377-023-01130-5 Abujetas D R, Barreda Á, Moreno F, et al (2021) High-Q transparency band in all-dielectric metasurfaces induced by a quasi bound state in the continuum Laser. Photon. Rev. 15:2000263. https://doi.org/10.1002/lpor.202000263 Liu C, Bai Y, Zhou J, et al (2021) Refractive index sensing by asymmetric dielectric gratings with both bound states in the continuum and guided mode resonances. Opt. Express 29:42978-42988. https://doi.org/10.1364/OE.446937 Koshelev K, Kruk S, Melik-Gaykazyan E, et al (2020) Melik-Gaykazyan, et al., Subwavelength dielectric resonators for nonlinear nanophotonics. Science 367:288-292. https://doi.org/ 10.1126/science.aaz3985 Koshelev K, Lepeshov S, Liu M, et al (2018) Asymmetric metasurfaces with high-Q resonances governed by bound states in the continuum. Phys. Rev. Lett. 121:193903. https://doi.org/10.1103/PhysRevLett.121.193903 Koshelev K, Bogdanov A, Kivshar Y (2019) Meta-optics and bound states in the continuum. Sci. Bull. 64:836-842. https://doi.org/10.1016/j.scib.2018.12.003 Koshelev K, Bogdanov A, Kivshar Y (2020) Engineering with bound states in the continuum. Optics & Photonics News 31:38-45. https://doi.org/10.1364/ OPN.31.1.000038 Sadrieva Z F, Sinev I S, Koshelev K L, et al (2017) Transition from optical bound states in the continuum to leaky resonances: role of substrate and roughness. ACS Photon. 4:723. https://doi.org/10.1021/acsphotonics.6b00860 Otsuji T, Popov V, Ryzhii V (2014) Active graphene plasmonics for terahertz device applications. J. Phys. D: Appl. Phys. 47:094006. https://doi.org/ 10.1088/0022- 3727/47/9/094006 Roy S, Mondal S, Debnath K (2022) Graphene-based chiral metasurface for generation of tunable circular dichroism-design and sensor applications. IEEE Sens. J. 23:285-292. https://doi.org/10.1109/JSEN.2022.3221437 He X, Huang Y, Yang X, et al (2017) Tunable electromagnetically induced transparency based on terahertz graphene metamaterial. RSC advances 7:40321-40326. https://doi.org/10.1039/C7RA06770D Wang W, Yan H, Xiong L, et al (2022) Dual quasi-bound states in the continuum in graphene/dielectric hybrid metamaterial from out-of-plane symmetry protection. Journal of Modern Optics 69:225-231. https://doi.org/10.1080/09500340.2021. 2024612 Li J, Li J, Zheng C, et al (2021) Free switch between bound states in the continuum (BIC) and quasi-BIC supported by graphene-metal terahertz metasurfaces. Carbon 182:506-515. https://doi.org/10.1016/j.carbon.2021.06.037 Wang W, You B (2022) Symmetry-protected bound states in the continuum in graphene nanoribbons. J. Opt. 24:055001. https://doi.org/10.1088/2040-8986/ac6194 Roy S, Mondal S, Debnath K (2023) Symmetric Bound States in the Continuum in an All Graphene Metasurface-Design and Sensor Applications. IEEE Sens. J. 23:8352-8359. https://doi.org/10.1109/JSEN.2023.3255223 Chen J, Zhang M, Gu P, et al (2022) Broadband, wide-incident-angle, and polarization-insensitive high-efficiency absorption of monolayer graphene with nearly 100% modulation depth at communication wavelength. Results Phys. 40 105833. https://doi.org/10.1016/j.rinp.2022.105833 Yan Z, Kong L, Tang C, et al (2022) Ultra-broadband and completely modulated absorption enhancement of monolayer graphene in near-infrared region. Opt. Express 30:34787-34796. https://doi.org/10.1364/OE.470792 Khardikov V V, Iarko E O, Prosvirnin S L (2010) Trapping of light by metal arrays. J. Opt. 12:045102. https://doi.org/10.1088/2040-8978/12/4/045102 Niu J, Zhai Y, Han Q, et al (2021) Resonance-trapped bound states in the continuum in metallic THz metasurfaces. Opt. Lett. 46:162-165. https://doi.org/10.1364 /OL.410791 Cong L, Singh R (2019) Symmetry-protected dual bound states in the continuum in metamaterials. Adv. Opt. Mater. 7:1900383. https://doi.org/10.1002/adom. 201900383 Guo B, Fang L, Zhang B, et al (2011) Graphene doping: a review. Insciences J. 1:80–89. https://doi.org/10.5640/insc.010280 Li Z, Xiang Y, Xu S, et al (2022) Ultrasensitive terahertz sensing in all-dielectric asymmetric metasurfaces based on quasi-BIC. J. Opt. Soc. Amer. B 39: 286-291. http://dx.doi.org/10.1364/josab.442660 Yin W, Shen Z, Li S, et al (2022) THz absorbers with an ultrahigh Q-factor empowered by the quasi-bound states in the continuum for sensing application. Opt. Express 30: 32162-32173. http://dx.doi.org/10.1364/OE.469962 Xie Q, Dong G X, Wang B X, et al (2018) High-Q Fano resonance in terahertz frequency based on an asymmetric metamaterial resonator. Nanoscale Res. Lett. 13. http://dx.doi.org/10.1186/s11671-018-2677-0 Hu S, Liu D, Yang H, et al (2019) Staggered H-shaped metamaterial based on electromagnetically induced transparency effect and its refractive index sensing performance. Opt. Commun. 450: 202-207. http://dx.doi.org/10. 1016/j.optcom.2019.05.038 Cen W, Lang T, Hong Z, et al (2022) Ultrasensitive flexible terahertz plasmonic metasurface sensor based on bound states in the continuum. IEEE Sens. J. 22: 12838-12845. http://dx.doi.org/10.1109/JSEN.2022.3177538 Sadafi M M, Karami H, Hosseini M (2021) A tunable hybrid graphene-metal metamaterial absorber for sensing in the THz regime. Curr. Appl. Phys. 31: 132-140. http://dx.doi.org/10.1016/j.cap.2021.07.020 Additional Declarations No competing interests reported. Cite Share Download PDF Status: Published Journal Publication published 14 Jul, 2023 Read the published version in Plasmonics → Version 1 posted Editorial decision: Accepted 06 Jul, 2023 Reviews received at journal 03 Jul, 2023 Reviewers agreed at journal 28 Jun, 2023 Reviewers invited by journal 27 Jun, 2023 Submission checks completed at journal 26 Jun, 2023 Editor assigned by journal 26 Jun, 2023 First submitted to journal 20 Jun, 2023 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-3087509","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":213367660,"identity":"a7b2b313-cb97-4917-b968-e8572995774a","order_by":0,"name":"Zhong Huang","email":"","orcid":"","institution":"Jiangsu Second Normal University","correspondingAuthor":false,"prefix":"","firstName":"Zhong","middleName":"","lastName":"Huang","suffix":""},{"id":213367661,"identity":"37e8ea2d-83be-411e-a14f-68cd0235465e","order_by":1,"name":"Qing Chen","email":"","orcid":"","institution":"Jiangsu Second Normal University","correspondingAuthor":false,"prefix":"","firstName":"Qing","middleName":"","lastName":"Chen","suffix":""},{"id":213367662,"identity":"ba3ce858-c206-4289-a343-2aec1623d6f3","order_by":2,"name":"Xiaoxiao Cao","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA10lEQVRIiWNgGAWjYBACefbm4x8+GLDxyMs/PkCcFsOeY2mMMyr45Awb0hKItOZGjhkzzxk5Y4YDOQbE6WBsSEt7zNtmltjYcObjjTcMdnK6DQS0sDMcPm44ty0tsZ2xd7PlHIZkY7MDhGxpbEuQeNt2LLGxmXebNA/DgcRthLQwHOYxkOBt+5/YcIznGZFajvGYSfKcYTNmOMPDRpwWwx62ZMMZFWxyhjPYjC3nGBDhF2AMHnwAjkoJ5oc33lTYyRHUggIkeIiMGmQtpOoYBaNgFIyCEQEAAdVEbQAUyQkAAAAASUVORK5CYII=","orcid":"","institution":"Jiangsu Second Normal University","correspondingAuthor":true,"prefix":"","firstName":"Xiaoxiao","middleName":"","lastName":"Cao","suffix":""},{"id":213367663,"identity":"863f13c3-9421-4949-be3e-7d7e28d437a6","order_by":3,"name":"Ji Ding","email":"","orcid":"","institution":"Jiangsu Second Normal University","correspondingAuthor":false,"prefix":"","firstName":"Ji","middleName":"","lastName":"Ding","suffix":""},{"id":213367664,"identity":"9156429b-9f70-42ea-9169-307a6673468f","order_by":4,"name":"Chaojun Tang","email":"","orcid":"","institution":"Zhejiang University of Technology","correspondingAuthor":false,"prefix":"","firstName":"Chaojun","middleName":"","lastName":"Tang","suffix":""},{"id":213367665,"identity":"35d6696b-12ef-427c-b789-c9933fab9b17","order_by":5,"name":"Zhendong Yan","email":"","orcid":"","institution":"Nanjing Forestry University","correspondingAuthor":false,"prefix":"","firstName":"Zhendong","middleName":"","lastName":"Yan","suffix":""},{"id":213367666,"identity":"f17c6538-9cae-45f4-92bb-11acb7da290e","order_by":6,"name":"Jing Chen","email":"","orcid":"","institution":"Nanjing University of Posts and Telecommunications","correspondingAuthor":false,"prefix":"","firstName":"Jing","middleName":"","lastName":"Chen","suffix":""}],"badges":[],"createdAt":"2023-06-20 13:14:21","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-3087509/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-3087509/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s11468-023-01942-6","type":"published","date":"2023-07-15T01:08:24+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":39328583,"identity":"fc0c0472-fa22-4bf1-8f2f-86cc5370407b","added_by":"auto","created_at":"2023-06-29 20:58:48","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":22820,"visible":true,"origin":"","legend":"\u003cp\u003e(a) Schematic diagram of proposed graphene dual Q-BICs metasurface. (b) Top view of a unit cell of the graphene metasurface.\u003c/p\u003e","description":"","filename":"1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-3087509/v1/8b0c8934e833a0f90d22a52a.jpg"},{"id":39328422,"identity":"e748e044-e8b1-44f9-800b-ee653aa80a1e","added_by":"auto","created_at":"2023-06-29 20:50:48","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":41837,"visible":true,"origin":"","legend":"\u003cp\u003e(a) Transmission spectra of the periodic graphene DSRRs with different asymmetric parameters \u003cem\u003eδ\u003c/em\u003e. (b) The relationship between the asymmetry degree \u003cem\u003eα\u003c/em\u003e and Q-factors of Q-BIC-I and II.\u003c/p\u003e","description":"","filename":"2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-3087509/v1/70500f8b1758061c9394bda0.jpg"},{"id":39327581,"identity":"c9654a23-34d8-48f9-a5cd-aca8e8442ff1","added_by":"auto","created_at":"2023-06-29 20:42:48","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":54404,"visible":true,"origin":"","legend":"\u003cp\u003e(a) Analysis between symmetric and asymmetric graphene DSRR (\u003cem\u003eδ\u003c/em\u003e = 0 and \u003cem\u003eδ\u003c/em\u003e = 4 μm) transmission spectra. (b) Fano fitting curves of Q-BIC-I and Q-BIC-II with asymmetric parameter \u003cem\u003eδ\u003c/em\u003e = 4 μm. (c) Electric field distribution and surface current flow distribution (the white arrows) at 0.59 THz with asymmetric parameter \u003cem\u003eδ\u003c/em\u003e= 0. (d) and (e) Electric field distribution and current flow distribution (the white arrows) of mode Q-BIC-I at 0.44 THz and Q-BIC-II at 1.19 THz with asymmetric parameter \u003cem\u003eδ\u003c/em\u003e= 4 μm, respectively.\u003c/p\u003e","description":"","filename":"3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-3087509/v1/6d56a0d5fb36602d433a281f.jpg"},{"id":39327582,"identity":"8e452755-2aef-44a4-acb8-42f160627de9","added_by":"auto","created_at":"2023-06-29 20:42:48","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":28910,"visible":true,"origin":"","legend":"\u003cp\u003eTransmission spectra of Q-BICs with various graphene’s chemical potentials from 0.4 eV to 0.55 eV.\u003c/p\u003e","description":"","filename":"4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-3087509/v1/e0d8f16f950262d9cfd9fd51.jpg"},{"id":39327577,"identity":"5488268e-3d01-418f-b787-dcea3e231e15","added_by":"auto","created_at":"2023-06-29 20:42:48","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":19797,"visible":true,"origin":"","legend":"\u003cp\u003eTransmission spectra as a function of the frequency and the incident angle.\u003c/p\u003e","description":"","filename":"5.jpg","url":"https://assets-eu.researchsquare.com/files/rs-3087509/v1/886247f5de7a14c8728b27f5.jpg"},{"id":39327580,"identity":"d94d305a-79de-4605-9955-70851e839a66","added_by":"auto","created_at":"2023-06-29 20:42:48","extension":"jpg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":37433,"visible":true,"origin":"","legend":"\u003cp\u003e(a) The phase of the Transmission. (b) The corresponding time delay.\u003c/p\u003e","description":"","filename":"6.jpg","url":"https://assets-eu.researchsquare.com/files/rs-3087509/v1/ca8965e778f37943253906f5.jpg"},{"id":39328423,"identity":"61d34085-4649-452c-a943-6958e1d071b6","added_by":"auto","created_at":"2023-06-29 20:50:48","extension":"jpg","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":77307,"visible":true,"origin":"","legend":"\u003cp\u003eTransmission spectra at different refractive indices of (a) mode Q-BIC-I and (c) mode Q-BIC-II, respectively. The relation between frequency shift and refractive index of (b) mode Q-BIC-I and (d) mode Q-BIC-II, respectively.\u003c/p\u003e","description":"","filename":"7.jpg","url":"https://assets-eu.researchsquare.com/files/rs-3087509/v1/7055d82c83f57e53c12389a2.jpg"},{"id":41736920,"identity":"7f939237-0faf-41dd-ac0d-bf8cb8d8c68b","added_by":"auto","created_at":"2023-08-18 04:13:28","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":466690,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-3087509/v1/0ef3b102-1dc5-4a27-8b19-e5d9a615b9c7.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"Excitation of tunable dual quasi-bound states in the continuum in graphene metasurface and terahertz sensing application","fulltext":[{"header":"1. Introduction","content":"\u003cp\u003eIn 1929, the idea of bound states in the continuum (BIC) was firstly reported by J. V. Neumann and E. Wigner in quantum mechanics during the process of unique Schr\u0026ouml;dinger equation calculation [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]. The BIC, however, is often brittle, necessitates uniquely tuned potentials, and is hence challenging to observe in the quantum experiments. Because the Maxwell equation and the Schrodinger equation are formal equals, the emergence of these exotic non-radiating states in optics, is strongly influenced by the quantum mechanical system [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]. Similar to quantum systems, the concept of BIC in optics is the result of two or more waves interfering destructively to reduce radiative loss, where radiation resonances can generally carry power from or to infinity. At the operating frequency of BIC, the embedded eigenvalues vanish with the leakage power and cause eigenstates having an infinitely high-quality factor (Q-factor) [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. Generally, there are two ways to excite BIC resonance. One is tuning the parameters of optic systems to make outgoing waves interfering destructively, and this type of BIC is named as accidental BIC [\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e]. The other way is breaking the in-plane symmetry to generate mismatch between the symmetry of outgoing wave and the resonance mode, and this type is symmetry protected BIC [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. Due to its novel ability to have the high Q-factor, BIC has been widely studied in various optical systems of quantum dots lasers [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e], filters [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e], sensors [\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e] and nonlinear optics [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eRecently, BIC has been introduced into the metasurfaces [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e] owing to the advantage of electromagnetic wave confinement in sub-wavelength, precise wavefront control, and the simplicity of fabrication in two-dimensional planes [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. The periodic arrangements of artificial sub-wavelength unit cells in metasurfaces result in some diffraction open channels. The symmetry protected BIC is easier to be existed in the absence of coupling to these open channels by engineering the periodic sub-wavelength structures [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e]. Interestingly, true optical BIC (ideal BIC) is a mathematically defined concept with the infinite value of Q-factor and its resonance linewidth infinitely tends towards 0. This is only possible in an ideal lossless infinite optical system with extreme structure parameters. In practice, suffering from material absorption, finite lateral size of structures and other leakage, the observable BIC often refers to quasi-BIC (Q-BIC) with high, yet finite Q-factor [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. Although Q-BIC has received extensive attention in metasurfaces, most of the devices reported in the literatures are lack of dynamic capability. Once the manufacturing parameters are fixed, the operating frequency of Q-BIC is synchronous fixed. So the ability to expand the Q-BIC\u0026rsquo;s operating frequency range without changing the geometrical parameters of the structure is very important and attractive.\u003c/p\u003e \u003cp\u003eCurrently, graphene has exhibited powerful tuning ability for active devices [\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e, \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e]. By controlling the chemical potential \u003cem\u003eE\u003c/em\u003e\u003csub\u003e\u003cem\u003ef\u003c/em\u003e\u003c/sub\u003e of the graphene, the target operating frequency region can be dynamically shifted [\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. However, the active tuning of Q-BIC reported in the literatures are mainly focused on the switching characteristics controlled by graphene absorption. For example, the authors in ref. [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e] proposed dual Q-BICs in a silicon stair-like metasurface, where the transmittance can be controlled with respect to a change in the chemical potential of graphene layer. Similarly, the authors in ref. [\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e] have reported graphene-metal hybrid terahertz metasurfaces with a similar kind of control. So far, to the best knowledge of us, only two groups have reported dynamic tuning the operating frequency of Q-BIC via the changing of the chemical potential \u003cem\u003eE\u003c/em\u003e\u003csub\u003e\u003cem\u003ef\u003c/em\u003e\u003c/sub\u003e of the graphene. Wang \u003cem\u003eet al.\u003c/em\u003e have investigated the symmetry-protected Q-BIC by using the structure of two contacting graphene nanoribbons [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e], and Roy \u003cem\u003eet al.\u003c/em\u003e have proposed rectangular dimer-shaped graphene metasurfaces that supported asymmetric Q-BIC [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e]. Therefore, the tunable aspects of graphene in Q-BIC are still lack of research. Besides that, tunable Q-BIC with only one operating band also faces time delay introduced by modulation, especially in sensing applications in detecting a variety of gas and chemical solvents.\u003c/p\u003e \u003cp\u003eIn this article, the exciting of symmetry-protected dual Q-BICs in an all graphene-based metasurface is theoretically demonstrated for the first time. The metasurface is made of a square lattice of which unit cell consists of a double-gap split ring resonator (DSRR). By displacing the two gaps from the center, the in-plane \u003cem\u003eC\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e rotational symmetry is broken, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e. This symmetry breaking provides an efficient way to transform BICs into Q-BICs, manifested as two new sharp dips in the transmission spectrum. The Q factor of dual Q-BICs can be well evaluated by the theoretical normalized Fano line shape equation, and follow the inverse square dependence on the asymmetry degree. The dual Q-BICs are further analyzed by the near-field electric field distributions and surface current flows, indicating the leakage of an ideal bound state can be caused by a slight perturbation on the effective electric dipole moments. We show the operating frequency of the dual Q-BICs can be actively tuned by varying the chemical potential \u003cem\u003eE\u003c/em\u003e\u003csub\u003e\u003cem\u003ef\u003c/em\u003e\u003c/sub\u003e of the graphene. Our Q-BIC metasurface is immune to the variation of incident angle, and has a significant slow light effect with time delay up to 23.8 ps. Moreover, our Q-BIC metasurface can operate as a refractive index sensor in the THz region. A maximum sensitivity of 267.5 GHz/RIU is obtained with a FOM of 24.08 RIU\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e. Our work shows an alternative way to design the class of tunable Q-BIC metasurface, which will provide a valuable reference for future dynamic sensor and other fields.\u003c/p\u003e"},{"header":"2. Design and simulation of dual Q-BICs metasurface","content":"\u003cp\u003e \u003c/p\u003e \u003cp\u003eThe schematic of the proposed Q-BIC metasurface is shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e(a). On a thick SiO\u003csub\u003e2\u003c/sub\u003e substrate, a monolayer graphene is patterned into periodic DSRRs along \u003cem\u003ex\u003c/em\u003e and \u003cem\u003ey\u003c/em\u003e directions. The relative permittivity of the THz transparent SiO\u003csub\u003e2\u003c/sub\u003e substrate is taken as 3.9. Figure\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e(b) shows the unit cell structure. Here, the period along \u003cem\u003ex\u003c/em\u003e and \u003cem\u003ey\u003c/em\u003e directions are \u003cem\u003eP\u003c/em\u003e\u003csub\u003e\u003cem\u003ex\u003c/em\u003e\u003c/sub\u003e and \u003cem\u003eP\u003c/em\u003e\u003csub\u003e\u003cem\u003ey\u003c/em\u003e\u003c/sub\u003e, and \u003cem\u003eP\u003c/em\u003e\u003csub\u003e\u003cem\u003ex\u003c/em\u003e\u003c/sub\u003e = \u003cem\u003eP\u003c/em\u003e\u003csub\u003e\u003cem\u003ey\u003c/em\u003e\u003c/sub\u003e = 50 \u0026micro;m. The arm length of the DSSR is \u003cem\u003el\u003c/em\u003e\u0026thinsp;=\u0026thinsp;31 \u0026micro;m, and the arm width is \u003cem\u003ew\u003c/em\u003e\u0026thinsp;=\u0026thinsp;4 \u0026micro;m. The gap length \u003cem\u003eg\u003c/em\u003e\u0026thinsp;=\u0026thinsp;2 \u0026micro;m. In order to break the in plane \u003cem\u003eC\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e rotational symmetry of the DSSR, we displace both two gaps from the center by a distance \u0026ldquo;\u003cem\u003eδ\u003c/em\u003e\u0026rdquo;. In this work, we use open souce software \u0026ldquo;East FDTD\u0026rdquo; [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e, \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e] based on the finite-difference time-domain method to design our structure and calculate the transmission spectrum numerically. Along the \u003cem\u003ex\u003c/em\u003e- and \u003cem\u003ey\u003c/em\u003e-axis, we choose periodic boundary conditions to simulate a system of infinite size. Along the \u003cem\u003ez\u003c/em\u003e-axis direction, the matesurface is excited by a normally incident broadband terahertz plane wave polarizing along \u003cem\u003ex\u003c/em\u003e-direction, and two perfectly matched layers (PMLs) behind Floquet periodic ports are used to absorb the scattered light. For the numerical calculations, graphene is treated as a surface conductive layer with the thickness of 1 nm, and its surface conductivity is defined by analytical Kubo\u0026rsquo;s formula [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e, \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]:\u003cdiv id=\"Equ1\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ1\" name=\"EquationSource\"\u003e\n$${\\sigma _g}=\\frac{{2{e^2}{k_B}T}}{{\\pi {\\hbar ^2}}}Ln\\left[ {2\\cosh (\\frac{{{E_f}}}{{2{k_B}T}})} \\right]\\frac{i}{{\\omega +i{\\tau ^{ - 1}}}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e1\u003c/div\u003e\u003c/div\u003e,\u003c/p\u003e \u003cp\u003ewhere \u003cem\u003eω\u003c/em\u003e represents the frequency of incident terahertz wave, \u003cem\u003eT\u003c/em\u003e represents the environmental temperature, \u003cem\u003ee\u003c/em\u003e represents an electron charge, \u003cem\u003ek\u003c/em\u003e\u003csub\u003e\u003cem\u003eB\u003c/em\u003e\u003c/sub\u003e represents the Boltzmann's constant, \u003cem\u003eħ\u003c/em\u003e represents the reduced Planck constant, and \u003cem\u003eE\u003c/em\u003e\u003csub\u003e\u003cem\u003ef\u003c/em\u003e\u003c/sub\u003e represents the chemical potential of graphene DSSRs. In our numerical calculations, room temperature is taken as \u003cem\u003eT\u003c/em\u003e\u0026thinsp;=\u0026thinsp;300 K, and the relaxation time of graphene DSSRs is taken as \u0026#120591; = 2 ps.\u003c/p\u003e"},{"header":"3. Results and discussions","content":"\u003cp\u003e \u003c/p\u003e \u003cp\u003eWe start with the investigation of transmission spectra from Q-BIC metasurface with in plane \u003cem\u003eC\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e rotational symmetry, which two split gaps of the graphene DSRR are located at the center of the left and right arms and corresponding to asymmetric parameter \u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0. One can see from the black line in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e(a) that the graphene DSRRs show an obvious plasmonic resonance at 0.59 THz when excited by the incident THz wave along \u003cem\u003ex\u003c/em\u003e-direction. The transmittance at this frequency point is 0.1. In this simulation, the chemical potential of the graphene is \u003cem\u003eE\u003c/em\u003e\u003csub\u003e\u003cem\u003ef\u003c/em\u003e\u003c/sub\u003e = 0.4 eV. The form of the electric field distribution at the frequency of the minimum value of transmittance demonstrates the trapped mode nature of the resonance, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e(c). The electric field of plasmonic resonance mode is mostly concentrated in the left and right branches of the graphene DSRR. There are the in-phased surface current flows parallel to the \u003cem\u003ex-\u003c/em\u003e and \u003cem\u003ey\u003c/em\u003e-axis of the DSRR elements (white arrows in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e(c)). The amplitude of the surface current flow parallel to the \u003cem\u003ex\u003c/em\u003e-axis is greater than that parallel to the \u003cem\u003ey\u003c/em\u003e-axis. This amplitudes distribution of currents provides a strong coupling between the trapped mode and the incident THz wave [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eWe then plot the transmission spectrum of the periodic graphene DSRRs with distinct asymmetric parameters \u003cem\u003eδ\u003c/em\u003e in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e(a). Besides the plasmonic resonance at 0.59 THz, there are two new transmission dips appear with the introduction of asymmetric parameter \u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.5 \u0026micro;m. The exhibited two new narrow resonances are strongly altered with increasing \u003cem\u003eδ\u003c/em\u003e from 0 to 4 \u0026micro;m. The two new transmission dips become prominent with the increase of \u003cem\u003eδ\u003c/em\u003e, and the linewidth gradually increase. Meantime, the two new transmission dips are slightly red shifted. As the non-zero value of asymmetric parameters \u003cem\u003eδ\u003c/em\u003e breaks the \u003cem\u003eC\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e rotational symmetry of the graphene DSRR, the leaky energy channel between the resonances and the free space continuum can be built. At \u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0, the two new transmission dips vanish in the transmission spectrum, which implies that the energy leakage channel is closed by the in plane \u003cem\u003eC\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e rotational symmetry. This is the ideal BIC point as marked by red dashed circles in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e(a). The two new resonances have finite linewidth induced by the leakage of ideal BICs are two quasi-BICs, and they are named as Q-BIC-I in the low frequency and Q-BIC-II in the high frequency, respectively. To further confirm, the relationship between Q-factors and asymmetry degree \u003cem\u003eα\u003c/em\u003e are discussed in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e(b). The dimensionless asymmetry degree is defined to be, \u003cem\u003eα\u003c/em\u003e\u0026thinsp;=\u0026thinsp;2\u003cem\u003eδ\u003c/em\u003e / (\u003cem\u003el\u003c/em\u003e\u003csub\u003e\u003cem\u003e1\u003c/em\u003e\u003c/sub\u003e\u0026thinsp;+\u0026thinsp;\u003cem\u003el\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e) \u0026times; 100%. Here, \u003cem\u003el\u003c/em\u003e\u003csub\u003e\u003cem\u003e1\u003c/em\u003e\u003c/sub\u003e and \u003cem\u003el\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e are the total length of upper and lower branches of the graphene DSRR, respectively. When \u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.5 \u0026micro;m, and \u003cem\u003el\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e\u0026thinsp;=\u0026thinsp;\u003cem\u003el\u003c/em\u003e\u003csub\u003e\u003cem\u003e1\u003c/em\u003e\u003c/sub\u003e \u0026minus;\u0026thinsp;2\u003cem\u003eδ\u0026thinsp;=\u003c/em\u003e\u0026thinsp;14 \u0026micro;m, \u003cem\u003eα\u003c/em\u003e equals to 3.45%. The value of Q-factors of Q-BIC-I and Q-BIC-II are 349.75 and 1456.46, respectively. With the decrease in the asymmetry degree \u003cem\u003eα\u003c/em\u003e, the increasing Q-factors of modes Q-BIC-I and Q-BIC-II are observed. The fitting solid lines in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e(b) indicate the Q-factors follows the asymmetry degree as Q\u0026prop;\u003cem\u003eα\u003c/em\u003e\u003csup\u003e\u0026minus;2\u003c/sup\u003e, and diverge at the ideal BIC points.\u003c/p\u003e \u003cp\u003eAll the Q-factor values of modes Q-BIC-I and Q-BIC-II discussed above are analyzed by using the following normalized Fano\u0026rsquo;s equation [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]:\u003cdiv id=\"Equ2\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ2\" name=\"EquationSource\"\u003e\n$${T_{Fano}}=a+\\frac{{{{(b+c)}^2}}}{{({b^2}+1)(1+{c^2})}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e2\u003c/div\u003e\u003c/div\u003e,\u003c/p\u003e \u003cp\u003ewhere \u003cem\u003ea\u003c/em\u003e donates the real Fano fitting number. \u003cem\u003eb\u003c/em\u003e\u0026thinsp;=\u0026thinsp;2(\u003cem\u003ef\u003c/em\u003e - \u003cem\u003ef\u003c/em\u003e\u003csub\u003e\u003cem\u003e0\u003c/em\u003e\u003c/sub\u003e)/\u003cem\u003eΓ\u003c/em\u003e\u003csub\u003e\u003cem\u003etot\u003c/em\u003e\u003c/sub\u003e, where \u003cem\u003ef\u003c/em\u003e donates the broadband frequency of incident THz wave, \u003cem\u003ef\u003c/em\u003e\u003csub\u003e\u003cem\u003e0\u003c/em\u003e\u003c/sub\u003e donates the central plasmonic resonant frequency, and \u003cem\u003eΓ\u003c/em\u003e\u003csub\u003e\u003cem\u003etot\u003c/em\u003e\u003c/sub\u003e is the total loss. \u003cem\u003ec\u003c/em\u003e is a non-zero Fano parameter for Q-BICs. Since the transmission dips are due to plasmonic Q-BICs, the Q-factor is given by [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]:\u003cdiv id=\"Equ3\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ3\" name=\"EquationSource\"\u003e\n$$Q=\\frac{{{f_0}}}{{{\\Gamma _{tot}}}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e3\u003c/div\u003e\u003c/div\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe match of fitting results via Eq.\u0026nbsp;\u003cspan refid=\"Equ2\" class=\"InternalRef\"\u003e2\u003c/span\u003e with the numerical spectra is shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e(b) by choosing the asymmetry parameter \u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;4 \u0026micro;m. The fitting spectral line (the red solid line) is consistent widely with the simulation curve (the black solid curve). The full transmission spectra of Q-BICs (\u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;4 \u0026micro;m) in comparison with ideal BIC (\u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0) is shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e(a). The spectral contrast ratio is defined as:\u003cdiv id=\"Equ4\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ4\" name=\"EquationSource\"\u003e\n$$\\frac{{{T_{\\hbox{max} }} - {T_{\\hbox{min} }}}}{{{T_{\\hbox{max} }}+{T_{\\hbox{min} }}}} \\times 100\\%$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e4\u003c/div\u003e\u003c/div\u003e,\u003c/p\u003e \u003cp\u003ewhere \u003cem\u003eT\u003c/em\u003e\u003csub\u003e\u003cem\u003emax\u003c/em\u003e\u003c/sub\u003e and \u003cem\u003eT\u003c/em\u003e\u003csub\u003e\u003cem\u003emin\u003c/em\u003e\u003c/sub\u003e donate the maximum and minimum values of transmittance. Mode Q-BIC-I at 0.44 THz has a spectral contrast ratio larger than 70%, and that of Q-BIC-II at 1.19 THz is close to 50%.\u003c/p\u003e \u003cp\u003eIn order to clarify the transition from ideal BICs to Q-BICs more clearly, the electric field distributions and surface current flow distributions at \u003cem\u003exoy\u003c/em\u003e sections at each Q-BIC resonance frequency are presented at Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e(d) and (e). Under linearly polarized THz plane wave incidence, for DSRRs with \u003cem\u003eC\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e rotational symmetry, only collective modes that have net effective electric dipole moments can be excited due to the symmetry incompatibility [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e]. Anti-phased and equal magnitude currents parallel to the \u003cem\u003ex\u003c/em\u003e- and \u003cem\u003ey\u003c/em\u003e-axis will result in the zero effective electric dipole moments, which cannot cause the coupling of ideal BIC to the incident THz wave. When the gap of the graphene DSRR is displaced (\u003cem\u003eδ\u003c/em\u003e\u0026thinsp;\u0026ne;\u0026thinsp;0), the symmetry of surface current flow is broken in terms of magnitude, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e(d). The currents at 0.44THz parallel to the \u003cem\u003ex\u003c/em\u003e- and \u003cem\u003ey\u003c/em\u003e-axis are anti-phased while the magnitude of current in lower branch is obviously higher than that of upper branch. Therefore, the symmetry of electric dipole moments will be broken not only in terms of magnitude, but also in terms of the orientation, and leads to non-zero net effective electric dipole moments which can cause the coupling of ideal BIC to the incident THz wave. As a result, the ideal BIC transforms into Q-BIC which reveals the characteristic of symmetry protection. Similar trend is observed for mode Q-BIC-II at 1.19 THz. Compared with mode Q-BIC-I which has two branches of current flows, there are four asymmetric branches of current flows are observed in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e(e). Here, more nodes of asymmetric electric field distribution and more branches of current distributions also explain why the mode Q-BIC-II can be excited at the higher frequency.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe most attractive property of graphene plasmonic resonance is its tunable chemical potential without structure reconfiguration via chemical doping, electrical gating, and ion gel as top gates [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e]. To further expand the Q-BIC\u0026rsquo;s operating frequency range, the transmission spectra of our Q-BIC metasurface with different chemical potentials are calculated in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e. For the graphene DSRRs with asymmetric parameter \u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;4 \u0026micro;m, we change the chemical potential from 0.4 eV to 0.55 eV. With the increase of chemical potential, there is a blueshift of both two Q-BICs resonances. The spectral contrast ratio of mode Q-BIC-I and Q-BIC-II are also increased. At higher chemical potential, graphene becomes more absorbing because of the band filling. Therefore, the total incident wave energy loss due to transmission is less.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eTo proceed, we perform an angle-dependent calculation for our proposed graphene metasurface. The asymmetric parameter is \u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;4 \u0026micro;m and the graphene\u0026rsquo;s chemical potential is \u003cem\u003eE\u003c/em\u003e\u003csub\u003e\u003cem\u003ef\u003c/em\u003e\u003c/sub\u003e = 0.4 eV. The 2D mesh plot in Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e shows that the mode Q-BIC-I is constantly maintained at 0.44 THz for the incident angle from 0\u0026deg; to 70\u0026deg;, with a spectral contrast ratio larger than 70%. The mode Q-BIC-II is constantly maintained at 0.19 THz as well when changing the incident angle. At higher angles\u0026thinsp;\u0026gt;\u0026thinsp;60\u0026deg;, its spectral contrast ratio begins to decrease, and the Q-BIC point is slight shift in the spectrum. The operating frequencies of Q-BICs are loosely dependent on incident angle because the geometric dimensions of the graphene DSRR are much smaller than the incident wavelength, thus the phase difference at neighboring graphene DSRR is rather independent of incident angle. Figure\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e shows our dual band Q-BIC metasurface are immune to the variation of incident angle.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eWe also perform the calculation of slow light effect for our proposed Q-BIC metasurface. In this simulation, the asymmetric parameter is \u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;4 \u0026micro;m and the graphene\u0026rsquo;s chemical potential is \u003cem\u003eE\u003c/em\u003e\u003csub\u003e\u003cem\u003ef\u003c/em\u003e\u003c/sub\u003e = 0.4 eV. The slow light effect can be evaluated by group time delay \u003cem\u003eτ\u003c/em\u003e\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e, which is given by:\u003cdiv id=\"Equ5\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ5\" name=\"EquationSource\"\u003e\n$${\\tau _g}= - \\frac{{d\\phi }}{{d\\omega }}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e5\u003c/div\u003e\u003c/div\u003e,\u003c/p\u003e \u003cp\u003ewhere \u003cem\u003eф\u003c/em\u003e donates the phase change of the transmission spectrum. In Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e(a), we show the phase oscillates rapidly from \u0026minus;\u0026thinsp;1.39 rad to -0.42 rad and from \u0026minus;\u0026thinsp;2.41 rad to -1.98 rad as close to the Q-BIC-I and Q-BIC-II region, respectively. The phase \u003cem\u003eф\u003c/em\u003e changes monotonously like usual plane wave while away from the Q-BIC region. Figure\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e(b) shows the largest positive time delay of 23.8 ps and 14.5 ps, which represent the remarkable slow light effect, can be achieved by the mode Q-BIC-I and Q-BIC-II, respectively. The very slow group velocity when light passing through the metasurface will extremely enhance the light-matter interaction for various applications.\u003c/p\u003e \u003cp\u003eThereafter, we investigate the sensing characteristics of our proposed Q-BIC metasurface. The asymmetric parameter and the graphene\u0026rsquo;s chemical potential are fixed to \u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;4 \u0026micro;m and \u003cem\u003eE\u003c/em\u003e\u003csub\u003e\u003cem\u003ef\u003c/em\u003e\u003c/sub\u003e = 0.4 eV, respectively. We calculate the change of transmission spectrum by varying the refractive index (\u003cem\u003en\u003c/em\u003e) of external liquid from 1.33 to 1.37. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e(a) and (c), both two Q-BICs\u0026rsquo; resonant frequency has a redshift when increasing the refractive index. As we know that, the important performance index of the refractive index sensor is its sensitivity, which is defined as [\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e]:\u003cdiv id=\"Equ6\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ6\" name=\"EquationSource\"\u003e\n$$S=\\frac{{df}}{{dn}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e6\u003c/div\u003e\u003c/div\u003e.\u003c/p\u003e \u003cp\u003eFigure\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e(b) and (d) show the frequency shifts of the mode Q-BIC-I and Q-BIC-II corresponding to different refractive indexes, respectively. For mode Q-BIC-I in Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e(b), the curve shows a relatively linear dependence, and the slope is 96 GHz/RIU. Similar trend is observed for mode Q-BIC-II in Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e(d), while the slope is relative larger and can reach up to 267.5 GHz/RIU. The measurable frequency offset caused by a slight change in the refractive index of external liquid indicates that our proposed Q-BIC metasurface can operate as a THz liquid sensor. The sensing characteristics can also be quantified by calculating the figure of merit (FOM). The FOM is defined as [\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e]:\u003cdiv id=\"Equ7\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ7\" name=\"EquationSource\"\u003e\n$$FOM=\\frac{S}{{{\\Gamma _{tot}}}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e7\u003c/div\u003e\u003c/div\u003e.\u003c/p\u003e \u003cp\u003eThe maximum FOM is calculated to be 5.87 RIU\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e and 24.08 RIU\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e for the mode Q-BIC-I and Q-BIC-II, respectively.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn addition, the sensing performance of the proposed graphene metasurface are compared with the other types of recently reported THz sensors as shown in Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e. The sensitivities and FOMs of several kinds of sensors are listed. Obviously, our work having high sensitivity is superior to the previous works [\u003cspan additionalcitationids=\"CR28 CR29 CR30 CR31\" citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e] in sensing applications. It is worth noting that the sensitivities of Ref. [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e] and [\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e] are higher than that of our proposed Q-BIC metasurface, however, their FOMs are much lower. When the analyte concentration of external liquid is very low, a sensing platform with high FOM is typically desired. Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e shows that our proposed Q-BIC metasurface exhibits improved sensing performance over most of the recently reported THz sensors. The simplicity of our tunable design will provide a valuable reference for future dynamic sensor.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eComparison of sensing performance of the proposed Q-BIC metasurface.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"3\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eReference\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSensitivity (GHz/RIU)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eFOM (RIU\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e[\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e]\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e170.58\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e[\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e]\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e187.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e285.8\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e[\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e]\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e105\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e7.501\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e[\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e]\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e533.7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e6.03\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e[\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e]\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e266\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e[\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e]\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e597\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e9.48\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eThis work\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e96\u003c/p\u003e \u003cp\u003e267.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e5.87\u003c/p\u003e \u003cp\u003e24.08\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e"},{"header":"4. Conclusions","content":"\u003cp\u003eIn summary, an all-graphene metasurface capable of simultaneously exciting dual Q-BICs is presented for the first time. The dual ideal BIC points are established by choosing periodic graphene DSRRs on a terahertz transparent SiO\u003csub\u003e2\u003c/sub\u003e substrate. This happens at the asymmetric parameter \u003cem\u003eδ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0. Dual Q-BICs are observed by varying the asymmetric parameter \u003cem\u003eδ\u003c/em\u003e to break the in-plane \u003cem\u003eC\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e rotational symmetry, which manifested as two new resonance dips in the transmission spectrum. We calculate the finite values of two Q-BICs\u0026rsquo; Q-factors, which follow the asymmetry degree as Q\u0026prop;\u003cem\u003eα\u003c/em\u003e\u003csup\u003e\u0026minus;2\u003c/sup\u003e and diverge at the ideal BIC points. Furthermore, the Q-BICs\u0026rsquo; operating frequency range can be actively expanded via tuning graphene\u0026rsquo;s chemical potential, and the dual Q-BICs are immune to the variation of incident angle. Meanwhile, a significant slow light effect is reported with time delay up to 23.8 ps. Our proposed Q-BIC metasurface can operate as a THz sensor. A maximum sensitivity of 267.5 GHz/RIU is obtained with a FOM of 24.08 RIU\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e. We envision that the proposed graphene metasurface empowered by Q-BIC resonances will exhibit applications in actively tunable slow light and sensing devices, and will promote the development of the class of graphene Q-BIC metasurfaces.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll authors have no competing interests\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026apos; contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eZhong Huang processed the data and wrote the manuscript. Qing chen and Ji Ding give the support of numerical calculation in software. Chaojun Tang, Zhendong Yan and Jing Chen helped with analysis and revised the manuscript. Xiaoxiao cao revised the manuscript. All authors reviewed the paper and approved the final version to be published.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work is financially supported by National Natural Science Foundation of China (11904139), and The Science and Technology Project of Jiangsu Education Department (ZD2020410).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll the data and materials are in the manuscript.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eVon Neumann J, Wigner E P (1929) \u0026Uuml;ber merkw\u0026uuml;rdige diskrete Eigenwerte. Phys. Z. 30:465-467.\u003c/li\u003e\n\u003cli\u003eHsu, C W, Zhen, B, Stone, A D, et al (2016) Bound states in the continuum. Nat. Rev. Mater. 1:1-13. https://doi.org/10.1038/natrevmats.2016.48\u003c/li\u003e\n\u003cli\u003eFan, K., Shadrivov, I V, Padilla, W J (2019) Dynamic bound states in the continuum. Optica 6:169-173. https://doi.org/10.1364/OPTICA.6.000169\u003c/li\u003e\n\u003cli\u003eAbujetas D R, Olmos-Trigo J, S\u0026aacute;nchez-Gil J A (2022) S\u0026aacute;nchez-Gil, Tailoring Accidental Double Bound States in the Continuum in All-Dielectric Metasurfaces. Adv. Opt. Mater. 10:2200301. https://doi.org/10.1002/adom.202200301\u003c/li\u003e\n\u003cli\u003eLi S, Zhou C, Liu T, et al (2019) Symmetry-protected bound states in the continuum supported by all-dielectric Metasurfaces. Phys. Rev. A 100:063803. https://doi.org/10.1103/PhysRevA.100.063803\u003c/li\u003e\n\u003cli\u003eZhong H, Yu Y, Zheng Z, et al (2023) Ultra-low threshold continuous-wave quantum dot mini-BIC lasers. Light: Science \u0026amp; Applications 12:100. https://doi.org/10.1038/s41377-023-01130-5\u003c/li\u003e\n\u003cli\u003eAbujetas D R, Barreda \u0026Aacute;, Moreno F, et al (2021) High-Q transparency band in all-dielectric metasurfaces induced by a quasi bound state in the continuum Laser. Photon. Rev. 15:2000263. https://doi.org/10.1002/lpor.202000263\u003c/li\u003e\n\u003cli\u003eLiu C, Bai Y, Zhou J, et al (2021) Refractive index sensing by asymmetric dielectric gratings with both bound states in the continuum and guided mode resonances. Opt. Express 29:42978-42988. https://doi.org/10.1364/OE.446937\u003c/li\u003e\n\u003cli\u003eKoshelev K, Kruk S, Melik-Gaykazyan E, et al (2020) Melik-Gaykazyan, et al., Subwavelength dielectric resonators for nonlinear nanophotonics. Science 367:288-292. https://doi.org/ 10.1126/science.aaz3985\u003c/li\u003e\n\u003cli\u003eKoshelev K, Lepeshov S, Liu M, et al (2018) Asymmetric metasurfaces with high-Q resonances governed by bound states in the continuum. Phys. Rev. Lett. 121:193903. https://doi.org/10.1103/PhysRevLett.121.193903\u003c/li\u003e\n\u003cli\u003eKoshelev K, Bogdanov A, Kivshar Y (2019) Meta-optics and bound states in the continuum. Sci. Bull. 64:836-842. https://doi.org/10.1016/j.scib.2018.12.003\u003c/li\u003e\n\u003cli\u003eKoshelev K, Bogdanov A, Kivshar Y (2020) Engineering with bound states in the continuum. Optics \u0026amp; Photonics News 31:38-45. https://doi.org/10.1364/ OPN.31.1.000038\u003c/li\u003e\n\u003cli\u003eSadrieva Z F, Sinev I S, Koshelev K L, et al (2017) Transition from optical bound states in the continuum to leaky resonances: role of substrate and roughness. ACS Photon. 4:723. https://doi.org/10.1021/acsphotonics.6b00860\u003c/li\u003e\n\u003cli\u003eOtsuji T, Popov V, Ryzhii V (2014) Active graphene plasmonics for terahertz device applications. J. Phys. D: Appl. Phys. 47:094006. https://doi.org/ 10.1088/0022- 3727/47/9/094006\u003c/li\u003e\n\u003cli\u003eRoy S, Mondal S, Debnath K (2022) Graphene-based chiral metasurface for generation of tunable circular dichroism-design and sensor applications. IEEE Sens. J. 23:285-292. https://doi.org/10.1109/JSEN.2022.3221437\u003c/li\u003e\n\u003cli\u003eHe X, Huang Y, Yang X, et al (2017) Tunable electromagnetically induced transparency based on terahertz graphene metamaterial. RSC advances 7:40321-40326. https://doi.org/10.1039/C7RA06770D\u003c/li\u003e\n\u003cli\u003eWang W, Yan H, Xiong L, et al (2022) Dual quasi-bound states in the continuum in graphene/dielectric hybrid metamaterial from out-of-plane symmetry protection. Journal of Modern Optics 69:225-231. https://doi.org/10.1080/09500340.2021. 2024612\u003c/li\u003e\n\u003cli\u003eLi J, Li J, Zheng C, et al (2021) Free switch between bound states in the continuum (BIC) and quasi-BIC supported by graphene-metal terahertz metasurfaces. Carbon 182:506-515. https://doi.org/10.1016/j.carbon.2021.06.037\u003c/li\u003e\n\u003cli\u003eWang W, You B (2022) Symmetry-protected bound states in the continuum in graphene nanoribbons. J. Opt. 24:055001. https://doi.org/10.1088/2040-8986/ac6194\u003c/li\u003e\n\u003cli\u003eRoy S, Mondal S, Debnath K (2023) Symmetric Bound States in the Continuum in an All Graphene Metasurface-Design and Sensor Applications. IEEE Sens. J. 23:8352-8359. https://doi.org/10.1109/JSEN.2023.3255223\u003c/li\u003e\n\u003cli\u003eChen J, Zhang M, Gu P, et al (2022) Broadband, wide-incident-angle, and polarization-insensitive high-efficiency absorption of monolayer graphene with nearly 100% modulation depth at communication wavelength. Results Phys. 40 105833. https://doi.org/10.1016/j.rinp.2022.105833\u003c/li\u003e\n\u003cli\u003eYan Z, Kong L, Tang C, et al (2022) Ultra-broadband and completely modulated absorption enhancement of monolayer graphene in near-infrared region. Opt. Express 30:34787-34796. https://doi.org/10.1364/OE.470792\u003c/li\u003e\n\u003cli\u003eKhardikov V V, Iarko E O, Prosvirnin S L (2010) Trapping of light by metal arrays. J. Opt. 12:045102. https://doi.org/10.1088/2040-8978/12/4/045102\u003c/li\u003e\n\u003cli\u003eNiu J, Zhai Y, Han Q, et al (2021) Resonance-trapped bound states in the continuum in metallic THz metasurfaces. Opt. Lett. 46:162-165. https://doi.org/10.1364 /OL.410791\u003c/li\u003e\n\u003cli\u003eCong L, Singh R (2019) Symmetry-protected dual bound states in the continuum in metamaterials. Adv. Opt. Mater. 7:1900383. https://doi.org/10.1002/adom. 201900383\u003c/li\u003e\n\u003cli\u003eGuo B, Fang L, Zhang B, et al (2011) Graphene doping: a review. Insciences J. 1:80\u0026ndash;89. https://doi.org/10.5640/insc.010280\u003c/li\u003e\n\u003cli\u003eLi Z, Xiang Y, Xu S, et al (2022) Ultrasensitive terahertz sensing in all-dielectric asymmetric metasurfaces based on quasi-BIC. J. Opt. Soc. Amer. B 39: 286-291. http://dx.doi.org/10.1364/josab.442660\u003c/li\u003e\n\u003cli\u003eYin W, Shen Z, Li S, et al (2022) THz absorbers with an ultrahigh Q-factor empowered by the quasi-bound states in the continuum for sensing application. Opt. Express 30: 32162-32173. http://dx.doi.org/10.1364/OE.469962\u003c/li\u003e\n\u003cli\u003eXie Q, Dong G X, Wang B X, et al (2018) High-Q Fano resonance in terahertz frequency based on an asymmetric metamaterial resonator. Nanoscale Res. Lett. 13. http://dx.doi.org/10.1186/s11671-018-2677-0\u003c/li\u003e\n\u003cli\u003eHu S, Liu D, Yang H, et al (2019) Staggered H-shaped metamaterial based on electromagnetically induced transparency effect and its refractive index sensing performance. Opt. Commun. 450: 202-207. http://dx.doi.org/10. 1016/j.optcom.2019.05.038\u003c/li\u003e\n\u003cli\u003eCen W, Lang T, Hong Z, et al (2022) Ultrasensitive flexible terahertz plasmonic metasurface sensor based on bound states in the continuum. IEEE Sens. J. 22: 12838-12845. http://dx.doi.org/10.1109/JSEN.2022.3177538\u003c/li\u003e\n\u003cli\u003eSadafi M M, Karami H, Hosseini M (2021) A tunable hybrid graphene-metal metamaterial absorber for sensing in the THz regime. Curr. Appl. Phys. 31: 132-140. http://dx.doi.org/10.1016/j.cap.2021.07.020\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"plasmonics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"plas","sideBox":"Learn more about [Plasmonics](https://www.springer.com/journal/11468)","snPcode":"11468","submissionUrl":"https://submission.nature.com/new-submission/11468/3","title":"Plasmonics","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"Graphene metasurface, Bound state in the continuum (BIC), Plasmonic quasi-BIC, Terahertz sensing, Tunable","lastPublishedDoi":"10.21203/rs.3.rs-3087509/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-3087509/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThis study proposes an all-graphene metasurface supporting dual band symmetric bound state in the continuum (BIC) in the terahertz (THz) range for the first time. The structure consists of a unit cell containing a double-gap split ring resonator periodically on a dielectric substrate. By introducing symmetry breaking, two plasmonic quasi-BICs (Q-BICs) transmission dips are observed with finite Q-factors and high modulation depth. Simulation and analysis results simultaneously exhibit that the Q-factors of Q-BICs follow an inverse square dependency with the asymmetry degree. Via changing the graphene\u0026rsquo;s chemical potential, the Q-BICs\u0026rsquo; operating frequency range can be actively expanded. The dual Q-BICs are immune to the variation of incident angle and has a significant slow light effect with time delay up to 23.8 ps. In addition, the sensing performance in THz region is investigated. A maximum sensitivity of 267.5 GHz/RIU is obtained with a FOM of 24.08 RIU\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e. Our work shows an alternative way to design the class of tunable Q-BIC metasurface, which will provide a valuable reference for future dynamic sensor and other fields.\u003c/p\u003e","manuscriptTitle":"Excitation of tunable dual quasi-bound states in the continuum in graphene metasurface and terahertz sensing application","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2023-06-29 20:42:43","doi":"10.21203/rs.3.rs-3087509/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Accepted","date":"2023-07-06T10:21:35+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2023-07-03T14:27:46+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"3b12eb3c-881d-47fa-9b79-1fee0976d085","date":"2023-06-28T07:23:32+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2023-06-27T17:02:51+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2023-06-26T16:12:01+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2023-06-26T16:12:01+00:00","index":"","fulltext":""},{"type":"submitted","content":"Plasmonics","date":"2023-06-20T13:05:10+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"plasmonics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"plas","sideBox":"Learn more about [Plasmonics](https://www.springer.com/journal/11468)","snPcode":"11468","submissionUrl":"https://submission.nature.com/new-submission/11468/3","title":"Plasmonics","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"f59a82d8-a962-404f-8bef-a98a4c2f9b7c","owner":[],"postedDate":"June 29th, 2023","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2023-08-18T04:01:33+00:00","versionOfRecord":{"articleIdentity":"rs-3087509","link":"https://doi.org/10.1007/s11468-023-01942-6","journal":{"identity":"plasmonics","isVorOnly":false,"title":"Plasmonics"},"publishedOn":"2023-07-15 01:08:24","publishedOnDateReadable":"July 15th, 2023"},"versionCreatedAt":"2023-06-29 20:42:43","video":"","vorDoi":"10.1007/s11468-023-01942-6","vorDoiUrl":"https://doi.org/10.1007/s11468-023-01942-6","workflowStages":[]},"version":"v1","identity":"rs-3087509","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-3087509","identity":"rs-3087509","version":["v1"]},"buildId":"_2-kVJe1T_tPrBINL-cwx","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}