Observation of spin-polarized directive coupling of light at bound states in the continuum

Spin-polarized directive coupling of light associated with the photonic quantum spin-Hall effect (QSHE) is a nanoscale phenomenon based on strong spin – orbit interaction that has recently attracted significant attention. Herein, we discuss the experimental manifestation of QSHE intrinsic in the Bloch waves associated with a bound state in the continuum (BIC) of a dielectric photonic crystal metasurface (PhCM). We show numerically that BICs in nanoscale PhCMs have photonic spin angular momentum density transverse to the orbital momentum not only at the interfaces but also inside the confining dielectric medium. Then, we experimentally demonstrate that the fundamental Bloch waves of the BIC mode, macroscopically amplified on resonance, propagate along the symmetry axes of the PhCM obeying spin-momentum locking also at normal incidence, i.e., with no symmetry breaking. This BIC-enhanced spin-directive coupling of light may enable versatile implementations of spin-optical structures, paving the way for novel photonic spin multiplatform devices.


INTRODUCTION
The concept of bound state in the continuum (BIC) was introduced by von Neumann and Wigner in quantum mechanics as unusual spatially localized states of electron waves [1]. Later, it was understood that BIC points can occur in the parametric space of all kinds of wave phenomena [2]. In the last years, optical BICs have been demonstrated in engineered photonic crystal metasurfaces (PhCMs) [3][4][5][6][7][8][9], either at specific conditions of field symmetry mismatch or topological constraints with respect to out-of-plane propagating waves in free space [10][11][12]. In particular, the BIC point is associated with a topological singularity in the momentum space [11][12][13][14][15]. Indeed, radiated fields with momentum close to the BIC point have a polarization winding in a vortex structure around the BIC core, as also recently experimentally verified [13,14]. A striking feature is that the electromagnetic field of a BIC is characterized by a diverging lifetime, or an infinite radiative Q-factor [3,12], from which light-matter interaction could be highly enhanced [16][17][18][19]. Reciprocity implies that far-field excitation of an ideal (infinite Q-factor) BIC mode is not possible. However, real structures can have finite, arbitrarily large Q-factors close enough to the BIC point, and partial coupling to the input light makes its spectral signature to appear as a special Fano-shaped resonance in quasi-BIC regime, which can be exploited for several applications [17][18][19][20]. In addition, of a certain interest is also the polarization structure of the quasi-BIC radiation [21,22]. As for instance, BIC tuning and structured polarization patterns of radiated light can be produced by suitable in-plane C 2 symmetry breaking [23,24]. Furthermore, unidirectional BICs produced by topological charge bouncing and tuned by broken C 2 z symmetry have also been recently demonstrated [25].
Herein, we experimentally find and discuss another phenomenon associated with the intrinsic surface-wave nature of BIC modes in subwavelength PhC thin slabs, the so-called quantum spin-Hall effect (QSHE) of light [26,27]. In solid-state physics, the QSHE is characterized by unique edge states of electrons [28] peculiar of a class of materials termed topological insulators [29]. For time-reversal invariant phases, so-called helical spin-polarized edge states are characterized by spin-momentum locking, i.e., matter waves with opposite spin propagate in opposite directions. Also photonic topological phases that are time-reversal invariant are characterized by two-way spin-polarized directive coupling of light [26,[30][31][32]. In particular, the hallmark of an intrinsic QSHE of light is the transverse spin angular momentum (SAM) density arising in inhomogeneous optical fields [27]. The simplest situation is that of an evanescent wave [26,33,34]. Several experiments have demonstrated the effect of spin-directive coupling in nanoscale dielectric [35][36][37] and plasmonic [38,39] metasurfaces. opposite spin by virtue of near-field coupling, through nanoparticles [35,39] or quantum emitters [36,37], travel in opposite directions at the interface between two media. This phenomenon of spin-polarized directive coupling [34] enables new quantum-information protocols [36,37,[40][41][42][43][44][45][46][47]. However, surface plasmon polaritons in metal are quickly extinguished by absorption and scattering losses. In addition, the efficiency of a nanoparticlemediated coupling, even in loss-free dielectrics, is limited by the small scattering cross section.
In this work, we discuss the observation of spin-polarized directive coupling of light, in a transparent subwavelength PhCM, achieved as a consequence of the surface wave nature of the BIC mode. We numerically show that the tight nanoscale confinement at the BIC point allows storing high levels of energy inside a highly inhomogenously structured optical field, thus carrying transverse SAM not only trivially in the evanescent regions but also inside the confining dielectric medium. Experimentally, resonantly enhanced spin-polarized directive coupling of light is demonstrated also in a symmetric excitation condition. The BIC field not only mediates the coupling to the waves characterized by spin-momentum locking but also resonantly amplifies them. This points out exciting opportunities with potential scalability in integrated photonic circuits and chiral quantum technology applications [35,36,[41][42][43][44][45][46][47]. Indeed, the nearly lossfree dielectric metasurface can be effective for a versatile device implementation. Furthermore, enhancing light-matter interaction at the BIC with concurrent polarization coupling may provide new solutions for spin-optical applications.

NUMERICAL MODELING AND RATIONALE OF THE PHENOMENON
At normal incidence, a square-lattice PhC thin slab (Fig. 1) has a C 4ν point group symmetry. For λ close to the first Bragg condition a ∼ λ∕n e , with a lattice constant and n e effective refractive index, four dominant guided modes exist for each polarization in correspondence to the high symmetry points in the first Brilluoin zone [16].
We considered a dielectric geometry consisting of a square lattice of cylindrical air holes etched in a thin film of silicon nitride (Si 3 N 4 ) [Figs. 1(a)-1(b)] (Supplement 1, Section 1). This material is transparent in both visible and infrared ranges of light and has a real-part refractive index of 2.15 at 532 nm [17]. The unit cell design does not break in-plane inversion symmetry (C 2 symmetry), as visible in Fig. 1(b). Imperfections due to border roughness are randomly oriented and not systematic (tolerance below 1%). These imperfections, together with lattice spacing imperfections, material density fluctuations and absorption loss, finite sample size, and also beam finite collimation eventually allow coupling far-field light with the quasi-BIC mode.
A scheme of the resonator is shown in Fig. 1(c), where the PhCM unit cell is shown with the bottom region made of quartz (SiO 2 , region 1), i.e., the coverslip experimentally supporting the silicon nitride layer, and top material of air (region 2).
When the excitation radiation tends to zero-transverse wavevector components, i.e., k i θ i 0 Γ 0, 0, k z , the dispersion curves ω n k i climb over the light line, and the confined, guided modes become leaky modes radiating in the far field [10]. However, symmetry incompatibility with free-space modes restores the perfect confinement in the PhCM close to θ i 0, approaching which the mode Q-factor diverges, as numerically calculated in Figs. 1(d) and 1(e) for our specific structure. Ideally, the leaky mode becomes a symmetry-protected BIC, of course at specific frequencies depending on the PhCM geometry and composition [3,4]. Numerical modeling was carried out using a rigorous coupled-wave approach (RCWA) [48] and finite element method (COMSOL Multiphysics 5.2a) as detailed in Supplement 1, Section 1. The angular frequency of the lower nontrivial dispersion curve at k 0 (Γ 2 point) is ω 2 a∕2πc a∕λ 0.67621225, as can be seen in the inset of Fig. 1(d). This value corresponds to a BIC point expected to occur at a free-space wavelength λ 1 541.4 nm for a 365 nm, in excellent agreement with our experimental findings. The associated mode has a diverging Q-factor [ Fig. 1(e)]. The BIC point occurs at a frequency below the frequencies a∕λ 1∕n 0.685 (n 1.46) and a∕λ 1 of the modes expected, respectively, in the quartz support and air. The Γ 2 wavevectors, parallel to the x, y plane, have all modulus β 0 2π∕a, i.e., larger than the wavevector modulus both in SiO 2 , given that 2π∕a > 2πn∕λ k SiO 2 , and air, given that 2π∕a > 2π∕λ k air . Since the mode is purely evanescent in regions 1 and 2, at normal incidence, the mode wavevector has a purely imaginary z component in both regions 1 and 2, k 1,2 z ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi k 2 air,SiO 2 − β 2 0 q ≐ iκ 1,2 and cannot couple to the far field. The mode is a Bloch wave exponentially decaying in both ẑ, with a larger decay length in the glass (region 1) since In the case of purely evanescent fields in regions 1 and 2, the allowed triads of transverse spin s and complex wavevector's real part k and imaginary part κ according to spin-momentum locking [34] are depicted in Fig. 1(c) [see Eq. S(2) in Supplement 1, Section 2)]. Figure 2 shows the solutions excited in the case of right-circularly polarized (RCP) and left-circularly polarized (LCP) input plane waves at normal incidence, where convenient finite approximation of the excitation frequency allows far-field source coupling [16]. The degenerate mode excited at normal incidence is a symmetry-protected BIC [4] (diverging lifetime, Q > 10 9 ) with hybrid characteristics of transverse magnetic and electric distributions, with main axial component E z . In Figs. 2(a)-2(b), the magnetic field intensity jH j 2 in the PhCM unit cell is superimposed to the arrow maps of the electric field E in the x, z plane of the cell border where the field is maximum. Top and bottom panels refer to RCP and LCP incident radiation, respectively. The gray region reproduces the Si 3 N 4 elementary volume. The jH j 2 profile in the z-cross section is that of evanescent waves in regions 1 and 2, with larger amplitude at the interface with the SiO 2 . The field structure E appears as a surface wave with strong longitudinal components in the meridional planes. It evolves with cycloidal rotations consistently with the existence of transverse SAM [49] not only trivially in the evanescent regions [30] but also inside the PhCM because of the nanoscale confinement. . The vector map of the field E (uniform magnitude for clarity of representation) shows its character of surface wave and reveals the intriguing behavior of the field excited by opposite circularly polarized input plane waves with opposite orientation of the electric field point by point in (a) and (b). The colormap is associated with magnetic field intensity jH j 2 . The tight confinement produces a significant transverse SAM not only at the interfaces but also inside the nanoscale slab, which locks the propagation direction of the wave along the x axis depending on the helicity of the input plane wave as can be seen by comparing top and bottom Poynting vector maps in the zx plane (c), (d) (arrow length is proportional to the vector magnitude). Because of symmetry, an analogous behavior of directional Poynting vector flow occurs also in the z cross section along y axis (zy plane). The colormap is associated with the energy density. In particular, the relative orientation of SAM and Poynting vector can be seen in detail in (e) and (f ) (arrow length imposed uniform for better visual inspection for both vectors). It is worth emphasizing that the energy density is much weaker at the interface of the PhC cell with air (top region) as shown in panels (c) and (d); thus, the SAM average z-axis projection (S z ) is that represented by the arrows in (e) and (f). The input spin (scheme on the left) matches the average S z of the mode. Since this is transverse to the Poynting vector, as visible in the insets, spin-momentum locking occurs.
The 3D time averaged Poynting vector density Π 1∕2ReE × H ⋆ shows reversed vectors orientation from RCP to LCP [Figs. 2(c) and 2(d)], point by point with a vortex structure. Following [33,49] for definition and calculation, we numerically evaluated the SAM density S Imϵ 0 E ⋆ × E μ 0 H ⋆ × H ∕4ω associated with this BIC mode, which resulted in having transverse components with respect to Π. Indeed, Figs. 2(e)-2(f ) show that the SAM always has components orthogonal to the linear momentum density Π point by point in the unit cell with reversed orientation at the two interfaces, respectively, with SiO 2 (region 1) and with air (region 2), consistent with the spin-momentum locking conditions expected for top and bottom evanescent waves having opposite imaginary wavevectors in the two regions. It is worth mentioning that while this locking mechanism occurs at the boundaries of any optically confining system [34], in the particular case of our nanoscale PhCM at the BIC, the boundaries are so close, with a distance inferior to the same optical wavelength, that the transverse SAM is large also inside the medium where it is resonantly enhanced. Back to Fig. 2, given that the x, y components of the field are comparable in amplitude to the z component, we can see that the spin density has strong components along z, which allows spin matching with the z-directed spin of the input excitation at normal incidence [Figs. 2(e)-2(f )]. Since the spin-momentum locking is a boundary condition that selects the waveguiding direction, of the four possible wave solutions at the boundaries, only two can be excited by a given input spin and correspond to counter-propagating surface waves at the two dielectric interfaces [ Fig. 1(c)]. The boundary spin structure is consistent with the spin structure of the mode inside the confining dielectric, which is on the other hand determined by the spin matching with input light. Importantly, we can easily see from the asymmetry of the mode profile along the z axis in Figs. 2(c)-2(d) that these two waves have different amplitudes at the interfaces. Along a given direction, say x axis in Fig. 2, in one way, there will be more light than the way back for a given input circular polarization, with reversed behavior for the opposite input helicity. Therefore, we can say that the input spin matching selects two surface waves that, if out-coupled again in free space, will be characterized by a different field intensity, depending on the input spin causing a breaking of the planar mirror symmetry of the system even at normal incidence, with a macroscopic spin-polarized directive coupling. Figures S2 and S3 in Supplement 1 show that transverse SAM exists also in the case of a quasi-BIC and in case of a nonsymmetry-protected BIC [5], existing for different values of h.

BIC SPIN-ORBIT ASYMMETRY PHENOMENOLOGY AND CHARACTERIZATION
Several 1 − mm 2 PhCMs of Si 3 N 4 were fabricated with various thicknesses (Supplement 1, Section 3), as also shown in Fig. 1(b). Herein, we will focus our attention on the case of h 144 nm, where the expected BIC mode in Fig. 2 emerges at 541.4 nm. A schematic layout of the optical characterization of the PhCM is shown in Fig. 3(a). The PhCM is placed at the center of the supporting glass coverslip and illuminated with a collimated supercontinuum laser source to measure the transmission spectrum and to identify the characteristic resonances of interest (Supplement 1, Section 3). The PhCM transmission spectrum collected at normal incidence θ i 0°is shown in Fig. 3(b). Two spectral dips are visible at λ 1 541.4 nm and λ 2 550.0 nm, corresponding, respectively, to the BIC point and to a lower Q-factor mode (leaky mode). It is worth noticing that the total Q-factor of the BIC is limited by finite sample size (Q H ), effective collimation, and scattering losses due to the sample imperfections (Q A ) (Supplement 1, Section 3) [50]. As such, we measured a total Q ∼ 10 3 corresponding to a linewidth ≃0.5 nm [ Fig. 3(b)], due mainly to finite lateral confinement and surface roughness, in agreement with previous papers [3,50].
We used an acousto-optical tunable filter to tune the incident wavelength to λ 1 [Fig. 3(c)] or λ 2 [ Fig. 3(d)]-with linewidth ∼4 nm. At λ 1 , the peak incident intensity I in of the transmitted beam through the sample decreased by about 80%, showing a clear dip with 0.5-nm width (Fig. S4). When the light source is resonant with this BIC, the light trapped inside the metasurface positioned at the center of the supporting glass [the green laser spot in Fig. 3(c)] generates four coherent collimated beams propagating along the four-fold symmetry directions of the PhCM. The representative image in Fig. 3(c) is obtained with LCP excitation. A significant fraction of the input intensity I in is conveyed to these beams, which we will term side waves, with typically 0.1I in per wave detected in the far field on the sides of the sample. We indicate the observed side waves by their wavevectors, as displayed in the scheme reported in Fig. 3(a).
Conversely, at λ 2 (leaky mode), no light was visible around the PhCM [center spot, Fig. 3(d)]. The intensity of the light possibly leaking on the sides of the sample was on the order of the detector noise level, i.e., at least four orders of magnitude less intense than the side waves observed at the BIC.
We interpreted these side waves generated at the BIC wavelength as the first four fundamental Bloch's waves Γ 2 composing the BIC mode, thus intrinsically amplified by the resonant mechanism (see additional text in Supplement 1, Section 2).
Given the geometry of the lattice and the input wavelength, if we examine the PhCM as a basic two-dimensional diffraction grating, the Bragg law implies that at normal incidence, the first orders of diffraction have complex wavevectors, since these firstorder waves are evanescent and have phase wavevectors parallel to the structure plane. When the input wavelength matches the singular point defined by the BIC, the resonant mechanism feeds such diffracted waves, which actually compose the mode, and makes them easily detectable even in the far field. Trivially, the rotation of the PhCM about the z axis of an angle ϕ produces an equal rotation of the side waves, which are therefore uniquely associated with the geometry of the system. Now, these side waves, which are collected back in free space and characterized [ Fig. 4(a)], have linear momenta in the x, y plane of the PhCM, and the optical field generating them is characterized by a transverse SAM (Fig. 2). As a consequence, the symmetry breaking due to spin-orbit, as discussed in the previous section, is expected to affect the intensity of the radiation (out-coupled along a specific wavevector) as a function of input polarization. Indeed, the surface wave at the air∕Si 3 N 4 boundary is more tightly confined and weaker than the wave at the Si 3 N 4 ∕SiO 2 interface, as can be seen in Figs. 2(c)-2(d). This last can propagate at the interface Si 3 N 4 ∕SiO 2 , since the Si 3 N 4 covers all quartz areas, and partly leaks to the quartz substrate where it is guided and then measured in the far field. A scheme is depicted in Fig. S5 (Supplement 1).  Table 1: S/P indicates that the values are, respectively, refereed to s-pol or p-pol of (b).
A similar spin-polarized directive coupling in a rotationally symmetric nanoscale structure has been experimentally observed by Petersen et al. [35] in the case of a two-fold symmetric optical nanofiber illuminated by circularly polarized light at normal incidence, indeed despite the rotational invariance of the configuration. In that case, there was an evanescent coupling provided by the light scattering of a gold nanoparticle. In our case, the scattering element is the cylindrical hole (coupled with the other ones in a lattice resonance). The resonant field has transverse SAM inside the confining region because of nanoscale localization. Thus, it is the very intrinsic nature of the eigenmode to determine the Poynting vector flow that obeys spin-momentum locking (see additional text in Supplement 1, Section 2).
Tunable waveplate retarders were used to set the input polarization E i with a polarization-state generator (PSG) consisting of a polarizer (axis //x), a half-wave plate (HWP), and a quarter-wave plate (QWP). The layout of the experiment is shown in Fig. 4(a). The Jones optical transformation of the input beam is detailed in Supplement 1, Section 4. The intensity I j α, β of the side waves (j ∈ fR, L, T , Bg) was measured as a function of the angle β (QWP) for fixed values of α (HWP) [see also Fig. 3(a)]. The measurements shown in Fig. 4(b) were carried out at θ i 0. In particular, we fixed the HWP angle α to two values, α S 0°( S-polarized input on the QWP) and α P 45°(P-polarized input on the QWP), whereas β was varied in the range (0,360°). This means that the handedness of the circular polarization generated by the QWP is reversed by switching α S to α P , which we used to check any possible spurious dependence due to optical misalignment.
Comparing the curves for the bottom side wave (k B ) in Fig. 4(b) that were obtained for α S (left panel) with that for α P (right panel), we can observe that when the input state is LCP, the intensity measured in the bottom side wave is always maximum. In other words, LCP light is coupled preferentially to the propagation direction k B , in agreement with the simulation for LCP input excitation in Figs. 2(b)-2(d) (k B −β 0x ). On the other hand, input RCP light is instead preferentially coupled to the propagation direction k T along positive x axis. Actually, this behavior can be verified also on the side waves R and L. Indeed, in Fig. 4(b), we can observe a clear dependence in Iβ on the handedness of the input polarization for all the side waves, with anticorrelated behaviors between homologous opposite beams, i.e., between waves R and L, and waves T and B.
When the exciting incident SAM (parallel toẑ) couples to the transverse spin of the BIC mode, it locks to the photon flow mediated by the BIC transverse SAM components (z-component on average), which induces opposite propagation directions in the x, y plane of the fundamental Bloch's waves composing the mode, as shown in the simulations in Fig. 2. The resonant mechanism of the BIC is then responsible for the macroscopic nature of this spin-polarized directive coupling.
It is worth stressing that our PhC unit cell, having a cylindrical hole, does not break in-plane C 2 inversion symmetry. Such a symmetry breaking would be responsible for splitting the BIC topological charge (destroying the BIC singularity) into a couple of circularly polarized states of half-integer charge, which can radiate in the far field along directions at opposite angles very close to the normal direction [24]. In our case, instead, the chiral behavior is observed in the light radiated from the Bloch's wave components of the BIC in the plane of the PhC, i.e., at 90°with respect to the normal incidence, a phenomenon that is related to the intrinsic nature of a surface wave of the eigenmode and its SAM orientation, transverse to the momentum.

PHENOMENOLOGICAL MODEL
Also in the case of linear input polarization (QWP angle β 0°), symmetry breaking is visible in the experimental data in Fig. 4(b) betweenx (T , B) andŷ (L, R) waves. When the input field is not spin polarized, i.e., the input is s or p polarized, the response of the system can be seen in terms of the coupling to a surface wave, which in the first place has an axial field component E z . Supposing an s-polarized input field E x , 0, 0 [ Fig. 3(a)], then the surface wave E x , H y , E z propagating along the x axis is favored, having a longitudinal component E x that is conserved by tangential component continuity. For p-polarized input 0, E y , 0, the waves directed along the y axis are similarly favored, again for matching of the field structure of a surface wave along that direction.
We built a phenomenological model to fit the intensity of the side waves, out-coupled in the far field, for all possible input states of polarization. We modeled the scattering of photons from the incidence directionẑ to the x, y plane, considering that redirection is accompanied by nonparaxial spin-to-angular momentum conversion [27,49,[51][52][53][54], as fully detailed in Supplement 1, Section 4.
The polarization E i incident on the sample depends on the wave plate retarder angles α and β, while the transformation U due to redirection introduces dependence on the orientation angles θ r π∕2 and ϕ in the x, y plane [ Fig. 3(a)]. In particular, we considered ϕ to assume only four discrete values determined by the four-fold symmetry of the PhCM lattice as verified experimentally [ Fig. 3(c)]. Compacting the variables of interest in the vector form v j ≡ α, β, ϕ j , with j ∈ fR, L, T , Bg (Supplement 1, Section 4), in the circular polarization basis with output field components fE w , E − w , E z w g, we can write Iv j c jE w v j j 2 c − jE − w v j j 2 c z jE z w v j j 2 , (1) with j ∈ fR, L, T , Bg, θ r π∕2 and ϕ T 0°, ϕ L 90°, ϕ B 180°, and ϕ R 270°. In the above expression, the output intensity of the field redirected along a PhCM symmetry axis is written in terms of its circular and longitudinal polarization components along the direction of propagation. Since the side waves possess a transverse SAM, which the direction of propagation is locked to, in the above expression, we balance the components with the coefficients c admitting values in the range (0,1) to account for a possible spin-polarized coupling of input light to the mode (Fig. 2). In the case of fc 1, c − 0, c z 0g, for instance, light would be right-handed circularly polarized. Similarly, fc 0, c − 1, c z 0g would correspond to an LCP state. These coefficients were used as parameters when fitting the experimental data in Figs. 4 and 5, quantifying the parameter of chiral directivity defined as η c −c − c c − [35]. In absence of spinpolarized coupling along a given observation direction ϕ j , j ∈ fR, L, T , Bg, since RCP and LCP states would have no preferential coupling, then c c − and η 0; thus, for each ϕ j , the intensity becomes Iβ Iβ mπ∕2, with m integer, i.e., periodic of period π∕2. Instead, experimentally, the curves I β show a periodic modulation of period π, indeed characterized by the presence of secondary maxima or minima [ Fig. 4(b)]. When we include spin-polarized coupling, c ≠ c − , the side wave intensity does depend on the handedness of the polarization,

Research Article
which explains the period π and fits the experimental data as a function of the direction, as shown in Fig. 4(b). In Fig. 4(c) ( Table 1), we report the fit parameters and the estimated chiral directivity η. The fit coefficient c z is below 1%, consistent with a far-field detection scheme for which no longitudinal components may exist along the propagation, and thus is not reported. The value jηj varies from 32% to 68% at θ i 0°.
In addition, in Fig. 4(b), once fixed, the side wave, by comparing left and right panels, we can see that by exchanging α from 0°to 45°, an intensity reversal occurs as a function of β. This is due to the geometric phase shift introduced by the HWP, which is well predicted by the model and further verified in the next data.
In Fig. 5, we analyze the case of a slight deviation from normal incidence, measuring the intensity of the side waves for θ i 0.03°. In this case, jηj reaches the largest value of 94% [ Table 2 in Fig. 5(b)]. We ascribe this increased spin-directive coupling to the breaking of the system symmetry due to a non-symmetric excitation condition, which favors the observed chiral coupling. Indeed, the BIC transverse SAM has components also in the PhCM plane (Fig. 2), which may thus influence the overall efficiency of light coupling. Further research is necessary to provide a general comprehension of the phenomenon.
As previously discussed, the side waves follow the axes of the PhCM. The second row in Fig. 5(a) shows the effect of rotating the PhCM about the z-axis of an angle ϕ 0 30°. The rotation with respect to the input polarizer, defining the reference, may be seen as the simple transformation of the geometric phase ϕ 0 j ϕ j ϕ 0 for each redirected side wave. As visible in the bottom panel of Fig. 5(a) (ϕ 0 30°), the intensity curves appear shifted in β because of the new ϕ 0 j and are characterized by a different overall shape with respect to the case ϕ 0 0°(top panel). This can be appreciated, for instance, by the different values measured at β 0°. The agreement with the model fit is also in this case very good, with an estimated directivity jηj 99.9% (Table 2).

CONCLUSION
Concluding, the optical field at the BIC point of a nanoscale PhCM is characterized by strong transverse SAM components that are induced by the tight confinement of the electromagnetic field and that can be used to devise macroscopic spin-orbit interactions. Experimentally, we have shown that the BIC resonantly feeds the fundamental Bloch's surface waves of the mode, providing macroscopic spin-momentum-locked propagation along the symmetry axes of the PhCM. Since a resonator in the BIC regime is characterized by strong near-field amplification, engineering such high Q-factor resonators may lead to metasurface multiplatform applications, also favored by the large-scale collective character of the BIC that makes the resonator intrinsically robust against external perturbations. Concurrently, other effects may be enhanced, such as unidirectional low-threshold lasing and strong-coupling regime phenomena, providing novel quantum chiral photonic applications.
Funding. Ministero dell'Istruzione, dell'Università e della Ricerca (PRIN 2017 0002501).  Table 2 reports the fit parameters with a chiral directivity η close to 1. S/P indicates that the values are respectively referred to s-pol or p-pol of (a).