Introduction

Chirality refers to the property of an object that cannot be superimposed with its mirror image after rotations or translations1,2,3. It is quite common in nature, e.g., in organic molecules4, quartz crystals5, and many others. The study of chirality is fundamentally important in various areas, including analytical chemistry, pharmaceutics, and even searching for extraterrestrial life. Light interactions with these geometries can induce chiroptical effects, including circular dichroism (CD) and optical activity, manifested by a difference in intensity and phase responses between left/right circularly polarized (LCP/RCP) light illuminations. Optical chirality was previously demonstrated in three-dimensional (3D) photonic structures6, such as the helices6,7, twisted cross structures8,9 and multi-layered structures10,11,12, all exhibiting strong CD and optical activity with preserved circular polarization, which however require very demanding 3D nanofabrication techniques. Recently, two-dimensional (2D) or planar structures have been shown to support both intrinsic planar chirality at normal incidence13,14,15,16 and extrinsic chirality at oblique incidences17,18,19,20,21. Being different from 3D chirality with broken mirror symmetry in the propagation direction, although even existing in some stacked structures15, planar chirality shows a circular polarization conversion between the output and incident light, and the cross-polarization behavior can be further exploited to modulate the geometric phase for arbitrary chiral wavefront shaping22,23,24,25,26. However, based on existing approaches it is challenging to achieve maximal planar chirality combined with ultra-high quality (Q-) factor due to absorption and scattering loss, hindering many applications that rely on strong chiral light-matter interactions.

Importantly, bound states in the continuum (BICs) can provide a feasible solution for significant problems of chiral photonics. BICs are identified as localized states coexisting with extended modes within the light cone, and they have attracted tremendous attention due to their unbounded Q-factors, which may boost light-matter interactions for applications such as surface-emitting lasing27,28,29, biomedical sensors30,31,32, and nonlinear frequency conventers33,34,35,36,37,38. BICs originate from destructive interference of several radiative channels, and hence they cannot be accessed externally. Quasi-BICs (q-BICs) with high-Q resonances are tailored perturbing the symmetry to couple out the resonant mode to free-space radiation39,40,41,42,43,44,45,46,47,48,49. Recently, chiroptical nanostructures mediated by BICs have been proposed, demonstrating perfect unitary chirality and ultrahigh Q-factors50,51, and greatly expanding the available platforms to achieve optical chirality51,52,53,54,55,56,57.

Here, we suggest and realize experimentally a design of planar chiral metasurfaces at optical frequencies with simultaneously mediated maximal chirality and ultrahigh Q-factors. A double-sided scythe (DSS-) shaped α-Si inclusion with in-plane inversion C2 symmetry but without in-plane mirror symmetry is employed to construct the BIC state. Such BIC state with strong birefringence supports a vortex polarization singularity (V point) surrounded by elliptical eigenstate polarizations with non-vanishing helicity, providing a convenient way to achieve maximum and high-Q chirality by slightly perturbing the inversion symmetry. With either breaking in-plane inversion symmetry or varying illumination symmetry, we realize planar chiral q-BIC states with strong intrinsic or extrinsic chirality, exceeding CD = 0.99 (in simulation) and CD = 0.93 (in experiment). For the intrinsic chirality, the Q-factor of the CD spectra has an inverse quadratic relation with respect to the geometrical asymmetry, with a nearly unchanged CD peak. For the extrinsic chirality, both value and sign of CD can be flexibly tuned by changing the incident angles encircling the Г point, which is connected with the intrinsic non-vanishing helicities of the eigen-polarizations near the BIC singularity58. It reveals that the maximal linear CD does not necessarily occur at the q-BIC with circular eigen-polarization due to the interference with background scattering channels; while the maximal nonlinear CD can be achieved exactly at the circularly polarized q-BIC with strong near-field enhancement CD based on the same planar chiral metasurface design. We believe that the demonstrated planar chiral metasurfaces governed by the BIC physics may find many applications in chiral lasers, nonlinear filters, and other active chiroptical devices.

Results

Most previous q-BIC responses in metasurfaces usually did not support polarization effects28,43,44 (Fig. 1a) or only supported linear polarization selectivity30,39,41. Very recently, 3D optical chiral q-BICs with versatile polarization responses (chirality selectivity) were theoretically proposed by breaking the out-of-plane symmetry with stereo-structures50,51,57 (Fig. 1b). Our present work shows that planar metasurfaces with only reduced in-plane symmetry can support 2D extrinsic and intrinsic chirality with q-BICs (Fig. 1c). Figure 1d shows the geometry of interest, a double-sided scythe (DSS) α-Si structure (refractive indexes refer to Supplementary Fig. S2) employed as the unit cell of our proposed planar metasurface. Such DSS structure (W1 = W2 and L1 = L2, but L1 ≠ W1) has in-plane inversion (C2) symmetry, but lacks any in-plane mirror symmetries. From the simulated band structure and quality factor as shown in Fig. 1e, f, we see that the DSS metasurface hosts a BIC at the Г point, characterized by a vertical magnetic dipole (MD) mode as shown in Fig. 1e, f. Different from the Г point BIC state supported by highly symmetric inclusions58, which is companied by a polarization singularity surrounded by vortex linear polarizations, the BIC state of the DSS structure supports a polarization singularity enclosed by a vortex elliptical eigen-polarization with non-vanishing helicity (see the inset of Fig. 1f and Supplementary Fig. S3). The planar chirality behavior can be predicted by temporal coupled mode theory (TCMT) (Supplementary Notes I & II), which suggests that a perfect planar chiral response can be obtained by tuning the local response to be birefringent using a geometry that has an inversion center, and then breaking inversion symmetry with a suitable perturbation to yield an intrinsic planar chiral q-BIC. Our proposed DSS inclusion with reduced symmetry provides a smart perturbation degree of freedom along the 45° direction, capable of achieving both maximal extrinsic and intrinsic chirality by perturbing a single parameter. As a general remark, in our systems the sign of the CD is determined by the sign of the helicity, while the maximum of CD does not necessarily coincide with circular eigen-polarization with the highest helicity, because there is non-zero background scattering accumulated by Jones matrix elements (see Supplementary Notes I and Supplementary Fig. S4).

Fig. 1: Schematics of different BIC states and planar chiral q-BIC metasurface with a polarization singularity enclosed by vortex elliptical eigen-polarizations with non-vanishing helicities.
figure 1

a–c Comparison of q-BICs (a) in planar metasurfaces without polarization effect, (b) in stereo-nanostructures with 3D chiral effect in reflection, and (c) in planar metasurfaces with 2D chiral effect in transmission. d Schematic of the unit cell of the planar chiral q-BIC metasurface formed by the amorphous silicon (α-Si) DSS structure on a SiO2 substrate with parameters P = 850 nm, R = 280 nm, L1 = L2 = 220 nm, W1 = W2 = 191 nm, and H = 350 nm. e Simulated band structure with real parts of the eigen-wavelength of the DSS metasurface near the Г point. Insets show the field patterns (color: magnetic distributions, arrow: electric vectors) of the three typical modes. f Simulated Q-factors of the eigenmode in the k-space in the vicinity of the Г point. Insets show the first Brillouin zone and eigen-polarization profiles. On the eigen-polarization profile, the polarization states are represented by ellipses of which the blue and red represent the eigen left-handed and right-handed states, respectively, and the black dot represents the V point (BIC).

Based on such a BIC state supported by DSS inclusion with reduced symmetry, one can expect both extrinsic and intrinsic planar chirality by slightly breaking the BIC symmetry. Figure 2a–d show the q-BIC with near-unity extrinsic chirality achieved by illumination symmetry breaking (varying incident angle θ and conical angle φ). Here, the CD is defined as the transmittance difference under right-handed polarization (RCP) and left-handed polarization (LCP) incidence [Eq. 1]:

$${{{{{\rm{CD}}}}}}=\frac{(T_{{rr}}+{T}_{{lr}})-(T_{{rl}}+{T}_{{ll}})}{(T_{{rr}}+{T}_{{lr}})+(T_{{rl}}+{T}_{{ll}})},$$
(1)

where Tij = |tij |2(i = r, l; j = r, l; r represents RCP, l represents LCP) is the transmittance of output polarization i from the input polarization j, and all element of tij construct the Jones matrix under circular polarization basis [Eq. 2]11,59

$${J}_{{circ}}=\,\left(\begin{array}{cc}{t}_{{rr}} & {t}_{{rl}}\\ {t}_{{lr}} & {t}_{{ll}}\end{array}\right).$$
(2)
Fig. 2: Proposed implementation of the DSS structure with both the extrinsic and intrinsic planar chirality by transferring the BIC to q-BICs.
figure 2

a, e Schematic of the symmetry-breaking processes that transfer the BIC to planar chiral q-BICs: (a) illumination symmetry breaking by varying incident angle θ, (e) in-plane geometry symmetry breaking by δ = W2-W1. b, f Simulated transmission Jones matrix spectra of Tll, Trr, Trl and Tlr as well as the CD spectrum of the metasurface with (b) the same parameters as Fig. 1d (δ = 0 nm) under oblique incidence (θ = 8°, φ = 90°) and with (f) an asymmetric structure parameter (δ = 40 nm) for normal incidences. c, g The evolution of CD spectra by continuous varying (c) incident angle θ along the φ = 90° direction and by (g) the geometry asymmetry parameter δ. d, h Dependence of the Q-factors of the planar chiral q-BIC mode on the relative asymmetry parameter (d) α1= sinθ, (h) α2= δ/W1 around the chiral q-BIC state. The solid line shows an inverse quadratic fitting.

Figure 2b shows the calculated transmission spectra of all Jones matrix elements and the CD spectrum of the DSS metasurface with C2 symmetry (W1= W2 and L1 = L2, but L1 ≠ W1) at oblique incidence (θ = 8o, φ = 90o). Tlr has a sharp peak >0.9, while all the other three Jones matrix elements (Trl, Tll and Trr) exhibit dips close to 0 at the resonant wavelength (1456 nm) of the q-BIC, resulting in an ultra-sharp CD spectrum with a maximum >0.95 (green curves in Fig. 2b). Figure 2c shows the evolution of CD spectra by continuously varying the incident angle θ along the y-direction (φ = 90°) (evolutions of transmission spectra of all Jones matrix elements are shown in Supplementary Fig. S5). The linewidth of the CD spectrum increases from 0 to finite values as θ increases, and the peak CD sustains near-unitary value with finite linewidths, manifesting the transfer process from a non-radiative BIC to a high-Q radiative q-BIC with strong extrinsic chirality. The Q-factor of the extrinsic chiral q-BIC shows the expected inverse quadratic law with the illumination asymmetry parameter α1= sinθ (Fig. 2d), which provides a convenient way to tailor the Q-factor of the chiral response. The evolutions of the extrinsic chiral q-BIC along the x-direction (φ = 0°) are shown in Supplementary Fig. S6, the variation trend is the same as the φ = 90° case, except that the sign of CD is flipped due to the opposite helicity of eigen-polarization supported by the DSS structure along x- and y-directions (the inset of Fig. 1f). To investigate the full CD flip picture, the evolution of CD spectra with varying conical angle φ at a fixed incident angle θ = 12o is shown in Supplementary Fig. S7. The opposite varying trend of Tlr and Trl can be observed as φ is varying, leading to the sign flip of CD at positions near (2n−1) × 45°, (n = 1,2,3,4), which is consistent with helicity sign flip shown in the inset of Fig. 1f.

Beyond the extrinsic chiral behavior, we can also access intrinsic high-Q chiral behavior under normal incidence. This can be achieved by introducing an in-plane geometrical asymmetry parameter (δ = W2−W1) as shown in Fig. 2e. As the DSS structure supports strong birefringence along 45°, an in-plane perturbation that breaks inversion symmetry will provide a degree of freedom to obtain planar chirality, as TCMT suggests in Supplementary Note II. Figure 2f depicts the transmission spectra of all Jones Matrix elements and the CD spectra of the metasurface with δ = 40 nm (W1 = 191 nm, W2 = 231 nm). The ultra-sharp resonant peak of Tlr and resonant dips of Trl, Tll and Trr are obtained at the wavelength of 1392 nm, resulting in a CD peak >0.99 and an ultra-narrow linewidth of 1.45 nm. We also confirm the planar chiral q-BIC modes by utilizing multipole expansions and studying the nearfield electromagnetic patterns of the modes (Supplementary Fig. S8), which further confirm that the high-Q planar chiral q-BIC states are dominated by the MD mode. Figure 2g, h show the CD evolution spectra as a function of the asymmetry parameter δ, the Q-factor of the CD spectrum also manifests an inverse quadratic relationship against the relative asymmetry parameter defined by α2 = δ/W1, which is the signature of symmetry-protected BICs (corresponding spectra for all components of the Jones matrix spectra refer to Supplementary Fig. S9). In addition, the circular eigen-polarization state gradually moves away from the Г point when the asymmetry parameter δ is increasing, as shown in Supplementary Fig. S10. Different from the intrinsic 3D chirality with preserved helicity between input and output polarizations, the intrinsic planar chirality is always accompanied by circular conversion dichroism (CCD) and asymmetry transmission (AT)2,13,60, as shown in Supplementary Fig. S11, due to the preserved symmetry in the propagation direction. The CCD is produced by the mutual orientation of chiral elements and light transmission direction, and the peak value of CCD is close to 1, which is slightly larger than CD (Supplementary Fig. S12). Perfect unitary circular conversion transmittance Tlr and CD can be achieved with a strictly symmetric background as shown in Supplementary Fig. S13. The existence of a practical SiO2 substrate introduces a perturbation of the background scattering, resulting in reduced transmittance and CD, as presented in Fig. 2b, f.

Due to the symmetry-protected origin of this phenomenon, our designed DSS metasurface has good tolerance against fabrication imperfections, as tested in Supplementary Fig. S14, which provides a convenient way for experimental realization. We fabricated the DSS metasurfaces with different in-plane asymmetry parameters (δ = 0, 20, 40, 60, and 80 nm), and employed a homemade optical setup (Fig. 3a) to measure the Jones matrix spectra of the metasurface. We first examine the extrinsic chiral q-BIC by measuring the DSS metasurface without in-plane geometry asymmetry (δ = 0) for varying incident angles (Fig. 3b). Figure 3c shows the comparison between simulated and measured transmission Jones matrix spectra at different incident angles (θ = 2°, 4°, 6°, 8°, 10°, 12°) along the φ = 90o direction. At normal incidence (middle panel of Fig. 3c), all transmission spectra are smooth without resonant features, because the BIC state at the Г point is decoupled from any scattering channels of the external environment. At other incident angles, sharp circular conversion peaks of Tlr and resonant dips of other three Jones matrix elements (Tll, Trl, Trr) emerge in both simulations and experiments. The measured absolute peak transmittances are >0.78 for all incident angles, which is a little below the simulated 0.9, mainly due to the practical scattering loss of the fabricated sample as well as a limited collection efficiency of the measurement setup. As expected, we can clearly observe the linewidth increment of the q-BIC peak with increasing incident angle in the experiment, consistent with our theory and simulations. From the measured Jones matrix spectra, we extracted the CD spectra for different incident angles corresponding to Eq. (1) in Fig. 3d. As we can see, sharp CD peaks appear at the q-BIC resonant wavelengths, the simulated/measured CD maxima reach 0.67/0.61, 0.80/0.72, 0.95/0.83, 0.95/0.82, 0.96/0.86 and 0.96/0.88 with the simulated/measured Q-factors of 1835/602, 616/258, 261/142, 143/96, 95/56 and 58/39 for incident angles θ = 2°, 4°, 6°, 6°, 10°, 12°, respectively. For normal incidence, the CD is always 0 as the BIC state is not accessed by external excitation. We also experimentally characterized the extrinsic chiral q-BIC along the φ = 0o direction in Fig. S15, where sharp transmission peaks of Trl are observed, achieving opposite signs of simulated/measured CD: −0.89/−0.82, −0.95/−0.88, and −0.97/−0.91 with the simulated/measured Q-factors of 798/312, 196/102, 85/55 at incident angles θ = ±4°, ±8°, ±12°, respectively, which further confirms the CD flip behavior.

Fig. 3: Experimental verification of the extrinsic planar chiral q-BIC with illumination symmetry breaking.
figure 3

a Experimental setup for the Jones matrix spectra measurement under the circular polarization basis. P1 and P2 represent the polarizers, QWP is the quarter-wave plate, L1, L2, and L3 are lenses. b Schematic of the metasurface with δ = 0 supporting chiral q-BICs through different incident angles along the φ = 90°direction. c Simulated and measured transmission Jones matrix spectra (Tlr, Trl, Trr and Tll) for different incident angles (θ = 0°, 2°, 4°, 6°, 8°, 10°, 12°) along the φ = 90° direction. d Simulated and measured CD spectra for different incident angles extracted from the Jones matrix spectra in (c).

To verify the intrinsic chiral q-BIC further experimentally, we characterized the DSS metasurfaces with different in-plane geometrical asymmetries, as shown in Fig. 4. Figure 4a shows the scanning electron microscope (SEM) images of fabricated samples with different asymmetry parameters δ. From both the top view and slant view, we see that the fabricated DSS silicon pillars are smooth and uniform, in agreement with the theoretical design. Figure 4b, c show the simulated and measured transmission Jones matrix and CD spectra of the five fabricated metasurfaces. When δ = 0 (lowest panel in Fig. 4b, c), there is no resonance with Trl = Tlr and Tll = Trr due to reciprocity, and CD = 0 crossing the whole frequency band. As the asymmetry parameter δ becomes non-zero, sharp circular conversion transmission peaks of Tlr and CD appear, and the bandwidth increases with the increase of the asymmetry parameter (upper four panels in Fig. 4b, c). The simulated/measured CD maxima for the intrinsic planar chirality are 0.98/0.88, 0.99/0.93, 0.99/0.92, and 0.99/0.91 with the simulated/measured Q-factors of 2457/390, 480/121, 198/62, and 128/38 for δ = 20, 40, 60, and 80 nm, respectively. We note that, the Q-factor can reach the order of 106 with still relatively high CD of 0.88 by further decreasing the asymmetry parameter to δ = 2 nm, as the simulation results in Supplementary Fig. S16 indicate. Such simultaneously achieved strong chirality and high Q-factors promise chiroptical functionalities with strong light-matter interactions.

Fig. 4: Experimental verification of the intrinsic planar chiral q-BIC with various in-plane geometry asymmetry parameters.
figure 4

a The top view of the SEM image of five fabricated DSS metasurfaces with asymmetry parameters δ = 0 nm, 20 nm, 40 nm, 60 nm, and 80 nm, respectively. The inset shows the slant view of the metasurface. Scale bar: 1 µm. b Simulated and measured transmission Jones matrix spectra (Tlr, Trl, Trr and Tll) of the five DSS metasurface samples with different δ. c Simulated and measured CD spectra for different δ extracted from the Jones matrix spectra in (b).

The planar chiral q-BIC metasurface can realize not only linear CD maximum, but also maximal nonlinear CD with promising active chiral applications. To pursue the active chiral behavior, a circularly polarized eigenstate whose excitation is independent of incident polarization is highly required. Based on the DSS-shaped design, we have found another set of parameters (P = 780 nm, R = 268 nm, L1 = 197 nm, L2 = 200 nm, W1 = 217 nm, W2 = 257 nm, and H = 654 nm) for generating the circular eigen-polarization state (Fig. 5). Figure 5a shows its transmission spectra under LCP/RCP normal incidences. The sharp resonance peak only exists for RCP incidence, while the overall spectrum is smooth for LCP incidence. It means that the q-BIC resonances could only be selectively excited by the RCP light, indicating a RCP eigenstate. Figure 5b showcases that the eigen-polarization at the Г point is indeed a circular eigen-polarization (C point). Although the linear CD is not high (0.6) at the C point, the near-field enhancement contrast between RCP/RCP incidences is very large (400:1) as shown in Fig. 5c. Such huge near-field enhancement CD is the key to achieving active chiral applications such as chiral lasing, nonlinear chiral emissions based on strong chiral light-matter interaction. We also confirm the planar chiral q-BIC modes by utilizing multipole expansions and study the near-field electromagnetic patterns of the modes for LCP and RCP incidences in Supplementary Fig. S17. Different from the largest linear CD caused by MD mode, the largest nonlinear CD is dominated by ED mode, which is selectively excited by the RCP incidence. To confirm the active chiral emission behavior, we examined the nonlinear CD in both simulation and experiment. Figure 5d shows the simulated results of third harmonic generation (THG). The THG intensity under RCP pumping is significantly enhanced at the q-BIC resonance frequency. While the THG signal under LCP pumping is negligible. Therefore, the extracted nonlinear CD reaches unity at the resonance frequency. Figure 5e, f show the experimental results for both linear CD and nonlinear CD. Figure 5e shows the measured transmission spectra and linear CD. The transmission spectra for LCP incidence are not ideally smooth, which can be attributed to the imperfect circular eigen-polarization caused by fabrication deviations. Figure 5f shows the experimentally measured THG intensity under RCP/LCP pumping. Clear THG emission contrast yields nonlinear CD as high as 0.81. The THG efficiency is significantly enhanced by the q-BIC, as compared to the THG emission intensity by a reference silicon thin film. Limited by the operating band and loss of the objective lens, the effect of THG visible to the naked eye can be realized when the incident power is large enough, as shown in the inset of Fig. 5f. Figure 5g illustrates the log-log plot of emission power as a function of pumping power showing a cubic power dependence, which is the signature of THG. TH conversion efficiency can be significantly enhanced by decreasing the asymmetry parameter of the chiral q-BIC as shown in Supplementary Fig. S18, which is attributed to the tunable and ultra-high Q-factor provided by the q-BIC. Therefore, our DSS structure design can satisfy not only the perfect linear chirality but also the nonlinear chirality. The maximal linear CD is not necessarily related to the circular eigen-polarization state due to the existence of background scattering. However, the nonlinear CD and active chiral behaviors heavily rely on the circular eigen-polarization, which could indeed be achieved by the planar chiral q-BIC metasurface, not necessarily relying on 3D chiral structures.

Fig. 5: Near-unity nonlinear CD with planar chiral q-BIC metasurface.
figure 5

a Simulated transmission spectra of Tlcp and Trcp as well as the linear CD spectra of the metasurface with parameters P = 780 nm, R = 268 nm, L1 = 197 nm, L2 = 200 nm, W1 = 217 nm, W2 = 257 nm, and H = 654 nm with an asymmetric structure parameter (δ = 40 nm) for normal incidences. Insets show the Structure diagram. b The eigen-polarization profile in the k-space. The red dot represents the circular eigen-polarization (C point). c Contrast of electric field enhancement |E | 2 of different circularly polarized light. d Simulated THG intensity under different circularly polarized light incidence and the corresponding nonlinear CD. The nonlinear CD is defined as the normalized THG difference between RCP and LCP lights. e Measured transmission spectra of Tlcp and Trcp as well as the linear CD spectrum. Insets show the SEM image of the metasurface sample. f Experimentally measured THG intensities under RCP, LCP incidences in planar chiral q-BIC metasurface, THG intensity in a reference Silicon thin film, as well as the measured nonlinear CD spectra. The inset shows a photographic image of the light spot of THG from the sample. g Power dependence of THG in logarithmic scale showing cubic power scaling law. The red dots show the measured data, the black dashed line is a fit to the data with a third-order power function.

Discussion

We have suggested and realized experimentally planar chiral metasurfaces with in-plane arrays of symmetry-reduced meta-atoms. Such metasurfaces are governed by the physics of bound states in the continuum with strong birefringence, identified as a vortex polarization singularity surrounded by ellipse eigenstate polarizations with non-vanishing helicity. Both extrinsic and intrinsic planar chirality can be realized by symmetry breaking. By illumination symmetry breaking, the BIC state can be transferred to q-BIC with strong extrinsic chirality accompanied by tunable linewidth and signs. By introducing in-plane structure asymmetry, the intrinsic planar chirality can be achieved under normal incidence with maximum CD of 0.99 (in theory) and 0.93 (in experiment) or ultra-high nonlinear CD of 0.99 (in theory) and 0.81 (in experiment) at optical frequencies. Beyond the chirality based on circular polarization, we can envision any elliptical chirality by varying simultaneously birefringence and axis angle of the meta-atoms. Due to the high Q-factor of the planar chiral q-BIC in metasurfaces with its accessibility and controllability, our results provide unique opportunities for many applications requiring ultrahigh Q-factors and chirality control, including chiral lasing and chiral nonlinear optics.

Methods

Simulation of reflection and transmission efficiency

Rigorous coupled wave analysis (S4)61 simulations were carried out to simulate the planar chiral q-BIC metasurface and to calculate efficiencies for the co-polarization handedness-preserved coefficients, and the circular cross-polarization conversion coefficients transmission spectra. The metasurfaces are designed by α-Si columns based on a SiO2 substrate.

Fabrication of samples

The metasurfaces were fabricated on a fused quartz substrate by utilizing electron beam lithography and including the processes such as deposition, patterning, lift-off, and etching. First, a 350 nm-thick amorphous silicon (α-Si) film was deposited by plasma-enhanced chemical vapor deposition method at 120 °C, during which the flow rates of SiH4 and Ar are 10 and 475 sccm. And the deposition pressure and RF power are 650 mTorr and 20 W, respectively. Then a PMMA film of 300 nm was spinning coated and covered by PEDOT: PSS film as a conducting layer. The desired structure was patterned by utilizing JEOL 6300FS EBL at a base dose of 1000 μC/cm2 with an accelerating voltage of 100 kV for 1 h. After the exposure process, the conducting layer was washed away and the resist was developed in 1:3 MIBK: IPA solution for 40 s and rinsed in IPA for 30 s successively, followed by a deposition of 80 nm Cr using electron beam evaporation deposition method. To realize the lift-off process, the sample was immersed in hot acetone of 75 °C and cleaned by ultrasonic. Finally, by using inductively coupled plasma (ICP) reactive ion etching (RIE) method with HBr at room temperature (RT) for 260 s (flow rate of 50 sccm, the pressure of 10 mTorr, RF and ICP power of 50 and 750 W, respectively), the desired structure was transferred from Cr to silicon and the residual Cr was removed by cerium (IV) ammonium nitrate. The five planar chiral q-BIC metasurfaces are composed of 500 × 500 periods (425 µm × 425 µm) with different cross-sections (δ = 0 nm, 20 nm, 40 nm, 60 nm, and 80 nm), and the final structures are shown in Fig. 4a.

Characterization of samples

A supercontinuum laser (Fianium-WL-SC480) was employed as the broadband light source for the measurement of the transmission spectra of all polarization components under circular polarization incidences. Incident light with a circular polarization was generated by cascading a broadband polarizer and a quarter waveplate from the supercontinuum laser. Then, the incident light was focused on the sample by a lens with a focal length of 5 cm. Subsequently, the polarization components are detected by a quarter waveplate and a polarizer. The transmittances were measured by means of an Ocean spectrometer (flame-NIR). Moreover, the sample was mounted on a rotation stage for varying incident angles in the measurement. In the nonlinear optical measurements, we use a spectrally tunable fs laser source (repetition frequency: 80 MHz, pulse duration: ∼146 fs). Incident light with the circular polarization was generated by cascading a broadband polarizer and a quarter waveplate. The laser beam was focused on the sample surface by an objective with effective NA of about 0.42. The transmitted light was collected on the other side by using another objective with NA = 0.95 and then was sent to an Ocean spectrometer (USB4000) for spectrum analysis.