Bound states in the continuum in anisotropic photonic crystal slabs

We investigate the bound states in the continuum (BICs) in photonic crystal slabs composed of alternating anisotropic and isotropic dielectric materials. According to the orientation of optical axis plane, three different configurations are proposed for analyzing various types of BICs, associated with extremely large quality factors and vanishing spectral linewidths. In particular, symmetry-protected (SP) BICs exist at the Brillouin zone center for zero rotation angle of the optical axis, which exhibit antisymmetric field patterns that are decoupled from the symmetric radiating fields. Accidental BICs and Friedrich-Wintgen (FW) BICs also occur at the Brillouin zone center for particular rotation angles of the optical axis. The former emerge on isolated bands with quasi-symmetric or quasi-antisymmetric field patterns, while the latter appear near the avoided crossing between two dispersion bands. At off the Brillouin zone center, SP BICs do not exist while accidental BICs and FW BICs appear at particular optical axis rotation angles, with similar features but somewhat more asymmetric field patterns than those at the Brillouin zone center.

tuned for destructive interference of the radiative waves to occur 45,46 .The ordinary and extraordinary refractive indices in the anisotropic material may alter the strength of electric field components.The orientation of optical axis can further rotate the field components [47][48][49][50][51] .As a consequence, the anisotropic structure acquire more system parameters in the tuning process to form the BICs.
In the present study, we investigate the BICs in photonic crystal slabs that compose of alternating anisotropic and isotropic dielectric materials.To analyze various types of BICs in the underlying structure, three configurations according to the orientation of optical axis plane are proposed.Basic features of BICs as extremely large quality factors and vanishing spectral linewidths are examined in the dispersion bands and transmittance diagram.At the Brillouin zone center ( k x = 0 ), SP BICs exist for zero rotation angle of the optical axis, which exhibit antisymmetric field patterns that are decoupled from the symmetric radiating fields.Accidental BICs and FW BICs occur at particular rotation angles of the optical axis.The former emerge at certain wave vector points on isolated bands with quasi-symmetric or quasi-antisymmetric field patterns, while the latter appear near the avoided crossing between two dispersion bands.At off the Brillouin zone center ( k x = 0 ), SP BICs do not exist while accidental BICs and FW BICs appear at particular optical axis rotation angles.Compared to BICs at k x = 0 , BICs at k x = 0 possess similar features but somewhat more asymmetric field patterns.

Results
Anisotropic photonic crystal slabs.Consider a photonic crystal slab composed of alternating anisotropic and isotropic dielectric materials, as schematically shown in Fig. 1, where a is the lattice constant, h is the slab height, [ε 1 ] is the dielectric tensor of the anisotropic material with the portion ρ in the unit cell, and ε 2 is the dielectric constant of the isotropic material.Assume that the anisotropic material is uniaxial, with ε o and ε e being the ordinary and extraordinary dielectric constants, respectively.According to the orientation of optical axis plane, there exist three different configurations for the anisotropic material.Here, I, II, III refer to the configurations where the optical axis vector c lies on the xy, yz, xz planes, with φ I , φ II , φ III being the rotation angles of the vector c with respect to the x, y, z axes, respectively.As the photonic crystal slabs are connected structures similar to the dielectric gratings, no substrate is present 25,43 in order to give more concise resonance features associated with BICs.In the presence of a substrate, the basic features of BICs will retain except that they are red-shifted and the associated resonances are slightly less significant.
To identify various types of the BICs, we compute the eigenfrequencies by the free and open-source software package MPB (MIT Photonic Bands) based on the plane wave expansion method and the supercell approach 52 , which have been extensively used in the study of photonic crystal slabs.To analyze the features of BICs, we also calculate the transmittance (ratio of the transmitted to incident power) based on the finite element method, which has been employed to study the transmission characteristics of grating structures and photonic crystal slabs [53][54][55] .The Bloch boundary condition is applied at the unit cell boundary and perfectly matched layers are applied on top and bottom of the computational domain 56 .An important factor for evaluating BICs is the quality factor (Q factor) based on the eigenfrequency , which is defined as 57 , where Re[•] and Im[•] denote the real and imaginary parts, respectively.Another definition of the quality factor is based on the spectral response as Q = f r /�f r , where f r is the resonance frequency and f r is the linewidth (full width at half maximum) of the transmittance (or reflectance).The above two definitions of the quality factor are nearly equivalent and amount to the ratio of energy stored in the resonator to the energy dissipated per cycle by damping processes 58 .A larger quality factor corresponds to a narrower linewidth with a smaller damping rate and a longer lifetime for the resonator.The higher the quality factor, the stronger the interaction between the confined light field and the material or structure.Theoretically, the linewidth of an ideal BIC becomes zero and the quality factor tends to infinity 22,59 .
In the present study, we report the results of BICs in the anisotropic photonic crystal slabs for transverse magnetic (TM) polarization.Here, TM refers to the case where the magnetic field is perpendicular to the longitudinal plane formed by the crystal lattice and the slab height, that is, the xy plane.The BICs in isotropic photonic crystal slabs for TE polarization, where the electric field is perpendicular to the longitudinal plane, can be found in Ref. 25 .
1. BICs at the Brillouin zone center ( k x = 0) Configuration I.In this configuration, the optical axis vector lies on the xy plane (cf.Fig. 1) and the permittivity tensor of the anisotropic medium in Cartesian coordinate system is given by 60 where [ ε I ] = diag[ε e , ε o , ε o ] and is the transformation matrix used to perform a rotation of angle φ I about the z axis 61 .Here, diag[•, •, •] denotes a 3 × 3 diagonal matrix.
Figure 2a shows the normalized frequencies ( a/ ) of four relevant dispersion bands for the anisotropic pho- tonic crystal slab in configuration I at the Brillouin zone center ( k x = 0 ) as functions of the rotation angle φ I .The frequencies are lowered as the rotation angle increases.In Fig. 2b, extremely large quality factors are attained at their peaks on band 1 (blue line) and band 3 (green line) for zero rotation angle ( φ I = 0 • ), which are labeled by (1) 1 and 2 , respectively.Because of symmetry in the ellipse representation of the anisotropic medium, the quality factor curves are symmetric about φ I = 90 • and therefore two similar peaks will occur at φ I = 180 • .In addition, the quality factor curves in the range of φ I = 0 • to 180 • coincide with those in the range of φ I = 180 • to 360 • .The band structure near the Brillouin zone center for φ I = 0 • is shown in Fig. 2c.The transverse magnetic fields ( Re[H z ] ) of the eigenmodes at 1 and 2 , as plot in Fig. 2f, display antisymmetric (TM 21 -like and TM 22 -like) patterns, which is characteristic of the SP BICs 25,62,63 .Here, the pattern symmetry is viewed along the periodic direction (x axis) 25 , and the magnetic fields are plot in one unit cell that consists of the anisotropic medium of width ρa in the middle and the isotropic medium of width (1 − ρ)a/2 on the left and right sides.The TM mn -like refers to the field pattern (in the unit cell) similar to the mode structure (in the cross section) in rectangular waveguides, where m and n are mode numbers in the horizontal and vertical directions, respectively 58 .The transverse magnetic fields of the BICs are well confined in the slab because of the dielectric contrast between inside and outside the slab, which are similar to those in metallic waveguides.
The SP BICs originates from the symmetry incompatibility and necessarily have an antisymmetric field pattern [ H z (x, y) = −H z (−x, y) ], so that they will not couple to the symmetric radiating fields.For this reason, the SP BICs always occur at k x = 0 and cannot exist at k x = 0 25 .Note that the SP BICs also occur at φ I = 90 • [cf.blue and green lines in Fig. 2b] with very similar quality factors and transverse magnetic field patterns, where the optical axis vector is perpendicular to that for φ I = 0 • .This feature has also been observed in semi-infinite 1D photonic crystals with anisotropic defect layers 48 .
The BICs at k x = 0 correspond to the bound states that occur at normal incidence.To evaluate the BICs from the aspect of electromagnetic radiation, the transmittance as a function of the frequency and rotation angle φ I is plot in Fig. 2i.The dispersion bands [cf.Fig. 2a] are overlaid in the same diagram to show consistent trends with the transmittance.Theoretically, the SP BICs are ideal bound states with infinite quality factors and zero spectral linewidths 22,59 .It is therefore hard to identify such BICs merely from the transmittance diagram without the help of dispersion bands.In fact, the ideal SP BICs cannot be excited from the far field radiation, as the former are decoupled from the latter.As a result, the corresponding transmittance is equal to unity, for no interactions occur between the incident field and the structure.In practice, a certain asymmetry either in the structure or incident angle allows for a small coupling between the bound state and incident field, which turns the ideal BICs into quasi-BICs.The transmittance will be somewhat less than unity and usually featured with the Fano resonance 49 .In this situation, the quality factor of a quasi-BIC is finite (though very large) and its spectral linewidth is not zero (though very small).This is also the circumstance when the quasi-BICs are used for practical applications, for instance, in sensors or filters.
In Fig. 2b, very large quality factors are attained at their peaks on band 1 (blue line) and band 4 (orange line) for particular rotation angles ( φ I = 49.8 • and 37.8 • ), which are labeled by 3 and 4 , respectively.The correspond- ing band structures near the Brillouin zone center are shown in Fig. 2d and e.Note that the eigenmodes with quality factor peaks occur on isolated bands 64 , which is characteristic of the accidental BICs 25,40,41 .This type of BICs are also named single-resonance parametric BICs 1,22 , which usually occur at oblique incidence 65 and are often referred to as off-Ŵ BICs 64 .It is, however, possible to create the accidental BICs at k x = 0 (or under normal incidence) because of the interference between two or more modes, which have been observed in 1D photonic crystal slabs 25 .In the present study, the accidental BICs also occur at k x = 0 for particular rotation angles of the anisotropic material.The transverse magnetic fields ( Re[H z ] ) of the eigenmodes at 3 and 4 are plot in Fig. 2g  and h, respectively, which display a quasi-antisymmetric (TM 21 -like) or quasi-symmetric (TM 32 -like) pattern.Because of the material anisotropy, no strict symmetry is presented in the field distribution, although the symmetric or antisymmetric pattern is almost, but not exactly, fulfilled.In photonic crystal slabs, the accidental BICs originate from the interference between two or more Bloch waves bouncing back and forth vertically inside the slab 25 .It can be examined in Fig. 2i that the accidental BICs correspond to the transmittance peaks with vanishing linewidths 29 .This feature is consistent with the fact that most electromagnetic fields of the BICs are confined in the structure [cf.Fig. 2f-h] and will not radiate to the far field.Equivalently, the BICs are decoupled from the far field radiation.
Configuration II.In this configuration, the optical axis vector lies on the yz plane (cf.Fig. 1) and the permittivity tensor of the anisotropic medium in Cartesian coordinate system is given by 60 where is the transformation matrix used to perform a rotation of angle φ II about the x axis 61 .This configuration has also been used in the study of BICs in semi-infinite 1D photonic crystals with anisotropic defect layers 48,49 .
Figure 3a shows the normalized frequencies ( a/ ) of four relevant dispersion bands for the anisotropic photonic crystal slab in configuration II at the Brillouin zone center ( k x = 0 ) as functions of the rotation angle φ II .The frequencies are raised as the rotation angle increases.In Fig. 3b, a very large quality factor is attained at its peak on band 3 (green line) for a particular rotation angle ( φ II = 35.7 • ), which is labeled by 1 .The corre- sponding band structure near the Brillouin zone center is shown in Fig. 3c.Note that the avoided crossing occurs between band 3 (green line) and band 4 (orange line) (with an enlarged view in the inset), which is characteristic of the FW BICs 29,63 .In particular, the quality factor for one of the two bands goes to a very large value near the avoided crossing point, while the quality factor for the other band is much smaller [cf.Fig. 3b].The transverse magnetic fields ( Re[H z ] ) of the eigenmodes near 1 with the quality factor peak, as plot in Fig. 3d, display a bound www.nature.com/scientificreports/state on one band (green line) and a leaky state on the other band (orange line) 63 .The former is well localized in the slab with a symmetric (TM 33 -like) pattern, while the latter is radiative outside the slab.Unlike the accidental BICs that are attributed to single resonances, the FW BICs are formed by tuning the interaction between two or more resonances 39 .This property can be understood from the coupled mode theory for a simple system with two resonances: i∂A/∂t = HA 1 , where A = (A 1 , A 2 ) T are the amplitudes of resonances and is the Hamiltonian of the system.Here, ω i and γ i ( i = 1, 2 ) are the resonant frequencies and damping rates, respectively, of the i-th resonance, and κ is the coupling strength between the two resonances.If the following condition is satisfied 1 : the eigensystem for the Hamiltonian H is solved to give a real eigenvalue ω r = (γ 1 ω 2 − γ 2 ω 1 )/(γ 1 − γ 2 ) and a complex eigenvalue The eigenvalue with a purely real part turns into a BIC and the other eigenvalue, with the imaginary part equal to the sum of two damping rates, becomes more lossy.This is the typical feature of FW BICs 6 .Note that the above solutions for the FW BICs are valid when ω 1 ≈ ω 2 and γ 1 ≈ γ 2 [cf.Eq. ( 6)], that is, the frequencies and damping rates of the two resonances are roughly equal.
Figure 3e shows the transmittance diagram as a function of the normalized frequency and rotation angle φ II .The dispersion bands [cf.Fig. 3a] are overlaid in the same diagram to show consistent trends with the transmittance.The feature of avoided crossing can be observed as two resonances approach each other at a particular rotation angle.To illustrate the FW BIC in a more clear manner, the transmittance curves as functions of the normalized frequency for selected rotation angles are plot in Fig. 3f.Note that the FW BIC [corresponding to the peak of quality factor in Fig. 3b] occurs near the transmittance peak, which is accompanied with a transmittance dip to form the feature of Fano resonance.The FW BIC can be located around or even right at the transmittance peak 24,39,63,64 .As the destructive interference in the FW BIC is hardly perfect in real structures, the corresponding transmittance may be slightly less than (though very close to) unity and the FW BIC has a finite yet very large quality factor.In isotropic photonic crystal slabs, the FW BICs occur at k x = 0 25 .In the present study, however, the FW BICs are observed at k x = 0 , as a consequence of additional degree of freedom in the system introduced by the anisotropic material parameters.
(5) www.nature.com/scientificreports/Configuration III.In this configuration, the optical axis vector lies on the xz plane (cf.Fig. 1) and the permittivity tensor of the anisotropic medium in Cartesian coordinate system is given by 60 where is the transformation matrix used to perform a rotation of angle φ III about the y axis 61 (cf.Fig. 1).
Figure 4a shows the normalized frequencies ( a/ ) of four relevant dispersion bands for the anisotropic photonic crystal slab in configuration III at the Brillouin zone center ( k x = 0 ) as functions of the rotation angle φ III .The frequencies do not change much as the rotation angle increases.In Fig. 4b, a very large quality factor is attained at its peak on band 4 (orange line) for a particular angle ( φ III = 73.3• ), which is labeled by 1 .The corresponding band structure near the Brillouin zone center is shown in Fig. 4c.Note that the eigenmode with the quality factor peak occurs on an isolated band, which is characteristic of the accidental BICs similar to those in configuration I [cf.Figs.2d and e].The transverse magnetic field ( Re[H z ] ) of the eigenmode at 1 , as plot in Fig. 4d, displays a quasi-symmetric (TM 32 -like) pattern, which is similar to that configuration I [cf. in Fig. 2h].
The transmittance diagram as a function of the normalized frequency and rotation angle φ III is plot in Fig. 4e.The dispersion bands [cf.Fig. 4a] are overlaid in same diagram to show consistent trends with the transmittance.To illustrate the accidental BIC in a more clear manner, the transmittance curves as functions of normalized frequency for selected rotation angles are plot in Fig. 4f.Note that the accidental BIC, which corresponds to the peak of quality factor in Fig. 4b, emerges at a frequency where the transmittance curve becomes flat and varies smoothly across the resonance without a noticeable change.The accidental BIC is associated with a high transmittance nearly equal to unity 29,40 , as a consequence of destructive interference in the formation of a BIC.
2. BICs off the Brillouin zone center (k x = 0) Configuration I. Figure 5a shows the band structure of the anisotropic photonic crystal slab in configuration I for zero rotation angle ( φ I = 0 • ).In Fig. 5b and c, very large quality factors are attained at their peaks on two isolated bands (blue and yellow lines) at particular wave numbers [ k x a/(2π) = 0.195, 0.085 ], which are labeled by 1 and 2 , respectively.At k x = 0 , the resonance modes do not present strict symmetry, even though they look approximately symmetric or antisymmetric, which will be in general coupled to the radiation in a continuous frequency spectrum.However, the radiative leakage may be exactly suppressed at certain wave vector components and disappears 'accidentally' because of the destructive interference between different channels, leading to the accident BICs 22 .These BICs are basically similar to those in Fig. 2, except that the former occur at a particular (7)   5b  and c].On the other hand, the avoided crossing occurs between two upper bands (orange and purple lines) at a particular wave number [ k x a/(2π) = 0.263 ], which is labeled by 3 .In Fig. 5e, a very large quality factor is attained at its peak on one of the two bands (orange line) that form the avoided crossing, which is characteristic of the FW BIC similar to the case for k x = 0 [cf.Fig. 3].The transverse magnetic fields ( Re[H z ] ) of the eigenmodes near 3 display a bound state with a quasi-antisymmetric (TM 21 -like) pattern and a leaky state that largely loses the symmetry [cf.insets in Fig. 5e].
The band structure of the anisotropic photonic crystal slab at a particular rotation angle ( φ I = 57 • ) is shown in Fig. 5d.The avoided crossing occurs between two upper bands (orange and purple lines) at a particular wave number [ k x a/(2π) = 0.107 ], which is labeled by 4 .Similar features of the FW BIC are shown in Fig. 5f, in which a very large quality factor is attained at its peak on one of the two bands (orange line), and the transverse magnetic fields ( Re[H z ] ) near the eigenmode at 4 display a bound state with a quasi-symmetric (TM 32 -like) pattern and a leaky state with a much less concentrated field distribution [cf.insets in Fig. 5f].
Configuration II. Figure 6a and b show the band structures of the anisotropic photonic crystal slab in configuration II for zero and a particular rotation angle ( φ II = 0 • and 45 • ), respectively.In Fig. 6c and d, very large quality factors are attained at their peaks on two isolated bands (blue and purple lines) at particular wave numbers [ k x a/(2π) = 0.0501, 0.138 ] for φ II = 0 • , which are labeled by 1 and 2 , respectively.The two eigen- modes have the characteristics of accidental BICs similar to those in configuration I.The transverse magnetic fields ( Re[H z ] ) of the eigenmodes at 1 and 2 display quasi-antisymmetric (TM 21 -like and TM 22 -like) patterns [cf.insets in Fig. 6c and d].
In Fig. 6e and f, very large quality factors are attained at their peaks on two isolated bands (blue and orange lines) at particular wave numbers [ k x a/(2π) = 0.087, 0.088 ] for φ II = 45 • , which are labeled by 3 and 4 , respectively.The two eigenmodes also have the characteristics of accidental BICs.The transverse magnetic fields ( Re[H z ] ) of the eigenmodes at 3 and 4 display quasi-antisymmetric (TM 21 -like and TM 22 -like) patterns [cf.insets in Fig. 6e and f].On the other hand, the avoided crossings occur between two pairs of two interacting bands (green and orange lines, purple and brown lines) at particular wave numbers [ k x a/(2π) = 0.175, 0.244 ], which are labeled by 5 and 6 , respectively.In Fig. 6 g and h, very large quality factors are attained at their peaks near the avoided crossing points, having the characteristics of FW BICs.The transverse magnetic fields ( Re[H z ] ) of the eigenmodes near 5 display a bound state with a quasi-antisymmetric (TM 21 -like) pattern and a leaky state, while the transverse magnetic fields for the eigenmodes near 6 display a bound state with a quasi-symmetric (TM 32 -like) pattern and a leaky state [cf.insets in Fig. 6g and h].
Configuration III. Figure 7a and b show the band structures of the anisotropic photonic crystal slab in configuration III for zero and a particular rotation angle ( φ III = 0 • and 50 • ), respectively.In Fig. 7c-e, very large quality factors are attained at their peaks on three isolated bands (blue, yellow, purple lines) at particular wave numbers [ k x a/(2π) = 0.194, 0.085, 0.0251 ] for φ III = 0 • , which are labeled by 1 , 2 , 3 , respectively.The three  7c-e].Note that the symmetry in higher order modes is less obvious, especially in the case of k x = 0 .On the other hand, the avoided crossing occurs between two upper bands (orange and purple lines) at a particular wave number [ k x a/(2π) = 0.262 ], which are labeled by 4 .In Fig. 7f, a 7g-i].On the other hand, the avoided crossing occurs between two upper bands (green and orange lines) at a particular wave number [ k x a/(2π) = 0.15 ], which is labeled by 8 .In Fig. 7j, a very large quality factor is attained at its peak near the avoided crossing point, having the characteristics of FW BIC similar to those in configurations I and II.The transverse magnetic fields ( Re[H z ] ) of the eigenmode near 8 display a bound state with a quasi-symmetric (TM 31 -like) pattern and a leaky state [cf.insets in Fig. 7j].
Based on the results reported in the above two subsections, various types of BICs in the photonic crystal slabs in the present study are summarized in Table 1 for k x = 0 and k x = 0 in three different configura- tions.Here, φ N = 0 • (N = I, II, III) refers to a particular angle of φ N not equal to zero degree.In the present study, the same height ratio h/a = 1.66 (see caption in Fig. 2) is used for different configurations, which is the smallest value that supports the occurrence of BICs for all three orientations of the optical axes for the anisotropic    material.Here, the ordinary and extraordinary dielectric constants of the anisotropic material do not correspond to real materials, for convenience of analyzing different types of BICs.They are, however, close to those used in Refs. 45,46,49.Finally, it is worthy of noting that some of the BICs demonstrated by numerical results without the support of theoretical models could be quasi-BICs even with extremely large quality factors.Nevertheless, it is not uncommon to use the terminology of BICs (might not be strictly true BICs) when the quality factors are very large 39,56,59,62 .
In conclusion, we have investigated the BICs in photonic crystal slabs composed of alternating anisotropic and isotropic dielectric materials.Three difference configurations according to the orientation of optical axis plane are proposed for analyzing the BICs based on the dispersion bands and transmittance diagram.Various types of BICs with very large quality factors and vanishing spectral linewidths are identified in the underlying structure.In particular, SP BICs associated with antisymmetric field patterns occur for zero rotation angle at the Brillouin zone center.Accidental BICs emerge on isolated bands while FW BICs appear near the avoided crossing between two interacting bands for zero or particular rotation angles at or off the Brillouin zone center.

Figure 1 .
Figure 1.Schematic diagram of the anisotropic photonic crystal slab, where the anisotropic and isotropic materials are shaded in green and gray colors, respectively.Insets are ellipse representations of the uniaxially anisotropic material for three different configurations.

1 Figure 2 .
Figure 2. BICs in the anisotropic photonic crystal slab at k x = 0 in configuration I, where ε o = 5.2 , ε e = 8.2 , ε 2 = 3.3 , ρ = 0.37 , and h/a = 1.66 . (a) Normalized frequencies and (b) quality factors as functions of the rotation angle φ I .(c)-(e) Band structures for φ I = 0 • , 49.8 • , 37.8 • , respectively.(f)-(h) Quality factors as functions of the wave number and transverse magnetic fields ( Re[H z ] ) of the SP BICs in (c) and the accidental BICs in (d), (e), respectively.Blue and red colors represent positive and negative values, respectively.(i) Transmittance diagram as a function of the normalized frequency and rotation angle.

Figure 3 .
Figure 3. BICs in the anisotropic photonic crystal slab at k x = 0 in configuration II, with the same geometry and material parameters as in Fig. 2. (a) Normalized frequencies and (b) quality factors as functions of the rotation angle φ II .(c) Band structure for φ II = 35.7 • .Inset is an enlarged view of two upper bands near k x = 0 , with the green dot corresponding to the FW BIC.(d) Quality factor as a function of the wave number and transverse magnetic field ( Re[H z ] ) of the FW BIC in (c).(e) Transmittance diagram as a function of the normalized frequency and rotation angle.(f) Transmittance curves as functions of the normalized frequency for selected rotation angles.

Figure 4 .
Figure 4. BICs in the anisotropic photonic crystal slab at k x = 0 in configuration III, with the same geometry and material parameters as in Fig. 2. (a) Normalized frequencies and (b) quality factors as functions of the rotation angle φ III .(c) Band structure for φ III = 73.3• .(d) Quality factor as a function of the wave number and transverse magnetic field ( Re[H z ] ) of the accidental BIC in (c).(e) Transmittance diagram as a function of the normalized frequency and rotation angle.(f) Transmittance curves as functions of the normalized frequency for selected rotation angles.
k x for φ I = 0 • , while the latter occur at k x = 0 for a particular φ I .The transverse magnetic fields ( Re[H z ] ) of the eigenmodes at 1 and 2 display quasi-antisymmetric (TM 21 -like and TM 22 -like) patterns [cf.insets in Fig.

Figure 5 .
Figure 5. BICs in the anisotropic photonic crystal slab at k x = 0 in configuration I, with the same geometry and material parameters as in Fig. 2. (a) and (d) Band structures for φ I = 0 • and 57 • , respectively.(b), (c) and (e) Quality factors as functions of the wave number and transverse magnetic fields ( Re[H z ] ) of the accidental BICs at k x a/(2π) = 0.195, 0.085 and the FW BIC at k x a/(2π) = 0.263 , respectively, for φ I = 0 • .(f) Quality factor and transverse magnetic field of the FW BIC at k x a/(2π) = 0.107 for φ I = 57 • .
very large quality factor is attained at its peak near the avoided crossing point, having the characteristics of FW BICs similar to those in configurations I and II.The transverse magnetic fields ( Re[H z ] ) of the eigenmodes near 4 display a bound state with a quasi-antisymmetric (TM 21 -like) pattern and a leaky state [cf.insets in Fig. 7f].Note also that the symmetry pattern is less obvious for k x = 0.In Fig. 7g-i, very large quality factors are attained at their peaks on three isolated bands (blue, green, brown lines) at particular wave numbers [ k x a/(2π) = 0.161, 0.0251, 0.0711 ] for φ III = 50 • , which are labeled by 5 , 6 , 7 , respectively.The three eigenmodes have the characteristics of accidental BICs similar to those in configura- tions I and II.The transverse magnetic fields ( Re[H z ] ) of the eigenmodes at 5 , 6 , 7 display quasi-antisymmetric (TM 21 -like, TM 22 -like, TM 23 -like) patterns [cf.insets in Fig.

Figure 7 .
Figure 7. BICs in the anisotropic photonic crystal slab at k x = 0 in configuration III, with the same geometry material parameters as in Fig. 2. (a) and (d) Band structures for φ III = 0 • and 50 • , respectively.(c-f) Quality factors and transverse magnetic fields ( Re[H z ] ) of the accidental BICs at k x a/(2π) = 0.194, 0.085, 0.0251 and the FW BIC at k x a/(2π) = 0.262 for φ III = 0 • .(g-j) Quality factors and transverse magnetic fields of the accidental FW BICs at k x a/(2π) = 0.161, 0.0251, 0.0711 and the FW BIC at k x a/(2π) = 0.15 , respectively, for φ III = 50 • .

Table 1 .
Summary of BICs in asymmetric photonic crystal slabs.