Polarization-selective modulation of supercavity resonances originating from bound states in the continuum

Bound states in the continuum (BICs) are widely studied for their ability to confine light, produce sharp resonances for sensing applications and serve as avenues for lasing action with topological characteristics. Primarily, the formation of BICs in periodic photonic band gap structures are driven by symmetry incompatibility; structural manipulation or variation of incidence angle from incoming light. In this work, we report two modalities for driving the formation of BICs in terahertz metasurfaces. At normal incidence, we experimentally confirm polarization driven symmetry-protected BICs by the variation of the linear polarization state of light. In addition, we demonstrate through strong coupling of two radiative modes the formation of capacitively-driven Freidrich-Wintgen BICs, exotic modes which occur in off-Γ points not accessible by symmetry-protected BICs. The capacitance-mediated strong coupling at 0° polarization is verified to have a normalized coupling strength ratio of 4.17% obtained by the Jaynes-Cummings model. Furthermore, when the polarization angle is varied from 0° to 90° (0° ≤ ϕ < 90°), the Freidrich-Wintgen BIC is modulated until it is completely switched off at 90°. Bound states in the continuum, offering long-lived, spatially-confined photon states, are typically realized via geometric manipulation of symmetry. Here, the formation of bound states in the continuum is instead driven in a THz metasurface by varying the polarization.

T he concept of bound states in the continuum, or BIC, was first proposed by von Neumann and Wigner for an electron in an artificial complex potential 1 . Thereafter, different types of BICs have been reported in quantum systems 2,3 , acoustic and water waves 4,5 and photonic systems [6][7][8][9][10][11][12] , where modes inside the radiation continuum above the light line are perfectly confined instead of radiating away. Due to the infinite quality (Q-) factor and zero linewidth of the BIC modes, only quasi-BICs with partial confinement and finite linewidth are observed experimentally. The partial confinement of waves has applications in lasing [13][14][15][16] , nonlinear phenomena [17][18][19][20] and highperformance sensing devices [21][22][23] .
In general, bound states in the continuum originate from two physical mechanisms. The more common type, symmetryprotected (SP-) BIC results from a symmetry mismatch between radiative modes and the mode profiles inside the Brillouin zone near the Γ point. The second classification, Friedrich-Wintgen (FW-) BICs (or accidental BICs), originates from destructive interference between two radiation modes modulated by a specific parameter 24 . Unlike SP-BICs, FW-BICs are observed at off-Γ points in regions with an accidental symmetry. These regions of the band dispersion result in extremely high Q factors around the FW-BIC, called near-BIC 25 or supercavity resonances 13,26 , and are attributed to the coupling strength of the two radiative modes that interfere destructively. The position and Q factor of the supercavity resonances can be extremely sensitive to structural parameters of the device 27 . Furthermore, in accordance with FW-BIC theory, changes in the coupling strength of the resonances result in a shift of the FW-BIC positions in the band dispersion 24,28 . Efficient control of the position and Q factor of the supercavity resonance is desirable for practical design of high Q-factor devices.
Limited experimental demonstrations of FW-BICs have been shown to be dependent on angles of incident radiation in the infrared regime 10,25,29 . At terahertz (THz) frequencies, SP-BICs [30][31][32] with few realizations of FW-BICs 26,33 have been reported. Recent investigations by Han et al. 26 reported resonance-trapped BICs with similar mechanism as FW-BICs for THz metasurfaces by tailoring the geometric lengths of the silicon resonators. They also reported frequency and Q-factor modulation of the supercavity resonance through optical pumping of the silicon resonators. Furthermore, Pankin et al. 34 reported FW-BICs at particular linear polarization angles in one-dimensional PhC structure with anisotropic defect layer (ADL). The FW-BIC in the structure is formed by complete destructive interference of the extraordinary and ordinary waves at the output from the ADL.
In this Article, we report the modulation of both types of BIC driven by two unique modalities: linear polarization state of light for SP-BIC and through capacitance-mediated strong coupling for FW-BIC. We then investigate the influence of polarization state of the light on FW-BIC formation. The results of this work provide an active methodology to manipulate the extremely high Q-factor resonances of BICs via manipulation of the incoming radiation or capacitive gaps of split-ring resonator (SRR) structure which could be advantageous for photonic, lasing, and sensing applications. The paper is organized as follows. First, we confirm that a quasi-BIC is formed by breaking the symmetry of the SRR unit cell. Next, we investigate the influence of polarization on the transmission spectra of the terahertz metasurface at highest asymmetry, where a SP-BIC is formed at 30°. Then, we demonstrate through modulation of the capactive gap, the phenomenon of FW-BIC in an asymmetric SRR from strong coupling between the radiative plasmonic dipole mode and higher order quadrupole mode. Finally, to investigate the polarization dependence of FW-BIC formation, we vary the polarization angle of incident light from 0°to 90°(0°≤ ϕ < 90°). Notably, the FW-BIC decreases gradually with increasing polarization angle until FW-BIC is completely switched off at 90°.

Results
SP-BICs can be verified by breaking the structural symmetry of the device. Here, we experimentally show the quasi-BICs originating from the translation of the top gap of the unit cell, confirming the results of Cong and Singh 30 . In general, any small perturbation of the BIC including off-normal incidence (Supplementary Note 1) gives rise to quasi-BICs that can be experimentally verified such that SP bound states can radiate through coupling with the incoming radiation.
The unit cell geometry of an asymmetric SRR with four capacitive gaps at δ = 72.5 μm is shown in Fig. 1a. The dimensions are length l of 250 μm, bracket width w of 35 μm and capacitive gap width c of 35 μm arranged in an array with 300 μm period. An asymmetry parameter, α has been defined as 100% × δ i /δ max where δ i is the distance between the gap center and the symmetry line and δ max is 72.5 μm. Experimentally, the SRR arms consist of 100 nm thick silver layer deposited by radiofrequency sputtering on polyimide substrate of 50.8 μm thickness using conventional lithography. An optical microscope image of the metasurface is also shown in Fig. 1b.
As shown in Fig. 1a, the top gap is translated laterally toward the right to break the C 4 symmetry and produce quasi-BICs in the form of a Fano resonance lineshape [35][36][37][38][39][40][41] . The asymmetry parameter α, defined previously, describes a general parameter common to all SP-BICs. The simulated and experimental transmission spectra are plotted when the incident radiation is horizontally polarized (x-pol) in Fig. 2a versus vertically polarized (y-pol) in Fig. 2b. The transition from a SP-BIC to a quasi-BIC, due to asymmetry, is observed for both experiment and simulation in x-and y-polarizations. The degree of asymmetry before quasi-BICs can be detected is larger for x-polarization where the transmission dip only becomes visible in simulations around 11% as compared to 3% asymmetry for y-polarization.
The extremely sharp quasi-BICs produced by very small perturbations are not resolved experimentally due to scattering losses from slight imperfections and defects in the sample. The quasi-BIC closest to SP-BIC modes are observed experimentally at α = 40% for x-polarization and at α = 20% for y-polarization. The inserts inside transmission graphs plotted in Fig. 2a, b show the quasi-BICs in more detail. The spectral feature of quasi-BIC also broadens with increasing asymmetry showing a strong dependence on asymmetric gap distance.
Further evidence of SP-BIC is represented by the disappearance of a resonance at the symmetry restoring region, where the linewidth becomes zero or quality factor tends to infinity 10,42 . The Q factors of the resonances are calculated as Q = ω 0 /2γ where ω 0 is center resonance frequency and γ is damping rate of the resonance extracted from the equation below 12 where a 1 , a 2, and b are constant real numbers. The calculated Q factors are then fitted by a Q ∝ α −2 curve 30 . In Fig. 2c, d, the calculated Q factors of quasi-BICs show a clear diverging trend when the structure approaches symmetry. However, the Q-factor diverges closer to the symmetric region for y-polarization. This shows the incident field is more sensitive to the asymmetry in the structure when the electric field is orientated perpendicular to the translated gap. The experimentally measured Q-factor values are resolved for higher degree of asymmetry with 23.7 for x-polarization and 17.2 for y-polarization at α = 40%. The extracted   To understand the dependence of incident polarization, the ultrahigh asymmetry sample (α = 100%) was measured for angle ϕ rotated counter-clockwise from 0°to 90°to mimic different linear polarization states; as illustrated in Fig. 3a. The resulting spectra for different rotation angles are plotted in Fig. 3b with the accompanying simulated spectra, where ϕ in that case corresponds to the polarization angle of the incoming THz wave. The measured Fano resonance disappears when the sample is rotated to ϕ = 30°and reappears at ϕ = 45°. This behavior is attributed to the phenomenon of SP-BIC. In a symmetric unit cell, the x and y components of the electric field cancel out such that the overlap integral is zero with respect to any polarization state of the incident beam and thus forming SP-BIC. Meanwhile for our asymmetric sample, the sample rotation alters the mode profile such that the overlap integral is 0 for a polarization angle of 30°r esulting in the disappearance of the Fano mode to form SP-BIC. A colormap for simulated amplitude as a function of incident polarization angle is presented in Fig. 3c. The Fano, dipole, and quadrupole modes are denoted as F, D, and Q respectively in the figure. Furthermore, the simulated transmission amplitude and Q factor of the Fano dip are plotted as a function of rotation angle in Fig. 3d. The transition from Fano to SP-BIC and back to a Fano-type resonance can be observed in the transmission amplitude with a peak at ϕ = 30°in Fig. 3d. This angle represents the SP-BIC point with vanishing linewidth as described by a cross in pictorial illustration of Fig. 3c. The Q factor also becomes infinite as the rotational angles approach the SP-BIC point near ϕ = 30°as seen in Fig. 3d.
In previous studies, SP-BICs occurred at nonzero incidence angles 30,32 and through geometric manipulation 30,31 . The SP-BIC experimentally observed in our system is driven by tuning the linear polarization states of electromagnetic field at normal incidence on the sample. This allows for a polarization selective method to produce the SP-BIC in an asymmetric sample through polarization engineering. The advantage of EM waves is the scalability of optical responses to different frequency bands. Therefore, we can take advantage of this result to induce the polarization selective SP-BICs in other deep subwavelength regions by tailoring the metasurface geometry.

Discussion
A noteworthy feature in the transmission spectra of Fig. 3b is the splitting of the dipole resonance into two distinct resonance dips. This occurs for a polarization angle θ between 45°and 60°a pproximately. Further numerical calculations and analysis have shown that this mode splitting is accompanied by a strong crosspolarized component t yx (Supplementary Note 2). At 0°incident polarization angle (Fig. 3b, ϕ = 0°) one can observe three main dips in the transmission spectrum around 0.4, 0.58, and 0.69 THz, respectively with a good agreement between simulation and experiment. In order to get an insight into the nature of the resonances, we simulate the surface current distribution at the aforementioned resonances (Fig. 4a-c). When the symmetry of the structure is broken with the increase of the parameter δ, the lengths of the right and left arms of the SRRs are no longer equal and the electric dipole moments of the coupled SRRs arms differ from each other. The coupling of these two different dipole modes leads to the emergence of two asymmetric Fano line shapes around 0.4 and 0.69 THz, respectively for the 100% of asymmetry case (i.e., α = 72.5 μm). For the sharp asymmetric resonance at 0.4 THz, denoted as F in Fig. 3b, we observe anti parallel currents in the right and left arcs (Fig. 4a). This resonance mode is similar in nature to the inductive capacitive (LC) resonance in a single gap SRR as both resonances result in current configurations that give rise to a magnetic dipole moment perpendicular to the metamaterial plane of the array.
The distinct surface current distribution profiles of the remaining resonances that appear around 0.58 and 0.69 THz suggest the excitation of dipole and quadrupole modes, denoted as D and Q, respectively in Fig. 3b. The dipole resonance at D 0.58 THz is attributed to the excitation of dipole-like plasmons in the vertical arms of the SRRs, where the currents are parallel to the polarization of the electric field. The current distribution as seen in Fig. 4b radiates strongly to the free space giving rise to a broader resonance in the transmission spectrum. The quadrupole resonance is observed at 0.69 THz. The four arrows show how the current distribution behaves in general. There are two sets of out of phase current distribution. Introducing the asymmetry leads to this kind of current distribution that scatters the electromagnetic field very weakly and dramatically reduces the coupling to the free space, thus causing a huge reduction in the radiation losses which eventually leads to very sharp resonance. The quadrupole resonance is inaccessible in the symmetrical structure 40 .
In contrast to SP-BIC, FW-BIC appears due to destructive interference of the resonances at the correct phase matching conditions 24 . This type of bound states is also sometimes referred to as "an accidental" BIC 25,46 . When two resonances approach each other as a function of a certain parameter, interferences cause an anticrossing of the two resonances in their frequency/ energy positions and for a certain value of the parameter, the width of one of the resonances may vanish, thus leading to the generation of a BIC. In this context, by simultaneously altering the width of the capacitive gaps of the structure with 100% of asymmetry, it is possible to strongly couple the radiative plasmonic dipole mode and the higher order quadrupole mode to form a BIC, as illustrated in Fig. 5.
The simulated transmission coefficients at different widths of the capacitive gaps are plotted in Fig. 5a illustrating an avoided crossing behavior. In the absence of strong coupling, the dispersion of the involved modes will merely cross each other at a specific width of the capacitive gaps, namely c = 78 μm. However, due to the strong coupling regime, an avoided resonances crossing is observed with a Rabi splitting Ω R and and a normalized coupling strength ratio g = Ω R /2ω d of about 4.17% (Fig. 5b). This hybrid dispersion state is accompanied by a vanishing of the spectral quadrupole line at c = 78 μm which is a distinct feature of the FW-BIC.
The Q factor of the FW-BIC is plotted in Fig. 5c. At the FW-BIC point, the Q-factor diverges to infinity due to absence of losses 25 . This is characteristic of a FW-BIC arising from destructive interference between radiative modes modulated by a specific parameter. The resonances with high Q-factor  around the FW-BIC region are termed as supercavity resonances.
The dipole-quadrupole strong coupling can be described by a full quantum mechanical Hamiltonian: H ¼ ω d ða y aÞ þ ω q ðb y bÞ þ gða y þ aÞðb y þ bÞ where the operators a † (a) and b † (b) represent the dipole and quadrupole creation (annihilation) terms, respectively. The first two terms of the Hamiltonian describe the dipole and quadrupole individual energies while the second two terms represent the interaction between dipole and quadrupole modes. The Jaynes-Cummings (JC) Hamiltonian also includes counterrotating interaction terms (CRTs), which can be neglected as an approximation to Eq. (2) (see Supplementary Note 3). Figure 5b presents the best fits to our data with the JC Hamiltonian. At the optimum fitting with Eq. (2), the normalized coupling strength ratio, g/ω d , was determined to be 0.0417 in agreement with the value for g/ω d obtained in Fig. 5. Rabi splitting from the effect of the hybridization of light-matter states can be demonstrated using either classical or quantum approaches and depending on coupling strength lead to very similar results. We emphasize that although our system is essentially classical, we choose to use the quantum approach because it is more inclusive. For example, it provides a unique and clear property of the ultrastrong coupling (USC) regime, such as the breakdown of the rotating waves approximation, which may lead to vacuum Bloch-Siegert shift. However, the spectral response of the system can be equally obtained from a classical description not involving any operator algebra, which is demonstrated above by the finite-differencetime-domain simulations.
Therefore, in a similar fashion to lattice-mediated strong coupling THz metasurfaces [47][48][49][50] , one can drive the formation of FW-BIC by capacitance-mediated strong coupling. The advantage of our system over conventional strong light-matter coupling systems is that our system achieves strong coupling with the use of a single structure. In commonly used hybrid light-matter coupling systems, the emitter and cavity are distinct and require a high degree of spatial overlap. As a result, the capacitancemediated strong coupling scheme not only drives the formation of FW-BIC but also provides flexibility in the engineering to study more exotic physics such as USC with normalized coupling strengths >10% 51 .
The possibility to tune the rotation angle ϕ in a continuous way allows to follow the evolution of the dispersion and evaluate the influence of the polarization state of the incident light on the FW-BIC formation. As the polarization angle ϕ is swept, we extract the positions of the minima of sample transmission and plot the dispersion curves for the eigenvalues as a function of the capacitive gap c (Fig. 6a-e). When ϕ is increased from 0°to 90°( 0°≤ ϕ < 90°), one can observe that the position of the FW-BIC (denoted by a white dashed line) decreases gradually. Indeed, when ϕ = 0°, 30°, 60°, and 75°, respectively, the corresponding widths of the capacitive gaps are 78, 72, 70, and 68 μm, respectively. Here, we attribute this FW-BIC to the degeneracy of the modal asymmetry during the tuning process, which is required by the classical anticrossing. When ϕ = 90°, there is no singularity in the dipole-quadrupole anticrossing since the structural asymmetry is quasi-preserved. In other words, the FW-BIC is completely switched off, which enables completing an ON/OFF cycle. This noticeable tuning of the position of the FW-BIC is accompanied by a linear enhancement of the normalized coupling strength ratio g/ω 0 , as shown in Fig. 6f. This polarization control of the FW-BIC can be strategically important in developing optically thin metadevices for quantum information processing in the USC regime.
In summary, we demonstrated two distinct methodologies to drive BICs in THz metasurfaces: (1) linear polarization state of light for SP-BIC and (2) capacitance-mediated strong coupling for FW-BIC. In the first method, SP-BIC is driven by linear polarization state in which the Fano due to asymmetry in α = 100% SRR sample is suppressed at polarization angle ϕ = 30°. For the second method, FW-BIC is driven by strong coupling between plasmonic dipole mode and higher order quadrupole mode in an asymmetric SRR (α = 100%) modulated by capacitance. We then demonstrated the polarization dependence of FW-BIC as confirmed by the gradually decreasing capacitive position of FW-BIC with increasing polarization angle ϕ. In addition, FW-BIC is completely switched off at ϕ = 90°which allows for ON/OFF cycling capability. The subwavelength nature of these metasurfaces allows for much smaller mode volume and lifts the dimensional limitations of dielectric photonic crystals 52 . The metadevices studied here have applications in lasing, nonlinear optics, and high-performance sensing devices.

Measurement.
A high-resolution continuous wave THz spectrometer (Teraview CW Spectra 400) that produces linearly polarized collimated terahertz beams was utilized to measure the transmission spectra of the aforementioned fabricated metasurfaces. In this system, the optical beat frequency of two distributed feedback near-IR diode lasers are tuned to produce coherent THz waves between 0.05 and 1.5 THz. The spectral resolution of around 100 MHz is enough to resolve relatively narrow spectral features, in contrast to the limited resolution of time-domain terahertz systems.
Samples were placed equidistant between the emitter and detector in ambient air conditions. After measurement, the transmission spectrum was calculated as T(f) = P M (f)/P sub (f), where P M (f) and P sub (f) are the filtered power spectra of the metasurface and substrate, respectively. For the polarization-dependence investigation, the metamaterial sample is rotated by a value θ to mimic the change of a polarization state of the incident beam 53 .
Simulation. For numerical simulations, we utilized a finite element method with periodic boundary conditions to simulate a 2D infinite array of unit cells. The length scale of the mesh was set to be less than or equal to λ 0 /10 throughout the simulation domain, where λ 0 is the central wavelength of the incident radiation. The input and output ports are located at 3λ 0 from the metasurface with open boundary conditions. The material parameters used for the structures are σ = 4.1 × 107 S/m for the metallic Ag layer and ϵ = 3.8 for the dielectric Kapton layer. For the polarization dependent study, each unit cell was excited with a THz plane wave of the correct linear polarization state as indicated by the polarization angle from 0°t o 90°.

Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.