• Photonics Research
  • Vol. 10, Issue 4, A57 (2022)
Carlos Pérez-Armenta1、*, Alejandro Ortega-Moñux1, José Manuel Luque-González1, Robert Halir1、2, Pedro J. Reyes-Iglesias1, Jens Schmid3, Pavel Cheben3, Íñigo Molina-Fernández1、2, and J. Gonzalo Wangüemert-Pérez1、2
Author Affiliations
  • 1Telecommunication Research Institute (TELMA), Universidad de Málaga, CEI Andalucía TECH, E.T.S.I. Telecomunicación, 29010 Málaga, Spain
  • 2Bionand Center for Nanomedicine and Biotechnology, Parque Tecnológico de Andalucía, Málaga 29590, Spain
  • 3National Research Council Canada, Ottawa, Ontario K1A 0R6, Canada
  • show less
    DOI: 10.1364/PRJ.446932 Cite this Article Set citation alerts
    Carlos Pérez-Armenta, Alejandro Ortega-Moñux, José Manuel Luque-González, Robert Halir, Pedro J. Reyes-Iglesias, Jens Schmid, Pavel Cheben, Íñigo Molina-Fernández, J. Gonzalo Wangüemert-Pérez, "Polarization-independent multimode interference coupler with anisotropy-engineered bricked metamaterial," Photonics Res. 10, A57 (2022) Copy Citation Text show less

    Abstract

    Many applications, including optical multiplexing, switching, and detection, call for low-cost and broadband photonic devices with polarization-independent operation. While the silicon-on-insulator platform is well positioned to fulfill most of these requirements, its strong birefringence hinders the development of polarization-agnostic devices. Here we leverage the recently proposed bricked metamaterial topology to design, for the first time, to our knowledge, a polarization-independent 2×2 multimode interference coupler using standard 220 nm silicon thickness. Our device can be fabricated with a single etch step and is optimized for the O-band, covering a wavelength range of 160 nm with excess loss, polarization-dependent loss, and imbalance below 1 dB and phase errors of less than 5°, as demonstrated with full three-dimensional finite-difference time-domain simulations.

    1. INTRODUCTION

    Silicon-on-insulator (SOI) is a prominent platform for the development of photonic integrated circuits (PICs). The high refractive index contrast between silicon and silicon dioxide enables the development of compact integrated photonic devices. Combined with CMOS-compatible fabrication, these benefits are driving the growth of the silicon photonic transceiver market [1] as well as research in biochemical sensing [2], lidars [3], quantum information processing [4], and other important applications. However, devices implemented in SOI typically suffer from high birefringence, especially for the widely used 220 nm waveguide thickness [1]. Thus, PICs often work for only one polarization at the operating wavelength.

    Since the light coming from an optical fiber is usually in a random state of polarization, some applications require the use of polarization-independent devices. This is the case, for example, when designing building blocks for passive optical networks (PONs). Furthermore, standardization groups from ITU-T and IEEE are working on increasing the PON line rate by operating in the O-band where the low chromatic dispersion of optical fibers reduces intersymbol interference [5]. The design of transceivers to operate in the O-band is also of great interest in metro and access networks to provide point-to-point interfaces with the high capacity and low latency required to support radio access networks for 5G and beyond [6]. On the other hand, a widely used solution to increase receiver sensitivity, and therefore to enhance the optical power budget link, is to properly combine the incoming signal with a local oscillator in an optical hybrid before beating them in photodiodes. As an alternative to the conventional intradyne coherent receiver assisted by digital signal processing, much simpler reception schemes for PON have been recently proposed [7,8]. In these approaches, the availability of integrated polarization-independent hybrids would enable the implementation of low-cost and compact transceivers.

    Approaches for implementing polarization-independent devices include using square waveguides [9], employing polarization diversity schemes at the cost of doubling the device footprint and additional insertion loss [10], or optimizing the device to match the behavior of both transverse electric (TE) and transverse magnetic (TM) polarizations [11]. Examples of the latter approach include power splitters using cascaded bent directional couplers [12], mode-evolution directional couplers [13], optimized tapered couplers [14], broadband power dividers with modal-engineered slot waveguides [15], or filters using Bragg gratings [16,17].

    Multimode interference (MMI) devices are key building blocks in PICs, as they can be used to build different components such as power splitters, Mach–Zehnder interferometers, or coherent receivers [18,19]. MMIs utilize wide and highly birefringent waveguides typically optimized for TE modes, making it difficult to achieve polarization-insensitive behavior. A strategy for making polarization-independent MMI devices is based on using a particular MMI width for which the TE and TM self-imaging distances are the same [20,21]. Nevertheless, this approach is not effective for the 220 nm thick SOI platform, since the necessary waveguide dimensions would be too small for proper multimode imaging. In Ref. [20], the minimum studied SOI thickness was 300 nm, for which the optimal width was found to be as small as 2 μm, hence difficult to scale down to 220 nm SOI. An additional silicon nitride overlay (SiN-on-SOI) [22] can be used to implement polarization-insensitive MMI couplers, taking advantage of the lower index contrast, at the cost of an increased fabrication complexity. Yet another strategy to achieve polarization-insensitive MMI couplers is using subwavelength grating (SWG) metamaterials. SWGs, since their first demonstration in silicon waveguides [2327], have been used as a powerful tool for overcoming performance limitations of conventional silicon-based integrated photonic devices [2830]. Various polarization-agnostic devices designed through SWG engineering have been demonstrated, including directional couplers [3133], grating couplers [34], and non-birefringent tilted SWG waveguides [35]. A 1×2 MMI power splitter using an SWG slot and a silicon layer thickness of 250 nm has also been proposed [36].

    Recently, a brick-patterned subwavelength metamaterial has been demonstrated [37], enabling engineering of the anisotropy of SWGs. Furthermore, it preserves a single etch step process and Manhattan-like geometry, thereby facilitating wafer-scale fabrication. In this work, we explore the properties of bricked SWG waveguides to design, for the first time, a polarization-insensitive 2×2  90° MMI coupler for the 220 nm SOI platform. We engineer the geometry of the waveguide to cover the wavelength ranges of 1.26–1.42 μm, with excess loss (EL), polarization-dependent loss (PDL), and imbalance (IB) below 1 dB and phase error (PE) lower than 5°, operating simultaneously for both TE and TM polarizations.

    The paper is organized as follows. In Section 2, we summarize the working principle of the proposed device. In Section 3, a two-dimensional (2D) anisotropic model for the bricked SWG waveguide is presented. In Section 4, the design methodology is described, and in Section 5, the simulation results are discussed. Finally, conclusions are drawn in Section 6.

    2. BRICKED METAMATERIAL MMI WORKING PRINCIPLE

    The proposed device is shown in Fig. 1(a). It is a 2×2 MMI with the multimode region comprising a bricked SWG metamaterial, including the corresponding transitions between the interconnecting conventional silicon (Si) waveguide of width win and the bricked SWG input/output access waveguides of width wa. As shown in Fig. 1(c), the brick pattern includes two periodicities: the first with a period Λz=az+bz and duty cycle (DC) DCz=az/Λz along the propagation direction (z), and the second with a period Λx=ax+bx and DCx=ax/Λx in the transversal direction (x). The longitudinal displacement dz shapes the bricked geometry and significantly contributes to the metamaterial anisotropy. For dz=0, the structure is a conventional SWG waveguide, whereas for DCz<dz/Λz<1DCz, the silicon segments are isolated from each other. The width and length of the multimode waveguide are wMMI and LMMI=NMMIΛz, respectively, where NMMI is the number of SWG periods. Wire-to-bricked SWG tapers, shown schematically in Fig. 1(b), are located at the input and output ports of the MMI. The tapering is used to widen the input access waveguides as in conventional MMI devices to ensure that only lower-order modes are excited, while at the same time to adapt the bricked SWG multimode waveguide to the interconnecting strip waveguides. The device is designed for the SOI platform, with a silicon core thickness of h=220  nm and silicon dioxide (SiO2) upper and lower cladding. The corresponding refractive indices are nSi=3.503 and nSiO2=1.466 at the wavelength of 1.31 μm.

    (a) Schematic of the proposed polarization-insensitive 2×2 MMI (SiO2 cladding is not shown for clarity). (b) Top view of the taper from input silicon wire to bricked SWG waveguide. (c) Top view of the bricked SWG multimode waveguide.

    Figure 1.(a) Schematic of the proposed polarization-insensitive 2×2 MMI (SiO2 cladding is not shown for clarity). (b) Top view of the taper from input silicon wire to bricked SWG waveguide. (c) Top view of the bricked SWG multimode waveguide.

    The operation of the MMI coupler is governed by the self-imaging effect [18]. For a general interference case, the distance at which the first two-fold image of the input field profile is formed is LMMI=32Lπ,where Lπ is the beat length of the two lowest-order modes given by Lπ=πβ0β1.

    In Eq. (2), β0 and β1 are the propagation constants of the two lowest-order modes of the multimode region for a given polarization, TE or TM. Ideally, when the MMI length is given by Eq. (1), a signal entering port 1 [Fig. 1(a)] should be equally split between ports 2 and 3 with a relative phase shift of 90°. For a conventional MMI, the beat length substantially differs for TE and TM polarizations. The primary objective of this work is to optimize the bricked SWG core to achieve identical beat lengths for both polarizations, Lπ(TE)=Lπ(TM), thus obtaining a polarization-insensitive device. For that purpose, the longitudinal displacement dz of the silicon blocks of size ax×az is used. This parameter is crucial for controlling the anisotropy of the bricked SWG waveguides. In Ref. [37], it was shown that TE modes are significantly more sensitive to longitudinal displacement dz than TM modes. Here we leverage this unique advantage of bricked SWG topology to design a polarization-insensitive MMI.

    3. 2D ANISOTROPIC MODEL OF THE BRICKED SWG MMI

    MMI couplers implemented in 220 nm SOI waveguides are single mode in the out-of-plane (y) direction. Therefore, an analytic formula for the beat length of a homogeneous MMI can be obtained based on paraxial approximation using equivalent 2D waveguides, e.g.,  by the effective index method (EIM) [18]. In Ref. [38], a generalization of the EIM was proposed to characterize anisotropic SWG multimode waveguides for TE modes. In this section, we generalize the procedure outlined in Ref. [38] by extending it to TM polarization, and apply this new model to the bricked SWG topology.

    The in-plane (TE) modes of a slab waveguide with an anisotropic core obey the dispersion relation [38] (kx,m(TE)nzz)2+(βm(TE)nxx)2=(2πλ)2,where kx,m(TE)=mπ/we,TE is the x component of the mth TE mode wave vector, we,TE is the effective width of the fundamental TE mode, βm(TE) is the propagation constant of the mth TE mode, nii2 is the i component of the electric permittivity tensor, and λ is the vacuum wavelength. The nxx and nzz terms can be calculated by solving the unit cell of Fig. 2(a), e.g.,  using the rigorous three-dimensional (3D) full-vectorial MIT Photonic Bands (MPB) tool [39]. As outlined in Fig. 2(a), to obtain nxx and nzz, the effective index of the fundamental TE mode is calculated when propagating in the z and x directions, respectively. Note that this structure is infinitely periodic in the x and z directions. Then, a first-order Taylor expansion in Eq. (3) yields the following approximate beat length: Lπ(TE)4we,TE23λnzz2nxx.

    (a) Procedure employed for modeling the bricked SWG. (b) Calculated refractive index components for a bricked SWG waveguide with h=220 nm, Λx=200 nm, DCx=50%, Λz=240 nm, and DCz=50% at the wavelength λ=1.31 μm. (c) Beat length for wMMI=3 μm using the 2D anisotropic model (solid curves) and full 3D simulations (dotted curves).

    Figure 2.(a) Procedure employed for modeling the bricked SWG. (b) Calculated refractive index components for a bricked SWG waveguide with h=220  nm, Λx=200  nm, DCx=50%, Λz=240  nm, and DCz=50% at the wavelength λ=1.31  μm. (c) Beat length for wMMI=3  μm using the 2D anisotropic model (solid curves) and full 3D simulations (dotted curves).

    Here we extend this 2D anisotropic model to TM polarization by incorporating two additional refractive index components, nyy and nyy. These are also calculated using MPB, as generally represented in Fig. 2(a). To find nyy and nyy, the fundamental TM mode effective index is calculated for propagation along the z and x directions, respectively. By analogy with Eq. (4), the beat length of TM modes is given by Lπ(TM)4we,TM23λ(nyy)2nyy.

    A polarization-independent operation of the MMI is achieved provided that the corresponding beat lengths of TE and TM polarizations are designed to be identical, i.e., we,TE2nzz2/nxx=we,TM2(nyy)2/nyy.

    In the design example of Fig. 2(b), the nii components of a bricked SWG waveguide with h=220  nm, Λx=200  nm, DCx=50%, Λz=240  nm, and DCz=50% are obtained for a central wavelength of λ=1.31  μm. Using these results, along with the effective width [18] corresponding to the waveguide width of wMMI=3  μm for both polarizations in Eqs. (4) and (5), yields the beat lengths shown in Fig. 2(c) (solid lines). Furthermore, a bricked SWG waveguide with the same parameters as the unit cell and the same width, wMMI=3  μm, is simulated in 3D with the MPB solver for the two first modes of each polarization. The calculated propagation constants are substituted in Eq. (2) to obtain the corresponding beat lengths. Results are shown in Fig. 2(c), as dotted curves. We observe that the 2D model correctly follows the trend confirmed by the full 3D simulations using MPB, which predicts the optimum shifting value dzΛz/2 to obtain Lπ(TE)=Lπ(TM). The accuracy of the simplified 2D model is highly dependent on the shifting and polarization, but it gives accurate results near the design point. The error in the 2D estimation of Lπ at dz=Λz/2 is negligible for TE polarization, and about 1 μm for TM polarization. It is also observed that for dz=0, which corresponds to a conventional SWG waveguide (without bricking), the difference between TE and TM beat lengths is maximum. A grid resolution of 128 pixels/μm and cell size of 6  μm×4  μm×Λz were used for the 3D waveguide simulation, while for the 2D anisotropic model, the cell size of Λx×4  μm×Λz was employed. The sweep of Fig. 2(c) lasted 2 h for the 3D simulation waveguide, while for 2D approximation, the permittivity tensor was calculated in less than 90 s (running on a single core of an Intel Xeon Gold 6242R at 3.10 GHz), yielding about two orders of magnitude difference in simulation time. Therefore, using our 2D anisotropic model, a good design starting point can be obtained, substantially reducing the optimization time for the full 3D structure.

    4. DESIGN OF THE BRICKED SWG MULTIMODE WAVEGUIDE

    The geometrical parameters of the central multimode waveguide are optimized using the procedure outlined in Fig. 3(a). The multimode waveguide width is initially chosen as wMMI=3  μm to ensure that it supports a sufficient number of modes without substantially increasing the waveguide length, as the latter increases with wMMI2. The DC in z direction is DCz=50%, and the transversal periodicity is set to Λx=200  nm and DCx=50%, to maintain a minimum feature size of 100 nm and work far from Bragg resonance. More details on this aspect are presented in Appendix A.

    Design methodology used in this work, comprising two main stages: (a) design of the bricked SWG multimode waveguide and (b) optimization of the complete device using a 3D-FDTD simulator.

    Figure 3.Design methodology used in this work, comprising two main stages: (a) design of the bricked SWG multimode waveguide and (b) optimization of the complete device using a 3D-FDTD simulator.

    To obtain the first design iteration, the SWG structure was simulated using the approximate 2D anisotropic model. Figure 4 shows a contour map of the relative difference between the beat lengths for TE and TM polarizations, defined as ΔLπ¯=ΔLπ/Lπ¯, where ΔLπ=Lπ(TE)Lπ(TM) and Lπ¯=(Lπ(TE)+Lπ(TM))/2. The parameter ΔLπ was calculated as a function of the longitudinal pitch Λz and the normalized shift dz/Λz, for λ=1310  nm. Some important observations on polarization effects are concluded from these results, as shown in Fig. 4: (i) for a given pitch Λx, there is a minimum Λz below which polarization insensitivity cannot be achieved, (ii) the optimal bricked shifting dz is approximately Λz/2, and (iii) longer periods restrict the selection of dz as the Bragg regime approaches the operation band. So, with a minimal computational effort, we choose Λz=240  nm and dz=Λz/2 as the initial design point, while the other structural parameters were set as stated above.

    Absolute value of the relative difference between TE and TM beat lengths |ΔLπ¯|=|ΔLπ/Lπ¯|, where ΔLπ=Lπ(TE)−Lπ(TM) and Lπ¯=(Lπ(TE)+Lπ(TM))/2, calculated using the anisotropic 2D model. Simulation parameters of the unit cell are h=220 nm, DCz=DCx=50%, Λx=200 nm, and λ=1.31 μm. The dashed line represents the points (dz,Λz) for which ΔLπ¯=0. The optimal parameters Λz=240 nm and dz=Λz/2 are marked with a green dot.

    Figure 4.Absolute value of the relative difference between TE and TM beat lengths |ΔLπ¯|=|ΔLπ/Lπ¯|, where ΔLπ=Lπ(TE)Lπ(TM) and Lπ¯=(Lπ(TE)+Lπ(TM))/2, calculated using the anisotropic 2D model. Simulation parameters of the unit cell are h=220  nm, DCz=DCx=50%, Λx=200  nm, and λ=1.31  μm. The dashed line represents the points (dz,Λz) for which ΔLπ¯=0. The optimal parameters Λz=240  nm and dz=Λz/2 are marked with a green dot.

    As the next step, we perform full-vectorial 3D simulations of the bricked SWG waveguide near this initial design point, to correct for the inaccuracy of the 2D model. Using MPB, we calculate the two lowest-order Floquet–Bloch modes of the multimode section. From the condition Lπ(TE)=Lπ(TM) at λ=1.31  μm, we obtain Λz=220  nm and dz=0.48Λz. The beat length is then calculated using Eq. (2) to be 18 μm. The calculated curves for TE and TM beat lengths as a function of wavelength are shown in Fig. 5.

    Beat length as a function of wavelength obtained from 3D Floquet–Bloch simulations of the (bricked) SWG waveguide. Solid lines correspond to dz=0, while dashed lines correspond to dz=0.48Λz. Other geometrical parameters are: wMMI=3 μm, h=220 nm, Λx=200 nm, DCx=50%, Λz=220 nm, and DCz=50%.

    Figure 5.Beat length as a function of wavelength obtained from 3D Floquet–Bloch simulations of the (bricked) SWG waveguide. Solid lines correspond to dz=0, while dashed lines correspond to dz=0.48Λz. Other geometrical parameters are: wMMI=3  μm, h=220  nm, Λx=200  nm, DCx=50%, Λz=220  nm, and DCz=50%.

    5. DESIGN OF THE COMPLETE DEVICE

    To complete the full design, the procedure outlined in Fig. 3(b) is applied. We first optimize the access ports and the transitions between the interconnecting waveguides and the bricked SWG waveguide. The width of the MMI access ports wa is designed to ensure that only lower-order modes of the multimode waveguide are excited [38]. Subsequently, the separation between the MMI ports ws is chosen to minimize the coupling between the access waveguides. The transition from the interconnecting single-mode waveguides of width win=400  nm to the access ports wa is shown in Fig. 1(b). It is divided into two sections of NT periods each, where the first one converts a photonic wire into an SWG waveguide and the second one progressively introduces the brick shifting. This design strategy effectively reduces Bragg reflections along the taper. Simulation of the full device is carried out using the 3D full-vectorial finite-difference time-domain (FDTD) open-source simulator Meep [40]. The simulation parameters used are a grid resolution of 128 pixels/μm and a simulation window size of 7  μm×5  μm×Lz, where the length Lz depends on the MMI and taper lengths.

    To assess the overall performance of the device, we calculated the EL [EL=10log(|S21|2+|S31|2)], IB (IB=10log|S31/S21|2), PE [PE=arg(S31/S21)90°], and PDL (PDL=|ELTEELTM|) from the S-parameters obtained by 3D FDTD simulations. Some of the device parameters (dz, Λz, NMMI, and wMMI) were fine-tuned to optimize the performance of the device in terms of EL, IB, PE, and PDL. The initial design point obtained from the first design stage [Fig. 3(a)] was good enough to require only minor adjustments in the second stage [Fig. 3(b)], so that no additional optimization routines to further optimize EL, IB, PE, and PDL were required. The optimized final design parameters are given in Table 1.

    Optimized Design Parameters of Polarization Insensitive 2×2 MMI Coupler

    ParameterwMMIΛz (DCz)Λx (DCx)dzNMMIwawsNT
    Value2.9 μm220 nm (50%)200 nm (50%)110 nm1211.2 μm0.8 μm21

    Final design performance of the optimized polarization-independent 2×2 MMI.

    Figure 6.Final design performance of the optimized polarization-independent 2×2 MMI.

    Regarding fabrication tolerance, the robustness of the nominal design to fabrication deviations has been studied for in-plane over-etching and under-etching errors up to 20 nm as illustrated in the inset of Fig. 7(b). Figure 7 shows the center wavelength and bandwidth of the device as a function of the fabrication error δ, where δ<0 means over-etching and δ>0 under-etching, in the chip plane. We define the bandwidth of the device as the wavelength range in which PDL, EL, and |IB| are simultaneously below 1 dB and |PE| is less than 5°. It can be seen that for |δ|<15  nm, the MMI has a bandwidth greater than 100 nm for both polarizations. In addition, the center wavelength has a small dependence on δ for TM polarization. Although TE center wavelength is heavily shifted when δ>0, it is robust to fabrication errors of δ=10  nm. This asymmetric behavior with respect to δ can be taken into account by introducing a small bias in the mask definition stage. We also studied a ±2% variation of the grating periods (Λx,Λz), which resulted in a ±30  nm detuning of the device wavelength response.

    Dependence of the bandwidth and center wavelength of the MMI on fabrication error δ. The bandwidth is defined as the wavelength range in which PDL, EL, and |IB| are simultaneously below 1 dB, and |PE| is less than 5°. The inset shows the effect on the in-plane geometry of positive and negative values of δ. The original geometry (top view, i.e., in the chip plane) is represented with a continuous line, while the over/under-etched geometry is represented with a dashed line.

    Figure 7.Dependence of the bandwidth and center wavelength of the MMI on fabrication error δ. The bandwidth is defined as the wavelength range in which PDL, EL, and |IB| are simultaneously below 1 dB, and |PE| is less than 5°. The inset shows the effect on the in-plane geometry of positive and negative values of δ. The original geometry (top view, i.e., in the chip plane) is represented with a continuous line, while the over/under-etched geometry is represented with a dashed line.

    6. CONCLUSION

    We have proposed the first polarization-independent 2×2 MMI coupler in a 220 nm SOI platform. Our design leverages the concept of bricked SWG metamaterial waveguide engineering. Our optimized MMI device design features sub-decibel level of insertion loss, PDL, and IB, and PEs under 5° in the wavelength range of 1.26–1.42 μm. This is an unprecedented design performance for a 2×2 MMI in a standard 220 nm SOI platform. The 2D model and design methodology proposed in this work is general and indeed can be readily applied to MMIs with larger numbers of ports. We believe that our design strategy based on anisotropy engineering using bricked SWG metamaterial waveguides opens promising prospects for development of polarization-insensitive integrated photonic components and devices.

    APPENDIX A: INFLUENCE OF BRICKED SWG PITCH AND DUTY CYCLE ON MMI PERFORMANCE

    In Section 4, we fixed three geometrical parameters of the bricked SWG core to Λx=200  nm, and DC=DCx=DCz=50%. Here we justify this choice by considering different periods and DCs.

    Figure 8 shows the contour maps of |ΔLπ¯| as a function of the longitudinal pitch size (Λz) and bricked shifting (dz) for different values of DCs (DCx=DCz=DC=40%,50%,60%) and transversal periods (Λx=200nm,250nm,300nm). As described in Section 4, such contour maps can be obtained following the procedure of Fig. 3(a). Figure 8(d) corresponds to the same case simulated in Fig. 4, where the chosen design point is marked with a green dot (Λx=200  nm, Λz=240  nm, DC=50%, dz=Λz/2).

    Calculation of |ΔLπ¯| as performed in Fig. 4 for different transversal periods Λx and duty cycles (DC=DCx=DCz) and the wavelength of λ=1310 nm. The dashed line corresponds to |ΔLπ|=0. The design point is marked with a green dot in (d). For Λx>200 nm or DC≠50%, the design space for |ΔLπ¯| is closer to the Bragg zone [see cases (a), (e), and (g)] or no solution can be found, as shown in (b), (c), (f), (h), and (i).

    Figure 8.Calculation of |ΔLπ¯| as performed in Fig. 4 for different transversal periods Λx and duty cycles (DC=DCx=DCz) and the wavelength of λ=1310  nm. The dashed line corresponds to |ΔLπ|=0. The design point is marked with a green dot in (d). For Λx>200  nm or DC50%, the design space for |ΔLπ¯| is closer to the Bragg zone [see cases (a), (e), and (g)] or no solution can be found, as shown in (b), (c), (f), (h), and (i).

    An analysis of the results of Fig. 8 yields the following conclusions: (i) the Bragg zone grows for DC increasing above 50%, resulting in a reduction of the design space, (ii) the possible design solutions (where ΔLπ=0) move towards the Bragg zone when DC is reduced below 50%, (iii) although increasing the period relaxes minimum feature size constraints, it also shifts the polarization-independent designs closer to the Bragg zone [see Fig. 8(e)] or even makes them impossible [see Figs. 8(b), 8(c), 8(f), 8(h), and 8(i)].

    Thus, increasing the grating pitch or shifting DC away from 50% will degrade the bandwidth of the device due to the proximity of the Bragg zone and the corresponding increase in backreflections. In addition, the proximity of the Bragg zone potentially increases the influence of disorder (jitter) for a given fabrication process [41]. Designing periodic devices with feature sizes greater than 100 nm indeed is more challenging at the O-band (1260–1360 nm) than at the C-band (1530–1565 nm), as the shorter wavelength requires the use of shorter periods to keep the Bragg resonance effect negligible. Therefore, the chosen preliminary design represents a good trade-off between keeping minimum feature size larger than 100 nm and ensuring that the Bragg zone is far from the operation band.

    References

    [1] A. Rahim, T. Spuesens, R. Baets, W. Bogaerts. Open-access silicon photonics: current status and emerging initiatives. Proc. IEEE, 106, 2313-2330(2018).

    [2] E. Luan, H. Shoman, D. M. Ratner, K. C. Cheung, L. Chrostowski. Silicon photonic biosensors using label-free detection. Sensors, 18, 3519(2018).

    [3] K. Y. Yang, J. Skarda, M. Cotrufo, A. Dutt, G. H. Ahn, M. Sawaby, D. Vercruysse, A. Arbabian, S. Fan, A. Alù, J. Vučković. Inverse-designed non-reciprocal pulse router for chip-based LiDAR. Nat. Photonics, 14, 369-374(2020).

    [4] A. W. Elshaari, W. Pernice, K. Srinivasan, O. Benson, V. Zwiller. Hybrid integrated quantum photonic circuits. Nat. Photonics, 14, 285-298(2020).

    [5] D. Zhang, D. Liu, X. Wu, D. Nesset. Progress of ITU-T higher speed passive optical network (50G-PON) standardization. J. Opt. Commun. Netw., 12, D99-D108(2020).

    [6] F. Saliou, P. Chanclou, L. A. Neto, G. Simon, J. Potet, M. Gay, L. Bramerie, H. Debregeas. Optical access network interfaces for 5G and beyond [Invited]. J. Opt. Commun. Netw., 13, D32-D42(2021).

    [7] M. S. Erkılınç, D. Lavery, K. Shi, B. C. Thomsen, P. Bayvel, R. I. Killey, S. J. Savory. Polarization-insensitive single-balanced photodiode coherent receiver for long-reach WDM-PONs. J. Lightwave Technol., 34, 2034-2041(2016).

    [8] M. S. Erkılınç, D. Lavery, K. Shi, B. C. Thomsen, R. I. Killey, S. J. Savory, P. Bayvel. Comparison of low complexity coherent receivers for UDWDM-PONs (λ-to-the-user). J. Lightwave Technol., 36, 3453-3464(2018).

    [9] S. T. Lim, C. E. Png, E. A. Ong, Y. L. Ang. Single mode, polarization-independent submicron silicon waveguides based on geometrical adjustments. Opt. Express, 15, 11061-11072(2007).

    [10] T. Barwicz, M. R. Watts, M. A. Popović, P. T. Rakich, L. Socci, F. X. Kärtner, E. P. Ippen, H. I. Smith. Polarization-transparent microphotonic devices in the strong confinement limit. Nat. Photonics, 1, 57-60(2007).

    [11] D. Dai, L. Liu, S. Gao, D.-X. Xu, S. He. Polarization management for silicon photonic integrated circuits. Laser Photon. Rev., 7, 303-328(2013).

    [12] X. Chen, W. Liu, Y. Zhang, Y. Shi. Polarization-insensitive broadband 2 × 2 3 dB power splitter based on silicon-bent directional couplers. Opt. Lett., 42, 3738-3740(2017).

    [13] Y. Wang, L. Xu, H. Yun, M. Ma, A. Kumar, E. El-Fiky, R. Li, N. Abadíacalvo, L. Chrostowski, N. A. F. Jaeger, D. V. Plant. Polarization-independent mode-evolution-based coupler for the silicon-on-insulator platform. IEEE Photon. J., 10, 4900410(2018).

    [14] H. Li, W. Chen, P. Wang, S. Dai, Y. Liu, Q. Fu, J. Li, Y. Li, T. Dai, H. Yu, J. Yang. Compact and low-loss 1 × 3 polarization-insensitive optical power splitter using cascaded tapered silicon waveguides. Opt. Lett., 45, 5596-5599(2020).

    [15] D. González-Andrade, C. Lafforgue, E. Durán-Valdeiglesias, X. Le Roux, M. Berciano, E. Cassan, D. Marris-Morini, A. V. Velasco, P. Cheben, L. Vivien, C. Alonso-Ramos. Polarization- and wavelength-agnostic nanophotonic beam splitter. Sci. Rep., 9, 3604(2019).

    [16] D. Liu, D. Dai. Silicon-based polarization-insensitive optical filter with dual-gratings. Opt. Express, 27, 20704-20710(2019).

    [17] H. Yun, L. Chrostowski, N. A. F. Jaeger. Narrow-band, polarization-independent, transmission filter in a silicon-on-insulator strip waveguide. Opt. Lett., 44, 847-850(2019).

    [18] L. Soldano, E. Pennings. Optical multi-mode interference devices based on self-imaging: principles and applications. J. Lightwave Technol., 13, 615-627(1995).

    [19] H. Guan, Y. Ma, R. Shi, X. Zhu, R. Younce, Y. Chen, J. Roman, N. Ophir, Y. Liu, R. Ding, T. Baehr-Jones, K. Bergman, M. Hochberg. Compact and low loss 90°; optical hybrid on a silicon-on-insulator platform. Opt. Express, 25, 28957-28968(2017).

    [20] D. Dai, S. He. Optimization of ultracompact polarization-insensitive multimode interference couplers based on Si nanowire waveguides. IEEE Photon. Technol. Lett., 18, 2017-2019(2006).

    [21] S. Wang, D. Dai. Polarization-insensitive 2 × 2 thermo-optic Mach-Zehnder switch on silicon. Opt. Lett., 43, 2531-2534(2018).

    [22] S. Y. Siew, B. Li, F. Gao, H. Y. Zheng, W. Zhang, P. Guo, S. W. Xie, A. Song, B. Dong, L. W. Luo, C. Li, X. Luo, G.-Q. Lo. Review of silicon photonics technology and platform development. J. Lightwave Technol., 39, 4374-4389(2021).

    [23] P. Cheben, D.-X. Xu, S. Janz, A. Densmore. Subwavelength waveguide grating for mode conversion and light coupling in integrated optics. Opt. Express, 14, 4695-4702(2006).

    [24] P. Cheben, S. Janz, D.-X. Xu, B. Lamontagne, A. Delage, S. Tanev. A broad-band waveguide grating coupler with a subwavelength grating mirror. IEEE Photon. Technol. Lett., 18, 13-15(2006).

    [25] J. H. Schmid, P. Cheben, S. Janz, J. Lapointe, E. Post, D.-X. Xu. Gradient-index antireflective subwavelength structures for planar waveguide facets. Opt. Lett., 32, 1794-1796(2007).

    [26] J. H. Schmid, P. Cheben, S. Janz, J. Lapointe, E. Post, A. Delâge, A. Densmore, B. Lamontagne, P. Waldron, D.-X. Xu. Subwavelength grating structures in silicon-on-insulator waveguides. Adv. Opt. Technol., 2008, 685489(2008).

    [27] P. J. Bock, P. Cheben, J. H. Schmid, A. Delâge, D.-X. Xu, S. Janz, T. J. Hall. Sub-wavelength grating mode transformers in silicon slab waveguides. Opt. Express, 17, 19120-19133(2009).

    [28] R. Halir, A. Ortega-Moñux, D. Benedikovic, G. Z. Mashanovich, J. G. Wangüemert-Pérez, J. H. Schmid, I. Molina-Fernández, P. Cheben. Subwavelength-grating metamaterial structures for silicon photonic devices. Proc. IEEE, 106, 2144-2157(2018).

    [29] P. Cheben, R. Halir, J. H. Schmid, H. A. Atwater, D. R. Smith. Subwavelength integrated photonics. Nature, 560, 565-572(2018).

    [30] J. M. Luque-González, A. Sánchez-Postigo, A. Hadij-ElHouati, A. Ortega-Moñux, J. G. Wangüemert-Pérez, J. H. Schmid, P. Cheben, I. Molina-Fernández, R. Halir. A review of silicon subwavelength gratings: building break-through devices with anisotropic metamaterials. Nanophotonics, 10, 2765-2797(2021).

    [31] N. Yang, J. Xiao. A compact silicon-based polarization-independent power splitter using a three-guide directional coupler with subwavelength gratings. Opt. Commun., 459, 125095(2020).

    [32] H. Xie, J. Zheng, P. Xu, J. Yao, J. Whitehead, A. Majumdar. Ultra-compact subwavelength-grating-assisted polarization-independent directional coupler. IEEE Photon. Technol. Lett., 31, 1538-1541(2019).

    [33] X. Jiang, D. Dai. Compact, broadband, and polarization-insensitive 3 dB power splitter on silicon. Asia Communications and Photonics Conference (ACP), 1-2(2018).

    [34] T. Hao, A. Sánchez-Postigo, P. Cheben, A. Ortega-Moñux, W. N. Ye. Dual-band polarization-independent subwavelength grating coupler for wavelength demultiplexing. IEEE Photon. Technol. Lett., 32, 1163-1166(2020).

    [35] A. Herrero-Bermello, J. M. Luque-González, R. Halir, P. Cheben, A. Ortega-Moñux, I. Molina-Fernández, A. V. Velasco. Zero-birefringence silicon waveguides based on tilted subwavelength metamaterials. IEEE Photon. J., 11, 2700308(2019).

    [36] W. Zhong, J. Xiao. Ultracompact polarization-insensitive power splitter using subwavelength-grating-based MMI couplers on an SOI platform. Appl. Opt., 59, 1991-1997(2020).

    [37] J. M. Luque-González, A. Ortega-Moñux, R. Halir, J. H. Schmid, P. Cheben, Í. Molina-Fernández, J. G. Wangüemert-Pérez. Bricked subwavelength gratings: a tailorable on-chip metamaterial topology. Laser Photon. Rev., 15, 2000478(2021).

    [38] R. Halir, P. Cheben, J. M. Luque-González, J. D. Sarmiento-Merenguel, J. H. Schmid, G. Wangüemert-Pérez, D.-X. Xu, S. Wang, A. Ortega-Moñux, Í. Molina-Fernández. Ultra-broadband nanophotonic beamsplitter using an anisotropic sub-wavelength metamaterial. Laser Photon. Rev., 10, 1039-1046(2016).

    [39] S. G. Johnson, J. D. Joannopoulos. Block-iterative frequency-domain methods for Maxwell’s equations in a planewave basis. Opt. Express, 8, 173-190(2001).

    [40] A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. Joannopoulos, S. G. Johnson. Meep: a flexible free-software package for electromagnetic simulations by the FDTD method. Comput. Phys. Commun., 181, 687-702(2010).

    [41] A. Ortega-Moñux, J. Čtyroký, P. Cheben, J. H. Schmid, S. Wang, I. Molina-Fernández, R. Halir. Disorder effects in subwavelength grating metamaterial waveguides. Opt. Express, 25, 12222-12236(2017).

    Carlos Pérez-Armenta, Alejandro Ortega-Moñux, José Manuel Luque-González, Robert Halir, Pedro J. Reyes-Iglesias, Jens Schmid, Pavel Cheben, Íñigo Molina-Fernández, J. Gonzalo Wangüemert-Pérez, "Polarization-independent multimode interference coupler with anisotropy-engineered bricked metamaterial," Photonics Res. 10, A57 (2022)
    Download Citation