Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Bi-layered composite gratings with high diffraction efficiency enabled by near-field coupling

Open Access Open Access

Abstract

In this paper, we present a design method for bi-layered composite gratings to achieve high diffraction efficiency. These composite gratings feature strong near-field coupling between their constituent dielectric subwavelength gratings, thus enabling high-efficiency first-order diffraction in the far-field. An intuitive explanation based on a wavevector matching condition for such high diffraction efficiency composite gratings is provided. According to theoretical analysis, a design strategy for the proposed composite gratings is developed and verified by numerical simulations with gratings working in both TE and TM modes. The proposed strategy could open door to develop bi-layered composite gratings for manipulating diffracted waves with high efficiency, thus may potentially enable new applications in photonic systems.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Flat optics [15] paves a new way to achieve high-efficiency wave manipulation within subwavelength thicknesses that are distinctly different from their conventional bulky counterparts [68]. Ease of fabrication, scalability of design, and convenience of integration make them a promising candidate, especially for long-wave applications. Among them, the subwavelength gratings (SWGs) have been extensively studied for a wide range of optoelectronic device applications, including lasers [9,10], filters [11,12], sensors [1315], and beam splitters [16,17], etc. SWGs are composed of subwavelength periodic structures that enable highly efficient wave manipulation. Through engineering the SWGs’ parameters, such as the duty cycle, grating depth, and tapered sidewalls, extraordinary characteristics including high-efficiency wideband reflection or transmission [18], narrowband filtering [19], and ultrathin photodetector [20] can be obtained.

In order to achieve high-efficiency wave manipulation, several working mechanisms utilizing guided-mode resonance (GMR), Fabry-Pérot (F-P) resonance, and metagratings have been adopted in SWGs. SWG featuring GMR is composed of a dielectric grating deposited on a thin film waveguide and a substrate [21,22]. By specifically designing the grating so that the diffraction orders of the grating can be coupled with the propagating waveguide modes. Some of these guided modes are diffracted out of the waveguide and coupled back to radiation, and subsequently interfere with the uncoupled waves, leading to high-efficiency reflection or transmission. SWG featuring F-P resonance is composed of a single grating or bi-layered gratings to form an F-P cavity [2326]. Due to the F-P interference between the waveguide array modes (WAMs) within the gratings or between the bi-layered gratings, high-efficiency reflection or transmission can occur. Metagrating [27,28] is a recently proposed concept to achieve high-efficiency wave manipulation. The electromagnetic response of each unit cell is specifically tailored to suppress undesired diffraction orders and reroute the incidence toward the desired diffraction order with high efficiency. Based on the concepts of GMR, F-P resonance, and metagratings, many dielectric optical components including wavelength-selective mirrors [29,30], bandpass filters [3133], and polarizers [34,35] have been developed.

SWG featuring near-field Moiré effect by superimposing two one-dimensional SWGs with different periods has also been reported [36]. By inserting a silver slab between the two SWGs, the evanescent field can be enhanced via surface plasmon polariton (SPP) excitation. However, this wave manipulation scheme suffers large Ohmic losses within near-infrared and visible wavelength ranges due to the metallic material used for near-field enhancement. This could limit the device’s efficiency and hinder its high optical power applications due to thermal heating. Therefore, all-dielectric bi-layered SWGs incorporating near-field coupling that can support distinct optical resonances and high-efficiency wave manipulation are much more desired.

In this paper, a new strategy is proposed for designing the bi-layered composite grating with high first-order diffraction efficiency (∼90%). The bi-layered composite grating is formed by superimposing two rectangular dielectric SWGs with different grating periods, which are placed close to their near-field region to incorporate the near-field coupling effect. Oblique incidence is coupled with WAM of the bottom grating via wavevector-matching enabled by the top grating. Then, the resonant WAM is out-coupled to the first-order far-field diffraction of the composite grating, resulting in high-efficiency diffraction in the reflection region. The coupling might be similar to that of Wood’s anomalies on gratings [37]. Based on this concept, an effective design method is developed to obtain bi-layered composite gratings with high-efficiency first-order diffraction that fits over a broad spectral range. Simulation results of the bi-layered composite grating designs operated in both TE and TM modes are presented to verify the proposed design method. Moreover, parametric dependencies of these designs are also investigated. The bi-layered composite gratings could provide tunable bandwidth and adjustable wavelength functionalities by controlling the gap spacing of the SWGs and the incident angle of the input light, thus showing great application potential in sensing and communication. We envision that the design method proposed here could advance the tunable diffractive optical elements (DOEs) devices and promote their applications.

2. Configuration of the bi-layered composite grating

A cross-section of a unit cell of the proposed bi-layered composite grating is schematically shown in Fig. 1(a). The bi-layered composite grating locates in the x-y plane and the grating bars lie along the y axis. Each unit cell consists of two sets of free-standing parallel silicon SWGs separated by an air gap with gap spacing g. The gap spacing is set smaller than the incident wavelength to ensure the near-field coupling effect. The SWGs in the upper and lower layer are denoted by SWG1 (gray) and SWG2 (blue), respectively. Each SWG has the same crystalline-silicon (c-Si) film thicknesses ti but different grating periods Λi, where i = 1, 2. The duty cycles of the gratings fi are expressed as aii with ai as their silicon beam widths. SWG1 and SWG2 form an equivalent composite grating with grating period Λ. The composite grating is polarization-sensitive by its nature of 1D periodicity. Incident beams with E-field polarization parallel and perpendicular to the grating bars are referred to transverse electric (TE) and transverse magnetic (TM) polarizations, respectively. With a proper oblique incidence of a free-space linearly polarized (TE or TM) input beam onto the composite grating, diffraction beams can be found in the reflection region. By specifically designing the grating structure, high-efficiency diffraction (∼90%) of the first-order can be achieved. The periods of the SWGs are carefully designed so that at the resonant wavelength λR, the WAM in SWG2 would be fully excited by the incident light via SWG1 through wavevector matching.

 figure: Fig. 1.

Fig. 1. (a) Cross-section of the unit cell of the proposed bi-layered composite grating. (b) Calculated diffraction efficiency spectrum of the 0th order and −1st order diffraction in the reflection region. Insets: y component electric field profile at $\bar{\lambda } = 1.775$ (left) and at $\bar{\lambda } = 1.785$ (right). The color scales are normalized by the modulus of the maximum electric field amplitude. The normalized structural parameters for this plot are: ${\bar{\Lambda }_1} = 1/2$, ${\bar{\Lambda }_2} = 1/3$, ${\bar{t}_1} = {\bar{t}_2} = 0.9$, $\bar{g} = 0.115$; f1 = f2 = 0.5, nsi = 3.47, θ = 63°, TE polarization of incidence.

Download Full Size | PDF

We employ the finite element method (FEM) method to numerically investigate the spectral response of the composite grating. In our simulations, the unit cell shown in Fig. 1(a) is built. To achieve high simulation accuracy, we use one unit cell of the composite grating with an equivalent grating period of Λ in the simulation. In the x- and y-directions of the simulation region, periodic boundary conditions are applied, while in the z-direction, perfectly matched layers (PMLs) are adopted. For PML regions a mapped mesh is defined, while for other regions a free triangular mesh is used. The numerical calculations are performed with a good convergence condition by setting the maximum mesh size to be λ/20. The unit cell is illuminated with a broadband linearly polarized (TE or TM) plane wave with an incident angle of θ. The dispersion effect of the c-Si film is modeled with its refractive index nsi measured by ellipsometry [38]. Without loss of generality, the structural parameters are normalized in the following manner: ${\bar{\Lambda }_i} = {\Lambda _i}/\Lambda $, ${\bar{t}_i} = {t_i}/\Lambda $, ${\bar{a}_i} = {a_i}/\Lambda $, $\bar{g} = g/\Lambda $ and the working wavelength is normalized to be $\bar{\lambda } = \lambda /\Lambda $. In this way, all parameters become dimensionless and can be applied to the whole spectrum of electromagnetic waves.

In Fig. 1(b), we present our numerical results on the far-field diffraction efficiency in the reflection region of the composite grating to illustrate the working mechanism. The composite grating structure is obtained by the design method we developed (see Section 4). The normalized device structural parameters are given in Fig. 1. A pronounced high-efficiency –1st order diffraction, manifesting a peak value of 89% with full width at half-maximum (FWHM) line width $\varDelta \bar{\lambda } = 0.01$ at $\bar{\lambda } = 1.785$(resonant wavelength) can be observed in Fig. 1(b). It can be seen from the right inset of Fig. 1(b) that at $\bar{\lambda } = 1.785$, TE0 WAM in SWG2 is strongly excited. On the contrary, the left inset shows weak field confinement in the grating layers at $\bar{\lambda } = 1.775$. Thus, the formation of the significant high-efficiency –1st order diffraction can be attributed to the coupling between the incident TE plane wave and TE0 WAM in SWG2.

3. Principle of high-efficiency diffraction by the bi-layered composite gratings

In this section, we explore the physics behind the high-efficiency first-order diffraction of the proposed bi-layered composite gratings. First, the oblique incidence is coupled to a WAM of SWG2 via wavevector-matching enabled by SWG1. Then, the resonant WAM of SWG2 is out-coupled to the first-order far-field diffraction via wavevector-matching by SWG1 for a second time. Meanwhile, all the other diffraction orders (0th order not included) are cut-off, resulting in high-efficiency first-order diffraction in the reflection region.

In the proposed bi-layered composite grating, the resonant SWG2 can be treated as a 1D photonic crystal operating at the bands below the light line [39]. It can be seen in the insets of Figs. 2(a) and 2(b) that SWG and 1D photonic crystal share similar TE0 and TM0 mode profiles, indicating that they can be unified with the same theoretical structure. By calculating the band diagram of SWG with a given grating profile, the normalized angular frequency ω of each WAM as functions of the normalized wavevector kx can be obtained [shown in Figs. 2(a) and 2(b)]. These bands are symmetric about kx = 0.5 due to Brillouin zone folding. We denote the WAMs as TE0, TE1, TE2, etc. for TE incidence; and TM0, TM1, TM2, etc. for TM incidence [40].

 figure: Fig. 2.

Fig. 2. Principle of high-efficiency diffraction of the bi-layered composite grating. (a) Calculated band diagram for the TE WAMs. Insets: normalized electric field profile exhibiting TE0 mode profile of (i) SWG and (ii) 1D photonic crystal. (b) Calculated band diagram for the TM WAMs. Insets: normalized magnetic field profile exhibiting TM0 mode profile of (i) SWG and (ii) 1D photonic crystal. The shaded regions represent radiated modes. Parameters for (a) and (b) are ${\bar{\Lambda }_2} = 1/3$, ${\bar{t}_2} = 0.9$, f2 = 0.5 and nsi = 3.47. (c) Schematic of the wavevector-matching enabled by the near-field coupling between SWG1 and SWG2. (d) Wavevector diagram sketching the wavevector-matching conditions. Parameters for (c) and (d) are ${\bar{\Lambda }_1} = 1/2$, ${\bar{\Lambda }_2} = 1/3$, $\bar{\lambda } = 1.782$, θ = 63°, θd = 63°. The color scales are normalized by the modulus of the maximum electric field amplitude.

Download Full Size | PDF

For a basic understanding of the physics behind high diffraction efficiency, the schematic shown in Fig. 2(c) provides an intuitive insight into the principle. The grating parameters are carefully designed such that TE0 WAM at $|{{k_x}} |= 0.5$ is coupled with the incident wave. When a linearly polarized (e.g., TE) plane wave obliquely incidents onto the bi-layered composite grating with an incident angle of $\theta$, the x-wavevector of the incoming light can be expressed as $k_{in}^\parallel{=} {k_0}\sin \theta$, where k0 = 2π/λ0 is the wavevector of the incident wave in free space with wavelength λ0. After being diffracted through SWG1, an additional wavevector component of 2πm11 will be added onto the incident x-wavevector, where m1 is the diffraction order of SWG1 and we choose m1 = −1 here. If the joint x-wavevector matches the wavevector of the WAM in SWG2 at kx = −0.5, which we call the in-couple condition:

$${k_0}\sin \theta - \frac{{2\pi }}{{{\Lambda _1}}} ={-} \frac{1}{2}\left( {\frac{{2\pi }}{{{\Lambda _2}}}} \right). $$

Light can be fully coupled to the WAM of SWG2 when the in-couple condition is satisfied. Then, the excited WAM in SWG2 will be further diffracted upwards or downwards by SWG2 and a wavevector component of 2πm22 will be added onto the joint x-wavevector, where m2 is the diffraction order of SWG2. The downward diffracted light cannot be further out-coupled to the far-field transmission region since the joint x-wavevector cannot match the wavevector in free space for any number of m2. On the contrary, as for the upward diffraction, the out-couple condition can be fulfilled when the joint x-wavevector picks up another x-wavevector component 2πm′11 provided by SWG1 where m′1 is the order of diffraction by SWG1 for a second time.

In order to develop an effective way to describe the wavevector-matching conditions, a conceptual bi-layered composite grating model based on a wavevector diagram is shown in Fig. 2(d). The circle in the diagram represents the allowed propagation wavevector in free space. The horizontal arrows represent the x-wavevectors of the incident and diffracted waves as well as the grating vectors. The wavevector-matching conditions can be described by consecutively adding the grating vectors to the x-wavevector of the incidence by the head-to-tail method. Firstly, the x-wavevector of the incident wave $k_{in}^\parallel$ (red dashed arrow) is diffracted by the –1st order of SWG1 (blue arrow), which fulfills the in-couple condition and leads to the excitation of a selected WAM in SWG2. Then, the wave is diffracted by the +1st order of SWG2 (green arrow). At this point, the joint x-wavevector is greater than the allowed propagation wavevector in free space thus leading to the cut-off of all diffractions in the transmission region except the 0th order. In the reflection region, on the other hand, the joint x-wavevector can pick up one additional –1st order grating vector from SWG1 thus leading to the final joint x-wavevector (yellow arrow) that can be coupled to free space thus leading to the high-efficiency –1st diffraction of the composite grating.

According to the above-analyzed wavevector-matching conditions, except for the –1st order diffraction in the reflection region, all the other diffraction orders (0th order not included) are cut-off. Hence, such high-efficiency –1st order diffraction is mainly attributed to the wavevector-matching enabled by near-field coupling between SWG1 and SWG2. From the wavevector diagram in Fig. 2(d), it is clear that the out-couple condition is fulfilled in the reflection region when m1 = −1, m2 = +1, and m′1 = −1, thus leading to:

$${k_0}\sin \theta - \frac{{2\pi }}{{{\Lambda _1}}} + \frac{{2\pi }}{{{\Lambda _2}}} - \frac{{2\pi }}{{{\Lambda _1}}} = {k_0}\sin \theta - \frac{{2\pi }}{\Lambda } = {k_0}\sin {\theta _d}. $$

The second half of Eq. (2) is the well-known grating equation representing the –1st order diffraction of the equivalent bi-layered composite grating with diffraction angle θd. Here, from the above equation, the equivalent period Λ of the composite grating is defined as:

$$\Lambda = \frac{{{\Lambda _1}{\Lambda _2}}}{{|{2{\Lambda _2} - {\Lambda _1}} |}}.$$

In the following design, the ratio of Λ12 is chosen as 3/2 to simplify the introduction of the design method. For other ratios of the SWG1 and SWG2 periods, high-efficiency –1st order diffraction can also be obtained by carefully selecting their grating periods following the wavevector-matching condition and optimizing their grating profiles.

4. Design method of the bi-layered composite gratings with high diffraction efficiency

To create a bi-layered composite grating with high-efficiency first-order diffraction, major structural parameters including the normalized grating periods, ${\bar{\Lambda }_i}$, the grating profile of SWG2 including t22 and a22, and the normalized gap spacing between the gratings $\bar{g}$ should be taken into consideration. Fig. 3 shows a flowchart of the proposed design approach.

 figure: Fig. 3.

Fig. 3. Flowchart of the proposed design procedure to obtain the bi-layered composite gratings with high-efficiency first-order diffraction.

Download Full Size | PDF

The design procedure is described as follows:

  • (1) Preset the working wavelength λ and resonant WAM in SWG2 for the desired application.
  • (2) Initialize the grating profile of SWG2. Since the high-efficiency performance of the composite grating originates from the WAM supported by SWG2, the grating profile of SWG2 including t22 and a22 is firstly initialized based on the available micro/nano fabrication capability.
  • (3) Determine the ratio of the grating periods. Choose a period ratio of Λ12 to form a composite grating with equivalent period Λ according to Eq. (3).
  • (4) Calculate the band diagram of SWG2. Given the determined grating profile of SWG2, its band diagram can be thus calculated. Consequently, the normalized angular frequency ω of the selected WAM and its corresponding normalized WAM wavevector, kx, can be obtained.
  • (5) Determine the SWG2 structural parameters according to the working wavelength λ. Calculate Λ2 according to the normalized angular frequency ω (ω = Λ2/λ) of the selected WAM determined in Step (4). Then, calculate t2 and a2 according to the grating profile of SWG2 (t22 and a22) initialized in Step (2).
  • (6) Calculate grating parameters of SWG1. Calculate Λ1 according to Λ12 determined in Step (3). Then, for the simplicity of the design, here we set the grating thicknesses and duty cycles of SWG1 to be t1 = t2 and f1 = f2. It is worth noting that SWG1 is used for wavevector-matching in the bi-layered composite grating structure. Therefore, when compared with SWG2, minor change of grating profile of SWG1 (t11 and a11) has little effect on the diffraction efficiency (see the simulation results in Section 5 for details).
  • (7) Sweep search for optimal gap spacing g. Since a change in the spacing between two layers of SWGs will affect the peak wavelength, efficiency, and bandwidth of the first-order diffraction, in order to determine the optimal spacing between the SWGs, it is necessary to consider specific application requirements including the desired diffraction efficiency and bandwidth, as well as the gap control accuracy in the actual device fabrication process.
  • (8) Output the design or repeat the design process. If an optimal gap spacing that meets the requirements is found, the loop is ended, and the final composite grating design is output. Otherwise, fine-tune the SWG2 profile and go back to Step (4), and repeat the process until a satisfactory design is found.

5. Results and discussions

To verify the proposed design method, two examples that employ respectively the TE0 and TM1 WAMs in the SWG2 are numerically investigated to achieve high-efficiency –1st order diffraction of the composite grating. We design the bi-layered composite gratings operate at wavelengths around 1550 nm. The refractive index of the c-Si measured by ellipsometry is 3.47 within this wavelength range, and the extinction coefficient of the c-Si is around 10−6, which is negligible in our numerical simulations. We then calculate the –1st order diffraction efficiency of the composite grating and run a parametric sweep to investigate the parametric dependence of both designs.

5.1 Design of the TE0 mode based bi-layered composite grating

We first design a TE0 mode based composite grating with its –1st order diffraction efficiency as high as 90% in the reflection region. In this design, the grating profile of SWG2 is set to be t22 = 0.897 and a22 = 0.5. For the simplicity of the design, the grating period ratio Λ12 is chosen at 3/2. As a result, the structural parameters of the bi-layered composite grating are determined to be Λ1 = 435 nm, Λ2 = 290 nm, Λ = 870 nm, t1 = t2 = 260 nm, and f1 = f2 = 0.5 according to the desired working wavelength around 1550 nm. Thus, the incident TE wave will be fully coupled with the TE0 mode of SWG2 leading to high-efficiency –1st order diffraction in the reflection region of the composite grating.

To make a systematic discussion of the TE0 mode based bi-layered composite grating, the dependences of the –1st order diffraction on various structural parameters are studied. Firstly, we investigate the influence of the variation of gap spacing on the –1st order diffraction spectrum of the composite grating while all the other parameters are fixed. As shown in Fig. 4(a), the TE0 mode supports a high-efficiency and narrow-bandpass performance. As the gap spacing increases from 50 nm to 300 nm when the incidence angle θ is 63° (corresponds to kx = 0.5), the diffraction efficiency decreases from 89% to 69% with a decreased FWHM (Δλ) from 6.4 nm to 0.1 nm.

 figure: Fig. 4.

Fig. 4. (a) Simulated –1st order diffraction efficiency spectrum of the composite grating in relation to the variation in gap spacing g when θ = 63° Insets: diffraction efficiency spectrum when g = 100 nm (top) and normalized electric field profile when g = 100 nm and λ = 1552 nm (bottom). The color scale is normalized by the modulus of the maximum electric field amplitude. (b) Simulated –1st order diffraction spectrum of the composite grating in relation to the variation in incident angle θ when g = 100 nm. Inset: a plot of diffraction efficiency as a function of θ at λ = 1552 nm (marked by dotted gray line). (c) Calculated resonant wavelength of TE0 mode when θ = 63° as a function of neff. (d) Calculated band diagram of the TE0 mode of SWG2. The range of kx when the incidence angle θ varying from 50° to 85° is marked by the black dotted lines.

Download Full Size | PDF

This phenomenon can be attributed to the near-field coupling of evanescent waves. Since the TE0 WAM is predominantly confined within SWG2 [see the bottom inset of Fig. 4(a)], an increase of gap spacing could weaken the coupling between SWG1 and SWG2, thus resulting in a decrease in diffraction efficiency and narrowing in bandwidth. Moreover, as g increases, a monotonous blueshift of peak wavelength can also be found. This characteristic can be roughly estimated by modeling the upper region of SWG2 as an isotropic medium with a refractive index of neff according to the effective medium theory (EMT) [41]. With an increase of g, the neff will be decreased, thus resulting in a blueshift of the resonant wavelength of the TE0 mode as shown in Fig. 4(c). In summary, by controlling the gap spacing, the diffraction efficiency and bandwidth of the composite grating can be adjusted.

Next, a calculated incident-angle-dependent –1st order diffraction spectrum is provided in Fig. 4(b). Setting at a gap spacing of 100 nm between two SWGs, when θ varies from 50° to 85°, the spectrum shows significant stability in the diffraction efficiency and bandwidth. It is interesting to observe from the inset of Fig. 4(b) that the spectrum exhibits certain tolerance against the incident angle variation. The calculated FWHM angular range (Δθ) is 17.6° at λ = 1552 nm with a maximum diffraction efficiency of 89%. Such an angular characteristic can be attributed to the flat angular dependence response of the TE0 mode at its band edge [seen from Fig. 4(d)]. This characteristic might be useful in applications requiring high angular stability.

Lastly, by scanning the duty cycles and thicknesses of the constituent SWGs while fixing all the other parameters, the influences of the SWG grating profiles on the –1st order diffraction efficiency of the bi-layered composite grating are studied. The incident angle θ is fixed at 63° and the gap spacing g between the two SWGs is kept at 100 nm. As shown in Figs. 5(a) and 5(c), the performance of the composite grating shows remarkable robustness against minor variations in the grating profile of SWG1. On the contrary, as shown in Figs. 5(b) and 5(d), the performance of the composite grating is quite sensitive to the variations in the grating profile of SWG2. The explanation to this phenomenon is as follows: since SWG1 in the bi-layered composite grating system serves as a wavevector matcher as can be seen from Eqs. (1) and (2), it is the grating period rather than the grating profile of SWG1 that affects the wavevector-matching conditions. In comparison, changes in the grating profile of SWG2, where the resonant TE0 WAM primarily resides, can affect the band structure and shift the resonant frequency of the WAM. As a result, the –1st order diffraction efficiency spectrum of the bi-layered composite grating is more sensitive to the variations of the SWG2 grating profile.

 figure: Fig. 5.

Fig. 5. Parametric dependences of the designed TE0 mode based bi-layered composite grating. Simulated –1st order diffraction efficiency spectra of the composite grating in relation to the variations in (a) and (b) grating thicknesses and (c) and (d) duty cycles of SWG1 and SWG2 at θ = 63° and g = 100 nm. Dotted white lines in (a)–(d) represent the designed composite grating parameters.

Download Full Size | PDF

5.2 Design of the TM1 mode based bi-layered composite grating

In the following design example, a high-efficiency TM1 mode based bi-layered composite grating with its –1st order diffraction efficiency as high as 90% is presented. In this design, the grating profile of SWG2 is set to t22 = 0.795 and a22 = 0.66. The grating period ratio Λ12 is also chosen to be 3/2. As a result, structural parameters of the bi-layered SWGs are determined to be Λ1 = 746 nm, Λ2 = 497 nm, Λ = 1492 nm, t1 = t2 = 395 nm, f1 = f2 = 0.66 according to our proposed design procedure with its working wavelength setting at 1550 nm.

We also investigate the dependence of the composite grating performance on the gap spacing and incident angle in a similar way as presented in Section 5.1. As shown in Fig. 6(a), when the gap spacing increases from 250 nm to 600 nm with the incident angle remained at 32° (corresponds to kx = 0.5), broadband diffraction will occur from 1.52 μm to 1.58 μm when g = 300 nm [see the bottom inset of Fig. 6(a)]. The bandwidth of the spectrum is gradually narrowed with the increase of g. This phenomenon can also be attributed to the near-field coupling of evanescent waves. As shown in the top inset of Fig. 6(a), the TM1 WAM decays relatively slowly into the near-field of SWG2. Therefore, the coupling between SWG1 and SWG2 in a TM1 mode composite grating is more tolerant of a larger gap spacing g when compared with a TE0 mode composite grating.

 figure: Fig. 6.

Fig. 6. (a) Simulated –1st order diffraction efficiency spectrum of the TM1 mode based composite grating in relation to the variation in gap spacing when θ = 32° Insets: normalized electric field profile when g = 300 nm and λ = 1544 nm (top) and diffraction efficiency spectrum when g = 300 nm (bottom). The color scale is normalized by the modulus of the maximum electric field amplitude. (b) Simulated –1st order diffraction spectrum of the composite grating in relation to the variation in incident angle when g = 300 nm. Inset: a plot of diffraction efficiency dependency on θ at λ = 1529 nm (marked by dotted gray line) showing high angular tolerance. (c) Calculated band diagram of the TM1 mode of SWG2. The range of kx when θ varying from 10° to 60° is marked by the black dotted lines.

Download Full Size | PDF

Next, a calculated incident-angle-dependent –1st order diffraction spectrum is provided in Fig. 6(b). When θ varies from 10° to 60° with the gap spacing g fixed at 300 nm, the high-efficiency region forms a curved strip as shown in the figure. The diffraction efficiency reaches its maximum (89.5%) at the center of the strip, where θ = 32° and λ = 1548 nm. This corresponds to the kx = 0.5 and ω = 0.33 point in the TM1 mode band diagram shown in Fig. 6(c). The band edge of the TM1 mode exhibits a concave curve when θ varying from 10° to 60° [marked by the black dotted lines in Fig. 6(c)], which represents the variation trajectory of the resonant frequency of the TM1 mode in the frequency domain. Therefore, the concave curve of TM1 mode in the frequency domain will translate into a convex strip region of high-efficiency –1st order diffraction in the wavelength domain. Moreover, the spectrum exhibits a high tolerance against the incident angle variation, and the calculated FWHM angular range (Δθ) is 42.5° at λ = 1529 nm with its maximum diffraction efficiency of 85% [see the inset of Fig. 6(b)].

Lastly, by scanning the duty cycles and thicknesses of the constituent SWGs while keeping all the other parameters unchanged, the influences of the SWG grating profiles on the –1st order diffraction efficiency of the bi-layered composite grating are studied when θ and g are kept at 32° and 300 nm, respectively. As shown in Figs. 7(a) and 7(c), the performance of the composite grating shows good robustness against minor variations in the grating profile of SWG1. On the contrary, as shown in Figs. 7(b) and 7(d), the –1st order diffraction efficiency of the composite grating is more sensitive to the variations in the grating profile of SWG2 as expected.

 figure: Fig. 7.

Fig. 7. Parametric dependence of the designed TM1 mode based composite grating. Simulated –1st order diffraction efficiency spectra of the composite grating in relation to the variations in (a) and (b) grating thicknesses and (c) and (d) duty cycles of SWG1 and SWG2 respectively when θ = 32° and g = 300 nm. Dotted black lines in (a)–(d) represent the designed composite grating parameters.

Download Full Size | PDF

5.3 Fabrication discussion

The proposed designs can be fabricated using CMOS-compatible micro/nano fabrication processes. First, a layer of the dielectric SWG pattern (SWG2) can be achieved on a silicon-on-insulator (SOI) wafer by lithography and dry etch. Next, a silicon dioxide (SiO2) layer can be deposited on the patterned SOI wafer, planarized using chemical mechanical polishing (CMP), and thinned down to the desired thickness. The SiO2 layer acts as a sacrificial layer to achieve the required gap spacing between the two SWGs. Next, a thin silicon layer is deposited on the SiO2 layer. With another lithography and dry etch process, another layer of SWG pattern (SWG1) can be obtained. Finally, the structure can be released through HF etching and CO2 critical drying to remove the sacrificial SiO2 layer and the buried oxide layer. Thus, the designed bi-layered composite grating structure can be achieved.

For the proposed bi-layered structure, the effects of alignment errors between the upper and lower layer are also investigated. The relative position of the bi-layered structure will be determined by the alignment of the two lithography steps during the device fabrication. Considering the matured micro/nano fabrication technology, alignment accuracy down to sub-micron scale can be easily achieved nowadays [42]. Misalignment occurring in the y direction [see Fig. 1(a)] makes no difference to the –1st order diffraction because the composite grating structure can be treated as sufficiently long in the y direction. On the other hand, for the misalignment occurring in the x-direction Δx, the influence of Δx on the –1st order diffraction efficiency of the designed TE0 and TM1 mode based composite gratings is investigated and the simulation results are provided in Figs. 8(a) and 8(b), respectively. It is obvious that both designs show significant robustness against the misalignment occurring in the x-direction.

 figure: Fig. 8.

Fig. 8. Simulated –1st order diffraction efficiency spectra of the designed (a) TE0 and (b) TM1 mode based composite gratings in relation to the variations in misalignment of the two grating layers in the x direction Δx. Parameters for plot (a) are Λ1 = 435 nm, Λ2 = 290 nm, Λ = 870 nm, t1 = t2 = 260 nm, f1 = f2 = 0.5, θ = 63° and g = 100 nm; Parameters for plot (b) are Λ1 = 746 nm, Λ2 = 497 nm, Λ = 1492 nm, t1 = t2 = 395 nm, f1 = f2 = 0.66, θ = 32° and g = 300 nm.

Download Full Size | PDF

6. Conclusion

In this letter, we propose a new strategy to design the bi-layered composite gratings with high-efficiency first-order diffraction in the reflection region. Although metasurface could also diffract the incident light to a desired order diffraction with high efficiency [4345], the proposed bi-layered composite gratings are different from metasurface in both design method and working mechanism. The proposed composite grating structure is composed of two different SWGs with thicknesses smaller than the operation wavelength, making the composite grating compact and integration-friendly. We build a simple and intuitive theoretical model to describe the bi-layered composite grating that yields high-efficiency first-order diffraction. Due to near-field coupling, oblique incident light is coupled with the WAM of the bottom grating by satisfying wavevector-matching conditions enabled by SWG1 and out-coupled to the first-order diffraction in the far-field reflection region. An effective design procedure based on the theoretical model is developed to obtain the desired bi-layered composite grating with high efficiency. In order to verify the proposed design method, two composite grating designs working on the TE0 or TM1 mode are demonstrated. Through controlling the gap spacing and the incident angle, both of the designed composite gratings could provide adjustable wavelength responses and show potentials in tunable bandwidth and diffraction efficiency. These results demonstrate the feasibility of the proposed design strategy for high-efficiency wave modulations. Such a bi-layered dielectric thin film structure can be effectively fabricated and characterized by several well-established techniques [46]. We believe that the rich properties of the bi-layered composite gratings could pave the way for the novel design of flat optical devices and aid in their future applications.

Funding

Ministry of Education - Singapore (MOE2019-T2-2-104).

Disclosures

The authors declare no conflict of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. N. Yu and F. Capasso, “Flat optics with designer metasurfaces,” Nat. Mater. 13, 139–150 (2014). [CrossRef]  

2. F. Capasso, “The future and promise of flat optics: a personal perspective,” Nanophotonics 7, 953–957 (2018). [CrossRef]  

3. N. Yu, P. Genevet, F. Aieta, M. A. Kats, R. Blanchard, G. Aoust, J.-P. Tetienne, Z. Gaburro, and F. Capasso, “Flat optics: Controlling wavefronts with optical antenna metasurfaces,” IEEE J. Sel. Top. Quantum Electron. 19, 4700423 (2013). [CrossRef]  

4. S. M. Kamali, E. Arbabi, A. Arbabi, and A. Faraon, “A review of dielectric optical metasurfaces for wavefront control,” Nanophotonics 7, 1041–1068 (2018). [CrossRef]  

5. H.-T. Chen, A. J. Taylor, and N. Yu, “A review of metasurfaces: Physics and applications,” Rep. Prog. Phys. 79, 076401 (2016). [CrossRef]  

6. C. L. Holloway, E. F. Kuester, J. A. Gordon, J. O’Hara, J. Booth, and D. R. Smith, “An overview of the theory and applications of metasurfaces: The two-dimensional equivalents of metamaterials,” IEEE Antennas Propag. Mag. 54, 10–35 (2012). [CrossRef]  

7. N. Meinzer, W. L. Barnes, and I. R. Hooper, “Plasmonic meta-atoms and metasurfaces,” Nat. Photonics 8, 889–898 (2014). [CrossRef]  

8. C. M. Soukoulis and M. Wegener, “Past achievements and future challenges in the development of three-dimensional photonic metamaterials,” Nat. Photonics 5, 523–530 (2011). [CrossRef]  

9. M. C. Huang, Y. Zhou, and C. J. Chang-Hasnain, “A surface-emitting laser incorporating a high-index-contrast subwavelength grating,” Nat. Photonics 1, 119–122 (2007). [CrossRef]  

10. B. Zhang, Z. Wang, S. Brodbeck, C. Schneider, M. Kamp, S. Höfling, and H. Deng, “Zero-dimensional polariton laser in a subwavelength grating-based vertical microcavity,” Light: Sci. Appl. 3, e135 (2014). [CrossRef]  

11. J. Wang, I. Glesk, and L. R. Chen, “Subwavelength grating filtering devices,” Opt. Express 22, 15335–15345 (2014). [CrossRef]  

12. C. Peng, L. Li, Q. Qiao, C. Yap, and G. Zhou, “Multifunctional silicon metagratings based on multiple periodicity design,” J. Opt. 22, 045103 (2020). [CrossRef]  

13. L. Zhang and D. Dai, “Silicon subwavelength-grating microdisks for optical sensing,” IEEE Photonics Technol. Lett. 31, 1209–1212 (2019). [CrossRef]  

14. L. Huang, H. Yan, X. Xu, S. Chakravarty, N. Tang, H. Tian, and R. T. Chen, “Improving the detection limit for on-chip photonic sensors based on subwavelength grating racetrack resonators,” Opt. Express 25, 10527–10535 (2017). [CrossRef]  

15. H. Yan, L. Huang, X. Xu, S. Chakravarty, N. Tang, H. Tian, and R. T. Chen, “Unique surface sensing property and enhanced sensitivity in microring resonator biosensors based on subwavelength grating waveguides,” Opt. Express 24, 29724–29733 (2016). [CrossRef]  

16. L. Liu, Q. Deng, and Z. Zhou, “Manipulation of beat length and wavelength dependence of a polarization beam splitter using a subwavelength grating,” Opt. Lett. 41, 5126–5129 (2016). [CrossRef]  

17. Y. Xu and J. Xiao, “Compact and high extinction ratio polarization beam splitter using subwavelength grating couplers,” Opt. Lett. 41, 773–776 (2016). [CrossRef]  

18. W. Yu, M. Ye, and Y. S. Yi, “Impacts of tapered sidewall profile on subwavelength grating wideband reflectors,” J. Nanophotonics 9, 093058 (2015). [CrossRef]  

19. J. Wang, I. Glesk, and L. Chen, “Subwavelength grating Bragg grating filters in silicon-on-insulator,” Electron. Lett. 51, 712–714 (2015). [CrossRef]  

20. M. Zohar, M. Auslender, and S. Hava, “Ultrathin high efficiency photodetectors based on subwavelength grating and near-field enhanced absorption,” Nanoscale 7, 5476–5479 (2015). [CrossRef]  

21. Z. Liu, S. Tibuleac, D. Shin, P. Young, and R. Magnusson, “High-efficiency guided-mode resonance filter,” Opt. Lett. 23, 1556–1558 (1998). [CrossRef]  

22. R. Magnusson and S. Wang, “Transmission bandpass guided-mode resonance filters,” Appl. Opt. 34, 8106–8109 (1995). [CrossRef]  

23. 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]  

24. H. Xu and Y. Shi, “Silicon-waveguide-integrated high-quality metagrating supporting bound state in the continuum,” Laser Photonics Rev. 14, 1900430 (2020). [CrossRef]  

25. K. Chen, X. Zhang, X. Chen, T. Wu, Q. Wang, Z. Zhang, Q. Xu, J. Han, and W. Zhang, “Active dielectric metasurfaces for switchable terahertz beam steering and focusing,” IEEE Photonics J. 13, 1–11 (2021). [CrossRef]  

26. W.-N. Liu, R. Chen, W.-Y. Shi, K.-B. Zeng, F.-L. Zhao, and J.-W. Dong, “Narrow-frequency sharp-angular filters using all-dielectric cascaded meta-gratings,” Nanophotonics 9, 3443–3450 (2020). [CrossRef]  

27. Y. Ra’di, D. L. Sounas, and A. Alù, “Metagratings: Beyond the limits of graded metasurfaces for wave front control,” Phys. Rev. Lett. 119, 067404 (2017). [CrossRef]  

28. D. Sell, J. Yang, S. Doshay, R. Yang, and J. A. Fan, “Large-angle, multifunctional metagratings based on freeform multimode geometries,” Nano Lett. 17, 3752–3757 (2017). [CrossRef]  

29. T. Kondo, S. Ura, and R. Magnusson, “Design of guided-mode resonance mirrors for short laser cavities,” J. Opt. Soc. Am. A 32, 1454–1458 (2015). [CrossRef]  

30. J. Inoue, T. Majima, K. Hatanaka, K. Kintaka, K. Nishio, Y. Awatsuji, and S. Ura, “Aperture miniaturization of guided-mode resonance filter by cavity resonator integration,” Appl. Phys. Express 5, 022201 (2012). [CrossRef]  

31. R. R. Boye and R. K. Kostuk, “Investigation of the effect of finite grating size on the performance of guided-mode resonance filters,” Appl. Opt. 39, 3649–3653 (2000). [CrossRef]  

32. W. Liu, Z. Lai, H. Guo, and Y. Liu, “Guided-mode resonance filters with shallow grating,” Opt. Lett. 35, 865–867 (2010). [CrossRef]  

33. K. Kintaka, T. Majima, J. Inoue, K. Hatanaka, J. Nishii, and S. Ura, “Cavity-resonator-integrated guided-mode resonance filter for aperture miniaturization,” Opt. Express 20, 1444–1449 (2012). [CrossRef]  

34. K. Lee, R. LaComb, B. Britton, M. Shokooh-Saremi, H. Silva, E. Donkor, Y. Ding, and R. Magnusson, “Silicon-layer guided-mode resonance polarizer with 40-nm bandwidth,” IEEE Photonics Technol. Lett. 20, 1857–1859 (2008). [CrossRef]  

35. K. J. Lee, J. Curzan, M. Shokooh-Saremi, and R. Magnusson, “Resonant wideband polarizer with single silicon layer,” Appl. Phys. Lett. 98, 211112 (2011). [CrossRef]  

36. Z. Liu, S. Durant, H. Lee, Y. Xiong, Y. Pikus, C. Sun, and X. Zhang, “Near-field Moire effect mediated by surface plasmon polariton excitation,” Opt. Lett. 32, 629–631 (2007). [CrossRef]  

37. A. Hessel and A. Oliner, “A new theory of Wood’s anomalies on optical gratings,” Appl. Opt. 4, 1275–1297 (1965). [CrossRef]  

38. D. F. Edwards, “Silicon (Si),” in Handbook of Optical Constants of Solids, E. D. Palik, ed. (Academic, 1985).

39. S. Fan, J. N. Winn, A. Devenyi, J. Chen, R. D. Meade, and J. Joannopoulos, “Guided and defect modes in periodic dielectric waveguides,” J. Opt. Soc. Am. B 12, 1267–1272 (1995). [CrossRef]  

40. C. J. Chang-Hasnain and W. Yang, “High-contrast gratings for integrated optoelectronics,” Adv. Opt. Photonics 4, 379–440 (2012). [CrossRef]  

41. P. Lalanne and J.-P. Hugonin, “High-order effective-medium theory of subwavelength gratings in classical mounting: application to volume holograms,” J. Opt. Soc. Am. A 15, 1843–1851 (1998). [CrossRef]  

42. R. S. Ghaida, G. Torres, and P. Gupta, “Single-mask double-patterning lithography for reduced cost and improved overlay control,” IEEE Trans. Semicond. Manufact. 24, 93–103 (2011). [CrossRef]  

43. K. E. Chong, L. Wang, I. Staude, A. R. James, J. Dominguez, S. Liu, G. S. Subramania, M. Decker, D. N. Neshev, and I. Brener, “Efficient polarization-insensitive complex wavefront control using Huygens’ metasurfaces based on dielectric resonant meta-atoms,” ACS Photonics 3, 514–519 (2016). [CrossRef]  

44. W. Zhao, H. Jiang, B. Liu, J. Song, Y. Jiang, C. Tang, and J. Li, “Dielectric Huygens’ metasurface for high-efficiency hologram operating in transmission mode,” Sci. Rep. 6, 1–7 (2016). [CrossRef]  

45. G. Zheng, H. Mühlenbernd, M. Kenney, G. Li, T. Zentgraf, and S. Zhang, “Metasurface holograms reaching 80% efficiency,” Nat. Nanotechnol. 10, 308–312 (2015). [CrossRef]  

46. W. Yang, T. Sun, Y. Rao, M. Megens, T. Chan, B.-W. Yoo, D. A. Horsley, M. C. Wu, and C. J. Chang-Hasnain, “High speed optical phased array using high contrast grating all-pass filters,” Opt. Express 22, 20038–20044 (2014). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1.
Fig. 1. (a) Cross-section of the unit cell of the proposed bi-layered composite grating. (b) Calculated diffraction efficiency spectrum of the 0th order and −1st order diffraction in the reflection region. Insets: y component electric field profile at $\bar{\lambda } = 1.775$ (left) and at $\bar{\lambda } = 1.785$ (right). The color scales are normalized by the modulus of the maximum electric field amplitude. The normalized structural parameters for this plot are: ${\bar{\Lambda }_1} = 1/2$ , ${\bar{\Lambda }_2} = 1/3$ , ${\bar{t}_1} = {\bar{t}_2} = 0.9$ , $\bar{g} = 0.115$ ; f1 = f2 = 0.5, nsi = 3.47, θ = 63°, TE polarization of incidence.
Fig. 2.
Fig. 2. Principle of high-efficiency diffraction of the bi-layered composite grating. (a) Calculated band diagram for the TE WAMs. Insets: normalized electric field profile exhibiting TE0 mode profile of (i) SWG and (ii) 1D photonic crystal. (b) Calculated band diagram for the TM WAMs. Insets: normalized magnetic field profile exhibiting TM0 mode profile of (i) SWG and (ii) 1D photonic crystal. The shaded regions represent radiated modes. Parameters for (a) and (b) are ${\bar{\Lambda }_2} = 1/3$ , ${\bar{t}_2} = 0.9$ , f2 = 0.5 and nsi = 3.47. (c) Schematic of the wavevector-matching enabled by the near-field coupling between SWG1 and SWG2. (d) Wavevector diagram sketching the wavevector-matching conditions. Parameters for (c) and (d) are ${\bar{\Lambda }_1} = 1/2$ , ${\bar{\Lambda }_2} = 1/3$ , $\bar{\lambda } = 1.782$ , θ = 63°, θd = 63°. The color scales are normalized by the modulus of the maximum electric field amplitude.
Fig. 3.
Fig. 3. Flowchart of the proposed design procedure to obtain the bi-layered composite gratings with high-efficiency first-order diffraction.
Fig. 4.
Fig. 4. (a) Simulated –1st order diffraction efficiency spectrum of the composite grating in relation to the variation in gap spacing g when θ = 63° Insets: diffraction efficiency spectrum when g = 100 nm (top) and normalized electric field profile when g = 100 nm and λ = 1552 nm (bottom). The color scale is normalized by the modulus of the maximum electric field amplitude. (b) Simulated –1st order diffraction spectrum of the composite grating in relation to the variation in incident angle θ when g = 100 nm. Inset: a plot of diffraction efficiency as a function of θ at λ = 1552 nm (marked by dotted gray line). (c) Calculated resonant wavelength of TE0 mode when θ = 63° as a function of neff. (d) Calculated band diagram of the TE0 mode of SWG2. The range of kx when the incidence angle θ varying from 50° to 85° is marked by the black dotted lines.
Fig. 5.
Fig. 5. Parametric dependences of the designed TE0 mode based bi-layered composite grating. Simulated –1st order diffraction efficiency spectra of the composite grating in relation to the variations in (a) and (b) grating thicknesses and (c) and (d) duty cycles of SWG1 and SWG2 at θ = 63° and g = 100 nm. Dotted white lines in (a)–(d) represent the designed composite grating parameters.
Fig. 6.
Fig. 6. (a) Simulated –1st order diffraction efficiency spectrum of the TM1 mode based composite grating in relation to the variation in gap spacing when θ = 32° Insets: normalized electric field profile when g = 300 nm and λ = 1544 nm (top) and diffraction efficiency spectrum when g = 300 nm (bottom). The color scale is normalized by the modulus of the maximum electric field amplitude. (b) Simulated –1st order diffraction spectrum of the composite grating in relation to the variation in incident angle when g = 300 nm. Inset: a plot of diffraction efficiency dependency on θ at λ = 1529 nm (marked by dotted gray line) showing high angular tolerance. (c) Calculated band diagram of the TM1 mode of SWG2. The range of kx when θ varying from 10° to 60° is marked by the black dotted lines.
Fig. 7.
Fig. 7. Parametric dependence of the designed TM1 mode based composite grating. Simulated –1st order diffraction efficiency spectra of the composite grating in relation to the variations in (a) and (b) grating thicknesses and (c) and (d) duty cycles of SWG1 and SWG2 respectively when θ = 32° and g = 300 nm. Dotted black lines in (a)–(d) represent the designed composite grating parameters.
Fig. 8.
Fig. 8. Simulated –1st order diffraction efficiency spectra of the designed (a) TE0 and (b) TM1 mode based composite gratings in relation to the variations in misalignment of the two grating layers in the x direction Δx. Parameters for plot (a) are Λ1 = 435 nm, Λ2 = 290 nm, Λ = 870 nm, t1 = t2 = 260 nm, f1 = f2 = 0.5, θ = 63° and g = 100 nm; Parameters for plot (b) are Λ1 = 746 nm, Λ2 = 497 nm, Λ = 1492 nm, t1 = t2 = 395 nm, f1 = f2 = 0.66, θ = 32° and g = 300 nm.

Equations (3)

Equations on this page are rendered with MathJax. Learn more.

k 0 sin θ 2 π Λ 1 = 1 2 ( 2 π Λ 2 ) .
k 0 sin θ 2 π Λ 1 + 2 π Λ 2 2 π Λ 1 = k 0 sin θ 2 π Λ = k 0 sin θ d .
Λ = Λ 1 Λ 2 | 2 Λ 2 Λ 1 | .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.