Photonics Research, 2019, 7 (11): 11001314, Published Online: Nov. 1, 2019   

Integrated flat-top reflection filters operating near bound states in the continuum Download: 507次

Author Affiliations
1 Image Processing Systems Institute—Branch of the Federal Scientific Research Centre “Crystallography and Photonics” of Russian Academy of Sciences, Samara 443001, Russia
2 Samara National Research University, Samara 443086, Russia
Abstract
We propose and theoretically and numerically investigate narrowband integrated filters consisting of identical resonant dielectric ridges on the surface of a single-mode dielectric slab waveguide. The proposed composite structures operate near a bound state in the continuum (BIC) and enable spectral filtering of transverse-electric-polarized guided modes propagating in the waveguide. We demonstrate that by proper choice of the distances between the ridges, flat-top reflectance profiles with steep slopes and virtually no sidelobes can be obtained using just a few ridges. In particular, the structure consisting of two ridges can optically implement the second-order Butterworth filter, whereas at a larger number of ridges, excellent approximations to higher-order Butterworth filters can be achieved. Owing to the BIC supported by the ridges constituting the composite structure, the flat-top reflection band can be made arbitrarily narrow without increasing structure size. In addition to the filtering properties, the investigated structures support another type of BIC—the Fabry–Perot BIC—arising when the distances between adjacent ridges meet the Fabry–Perot resonance condition. In the vicinity of the Fabry–Perot BIC, an effect similar to electromagnetically induced transparency is observed, namely, sharp transmittance peaks against the background of a wide transmittance dip.

1. INTRODUCTION

Spectral filters that selectively reflect or transmit incident light are indispensable in various optical devices, including spectrum analyzers, image sensors, and wavelength demultiplexers. In a wide class of planar (integrated) optoelectronic systems, spectral or spatial filtering is performed in a slab waveguide [17" target="_self" style="display: inline;">7]. Such geometry corresponds to the “insulator-on-insulator” platform and is suitable for the creation of fully integrated optical devices.

In Refs. [1,812" target="_self" style="display: inline;">12], planar Bragg gratings (BGs) and phase-shifted BGs (PSBGs) for spectral filtering of optical radiation propagating in the waveguide were proposed. A planar BG corresponds to a single-mode slab waveguide with a periodically corrugated surface and acts as a reflection filter centered at the Bragg wavelength. A change in the waveguide thickness leads to a change in the effective refractive index of the guided mode, which makes it possible to encode the required refractive index distribution analogous to a free-space BG. A planar PSBG consists of two or more symmetric BGs separated by phase-shift regions and enables obtaining a narrow transmittance peak at the center of the stopband [1,11,12].

Recently, a promising alternative to a narrowband planar BG-based filter was proposed [1317" target="_self" style="display: inline;">17]. It was shown that an extremely simple resonant structure consisting of a single dielectric ridge on the surface of a single-mode slab waveguide supports high-Q resonances associated with the excitation of the TM-polarized modes of the ridge by the obliquely incident TE-polarized mode of the waveguide. Moreover, at properly chosen parameters, the ridge supports bound states in the continuum (BICs).

Photonic BICs are nonradiating eigenmodes supported by structures having open scattering channels [1820" target="_self" style="display: inline;">–20]. The leakage of BIC energy into the open channels can be prevented by means of different mechanisms, including symmetry protection, interaction of several resonators, and interference of several resonances in the same cavity [18]. Various photonic structures were shown to support BICs, including periodic structures (diffraction gratings, photonic crystal slabs, and arrays of rods or spheres) [1927" target="_self" style="display: inline;">27] and waveguide structures of different configurations [1317" target="_self" style="display: inline;">17,28]. In the vicinity of a BIC, the width of the resonance can be made arbitrarily small. Therefore, the photonic structures supporting BICs are not only of theoretical, but also of great practical interest due to their potential applications in lasing, filtering, and sensing [18,26,27].

Near a BIC supported by the ridge structure considered in Refs. [1317" target="_self" style="display: inline;">17], the reflectance spectrum has a symmetric Lorentzian line shape, which makes it possible to use the ridge structure as a narrowband filter operating in reflection. At the same time, the Lorentzian line shape of the reflectance peak does not always fit the design requirements. For many practical applications, achieving a nearly rectangular reflectance spectrum with a flat top and steep slopes is highly desirable [11,29].

In the present work, we investigate composite resonant structures comprising several identical ridges on the surface of a slab waveguide separated by phase-shift regions. We show that by proper choice of the widths of the phase-shift regions, flat-top reflectance profiles with steep slopes can be easily obtained using just a few ridges. Moreover, the width of the reflection band near the BICs can be arbitrarily small by choosing the ridge width. To the best of our knowledge, this work presents the first theoretical and numerical demonstration of integrated flat-top filters operating in the near-BIC regime. In addition to applications in spectral filtering, the proposed composite structures have another interesting property. At the distances between the ridges, which satisfy the condition of the Fabry–Perot resonance, (N1)-degenerate BICs of Fabry–Pérot type arise in the composite structure comprising N ridges. In the vicinity of these BICs, an effect similar to electromagnetically induced transparency (EIT) is observed.

2. GEOMETRY AND REFLECTANCE SPECTRA OF A SINGLE RIDGE ON THE SURFACE OF A SLAB WAVEGUIDE

One of the simplest integrated structures possessing remarkable resonant properties is a dielectric ridge located on the surface of a single-mode slab waveguide [Fig. 1(a)]. In our recent works [13,14], it is shown that in the case of diffraction of an obliquely incident TE-polarized mode of the waveguide, the ridge exhibits BICs and high-Q resonances associated with the excitation of the cross-polarized modes of the ridge.

Fig. 1. Geometry of a ridge on a (a) waveguide layer and of a (b) composite structure consisting of three ridges separated by phase-shift regions. Red arrows indicate the propagation directions of the incident wave I, reflected wave R, and transmitted wave T.

下载图片 查看所有图片

As an example, Fig. 2 shows the TE-polarized mode reflectance versus the ridge width w and angle of incidence θ. The plot was calculated using an efficient in-house implementation of the aperiodic rigorous coupled-wave analysis (aRCWA) technique [3032" target="_self" style="display: inline;">–32]. The simulations were carried out with the following parameters: free-space wavelength λ0=630  nm; refractive index of the waveguide core layer nf=3.3212 (GaP), refractive indices of the superstrate and substrate nu=1 and nl=1.45 (fused silica), respectively; and thickness of the waveguide core layer hf=80  nm and thickness of the waveguide in the ridge region hr=110  nm. At these parameters, the waveguide is single-mode both outside and inside the ridge region. Effective refractive indices of the TE- and TM-polarized modes supported by the waveguide amount to nwg,TE=2.5913 and nwg,TM=1.6327, respectively. In the ridge region (i.e., at the waveguide thickness hr=110  nm), effective refractive indices equal nr,TE=2.8192 and nr,TM=2.1867. Let us note that in the rest of the present work, the chosen free-space wavelength λ0=630  nm is used as the “central” wavelength for all the designed filters.

Fig. 2. Reflectance of an obliquely incident TE-polarized guided mode with free-space wavelength λ0=630  nm from the ridge versus ridge width w and angle of incidence θ. The white circle indicates the BIC position. White dashed lines indicate the ridge widths, at which the angular spectra in Fig. 3(a) are plotted.

下载图片 查看所有图片

The reflectance spectrum shown in Fig. 2 has a pronounced resonant maximum. As it was shown in Ref. [13], this resonance is caused by the excitation of a cross-polarized (“TM-like”) mode of the ridge. It is important to note that in the chosen range of angles of incidence θ[51°,55°], no out-of-plane scattering of the incident mode on the ridge occurs. Moreover, no outgoing (reflected and transmitted) cross-polarized modes propagate. A detailed description of this scattering cancellation mechanism is presented in Refs. [13,14]. In this case, at maxima, the reflectance reaches unity. The angular width of the resonant maximum varies from approximately 0.1° at the boundaries of the considered interval of the ridge width w[300,400]nm to zero at w=344  nm. The point at which the resonance disappears (w=344  nm, θ=53.06°) corresponds to a BIC.

Angular and wavelength reflectance spectra in the vicinity of a BIC have Lorentzian line shapes [13,14]. As an example, Fig. 3 shows the reflectance spectra at w=355  nm, w=360  nm, and w=380  nm. Let us note that these values are not unique and were chosen to illustrate the possibility of obtaining resonant peaks with different widths by changing the ridge width. The angular spectra at the design wavelength 630 nm, which are the cross sections of the reflectance distribution shown in Fig. 2 along the white dashed lines, are shown in Fig. 3(a). The maximum of the reflectance is reached at the angles of incidence θ1=53.28° (w=355  nm), θ2=53.37° (w=360  nm), and θ3=53.72° (w=380  nm); the full widths at half-maximum (FWHMs) of the peaks amount to 0.0045°, 0.009°, and 0.04°, respectively. Figure 3(b) shows the wavelength spectra calculated at the corresponding angles of incidence θi, i=1,2,3. The “spectral” FWHM values of the reflectance peaks centered at λ=λ0=630  nm equal 0.04 nm (w=355  nm), 0.084 nm (w=360  nm), and 0.367 nm (w=380  nm). It is evident from Fig. 3 that the width of the reflectance peak decreases when approaching the BIC. Theoretically, in the vicinity of a BIC, it is possible to obtain a reflectance peak having an arbitrarily small width. In practical implementation, the minimum achievable peak width will be limited by structure finiteness and by material and technological imperfections [17]. Let us mention that the possibility of obtaining near-BIC transmittance peaks using high-contrast dielectric gratings was also recently demonstrated [26,27]. However, the structures considered in Refs. [26,27] are non-integrated and proposed for narrowband filtering of free-space radiation.

Fig. 3. (a) Angular and (b) wavelength ridge reflectance spectra at w=355  nm (solid blue lines), w=360  nm (dashed red lines), and w=380  nm (dotted yellow lines). The angular spectra are calculated at the fixed free-space wavelength of 630 nm and are the cross sections of the reflectance distribution in Fig. 2 along the white dashed lines. The wavelength spectra are calculated at the angles of incidence θ1=53.28° (w=355  nm), θ2=53.37° (w=360  nm), and θ3=53.72° (w=380  nm), at which the resonant peaks are centered at λ=λ0=630  nm.

下载图片 查看所有图片

3. THEORETICAL MODEL DESCRIBING THE SPECTRA OF A COMPOSITE STRUCTURE COMPRISING SEVERAL RIDGES

Lorentzian shape of the reflectance spectra shown in Fig. 3 makes it possible to use the ridge as an angular or as a frequency filter. However, a rectangular reflectance peak with a flat top, steep slopes, and low sidebands would be more suitable for many practical applications [11,29]. In what follows, we demonstrate that composite structures comprising several ridges on the surface of a slab waveguide, which operate in the near-BIC regime, enable obtaining narrow flat-top reflectance peaks with steep slopes. In this case, the BICs supported by the ridges allow one to make the width of the flat-top reflection band arbitrarily small by choosing the ridge width.

Geometry of the proposed composite structure is shown in Fig. 1(b). For illustrative purposes, the structure consisting of three ridges is depicted. We assume that the ridges constituting the composite structure are identical, whereas the widths of the phase-shift regions separating the ridges may, in general, be different.

To analyze the composite structures, it is convenient to describe the optical properties of a single ridge at a fixed angle of incidence θ by the scattering matrix [3234" target="_self" style="display: inline;">–34] where r1(λ) and t1(λ) are the complex reflection and transmission coefficients, respectively, which are considered as functions of the wavelength. For Lorentzian resonances (see Fig. 3), these coefficients can be approximated by the following expressions: where the pole λp is the complex wavelength of the eigenmode supported by the ridge. This form of the reflection and transmission coefficients ensures the unitarity of the scattering matrix of Eq. (1). According to Eq. (2), the reflectance R1(λ)=|r1(λ)|2=1|t1(λ)|2 reaches unity at λ=λ0=Reλp. The width of the resonance is determined by the imaginary part of the pole λp. Indeed, from Eq. (2), it is easy to obtain that the FWHM of the resonant reflectance peak R(λ) (and of the resonant transmittance dip T(λ)=|t1(λ)|2) amounts to Δ=2|Imλp|. Note that Eq. (2) remains valid in the vicinity of a BIC. If the BIC condition is strictly satisfied (i.e., Imλp=0), the resonance disappears (see Fig. 2), the reflection coefficient r1(λ) vanishes, and the magnitude of the transmission coefficient becomes unity: |t1(λ)|=1.

3.4 A. Two-Ridge Composite Structures

Let us now consider a composite structure consisting of two ridges described by identical scattering matrices defined by Eq. (1) and separated by a phase-shift region with the width l1. At a fixed angle of incidence θ0, the scattering matrix of the composite structure can be expressed through the matrix S1(λ) as [32,35] S2(λ)=S1(λ)*L(l1)*S1(λ),where the symbol * denotes the Redheffer star product [32], and L(l1) is the scattering matrix describing the phase-shift region: L(l1)=exp[iψ(l1)]I.Here, ψ(l1) is the phase shift acquired by the mode upon propagation through the phase-shift region, and I is the 2×2 identity matrix. For simplicity, in Eq. (4), we assume that the dependence of the phase shift on the wavelength can be neglected, so that this phase shift can be approximated as ψ(l1)=kxl1=k0nwg,TEcosθ0·l1,where kx=k0nwg,TEcosθ0 is the x component of the wave vector of the incident TE-polarized guided mode, k0=2π/λ0 is the wavenumber, and nwg,TE is the effective refractive index of the incident TE-polarized mode. Note that the made assumption is valid when we consider a narrow wavelength range around the resonant wavelength λ0=Reλp.

Using the definition of the Redheffer star product [32], one can easily obtain the scattering matrix of the composite structure in the following form: S2(λ)=(t2(λ)r2(λ)r2(λ)t2(λ))=11exp(2iψ)r12(exp(iψ)t12r1[1exp(2iψ)(r12t12)]r1[1exp(2iψ)(r12t12)]exp(iψ)t12).By substituting Eq. (2) into Eq. (6), we obtain the following resonant approximations of the reflection and transmission coefficients of the composite structure: r2(λ)=γr,2λλz,1(λλp,1)(λλp,2),t2(λ)=γt,2(λReλp)2(λλp,1)(λλp,2),where λz,1=Reλp+tan[ψ(l1)+φ]Imλp,λp,1,2=Reλp+i(1±σ)Imλp,where γt,2=exp[iψ(l1)+2iφ], γr,2=2iγt,2cos[ψ(l1)+φ]Imλp, and σ=exp[iψ(l1)+iφ]. According to Eqs. (7)–(9), the composite structure consisting of two ridges supports two eigenmodes. These eigenmodes have the complex wavelengths defined by Eq. (9), which are the poles of the reflection and transmission coefficients. Also, the reflection coefficient has a real-valued zero λ=λz,1 defined by Eq. (8), whereas the transmission coefficient has a real-valued second-order zero at λ=Reλp.

It is interesting to discuss the behavior of Eq. (7) in the case when the two ridges constituting the composite structure support a BIC (λpR). In this case, the imaginary part Imλp in Eqs. (8) and (9) vanishes, and the following equalities take place: λp,1=λp,2=λz,1=Reλp=λp. Therefore, taking into account the expression for γr,2, the two poles in both fractions in Eq. (7) are canceled out by the zeros. Therefore, the considered composite structure supports two BICs having the same wavelength λp, or, in other words, a doubly degenerate BIC.

It is evident from Eq. (7) that in the general case (Imλp>0), the expressions for the reflection and transmission coefficients r2(λ) and t2(λ) have a significantly more complex form than the corresponding coefficients r1(λ) and t1(λ) of a single ridge. This gives more flexibility for the control of the resonant peak shape. In particular, under the condition ψ(l1)+φ=π(m1/2),mN,where the phase shift ψ(l1) is defined by Eq. (5) and φ is the phase of the reflection coefficient r1(λ0) [Eq. (2)], the coefficients r2(λ) and t2(λ) take the form r2(λ)=2exp(iφ)1[(λλp)/Imλp]2,t2(λ)=(1)miexp(iφ)(λReλp)/Imλp1[(λλp)/Imλp]2.It is important to note that the expression for r2(λ) in Eq. (11) coincides with the transfer function of the second-order Butterworth filter and, compared to r1(λ), provides a much more rectangular shape of the reflectance peak.

3.5 B. Multiple-Ridge Composite Structures

Similarly, one can obtain the reflection and transmission coefficients of the composite structure consisting of three equally spaced ridges. Indeed, by calculating S3(λ)=S2(λ)*L(l1)*S1(λ), we obtain the reflection and transmission coefficients of the composite structure in the form r3(λ)=γr,3(λλz,1)(λλz,2)(λλp,1)(λλp,2)(λλp,3),t3(λ)=γt,3(λReλp)3(λλp,1)(λλp,2)(λλp,3),where λz,1,2=Reλp+i1σ21±σ+σ2Imλp,λp,1=Reλp+i(1σ2)Imλp,λp,2,3=λp+iσ2(σ±8+σ2)Imλp,where σ=exp[iψ(l1)+iφ], γr,3=iexp(iφ){1+2σ3cos[ψ(l1)+φ]}Imλp, and γt,3=exp[2iψ(l1)+3iφ]. According to Eqs. (12)–(14), the reflection and transmission coefficients of the composite structure comprising three ridges have three poles [three eigenmodes with the complex wavelengths defined by Eq. (14)]. The reflection coefficient has two complex-valued zeros defined by Eq. (13), whereas the transmission coefficient has a third-order zero λ=Reλp.

In the general case of a composite structure consisting of N ridges, the reflection and transmission coefficients of such a composite structure have N poles (N eigenmodes with the complex wavelengths λp,m) and can be described by the following expressions: rN(λ)=γr,Nm=1N1(λλz,m)m=1N(λλp,m),tN(λ)=γt,N(λReλp)Nm=1N(λλp,m).Let us note that the reflection coefficient rN(λ) has N1 complex-valued zeros λz,m, whereas the transmission coefficient tN(λ) has a real-valued zero λ=Reλp of the order N. If the ridges constituting the composite structure support a BIC, λp is real, and, as one can show similarly to the two-ridge structure discussed above, λp,m=λz,m=λp, therefore, all the poles in Eq. (15) get canceled out by the zeros. This means that N independent BICs having the same frequency coexist in the structure comprising N ridges. Such BICs can be referred to as N-degenerate. Detuning from the BIC condition lifts the degeneracy, which results in N distinct poles (resonances) and N1 complex-valued zeros λz,m of the reflection coefficient. These zeros and poles of the reflection and transmission coefficients are functions of the widths li, i=1,,N of the phase-shift regions separating the ridges. In the following section, we will show that optimization with respect to these parameters enables obtaining a nearly rectangular reflectance peak with a flat top, steep slopes, and no sidelobes.

Let us also emphasize that the theoretical results obtained in this section are not specific for the considered ridge structures and can be applied to any composite structure comprising Lorentzian-line-shaped resonators separated by phase-shift layers.

4. NUMERICAL INVESTIGATION AND OPTIMIZATION OF COMPOSITE STRUCTURES

In Section 3, it is shown that if the condition of Eq. (10) is met, the complex reflection coefficient of the composite structure consisting of two ridges coincides with the transfer function of the second-order Butterworth filter, which provides a more rectangular shape of the resonance than that of a single ridge. Unfortunately, this effect is not preserved at a larger number of ridges N>2. For example, Fig. 4 shows the rigorously calculated transmittance and reflectance spectra of the composite structures consisting of N=2, N=4, and N=6 ridges. The spectra were calculated at the angle of incidence θ=53.72°, ridge width w=380  nm, and a fixed distance between the ridges l=948  nm satisfying the condition of Eq. (10) at m=3. For comparison, the dashed lines in Fig. 4 show the spectra of a single ridge. At N=2 [Fig. 4(a)], the spectra are with high accuracy described by the expressions R2(λ)=|r2(λ)|2 and T2(λ)=1|r2(λ)|2, where r2(λ) is the transfer function of the second-order Butterworth filter defined by Eq. (11). With an increase in N, the central parts (reflectance peak and transmittance dip) of the spectra become closer to a rectangle, but the sidelobes (two at N=4 and four at N=6) arise. These sidelobes emerge due to the appearance of additional zeros of the reflection and transmission coefficients.

Fig. 4. Transmittance (black) and reflectance (red) spectra of composite structures consisting of (a) N=2, (b) N=4, and (c) N=6 ridges at the ridge width w=380  nm and the distance between the ridges l=948  nm (solid lines). Dashed lines show the spectra of a single ridge.

下载图片 查看所有图片

The shape of the spectra can be altered by changing the widths lj of the phase-shift regions separating the ridges. Indeed, in the general case, the zeros and poles in Eq. (15) describing the reflection and transmission coefficients of a composite structure comprising N ridges are functions of the variables lj,j=1,,N1. Considering lj as optimization parameters, one can try to obtain the reflectance and transmittance spectra with a required shape. Let us define the desired shape of the reflectance spectrum RN(ω)=|rN(ω)|2 by the function which describes a peak with a smoothed rectangular shape having the width 2σ and no sidelobes. Note that Eq. (16) generalizes Eq. (11) and corresponds to the squared modulus of the transfer function of the Butterworth filter of the order N [36]. The function defined by Eq. (16) can be referred to as Butterworth line shape, which becomes the conventional Lorentzian line shape at N=1.

Let us demonstrate the possibility of tailoring the shape of the resonant reflectance peak (transmittance dip) of the composite structures at N=4 and N=6. The distances lj were optimized using the conjugate gradient method in order to obtain reflectance spectra described by the functions RBW,4(λ) (at N=4) and RBW,6(λ) (at N=6) at 2σ=1  nm. The used σ value corresponds to the reflectance peak width of a single ridge at the 0.1 level (Fig. 3). As a result of the optimization, the following widths of the phase-shift regions were found: l1=l3=1278  nm and l2=1395  nm at N=4 and l1=l5=1486  nm, l2=l4=1383  nm, and l3=1485  nm at N=6. The reflectance and transmittance spectra of the composite structures calculated at the obtained values of lj are shown in Fig. 5. The calculated spectra are very close to the functions TBW,N(λ) and RBW,N(λ) at both N=4 [Fig. 5(a)] and N=6 [Fig. 5(b)]. Let us note that by further increasing the number of ridges constituting the composite structure, one can obtain steeper reflectance spectra corresponding to the Butterworth line shapes at greater N values. However, here we restrict our consideration to structures consisting of no more than six ridges, since this value provides a good peak quality [Fig. 5(b)] while maintaining a reasonable lateral footprint of the structure (approximately 9.5 μm).

Fig. 5. Transmittance (black) and reflectance (red) spectra of the optimized composite structures consisting of (a) N=4 and (b) N=6 ridges at w=380  nm (solid lines). Dashed lines show the spectra of a single ridge.

下载图片 查看所有图片

Let us now consider an important advantage of the proposed composite filters operating in the near-BIC regime, which consists of the possibility of generating a near-rectangular reflectance peak with an essentially subnanometer width. The width of the resonance of a single ridge decreases when approaching a BIC and finally vanishes (Fig. 3). Therefore, in the vicinity of a BIC, it is theoretically possible to obtain a Lorentzian reflectance peak with an arbitrarily small width. Let us recall that the width of the flat-top peak of the composite structure is very close to the width of the initial Lorentzian peak (at the 0.1 level). That means that by tuning the width of the ridges constituting the composite structure, one can engineer a very narrow peak having a flat-top (rectangular) shape described by Eq. (16). Let us demonstrate this possibility. In the example considered above (θ=θ3=53.72°, w=380  nm), the width of the initial Lorentzian-shaped reflectance peak at the 0.1 level amounts to Δ1  nm (Figs. 3 and 5). At the angle of incidence θ=θ2=53.37° and the ridge width w=360  nm, the width of the peak (Fig. 3) decreases by approximately five times to Δ0.2  nm. Figure 6 shows the spectra of the composite structures constructed using this “high-Q” ridge with the width w=360  nm and consisting of N=2, N=4, and N=6 ridges. In the structure with N=2, the width of the phase-shift region between the ridges l=948  nm satisfies the condition of Eq. (10). In this case, the reflectance spectrum is very close to the squared modulus of the second-order Butterworth filter. In the cases of N=4 and N=6, the widths lj of the phase-shift regions were optimized in order to obtain the spectra with the Butterworth line shapes described by the functions RBW,4(λ) and RBW,6(λ), respectively. The following lj values were obtained: l1=l3=1285  nm and l2=1410  nm at N=4 and l1=l5=1286  nm, l2=l4=1405  nm, and l3=1282  nm at N=6. Figure 6 shows that the obtained spectra indeed contain smoothed rectangular reflectance peaks with the width of about 0.2 nm, i.e., approximately five times narrower than in Fig. 5 (note the different wavelength ranges shown in Figs. 5 and 6). By further approaching the BIC located at w=344  nm, θ=53.06°, one can further decrease the width of the resonant peak provided by the resonant structure. It is important to note that, in contrast to many conventional structures providing narrow reflectance or transmittance peaks [1,7,8,11], a decrease in the peak width does not lead to an increase in the structure size.

Fig. 6. Transmittance (black) and reflectance (red) spectra of the optimized composite structures consisting of (a) N=2, (b) N=4, and (c) N=6 ridges at w=360  nm (solid lines). Dashed lines show the spectra of a single ridge.

下载图片 查看所有图片

5. FABRY–PEROT BOUND STATES IN THE CONTINUUM IN THE COMPOSITE STRUCTURE

In addition to the possibility of obtaining a nearly rectangular reflectance peak with a flat top, steep slopes, and virtually no sidebands, the studied composite structures possess one more remarkable property. Before discussing it, let us recall that the single ridge on the surface of the waveguide used as the building block of the proposed composite structures supports BIC, which is caused by the interaction of TE- and TM-polarized modes in the ridge region [13]. Surprisingly, even when the width of the ridges constituting the considered composite filters is detuned from the BIC condition, the composite structure can by itself support BICs, albeit of a different type. Indeed, let us show that if the distances between the ridges satisfy the Fabry–Perot resonance condition the so-called Fabry–Perot BICs are formed in the composite structure [18,19,37]. Let us recall that in Eq. (17), ψ(li) is the phase shift acquired by the mode upon propagation through the phase-shift region and defined by Eq. (5), and φ is the phase of the reflection coefficient r1(λ0) in Eq. (2).

First, let us consider the structure comprising two ridges. The reflection and transmission coefficients of this structure are described by Eqs. (7)–(9). If the condition of Eq. (17) is met, one of the poles in Eq. (9) becomes real and coincides with the real-valued transmittance zero. To be specific, let Imλp,1=0. In this case, the integer m in Eq. (17) is even, and λp,2=Reλp+2iImλp. The eigenmode of the composite structure with the wavelength λ=Reλp is a BIC [18,19,37]. The considered structure consisting of two ridges has only two scattering channels corresponding to the reflected and transmitted TE-polarized modes. In the case of a BIC, the leakage to these channels is prevented due to the Fabry–Perot resonance formed between the ridges under the condition of Eq. (17).

Next, let us consider a composite structure comprising three ridges. The spectra of this structure are described by Eqs. (12)–(14). Under the condition of Eq. (17), the poles defined by Eq. (14) (the complex wavelengths of the eigenmodes of the structure) take the form Therefore, if the condition of Eq. (17) is met, the poles λp,1 and λp,2 become real and coincide with the real-valued zeros of the reflection and transmission coefficients. This means that two equal-wavelength BICs (or, in other words, a double-degenerate BIC) are formed in the composite structure containing three ridges.

In a composite structure consisting of N ridges, the reflection and transmission coefficients defined by Eq. (15) have N poles. By induction, it can be shown that, under the condition of Eq. (17), N1 poles become real, coinciding with the real-valued zeros of the reflection and transmission coefficients. In this case, (N1)-degenerate BICs are formed in the composite structure. One can show that the only remaining pole will have an imaginary part N times greater than the imaginary part of the pole λp of the reflection and transmission coefficients of the single ridge: λp,N=Reλp+NiImλp. Let us note that the field of the N1 degenerate BICs is localized between each pair of adjacent ridges, as opposed to the N degenerate BICs discussed in Section 3, the field of which is localized inside each ridge.

Let us note that the formation of the Fabry–Perot BICs in a composite structure consisting of two resonators is well known [18,19,37]. At the same time, to the best of our knowledge, the formation of (N1)-degenerate BICs in a composite structure consisting of N resonant ridges (or other resonators) has not been previously studied.

In order to confirm the formation of the Fabry–Perot BICs, we calculated the reflectance RN(l,λ)=|rN(l,λ)|2 and the transmittance TN(l,λ)=|tN(l,λ)|2 of several composite structures versus the wavelength and the distance l between the ridges. The calculation was carried out at the fixed angle of incidence θ=53.72° and the ridge width w=380  nm. Figure 7 shows the reflectance and transmittance of the composite structures consisting of two, three, and four ridges. Vertical dashed lines show the width of the phase-shift regions l=lFP=970.2  nm, which satisfies the Fabry–Perot resonance condition of Eq. (17) at m=4. Horizontal dashed lines show the resonant wavelength of the single ridge λ0=630  nm. In the vicinity of the intersection of these lines, resonant features (reflectance minima and transmittance maxima) are clearly visible. When approaching the intersection, the width of the resonant features decreases (the quality factor of the resonances increases). It is important to note that N1 sharp transmittance peaks occur against the background of a relatively wide transmittance dip caused by the multiple zeros of the transmission coefficient tN(λ) defined by Eq. (15). This effect is similar to the EIT effect and, in the considered case, is associated with the presence of N1 complex-valued zeros in the reflection coefficient rN(λ). At the vertical dashed lines, the resonant features vanish, which confirms the formation of the BICs. Let us note that according to the presented theoretical description, the number of the vanishing resonant features (and the number of the emerging EIT peaks) is one less than the number of the ridges N, which constitute the composite structure.

Fig. 7. Reflectance RN(l,λ)=|rN(l,λ)|2 (top) and transmittance TN(l,λ)=|tN(l,λ)|2 (bottom) of the composite structures consisting of (a), (d) N=2, (b), (e) N=3, and (c), (f) N=4 ridges versus the distance between the ridges l and free-space wavelength. Vertical dashed lines show the distance lFP=970.2  nm corresponding to the Fabry–Perot resonance.

下载图片 查看所有图片

As an additional illustration of the BIC formation and the EIT effect, let us consider the transmittance spectra of a composite structure consisting of three ridges (Fig. 8). The spectra were calculated at the phase-shift region width l=lFP=970.2  nm corresponding to the BIC condition of Eq. (17), and at the “offset” widths l=lFP+5  nm=975.2  nm and l=lFP+10  nm=980.2  nm. At l=lFP, the transmittance spectrum has a smooth dip centered at λ=λ0=630  nm. This agrees with Eq. (18), according to which the transmission coefficient has a single pole λp,3=Reλp+3iImλp. If we violate the BIC condition by increasing the distance between the ridges by, e.g., δ=5  nm, two sharp peaks occur in the transmittance spectra at the wavelengths greater than the central wavelength λ0 (see the red curve in Fig. 8). These peaks correspond to the EIT effect and are associated with the zeros of the reflection coefficient defined by Eq. (13). At δ=10  nm, the EIT peaks become wider and experience a further red shift (see the yellow curve in Fig. 8). It is interesting to note that the different widths of the EIT peaks follow from the expression for the poles of the transmission coefficient [Eq. (14)]. Indeed, the widths of the EIT peaks are determined by the magnitudes of the imaginary parts of the poles λp,i,i=1,2. From Eq. (14), it is easy to obtain that near the BIC [in the vicinity of the point l=lFP satisfying the condition of Eq. (17)], the imaginary parts of the poles λp,i are quadratic with respect to the “offset” magnitude δ: where α1=2 and α2=2/27. Since α1/α2=27, the EIT peak corresponding to the pole λp,2 should have an order-of-magnitude smaller width, which agrees with the spectra shown in Figs. 7(e) and 8.

Fig. 8. Transmittance spectra of the composite structure consisting of N=3 ridges calculated at the width of the phase-shift region between the ridges l=lFP=970.2  nm corresponding to the BIC condition (solid blue line), and at the widths l=lFP+5  nm=975.2  nm (dashed red line) and l=lFP+10  nm=980.2  nm (dotted yellow line).

下载图片 查看所有图片

6. CONCLUSION

In the present work, we investigated resonant optical properties of composite structures consisting of N identical ridges on the surface of a slab waveguide separated by phase-shift regions. The single ridge used as the building block of the proposed composite structures supports BICs and high-Q resonances having a Lorentzian line shape. Using the scattering matrix formalism, we obtained resonant approximations of the reflection and transmission coefficients of the composite structure. According to the derived expressions, in the case when the ridges constituting the composite structure support BICs, the composite structure supports N-degenerate BICs, whereas the detuning from the BIC condition allows one to obtain flat-top reflection filters. In particular, we demonstrated that at N=2, the reflection coefficient coincides with the transfer function of the second-order Butterworth filter. It was shown that by optimizing the widths of the phase-shift regions at N>2, one can obtain a flat-top reflectance peak with steep slopes using just a few ridges. In our opinion, the proposed composite structure has a simpler geometry and is more compact than the planar filters based on integrated BGs.

To the best of our knowledge, this work presents the first theoretical and numerical demonstration of integrated flat-top filters operating in the near-BIC regime. An important advantage of the proposed structure consists of the fact that the flat-top reflectance peak can be made arbitrarily narrow by choosing the width of the ridges in the vicinity of a BIC supported by a single ridge. In this case, the corresponding composite structure can be considered as a perturbation of the composite structure supporting an N-degenerate BIC. In contrast to many conventional structures providing narrow reflectance or transmittance peaks, a decrease in the peak width does not lead to an increase in the structure size.

In addition to the spectral filtering applications, the considered composite structures possess one more remarkable optical property, namely, the formation of BICs having a different origin than that of the BICs supported by the initial ridge. At the distances between the adjacent ridges satisfying the Fabry–Perot condition, (N1)-degenerate Fabry–Perot BICs arise in the composite structure comprising N ridges. In the vicinity of the formed BICs, an effect very similar to the EIT effect is observed. In the present case, the EIT peaks can be explained by the presence of N1 zeros of the reflection coefficient.

The obtained results may find application in the design of optical filters, sensors, and the devices for transformation of optical signals. In particular, the proposed composite structures can be used for higher-order differentiation of optical signals propagating in the slab waveguide. In addition, the theoretical results presented in Section 3 and concerning the description of the resonant properties of the composite structures are not specific for the considered ridge structures and can be extended to other composite structures comprising Lorentzian-line-shaped resonators separated by phase-shift layers.

7 Acknowledgment

Acknowledgment. This work was funded by the Russian Foundation for Basic Research (design and investigation of the composite filters, Sections 3 and 4), by the Russian Federation Ministry of Science and Higher Education (simulation of the single-ridge structure, Section 2), and by the Russian Science Foundation (investigation of bound states in the continuum and electromagnetically induced transparency in the composite structures, Section 5).

References

[1] HausH. A., Waves and Fields in Optoelectronics (Prentice-Hall, 1984).

[2] T. Mossberg. Planar holographic optical processing devices. Opt. Lett., 2001, 26: 414-416.

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

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

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

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

[7] L. L. Doskolovich, E. A. Bezus, D. A. Bykov. Two-groove narrowband transmission filter integrated into a slab waveguide. Photon. Res., 2018, 6: 61-65.

[8] R. V. Schmidt, D. C. Flanders, C. V. Shank, R. D. Standley. Narrow-band grating filters for thin-film optical waveguides. Appl. Phys. Lett., 1974, 25: 651-652.

[9] C. S. Hong, J. B. Shellan, A. C. Livanos, A. Yariv, A. Katzir. Broad-band grating filters for thin film optical waveguides. Appl. Phys. Lett., 1977, 31: 276-278.

[10] L. A. Weller-Brophy, D. G. Hall. Analysis of waveguide gratings: application of Rouard’s method. J. Opt. Soc. Am. A, 1985, 2: 863-871.

[11] R. Zengerle, O. Leminger. Phase-shifted Bragg grating filters with improved transmission characteristics. J. Lightwave Technol., 1995, 13: 2354-2358.

[12] J. N. Damask, H. A. Haus. Wavelength-division multiplexing using channel-dropping filters. J. Lightwave Technol., 1993, 11: 424-428.

[13] E. A. Bezus, D. A. Bykov, L. L. Doskolovich. Bound states in the continuum and high-Q resonances supported by a dielectric ridge on a slab waveguide. Photon. Res., 2018, 6: 1084-1093.

[14] E. A. Bezus, L. L. Doskolovich, D. A. Bykov, V. A. Soifer. Spatial integration and differentiation of optical beams in a slab waveguide by a dielectric ridge supporting high-Q resonances. Opt. Express, 2018, 26: 25156-25165.

[15] C.-L. Zou, J.-M. Cui, F.-W. Sun, X. Xiong, X.-B. Zou, Z.-F. Han, G.-C. Guo. Guiding light through optical bound states in the continuum for ultrahigh-Q microresonators. Laser Photon. Rev., 2015, 9: 114-119.

[16] A. P. Hope, T. G. Nguyen, A. Mitchell, W. Bogaerts. Quantitative analysis of TM lateral leakage in foundry fabricated silicon rib waveguides. IEEE Photon. Technol. Lett., 2016, 28: 493-496.

[17] NguyenT. G.RenG.SchoenhardtS.KnoerzerM.BoesA.MitchellA., “Ridge resonance: a new resonance phenomenon for silicon photonics harnessing bound states in the continuum,” arXiv:1810.06734 (2018).

[18] C. W. Hsu, B. Zhen, A. D. Stone, J. D. Joannopoulos, M. Soljačić. Bound states in the continuum. Nat. Rev. Mater., 2016, 1: 16048.

[19] D. C. Marinica, A. G. Borisov, S. V. Shabanov. Bound states in the continuum in photonics. Phys. Rev. Lett., 2008, 100: 183902.

[20] C. W. Hsu, B. Zhen, J. Lee, S. L. Chua, S. G. Johnson, J. D. Joannopoulos, M. Soljačić. Observation of trapped light within the radiation continuum. Nature, 2013, 499: 188-191.

[21] J. W. Yoon, S. H. Song, R. Magnusson. Critical field enhancement of asymptotic optical bound states in the continuum. Sci. Rep., 2015, 5: 18301.

[22] Y. Wang, J. Song, L. Dong, M. Lu. Optical bound states in slotted high-contrast gratings. J. Opt. Soc. Am. B, 2016, 33: 2472-2479.

[23] C. Blanchard, J.-P. Hugonin, C. Sauvan. Fano resonances in photonic crystal slabs near optical bound states in the continuum. Phys. Rev. B, 2016, 94: 155303.

[24] D. A. Bykov, E. A. Bezus, L. L. Doskolovich. Coupled-wave formalism for bound states in the continuum in guided-mode resonant gratings. Phys. Rev. A, 2019, 99: 063805.

[25] E. N. Bulgakov, A. F. Sadreev. Bloch bound states in the radiation continuum in a periodic array of dielectric rods. Phys. Rev. A, 2014, 90: 053801.

[26] J. M. Foley, S. M. Young, J. D. Phillips. Symmetry-protected mode coupling near normal incidence for narrow-band transmission filtering in a dielectric grating. Phys. Rev. B, 2014, 89: 165111.

[27] J. M. Foley, J. D. Phillips. Normal incidence narrowband transmission filtering capabilities using symmetry-protected modes of a subwavelength, dielectric grating. Opt. Lett., 2015, 40: 2637-2640.

[28] E. N. Bulgakov, A. F. Sadreev. Bound states in the continuum in photonic waveguides inspired by defects. Phys. Rev. B, 2008, 78: 075105.

[29] Y. H. Ko, R. Magnusson. Flat-top bandpass filters enabled by cascaded resonant gratings. Opt. Lett., 2016, 41: 4704-4707.

[30] E. Silberstein, P. Lalanne, J.-P. Hugonin, Q. Cao. Use of grating theories in integrated optics. J. Opt. Soc. Am. A, 2001, 18: 2865-2875.

[31] J. P. Hugonin, P. Lalanne. Perfectly matched layers as nonlinear coordinate transforms: a generalized formalization. J. Opt. Soc. Am. A, 2005, 22: 1844-1849.

[32] L. Li. Formulation and comparison of two recursive matrix algorithms for modeling layered diffraction gratings. J. Opt. Soc. Am. A, 1996, 13: 1024-1035.

[33] N. A. Gippius, S. G. Tikhodeev, T. Ishihara. Optical properties of photonic crystal slabs with an asymmetrical unit cell. Phys. Rev. B, 2005, 72: 045138.

[34] D. A. Bykov, L. L. Doskolovich, N. V. Golovastikov, V. A. Soifer. Time-domain differentiation of optical pulses in reflection and in transmission using the same resonant grating. J. Opt., 2013, 15: 105703.

[35] N. V. Golovastikov, D. A. Bykov, L. L. Doskolovich. Temporal differentiation and integration of 3D optical pulses using phase-shifted Bragg gratings. Comput. Opt., 2017, 41: 13-21.

[36] H.-C. Liu, A. Yariv. Synthesis of high-order bandpass filters based on coupled-resonator optical waveguides (CROWs). Opt. Express, 2011, 19: 17653-17668.

[37] V. Liu, M. Povinelli, S. Fan. Resonance-enhanced optical forces between coupled photonic crystal slabs. Opt. Express, 2009, 17: 21897-21909.

Leonid L. Doskolovich, Evgeni A. Bezus, Dmitry A. Bykov. Integrated flat-top reflection filters operating near bound states in the continuum[J]. Photonics Research, 2019, 7(11): 11001314.

本文已被 1 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!