## Abstract

We investigate the diffraction of the guided modes of a dielectric slab waveguide on a simple integrated structure consisting of a single dielectric ridge on the surface of the waveguide. Numerical simulations based on aperiodic rigorous coupled-wave analysis demonstrate the existence of sharp resonant features and bound states in the continuum (BICs) in the reflectance and transmittance spectra occurring at the oblique incidence of a transverse-electric (TE)-polarized guided mode on the ridge. Using the effective index method, we explain the resonances by the excitation of cross-polarized modes of the ridge. Formation of the BICs are confirmed using a theoretical model based on coupled-wave theory. The model suggests that the BICs occur due to the coupling of quasi-TE and quasi-transverse-magnetic modes of the structure. Simple analytical expressions for the angle of incidence and the ridge width predicting the location of the BICs are obtained. The existence of high-$Q$ resonances and BICs enables using the considered integrated structure for sensing, transformation of optical signals, and enhancing nonlinear light–matter interactions. Due to the Lorentzian line shape of the resonances near the BICs, the structure is also promising for filtering applications.

© 2018 Chinese Laser Press

## 1. INTRODUCTION

In a wide class of planar (integrated) optoelectronic systems, processing of the optical signal is performed in a slab waveguide [1–11]. Such a geometry corresponds to the “insulator-on-insulator” platform and is suitable for the creation of fully integrated optical devices. In this case, the processed signal corresponds to a superposition of slab waveguide modes having different propagation directions (in the case of spatial processing) or different frequencies (in the case of spectral processing). In this regard, the design of integrated resonant photonic structures is of great interest [12,13].

In the last few years, much attention has been paid to the investigation of so-called bound states in the continuum (BICs) in resonant photonic structures (see the recent review [14] and references therein). BICs are eigenmodes that have an infinite lifetime (and an infinitely high quality factor), although they are supported by a structure having open scattering channels. The mode leakage of these channels is canceled either due to symmetry reasons or by means of parameter tuning [14]. In photonics, BICs were studied in periodic structures (diffraction gratings, photonic crystal slabs, and infinite arrays of dielectric rods or spheres) [15–23], defects and interfaces of photonic crystals [24–26], arrays of optical waveguides [27–29], and leaky optical waveguides [30], among others [14]. The structures supporting BICs can be used as “unconventional” narrowband waveguides [30]. A slight deviation from the BIC condition enables obtaining very high-$Q$ resonators, which leads to various potential applications including lasers, sensors, and filters.

In this work, we investigate the resonant optical properties of a very simple structure consisting of a single subwavelength or near-wavelength ridge on the surface of a single-mode dielectric slab waveguide. We analytically and numerically demonstrate that, in the case of diffraction of a transverse-electric (TE)-polarized mode of the waveguide, the studied structure exhibits BICs and high-$Q$ resonances associated with the excitation of cross-polarized modes of the ridge. Comparing to Ref. [30], we present a theoretical model that accurately describes the optical properties of the studied structure. In particular, we prove the existence of robust BICs and obtain simple closed-form expressions predicting their locations. The discovered high-$Q$ resonances and BICs supported by the presented planar structure make it promising for various applications, including the transformation of optical signals, sensing, enhancing nonlinear light–matter interactions, and spatial (angular) and spectral filtering. The latter is possible due to the Lorentzian line shape of the resonances.

The paper is organized as follows. In Section 2, we define the geometry and discuss the scattering channels of the considered integrated structure. In Section 3, we present full-wave numerical simulation results demonstrating the existence of high-$Q$ resonances in the reflectance and transmittance spectra. Using the effective index method, we qualitatively explain the formation of the resonances. By rigorously calculating the quality factor of the eigenmodes of the structure, we show that it supports BICs. Section 4 is dedicated to the derivation of a theoretical model describing the optical properties of the structure and explaining the BIC formation mechanism. Section 5 concludes the paper.

## 2. GEOMETRY AND SCATTERING CHANNELS

Let us consider a simple integrated structure shown in Fig. 1(a). The structure consists of a dielectric ridge with height ${h}_{r}>0$ and width $w$ located on the surface of a dielectric slab waveguide with thickness ${h}_{\mathrm{wg}}$. For simplicity, we consider the case when the waveguide core and the ridge on its surface are made of the same material with dielectric permittivity ${\u03f5}_{\mathrm{wg}}$. Additionally, we assume that the substrate is optically denser than the superstrate: ${\u03f5}_{\mathrm{sub}}>{\u03f5}_{\mathrm{sup}}$. In the present work, we restrict our consideration to the case when the waveguides with thicknesses ${h}_{\mathrm{wg}}$ and ${h}_{\mathrm{wg}}+{h}_{r}$ (corresponding to the regions outside and inside the ridge, respectively) are single-mode for both TE and transverse-magnetic (TM) polarizations. However, one can show that optical effects similar to the ones described below arise also in multimode slab waveguides with properly chosen parameters.

We study the diffraction of a TE-polarized guided mode with an effective refractive index ${n}_{\mathrm{TE}}$ by an “inclined” ridge, as shown in Fig. 1(a). The angle between the ridge and the wave front of the incident mode is denoted by $\theta $ and referred to as the angle of incidence.

In the general case, reflected and transmitted TE- and TM-polarized guided modes are generated upon diffraction of the incident mode on the ridge, as well as a continuum of non-guided waves in the superstrate and substrate propagating away from the waveguide core layer. However, it was recently demonstrated that, under certain conditions, polarization conversion and out-of-plane scattering in such “integrated” diffraction problems can be completely eliminated [10,11,31].

The considered structure is invariant to translation in the direction parallel to the ridge interfaces (in the $y$ direction). According to Maxwell’s equations, the corresponding tangential wave vector component of the incident wave ${k}_{y}={k}_{0}{n}_{\mathrm{TE}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta $ has to be conserved for all outgoing waves. Here, ${k}_{0}=2\pi /\lambda $ is the wave number, and $\lambda $ is the free-space wavelength. Therefore, if the incidence angle is chosen so that ${k}_{y}^{2}>{k}_{0}^{2}{\u03f5}_{\mathrm{sub}}$, this tangential component has greater magnitude than that of the wave vector of the propagating plane waves in the substrate. It means that at angles of incidence greater than the “critical” angle

all the waves scattered to the substrate are evanescent and do not carry energy. A similar expression can be written for the critical angle ${\theta}_{\mathrm{sup}}$ describing the scattering to the superstrate. The assumed inequality ${\u03f5}_{\mathrm{sub}}>{\u03f5}_{\mathrm{sup}}$ implies that ${\theta}_{\mathrm{sub}}>{\theta}_{\mathrm{sup}}$, and therefore no out-of-plane scattering to the superstrate occurs once $\theta >{\theta}_{\mathrm{sub}}$.A similar approach can be applied to eliminate TM-polarized reflected and transmitted modes (with effective refractive index ${n}_{\mathrm{TM}}$) generated upon diffraction of the incident TE-polarized mode. Indeed, since in a single-mode slab waveguide ${n}_{\mathrm{TE}}>{n}_{\mathrm{TM}}$ [32], a cutoff angle for the TM-polarized modes exists. Similar to the derivation of Eq. (1), one can show that at angles of incidence greater than

no propagating outgoing TM-polarized modes are generated.Note that ${n}_{\mathrm{TM}}>\sqrt{{\u03f5}_{\mathrm{sub}}}$, therefore, ${\theta}_{\mathrm{TM}}>{\theta}_{\mathrm{sub}}$. Hence, at $\theta >{\theta}_{\mathrm{TM}}$, only two scattering channels remain open in the considered structure, namely, reflected and transmitted TE-polarized modes.

We obtained the cutoff conditions for the waveguide of thickness ${h}_{\mathrm{wg}}$ concerning outgoing waves in the considered diffraction problem. In a similar way, the ridge region can be considered as a segment of a slab waveguide with greater thickness ${h}_{\mathrm{wg}}+{h}_{r}$. Denoting the effective index of the TM-polarized mode of this slab waveguide by ${n}_{\mathrm{TM},r}$, we can write the corresponding cutoff angle as

Therefore, in the angular range ${\theta}_{\mathrm{TM}}<\theta <{\theta}_{\mathrm{TM},r}$, the ridge region supports both TE- and TM-polarized modes, while the scattered field contains only two TE-polarized modes. As we show below, it is in this region that the structure exhibits remarkable optical properties.## 3. NUMERICAL SIMULATIONS: HIGH-$Q$ RESONANCES AND BICs

In this section, we numerically study the optical properties of a dielectric ridge on a slab waveguide. Simulations were carried out for the following parameters: free-space wavelength $\lambda =630\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm}$, dielectric permittivities of the superstrate, waveguide layer, and substrate ${\u03f5}_{\mathrm{sup}}=1$, ${\u03f5}_{\mathrm{wg}}={3.3212}^{2}$ (GaP), and ${\u03f5}_{\mathrm{sub}}={1.45}^{2}$, respectively. For these parameters, a slab waveguide supports only the fundamental TE- and TM-polarized modes if its thickness is in the range from 44 nm to 117 nm. In what follows, as a matter of example, we consider waveguide thickness ${h}_{\mathrm{wg}}=80\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm}$ and ridge height ${h}_{r}=30\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm}$. Thus, the slab waveguides with thicknesses ${h}_{\mathrm{wg}}$ and ${h}_{\mathrm{wg}}+{h}_{r}$ are both single-mode, and the assumptions in Section 2 are fulfilled.

By solving the dispersion equation of a slab waveguide [33], let us first calculate the effective refractive indices and cutoff angles [Eqs. (1)–(3)] discussed in the previous section. At the considered parameters, the slab waveguide outside the ridge (thickness ${h}_{\mathrm{wg}}$) supports a TE-polarized mode with effective refractive index ${n}_{\mathrm{TE}}=2.5913$ and a TM-polarized mode with effective refractive index ${n}_{\mathrm{TM}}=1.6327$. The slab waveguide corresponding to the ridge region (thickness ${h}_{\mathrm{wg}}+{h}_{r}$) supports TE- and TM-polarized modes with effective refractive indices ${n}_{\mathrm{TE},r}=2.8192$ and ${n}_{\mathrm{TM},r}=2.1867$, respectively. Critical angles of incidence corresponding to the closing of different scattering channels discussed above amount to ${\theta}_{\mathrm{sup}}=22.7\xb0$ (out-of-plane scattering to the superstrate), ${\theta}_{\mathrm{sub}}=34.0\xb0$ (out-of-plane scattering to the substrate), and ${\theta}_{\mathrm{TM}}=39.1\xb0$ (scattering to the TM-polarized guided mode). The cutoff angle of the TM-polarized mode in the ridge region equals ${\theta}_{\mathrm{TM},r}=57.6\xb0$.

#### A. Transmittance and Reflectance Spectra

In this subsection, we present rigorous simulation results for the investigated structure. Numerical simulations were performed using an efficient in-house implementation of the conical diffraction formulation of the rigorous coupled-wave analysis (RCWA) method [34,35] extended to integrated optics problems [36,37]. The RCWA, also called the Fourier modal method, is an established numerical technique for solving Maxwell’s equations. For simulation of the considered integrated structure within this approach, it was artificially periodized in the $z$ direction [36]. To eliminate optical interaction between adjacent periods, perfectly matched layers (PMLs) implemented as nonlinear coordinate transforms [37] were utilized. The simulations were performed with $2\times 350+1$ Fourier harmonics and with the thickness of the PML equal to 600 nm.

Figures 2(a) and 2(b) show the reflectance and transmittance of a TE-polarized mode upon diffraction by a ridge on the waveguide surface versus angle of incidence $\theta $ and ridge width $w$. For comparison, Figs. 2(c) and 2(d) show the reflectance and transmittance of a plane wave upon refraction by a thin film with thickness equal to ridge width $w$. In this plane-wave diffraction problem, effective refractive indices of the TE-polarized modes in the corresponding regions of the integrated structure are used as refractive indices of the superstrate and substrate (${n}_{\mathrm{TE}}$) and of the film (${n}_{\mathrm{TE},r}$) [see Fig. 1(b)]. The polarization of the incident plane wave coincides with the polarization of the incident mode in the “integrated” problem (i.e., the plane wave is polarized in the plane of incidence). The model diffraction problem in Fig. 1(b) is very similar to the effective index method (EIM), which is widely used for quick estimation of the effective refractive indices of the modes of dielectric photonic wire and rib waveguides [33,38].

It is evident in Fig. 2 that when $\theta >{\theta}_{\mathrm{TM},r}$, the reflectance and transmittance of the guided mode (obtained using RCWA) are very close to the reflectance and transmittance of the plane wave (obtained within the EIM approach). The periodically arranged reflectance and transmittance extrema visible in this region are due to Fabry–Perot resonances of the TE-polarized mode inside the ridge. In the region $\theta <{\theta}_{\mathrm{TM}}$, some discrepancies in transmittance are present, which are caused by scattering to the reflected and transmitted cross-polarized (TM-polarized) modes and to non-guided waves in the substrate (at $\theta <{\theta}_{\mathrm{sub}}$) and in the superstrate (at $\theta <{\theta}_{\mathrm{sup}}$). Nevertheless, in the considered angular range ($\theta <{\theta}_{\mathrm{TM}}$ or $\theta >{\theta}_{\mathrm{TM},r}$), the EIM can be used to obtain reasonable estimates of the solution of the “integrated” diffraction problem.

At angles of incidence ${\theta}_{\mathrm{TM}}<\theta <{\theta}_{\mathrm{TM},r}$, the situation is drastically different. In this angular range, sharp resonant peaks and dips, which are not explained by the plane-wave model, are present in the reflectance and transmittance spectra shown in Figs. 2(a) and 2(b), respectively.

#### B. High-$Q$ Resonances

Starting from this subsection, we focus on the angular range ${\theta}_{\mathrm{TM}}<\theta <{\theta}_{\mathrm{TM},r}$, where the structure exhibits the most interesting optical properties, namely, narrow resonances with high quality factor.

To reveal the nature of these high-$Q$ resonances, let us investigate field distribution in the structure at resonance conditions. Electric field distributions in the integrated structure corresponding to two points depicted with white asterisks in Fig. 2(a) are shown in Figs. 3(a) and 3(b). It is evident in Fig. 3 that at resonances, the field distributions possess a strong ${E}_{z}$ component, which is absent in the incident TE-polarized mode. This, along with the fact that the resonance region is bounded by the cutoff angles of TM-polarized modes outside and inside the ridge, suggests that the resonances are associated with excitation of cross-polarized (quasi-TM) modes of the ridge, which, in this case, acts as a leaky rib waveguide [39,40].

To verify this assumption, let us estimate the effective refractive indices of the quasi-TM modes of the ridge using the EIM [see Fig. 1(b); the used refractive indices are shown in parentheses] [33,38]. Under this approach, the estimates of the wave numbers ${k}_{y,\mathrm{TM},\mathrm{rib}}$ can be found by solving the dispersion relation of the symmetric three-layer slab waveguide [33]

where $k=\sqrt{{k}_{0}^{2}{n}_{\mathrm{TM},r}^{2}-{k}_{y,\mathrm{TM},\mathrm{rib}}^{2}}$ and $\gamma =\sqrt{{k}_{y,\mathrm{TM},\mathrm{rib}}^{2}-{k}_{0}^{2}{n}_{\mathrm{TM}}^{2}}$. The excitation condition of the found rib waveguide modes, which relates $w$ and $\theta $, reads as ${k}_{y,\mathrm{TM},\mathrm{rib}}={k}_{0}{n}_{\mathrm{TE}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta $. The values $(w,\theta )$ satisfying this condition are shown with dashed red curves in Fig. 2(b), which are close to the location of the resonances. Thus, the observed high-$Q$ resonances are indeed associated with excitation of cross-polarized modes of the ridge.#### C. Bound States in the Continuum

Figure 4(a) shows a magnified fragment of Fig. 2(a) corresponding to the angular range ${\theta}_{\mathrm{TM}}<\theta <{\theta}_{\mathrm{TM},r}$. It is evident that the angular width (and, consequently, the quality factor) of the resonances strongly varies along the dispersion curves, apparently vanishing at certain $(w,\theta )$ values marked with white circles in Fig. 4(a). This phenomenon indicates the potential presence of BICs in the considered structure. Let us recall that BICs are eigenmodes of the structures with open diffraction channels, which, however, have an infinite quality factor and a real frequency due to symmetry reasons or parameter tuning. In the considered structure, the bound states arise in the continuum of guided modes supported by the slab waveguide with thickness ${h}_{\mathrm{wg}}$ [see Fig. 1(a)].

In order to determine whether the structure actually supports BICs, let us investigate the quality factor of the modes in the region shown by the white rectangle in Fig. 4(a). In Figs. 5(a) and 5(b), we present resonance angle $\theta $ and quality factor $Q$ of the eigenmodes of the ridge, which were calculated using the RCWA-based approach proposed by some of the present authors in Refs. [41,42] (see the details in Appendix A). It is evident in Fig. 5 that at $w=344\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm}$ and $\theta =53.06\xb0$, the quality factor diverges. This indicates the presence of a BIC.

The BICs in the considered structure are robust, i.e., when the geometrical parameters of the structure are slightly changed, the BICs continue to exist at a slightly different wavelength $\lambda $ and angle $\theta $. This is confirmed by the rigorous numerical simulation results. For example, if we change the ridge width from $w=344\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm}$ to $w=343\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm}$, the BIC position will shift from $(\lambda ,\theta )=(630\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm},53.06\xb0)$ to $(\lambda ,\theta )=(625\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm},53.58\xb0)$. Similarly, at a larger ridge width $w=345\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm}$, the BIC emerges at $(\lambda ,\theta )=(635\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm},52.55\xb0)$.

Figure 5(c) shows the angular reflectance spectra of the structure at different ridge widths. According to Figs. 5(b) and 5(c), the quality factor of the resonance increases when approaching the BIC, which results in narrowing of the resonant reflectance peak. The full width at half maximum (FWHM) of the reflectance peak in Fig. 5(c) varies from 0° (at the BIC) to 0.1°.

It is important to note that the resonances in Fig. 5(c) have symmetric Lorentzian line shape in contrast to the general, essentially asymmetric, Fano line shape of the resonances that usually occur near the BICs [22,23]. Lorentzian line shape emerges when the non-resonant scattering pathway is negligible, i.e., no reflection takes place in non-resonant conditions. Indeed, this is the case, since the contrast between the effective refractive indices inside and outside the ridge (${n}_{\mathrm{TE}}$ and ${n}_{\mathrm{TE},r}$) is relatively low, and therefore almost no reflection occurs away from the resonances. Symmetric line shape on a zero background is essential when designing optical filters. The plots in Fig. 5(c) suggest that the proposed structure can be used as a narrowband spatial (angular) filter. Further calculations not presented here show that the same structure can also be utilized as a spectral (frequency) filter.

In contrast to the positions of the resonances, the variation of the quality factor and the appearance of the BICs cannot be explained by an EIM-based approach. Let us note that the existence of BICs in a similar structure was predicted in Ref. [30], where the ridge was considered as a photonic rib waveguide with an unconventional (BIC-based) guiding mechanism. In contrast to Ref. [30], we consider this structure as an integrated optical element for slab waveguide modes and, in the next section, derive a simple and accurate coupled-wave model describing its reflectance and transmittance spectra, predicting the high-$Q$ resonances in the ridge, and confirming the existence of BICs.

## 4. COUPLED-WAVE ANALYSIS OF THE RESONANCES

In this section, we present a simple model based on coupled-wave equations, which explains the formation of high-$Q$ resonances and predicts the location of BICs supported by the considered structure.

#### A. Coupled-Wave Equations

In order to obtain coupled-wave equations, let us revisit the effective index model in Fig. 1(b), where the integrated structure shown in Fig. 1(a) is replaced with a $z$-invariant slab. In contrast to this model considered above in Subsection 3.A, here we will take into account the cross-polarization mode coupling that takes place in the considered integrated structure. To do this, we represent the field inside the slab as TE and TM plane waves, which are coupled at the interfaces. The field over and under the slab contains the incident, reflected, and transmitted TE waves only (see Fig. 6). The wave numbers of the plane waves are defined by the effective refractive indices: ${n}_{\mathrm{TE},r}$ and ${n}_{\mathrm{TM},r}$ for the waves inside the slab; ${n}_{\mathrm{TE}}$ for the waves outside the slab.

According to Fig. 6, let us introduce the following plane waves: the incident wave (with complex amplitude $I$); reflected wave ($R$); TE and TM waves propagating downwards inside the slab (${U}_{1}$ and ${V}_{1}$); and upward-propagating TE and TM waves inside the slab (${U}_{2}$ and ${V}_{2}$). Amplitudes ${U}_{1}$ and ${V}_{1}$ are defined at the lower interface of the slab; amplitudes of the same waves at the upper interface equal ${U}_{1}{e}^{-i\varphi}$ and ${V}_{1}{e}^{-i\psi}$. Similarly, amplitudes ${U}_{2}$ and ${V}_{2}$ are defined at the upper interface of the slab; amplitudes of these waves at the lower interface are ${U}_{2}{e}^{-i\varphi}$ and ${V}_{2}{e}^{-i\psi}$. Phases $\varphi $ and $\psi $ are defined through the effective refractive indices of the TE and TM waves inside the slab, respectively

At the upper interface, the incident and scattered plane waves are related by a $3\times 3$ scattering matrix

By applying the same scattering matrix at the lower slab interface, we obtain the following set of coupled-wave equations:

By eliminating ${U}_{1}$, ${U}_{2}$, ${V}_{1}$, and ${V}_{2}$ from Eq. (7), we represent the reflection coefficient $R$ in the following form:

Here, the denominator is the product of the following two terms: whereas the numerator $\mathcal{N}$ reads asThe modes of the structure correspond to the poles of the reflection coefficient [41,43]. Therefore, we can obtain the dispersion equation of the modes by equating the denominator in Eq. (8) to zero:

The fact that the denominator is a product of two terms means that the considered structure supports two uncoupled sets of modes, namely, the symmetric and antisymmetric modes of the structure. One can calculate these modes independently by solving the equations ${\mathcal{D}}_{1}=0$ and ${\mathcal{D}}_{2}=0$.Let us consider the antisymmetric modes defined by equation ${\mathcal{D}}_{1}=0$. If no cross-polarization coupling occurs (if ${r}_{c}=0$), this equation has two independent solutions, ${e}^{-i\varphi}+{r}_{1}=0$ and ${e}^{-i\psi}+{r}_{2}=0$, which correspond to the antisymmetric TE and the antisymmetric TM modes, respectively. If ${r}_{c}\ne 0$, these modes are coupled and form antisymmetric modes with a more complex polarization state. Similarly, the equation ${\mathcal{D}}_{2}=0$ describes the coupling of symmetric TE modes with symmetric TM modes.

It is noteworthy that the equation ${\mathcal{D}}_{1}{\mathcal{D}}_{2}=0$ can be written in the following matrix form:

#### B. Bound States in the Continuum

Here, we show that the coupled-wave model derived in the previous subsection predicts BICs in the considered structure. In what follows, we derive explicit expressions for frequency $\omega $ and wave number ${k}_{y}$ of the BIC, as well as expressions for ridge width $w$ and angle of incidence $\theta $ describing the BIC locations in Fig. 4(a).

The modes of the structure are described by the zeros of the denominator in Eq. (8). Therefore, the BICs are real-$\omega $, real-${k}_{y}$ solutions of one of the following equations: ${\mathcal{D}}_{1}=0$ or ${\mathcal{D}}_{2}=0$. Despite the simple form of these equations, it is difficult to solve them for real $\omega $ and ${k}_{y}$. Instead of doing this, we will use an alternative approach.

Assume that $\omega $ is the real frequency of a certain BIC. Since the modes are described by the zeros of the denominator in Eq. (8), the reflection coefficient in Eq. (8) at frequency $\omega $ tends to infinity. This apparently violates the energy conservation condition. To overcome this contradiction, $\omega $ should also be a zero of the numerator $\mathcal{N}$ in Eq. (8). In this case, the pole of $R$, which is associated with the zero of the denominator, will be compensated by the zero of the numerator. We will use this fact to find the BICs.

Let us first consider the modes defined by the equation ${\mathcal{D}}_{1}=0$. To find the corresponding BICs, we solve the following system with respect to the exponents ${e}^{i\varphi}$ and ${e}^{i\psi}$:

After some simple transformations, we find two different solutions of Eq. (13). The first solution reads asNow, let us obtain explicit expressions for phases $\varphi $ and $\psi $. To do this, we equate the arguments of the left- and right-hand sides of Eq. (14)

Having obtained the expressions for $\varphi $ and $\psi $, we can finally solve Eq. (5) for $\omega $ and ${k}_{y}$

Let us also note that since the BICs are simply real-$\omega $ solutions of the equation ${\mathcal{D}}_{1}=0$ (or ${\mathcal{D}}_{2}=0$), BIC conditions are uniquely defined by reflection coefficients ${r}_{1}$, ${r}_{2}$, and ${r}_{c}$ used in Eqs. (9) and (10). However, to obtain simple closed-form expressions for the BIC condition, we had to use the values of $t$ and ${t}_{c}$ in Eq. (16).

#### C. Comparison with the Rigorous Simulations

To verify the presented theoretical model, we used it to calculate the reflectance spectrum, mode dispersion, and locations of the BICs. The obtained results were compared with RCWA simulation results.

Figure 4(b) shows the reflectance spectrum calculated using Eqs. (8)–(11). This model spectrum is in excellent agreement with the rigorously calculated one shown in Fig. 4(a). Slight discrepancies between the two spectra can be seen only in the vicinity of $w=0$. These discrepancies are due to the near-field interactions in a narrow ridge, which are neglected in the presented theoretical model.

Dispersion of the modes was calculated by numerically solving the equation ${\mathcal{D}}_{1}{\mathcal{D}}_{2}=0$, where ${\mathcal{D}}_{1}$ and ${\mathcal{D}}_{2}$ are given by Eqs. (9) and (10). For the matter of illustration, we show the dispersion not in Fig. 4(b) but in the rigorously calculated Fig. 4(a). The presented analytically calculated dispersion curves are in perfect agreement with the resonances in the rigorously calculated spectrum. Using the model, we also calculated the quality factors of the modes; however, they are visually identical to the rigorously calculated curve in Fig. 5 and thus are not shown.

The white circles depicting the BICs in Fig. 4(a) were calculated using Eqs. (18) and (16). It is evident from the figure that the presented coupled-wave model provides an accurate estimate of the BIC positions. As mentioned in the previous subsection, these white circles were obtained when $m$ and $l$ in Eq. (16) have the same parity. It is interesting to note that when $m$ and $l$ are of different parity, the corresponding $(w,\theta )$ points obtained from Eq. (18) lie on the “resonant” curves $|R|=1$; this can be shown by direct substitution into Eq. (8). These points are shown with black crosses in Fig. 4(a). In contrast to the BICs (infinite-$Q$ resonances shown with white circles), the black crosses describe the relatively low-$Q$ resonances of the structure.

In Section 3, using the EIM approach, we argued that the BICs appear on the quasi-TM dispersion curves of the ridge. The derived theoretical model provides us with deeper insight into the BIC formation mechanism. According to the model, the BICs in the considered structure emerge due to the interaction of the quasi-TE and quasi-TM modes. If we exclude the TE waves propagating inside the ridge from the presented coupled-wave model, the BICs will disappear from the model spectrum.

## 5. CONCLUSION

In the present work, we demonstrated that in the case of diffraction of a TE-polarized mode of a slab waveguide by a dielectric ridge located on the waveguide surface, pronounced resonances in the reflectance and transmittance spectra occur. These resonances take place at the oblique incidence of the mode when only two scattering channels are open, which correspond to the reflected and transmitted modes of the same polarization. The considered structure supports resonances of arbitrarily high-quality factor, as well as modes with an infinite $Q$-factor, i.e., BICs.

We developed a coupled-wave model explaining the optical properties of the structure. By considering the interaction of the TE and TM modes inside the ridge, we derived an accurate expression for the complex reflection coefficient of the structure. Using the proposed model, we analytically confirmed that the considered structure indeed supports robust BICs, which can be obtained by tuning two parameters: ridge width and angle of incidence. We also obtained simple closed-form expressions for these parameters providing the BICs. The developed model is in perfect agreement with the presented full wave simulation results.

The authors believe that the investigated planar structure exhibiting high-$Q$ resonances with Lorentzian line shape can be used as a narrowband integrated optical spatial (angular) or spectral filter. The existence of these high-$Q$ resonances and BICs makes the structure promising also for sensing, enhancing nonlinear light–matter interactions, and analog optical signal processing. In particular, since the transmission coefficient of the ridge strictly vanishes at the resonances, the structure enables performing spatial optical differentiation in transmission as well as spatial optical integration in reflection. This application of the structure will be the subject of a separate work.

## Appendix A. Details of the eigenmode calculation

In this appendix, we present the details of the numerical calculation of the eigenmodes of the structure. In particular, we are interested in the quality factor of the modes, which is defined as $Q=\mathrm{Re}\text{\hspace{0.17em}}\omega /(-2\text{\hspace{0.17em}}\mathrm{Im}\text{\hspace{0.17em}}\omega )$, where $\omega $ is the complex eigenfrequency.

The complex eigenfrequency $\omega $ can be found as the pole of the scattering matrix $\mathbf{S}(\omega ,{k}_{y})$ [41], where ${k}_{y}$ is the wave number of the mode. However, at a fixed angle of incidence $\theta $, the wave number itself depends on angular frequency as ${k}_{y}=(\omega /c){n}_{\mathrm{TE}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta $. Here, ${n}_{\mathrm{TE}}$ is the effective refractive index of the incident mode, which also depends on frequency $\omega $ due to material and structural dispersion. Therefore, to calculate the complex frequency of the mode, we numerically find the complex pole of the matrix $\mathbf{S}[\omega ,(\omega /c)\xb7{n}_{\mathrm{TE}}(\omega )\mathrm{sin}\text{\hspace{0.17em}}\theta ]$, which is considered as a function of $\omega $. Note that the calculated mode will have both a complex frequency and a complex wave number. Therefore, it will decay in both time and space. The presented approach guarantees that the imaginary part of the complex frequency is accurately calculated and describes the line width of the resonance in the frequency spectrum.

It is worth noting that the calculation of the scattering matrix requires defining the basis of incident and scattered waves. This requires special attention in the case of complex frequencies, as discussed in Ref. [42].

Having a method for calculating the complex frequencies, and, hence, the complex wavelengths of the modes, we can numerically solve the equation $\mathrm{Re}\text{\hspace{0.17em}}\lambda =630\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{nm}$ for real $\theta $, which gives us the method for calculating angle of incidence $\theta $, at which the mode is excited exactly at the given wavelength. Moreover, since we calculate the complex $\omega $ as the pole of the scattering matrix, we can also calculate the quality factor. This is the approach we used to obtain Fig. 5.

## Appendix B. Unity amplitudes of the right-hand sides of Eq. (14)

Here, we show that the right-hand sides of Eq. (14) have unity amplitude. To do this, we consider the scattering matrix $\mathbf{S}$ used in Eq. (6). We write the inverse of $\mathbf{S}$ in two ways: using the unitarity property (${\mathbf{S}}^{-1}={\mathbf{S}}^{*}$) and in terms of the adjugate matrix. By equating these two matrices, we obtain

## Funding

Russian Foundation for Basic Research (RFBR) (16-29-11683, 17-47-630323); Ministry of Science and Higher Education of the Russian Federation (State assignment to the FSRC “Crystallography and Photonics” RAS).

## Acknowledgment

Numerical investigation of the optical properties of the structure was supported under Russian Foundation for Basic Research project 16-29-11683, the analysis of the resonances based on the effective index method was supported under Russian Foundation for Basic Research project 17-47-630323, and the derivation of the coupled-wave model and analysis of the BIC formation mechanism were supported by Ministry of Science and Higher Education of the Russian Federation.

## REFERENCES

**1. **H. A. Haus, *Waves and Fields in Optoelectronics* (Prentice-Hall, 1984).

**2. **T. Mossberg, “Planar holographic optical processing devices,” Opt. Lett. **26**, 414–416 (2001). [CrossRef]

**3. **G. Calafiore, A. Koshelev, S. Dhuey, A. Goltsov, P. Sasorov, S. Babin, V. Yankov, S. Cabrini, and C. Peroz, “Holographic planar lightwave circuit for on-chip spectroscopy,” Light Sci. Appl. **3**, e203 (2014). [CrossRef]

**4. **S. Babin, A. Bugrov, S. Cabrini, S. Dhuey, A. Goltsov, I. Ivonin, E.-B. Kley, C. Peroz, H. Schmidt, and V. Yankov, “Digital optical spectrometer-on-chip,” Appl. Phys. Lett. **95**, 041105 (2009). [CrossRef]

**5. **C. Peroz, C. Calo, A. Goltsov, S. Dhuey, A. Koshelev, P. Sasorov, I. Ivonin, S. Babin, S. Cabrini, and V. Yankov, “Multiband wavelength demultiplexer based on digital planar holography for on-chip spectroscopy applications,” Opt. Lett. **37**, 695–697 (2012). [CrossRef]

**6. **C. Peroz, A. Goltsov, S. Dhuey, P. Sasorov, B. Harteneck, I. Ivonin, S. Kopyatev, S. Cabrini, S. Babin, and V. Yankov, “High-resolution spectrometer-on-chip based on digital planar holography,” IEEE Photon. J. **3**, 888–896 (2011). [CrossRef]

**7. **X. Ma, M. Li, and J. J. He, “CMOS-compatible integrated spectrometer based on echelle diffraction grating and MSM photodetector array,” IEEE Photon. J. **5**, 7101307 (2013). [CrossRef]

**8. **L. L. Doskolovich, E. A. Bezus, and D. A. Bykov, “Two-groove narrowband transmission filter integrated into a slab waveguide,” Photon. Res. **6**, 61–65 (2018). [CrossRef]

**9. **L. L. Doskolovich, E. A. Bezus, N. V. Golovastikov, D. A. Bykov, and V. A. Soifer, “Planar two-groove optical differentiator in a slab waveguide,” Opt. Express **25**, 22328–22340 (2017). [CrossRef]

**10. **M. Hammer, A. Hildebrandt, and J. Förstner, “How planar optical waves can be made to climb dielectric steps,” Opt. Lett. **40**, 3711–3714 (2015). [CrossRef]

**11. **M. Hammer, A. Hildebrandt, and J. Förstner, “Full resonant transmission of semiguided planar waves through slab waveguide steps at oblique incidence,” J. Lightwave Technol. **34**, 997–1005 (2016). [CrossRef]

**12. **A. E. Miroshnichenko, S. Flach, and Y. S. Kivshar, “Fano resonances in nanoscale structures,” Rev. Mod. Phys. **82**, 2257–2298 (2010). [CrossRef]

**13. **W. Zhou, D. Zhao, Y.-C. Shuai, H. Yang, S. Chuwongin, A. Chadha, J.-H. Seo, K. X. Wang, V. Liu, Z. Ma, and S. Fan, “Progress in 2D photonic crystal Fano resonance photonics,” Prog. Quantum Electron. **38**, 1–74 (2014). [CrossRef]

**14. **C. W. Hsu, B. Zhen, A. D. Stone, J. D. Joannopoulos, and M. Soljačić, “Bound states in the continuum,” Nat. Rev. Mater. **1**, 16048 (2016). [CrossRef]

**15. **Z. F. Sadrieva and A. A. Bogdanov, “Bound state in the continuum in the one-dimensional photonic crystal slab,” J. Phys. Conf. Ser. **741**, 012122 (2016).

**16. **Z. F. Sadrieva, I. S. Sinev, K. L. Koshelev, A. Samusev, I. V. Iorsh, O. Takayama, R. Malureanu, A. A. Bogdanov, and A. V. Lavrinenko, “Transition from optical bound states in the continuum to leaky resonances: role of substrate and roughness,” ACS Photon. **4**, 723–727 (2017). [CrossRef]

**17. **D. C. Marinica, A. G. Borisov, and S. V. Shabanov, “Bound states in the continuum in photonics,” Phys. Rev. Lett. **100**, 183902 (2008). [CrossRef]

**18. **C. W. Hsu, B. Zhen, J. Lee, S. L. Chua, S. G. Johnson, J. D. Joannopoulos, and M. Soljačić, “Observation of trapped light within the radiation continuum,” Nature **499**, 188–191 (2013). [CrossRef]

**19. **E. N. Bulgakov and A. F. Sadreev, “Bloch bound states in the radiation continuum in a periodic array of dielectric rods,” Phys. Rev. A **90**, 053801 (2014). [CrossRef]

**20. **F. Monticone and A. Alù, “Bound states within the radiation continuum in diffraction gratings and the role of leaky modes,” New J. Phys. **19**, 093011 (2017). [CrossRef]

**21. **C. Blanchard, J.-P. Hugonin, and C. Sauvan, “Fano resonances in photonic crystal slabs near optical bound states in the continuum,” Phys. Rev. B **94**, 155303 (2016). [CrossRef]

**22. **S. P. Shipman and S. Venakides, “Resonant transmission near nonrobust periodic slab modes,” Phys. Rev. E **71**, 026611 (2005). [CrossRef]

**23. **L. Yuan and Y. Y. Lu, “Propagating Bloch modes above the lightline on a periodic array of cylinders,” J. Phys. B **50**, 05LT01 (2017). [CrossRef]

**24. **C. W. Hsu, B. Zhen, S.-L. Chua, S. G. Johnson, J. D. Joannopoulos, and M. Soljačić, “Bloch surface eigenstates within the radiation continuum,” Light Sci. Appl. **2**, e84 (2013). [CrossRef]

**25. **E. N. Bulgakov and A. F. Sadreev, “Bound states in the continuum in photonic waveguides inspired by defects,” Phys. Rev. B **78**, 075105 (2008). [CrossRef]

**26. **I. V. Timofeev, D. N. Maksimov, and A. F. Sadreev, “Optical defect mode with tunable *Q* factor in a one-dimensional anisotropic photonic crystal,” Phys. Rev. B **97**, 024308 (2018). [CrossRef]

**27. **Y. Plotnik, O. Peleg, F. Dreisow, M. Heinrich, S. Nolte, A. Szameit, and M. Segev, “Experimental observation of optical bound states in the continuum,” Phys. Rev. Lett. **107**, 183901 (2011). [CrossRef]

**28. **M. I. Molina, A. E. Miroshnichenko, and Y. S. Kivshar, “Surface bound states in the continuum,” Phys. Rev. Lett. **108**, 070401 (2012). [CrossRef]

**29. **S. Weimann, Y. Xu, R. Keil, A. E. Miroshnichenko, A. Tünnermann, S. Nolte, A. A. Sukhorukov, A. Szameit, and Y. S. Kivshar, “Compact surface Fano states embedded in the continuum of waveguide arrays,” Phys. Rev. Lett. **111**, 240403 (2013). [CrossRef]

**30. **C. L. Zou, J. M. Cui, F. W. Sun, X. Xiong, X. B. Zou, Z. F. Han, and G. C. Guo, “Guiding light through optical bound states in the continuum for ultrahigh-*Q* microresonators,” Laser Photon. Rev. **9**, 114–119 (2015). [CrossRef]

**31. **E. A. Bezus, L. L. Doskolovich, and N. L. Kazanskiy, “Low-scattering surface plasmon refraction with isotropic materials,” Opt. Express **22**, 13547–13554 (2014). [CrossRef]

**32. **G. Lifante, *Integrated Photonics: Fundamentals* (Wiley, 2003).

**33. **R. D. Kekatpure, A. C. Hryciw, E. S. Barnard, and M. L. Brongersma, “Solving dielectric and plasmonic waveguide dispersion relations on a pocket calculator,” Opt. Express **17**, 24112–24129 (2009). [CrossRef]

**34. **M. G. Moharam, E. B. Grann, D. A. Pommet, and T. K. Gaylord, “Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings,” J. Opt. Soc. Am. A **12**, 1068–1076 (1995). [CrossRef]

**35. **L. Li, “Use of Fourier series in the analysis of discontinuous periodic structures,” J. Opt. Soc. Am. A **13**, 1870–1876 (1996). [CrossRef]

**36. **E. Silberstein, P. Lalanne, J.-P. Hugonin, and Q. Cao, “Use of grating theories in integrated optics,” J. Opt. Soc. Am. A **18**, 2865–2875 (2001). [CrossRef]

**37. **J. P. Hugonin and P. Lalanne, “Perfectly matched layers as nonlinear coordinate transforms: a generalized formalization,” J. Opt. Soc. Am. A **22**, 1844–1849 (2005). [CrossRef]

**38. **C. R. Pollock, *Fundamentals of Optoelectronics* (McGraw-Hill, 2003).

**39. **S.-T. Peng and A. A. Oliner, “Guidance and leakage properties of a class of open dielectric waveguides: part I — mathematical formulations,” IEEE Trans. Microw. Theory Tech. **29**, 843–855 (1981). [CrossRef]

**40. **A. A. Oliner, S.-T. Peng, T. I. Hsu, and A. Sanchez, “Guidance and leakage properties of a class of open dielectric waveguides: part II—new physical effects,” IEEE Trans. Microw. Theory Tech. **29**, 855–869 (1981). [CrossRef]

**41. **D. A. Bykov and L. L. Doskolovich, “Numerical methods for calculating poles of the scattering matrix with applications in grating theory,” J. Lightwave Technol. **31**, 793–801 (2013). [CrossRef]

**42. **D. A. Bykov and L. L. Doskolovich, “On the use of the Fourier modal method for calculation of localized eigen modes of integrated optical resonators,” Comput. Opt. **39**, 663–673 (2015). [CrossRef]

**43. **N. A. Gippius and S. G. Tikhodeev, “Application of the scattering matrix method for calculating the optical properties of metamaterials,” Phys. Usp. **52**, 967–971 (2009). [CrossRef]

**44. **V. Karagodsky, C. Chase, and C. J. Chang-Hasnain, “Matrix Fabry–Perot resonance mechanism in high-contrast gratings,” Opt. Lett. **36**, 1704–1706 (2011). [CrossRef]

**45. **R. Orta, A. Tibaldi, and P. Debernardi, “Bimodal resonance phenomena — part II: high/low-contrast grating resonators,” IEEE J. Quantum Electron. **52**, 6600409 (2016). [CrossRef]

**46. **R. Orta, A. Tibaldi, and P. Debernardi, “Bimodal resonance phenomena — part I: generalized Fabry–Pérot interferometers,” IEEE J. Quantum Electron. **52**, 6100508 (2016). [CrossRef]