• Photonics Research
  • Vol. 8, Issue 12, 1843 (2020)
Shuyi Li1、†, Lifeng Cai1、†, Dingshan Gao1、*, Jianji Dong1, Jin Hou2, Chunyong Yang2, Shaoping Chen2, and Xinliang Zhang1
Author Affiliations
  • 1Wuhan National Laboratory for Optoelectronics & School of Optical and Electronic Information, Huazhong University of Science and Technology, Wuhan 430074, China
  • 2Hubei Key Laboratory of Intelligent Wireless Communications, College of Electronics and Information Engineering, South-Central University for Nationalities, Wuhan 430074, China
  • show less
    DOI: 10.1364/PRJ.403818 Cite this Article Set citation alerts
    Shuyi Li, Lifeng Cai, Dingshan Gao, Jianji Dong, Jin Hou, Chunyong Yang, Shaoping Chen, Xinliang Zhang, "Compact and broadband multimode waveguide bend by shape-optimizing with transformation optics," Photonics Res. 8, 1843 (2020) Copy Citation Text show less

    Abstract

    Multimode waveguide bend is one of the key components for realizing high-density mode-division multiplexing systems on chip. However, the reported multimode waveguide bends are either large, bandwidth-limited or fabrication-complicated, which hinders their applications in future high-density multimode photonic circuits. Here we propose a compact multimode waveguide bend supporting four TE modes simply by shape-optimizing with transformation optics. The shape of the waveguide is optimized in the virtual space with gradient distribution of the refractive index, so that the scattering loss and intermode cross talk are well suppressed. After conformal mapping back into the physical space, a compact (effective radius of 17 μm) multimode bending waveguide is obtained. Simulations show that the proposed multimode waveguide bend has little loss (<0.1 dB) and low cross talk (<-20 dB) throughout an ultrabroad wavelength range of 1.16–1.66 μm. We also fabricated the shape-optimized multimode bending waveguide on a silicon-on-insulator wafer. At 1550 nm wavelength, the measured excess losses for the four lowest-order TE modes are less than 0.6 dB, and the intermode cross talks are all below -17 dB. Our study paves the way for realizing high-density and large-scale multimode integrated optical circuits for optical interconnect.

    1. INTRODUCTION

    Mode-division multiplexing (MDM) can efficiently increase the communication capability on a photonic chip by utilizing several waveguide modes as communication channels [16]. And it is also compatible with other multiplexing methods, such as wavelength-division multiplexing (WDM) [79] and polarization-division multiplexing (PDM) [10]. For on-chip MDM photonic circuits, the multimode waveguide bend (MWB) is an essential component [11,12]. However, due to the mode mismatch between the straight waveguide and bent waveguide, the bending radius of a traditional multimode arc bend needs to be larger than 500 μm to support four TE modes [13]. Otherwise, it will cause serious loss and cross talk between the different waveguide modes. Therefore, designing a compact and high-performance MWB is still a big challenge.

    Previously, some design methods of MWB have been proposed. These methods can be classified into four types. The first is based on some gradient index profiles. In 2012, Gabrielli et al. proposed a novel MWB designed by transformation optics (TO) [14]. The MWB has a gradient index profile and theoretically supports three TE modes with the effective radius of 78.8 μm. But its fabrication needs complicated gray-scale lithography. In 2019, Wu et al. used shallowly etched subwavelength grating to realize a linearly varied gradient index along the radius direction [15]. It can support four TE modes with the effective radius of 20 μm. But it needs two electron-beam lithography (EBL) processes and precise etch depth control. Its bandwidth is also limited by the subwavelength grating. The second is to introduce the nonarc bending shape. In 2018, Jiang et al. presented that modified Euler curves could realize an adiabatic mode conversion between straight waveguide mode and bend waveguide mode because its curvature radius changes gradually [13]. It can support four TM modes with low loss and low cross talk, but the effective radius (45 μm) is large. The third is to use a mode converter to eliminate the mode mismatch between the straight and bent multimode waveguides. In 2017, Sun et al. used a step-tapered mode converter to convert two TE straight waveguide modes to corresponding bend waveguide modes [16]. Unfortunately, it can only support two waveguide modes. In 2018, Xu et al. used the straight waveguide with subwavelength-structured PMMA upper cladding as a mode converter for four TM modes [17]. The fabrication is complex, and the effective radius is 45.8 μm (including the length of the mode converters). The fourth is to introduce a total internal reflection (TIR) corner at the joint of two multimode straight waveguides [18]. It can support more than four modes, but needs a large waveguide width to guarantee similar TIR performances for different waveguide modes. For example, the TIR of two TE modes needs a 7 μm wide waveguide, which means a long multimode taper is required to connect the wide TIR corner to the narrow multimode waveguide. In summary, the reported MWBs are either large, bandwidth-limited, hard to extend to more waveguide modes, or fabrication-complicated.

    In this paper, we propose a compact and broadband MWB supporting four TE modes, which is simply designed by optimizing the bent waveguide boundary shape with TO. A traditional 90° arc-bend waveguide in the physical space is mapped to a distorted graded index waveguide in the virtual space. The boundary shape of the graded index waveguide is optimized to remove the scattering loss and intermode cross talk. When mapping back into the physical space, a multimode bending waveguide with uniform index and optimized boundary shape is formed. Unlike previous gradient index MWBs [14,15], our MWB avoids the introduction of complex gradient index distribution and can be easily fabricated in one-step lithography. Its effective radius is only 17 μm. The simulation results at 1550 nm show that the excess losses for TE0, TE1, TE2, and TE3 modes are 0.010, 0.010, 0.018, and 0.022 dB, respectively, and the intermode cross talks are all below 24  dB. The bandwidth (below 0.1 dB loss and 20  dB cross talk) is ultrawide from 1.16 to 1.66 μm, and the theoretical permitted width fabrication deviation of our MWB can be as large as ±50  nm. We also successfully fabricated this MWB on a commercial silicon-on-insulator (SOI) wafer with SiO2 upper cladding. For the experimental results at 1550 nm, the excess losses for TE0, TE1, TE2, and TE3 modes are 0.232, 0.327, 0.477, and 0.546 dB, respectively, and intermode cross talks are all below 17  dB.

    2. DESIGN AND SIMULATIONS

    TO is a powerful and flexible design method proposed by Pendry et al. [19] and Leonhardt [20], which uses spatial coordinate transformations to manipulate light propagation at will. Based on the TO theory, it often requires the permittivity ε and permeability μ to be anisotropic tensors [21]. But when the coordinate transformation is conformal mapping (CM), the ε and μ of the transformed device can be simplified to a gradient index distribution [22]. Many marvelous integrated optical devices based on CM have been designed and fabricated, such as invisibility carpets [23,24], optical Janus [25], MWBs [14], and multimode waveguide crossings [2628]. Generally, they need to be transformed from the original devices. Thanks to the form invariance of the Maxwell equations, the transformed device will inherit the same optical performance from the original devices. Obviously, a perfect MWB with graded index distribution can be transformed from a straight multimode waveguide [14]. However, its gradient index will lead to tough fabrication difficulty, especially when the bending radius is small.

    For better understanding of the multimode propagation in conventional waveguide arc-bend, we transform the arc-bend back to a rectangle by CM in reverse. The coordinates (x,y) of the arc-bend are commonly called physical space, while the coordinates (u,v) of the rectangle are called virtual space, which are shown in Figs. 1(a) and 1(b). The arc-bend in physical space includes both the waveguide core and cladding areas, which is composed by two additional rectangles with length La=8  μm and width Ws=6  μm and a 90° arc with radius R=9  μm, while the rectangle in virtual space has length Lt and width Ws=6  μm.

    CM from the arc-bend in the physical space to a rectangle in the virtual space. (a) The arc-bend is composed by two additional rectangles with length La=8 μm and width Ws=6 μm and a 90° arc with radius R=9 μm. (b) The rectangle in the virtual space has length Lt and width Ws=6 μm. Lt is optimized to ensure u and v satisfy the Cauchy–Riemann condition of CM. The optimal Lt is 29.683 μm. All the grids represent coordinates (u,v) in subplots (a) and (b).

    Figure 1.CM from the arc-bend in the physical space to a rectangle in the virtual space. (a) The arc-bend is composed by two additional rectangles with length La=8  μm and width Ws=6  μm and a 90° arc with radius R=9  μm. (b) The rectangle in the virtual space has length Lt and width Ws=6  μm. Lt is optimized to ensure u and v satisfy the Cauchy–Riemann condition of CM. The optimal Lt is 29.683 μm. All the grids represent coordinates (u,v) in subplots (a) and (b).

    The CM from physical space to virtual space can be easily obtained by solving Laplace’s equation on the arc-bend in (x,y) coordinates, with Dirichlet and Neuman boundary conditions [29,30]. We set the inner bend curve as v=Ws and outer bend curve as v=0, respectively, which correspond to the upper and bottom boundaries of the rectangle in Fig. 1(b). And we set the boundary conditions u=Lt/2 and u=Lt/2 for the left and upper straight boundaries of the arc-bend, respectively, which correspond to the left and right boundaries of the rectangle in Fig. 1(b). To ensure u and v satisfy the Cauchy–Riemann condition of CM, the parameter Lt needs to be optimized. The ratio between u/x and v/y is the optimization objective, which should get close to 1. The Laplace equation-solving and Lt optimization are realized by combining the Laplace and optimization modules in COMSOL Multiphysics software (http://www.comsol.com). Here the optimal Lt is 29.683 μm.

    Figure 2(a) shows the refractive index distribution n of the arc-bend in the physical space. The width of waveguide is 2 μm, which can support the four lowest-order TE modes. The waveguide core index is 2.84, which is the TE0 mode effective index of 220 nm thick silicon slab waveguide with SiO2 upper cladding, while the index distribution n of the rectangle in the virtual space is shown in Fig. 2(b), which is calculated by the equation [22] n=n/det(P),where P is the Jacobian matrix of CM. From Fig. 2(b), we can find that the original waveguide in virtual space is not straight and has a gradient index distribution. The slight bending of the original waveguide is caused by the direct connection between the waveguide arc-bend and the straight waveguide, where the curvature suddenly changes from 0 to 1/R at their interfaces. In contrast, the curvature of CM grids changes gradually, since the conformal coordinates (u,v) are second-order derivable. And the gradient index profile along the cross section of the original waveguide will cause wavefront propagating inconsistently between the inner and outer sides of the arc-bend. These two negative factors cause the big loss and intermode cross talk of the traditional multimode waveguide arc-bend.

    (a) Index distribution of the waveguide arc-bend in physical space. The width of waveguide arc-bend is 2 μm, which can support four TE modes. The waveguide core index is 2.84, which is the TE0 mode effective index of a 220 nm thick silicon slab waveguide with SiO2 upper cladding. (b) Index distribution of the original waveguide in virtual space. It is obtained by Eq. (1) of CM. This original waveguide is not straight and has a gradient index distribution. These two negative factors cause the big loss and intermode cross talk of the traditional multimode waveguide arc-bend.

    Figure 2.(a) Index distribution of the waveguide arc-bend in physical space. The width of waveguide arc-bend is 2 μm, which can support four TE modes. The waveguide core index is 2.84, which is the TE0 mode effective index of a 220 nm thick silicon slab waveguide with SiO2 upper cladding. (b) Index distribution of the original waveguide in virtual space. It is obtained by Eq. (1) of CM. This original waveguide is not straight and has a gradient index distribution. These two negative factors cause the big loss and intermode cross talk of the traditional multimode waveguide arc-bend.

    According to the analysis of the waveguide arc-bend above, we propose to modify and optimize the original waveguide for multimode propagation in virtual space. This can be simply realized by the waveguide shape-optimizing. Different from previous TO works [14], the optimization and simulation are all implemented in virtual space. Then after a CM, the corresponding MWB in the physical space will be only shape-optimized but without introducing any gradient index distribution. Considering the mirror symmetry along u=0 of the original waveguide in Fig. 2(b), the waveguide boundary curves can be expressed in a Fourier cosine series. The expressions of the upper and lower boundary curves v1 and v2 are v1=a0+a1cos(wu)+a2cos(2wu)+a3cos(3wu),v2=b0+b1cos(wu)+b2cos(2wu)+b3cos(3wu),where w=2π/Lt. If more cosine terms are introduced, better performance will be obtained, but the optimization will take longer. Here we introduce three cosine terms to achieve a balance between performance and optimization time. The parameters a1, a2, a3, b1, b2, and b3 are the optimization variables. Since the MWB needs to connect smoothly to the input and output multimode straight waveguides, a0 and b0 can be determined by the six optimization variables and the waveguide width. We set the transmission efficiencies of TE0TE3 modes as the optimization objectives and the intermode cross talks as the optimization constraints. We use the “Bound Optimization BY Quadratic Approximation” (BOBYQA) method in the COMSOL optimization module and couple it with the 2D wave-optics module. Once the optimization variables vary, the program will input TE0TE3 modes in turn. Then the transmission efficiencies of every mode and intermode cross talks will be calculated. The optimization module will estimate the next values of the optimization variables based on each mode transmission efficiency. But if one intermode cross talk is higher than the constraint of 20  dB, the optimization variables of this cycle will be abandoned. When the optimization objectives are nearly achieved within a predefined threshold, the optimization process converges and stops. The optimal parameters are a0=0.946087, a1=0.11757, a2=0.079818, a3=0.021599, b0=0.7622125, b1=0.19865, b2=0.0042085, and b3=0.0042154. Figure 3(a) shows the index distribution of the optimal original waveguide in virtual space, while Fig. 3(b) shows the index distribution of the corresponding MWB in physical space. This MWB has a uniform core index, which makes fabrication easier. Different from the modified Euler curves in Ref. [13], the waveguide width of our optimized MWB also changes along the waveguide bend, and the effective radius of our MWB is 17 μm, which is much smaller than that of the modified Euler curves of 45 μm.

    (a) Index distribution of the shape optimized original waveguide in virtual space; (b) index distribution of the MWB in physical space. This MWB has a uniform core index, which makes fabrication easier.

    Figure 3.(a) Index distribution of the shape optimized original waveguide in virtual space; (b) index distribution of the MWB in physical space. This MWB has a uniform core index, which makes fabrication easier.

    The multimode propagation performance of the shape-optimized original waveguide in virtual space is simulated by the 2D wave-optics module in COMSOL and presented in Fig. 4. Figures 4(a)–4(d) show the profiles of the Hz field component for TE0, TE1, TE2, and TE3 modes at 1550 nm wavelength, respectively. All the modes transmit adiabatically in the optimized original waveguide, and no obvious scattering occurs. The simulation results of the shape-optimized MWB in virtual space at 1550 nm wavelength show that the excess losses for the four lowest-order TE modes are less than 0.01 dB and the intermode cross talks are all below 25  dB.

    Multimode propagation performance of the shape-optimized original waveguide in virtual space. (a)–(d) The profiles of the Hz field component for (a) TE0, (b) TE1, (c) TE2, and (d) TE3 modes at 1550 nm wavelength, respectively, are shown.

    Figure 4.Multimode propagation performance of the shape-optimized original waveguide in virtual space. (a)–(d) The profiles of the Hz field component for (a) TE0, (b) TE1, (c) TE2, and (d) TE3 modes at 1550 nm wavelength, respectively, are shown.

    We also use rigorous 3D finite-difference time-domain (FDTD) simulations to verify the performance of our shape-optimized MWB in physical space. Figures 5(a)–5(d) show the profiles of the Hz field component for TE0, TE1, TE2, and TE3 modes at 1550 nm wavelength, respectively. The insets show the input/output Hz field profiles. Meanwhile, we have simulated the transmission spectra of the four TE modes in the MWB in the wavelength range of 1.1–2.0 μm. To reduce the simulation cost, we divide this broad wavelength range into nine equal subranges to simulate them separately, each with a bandwidth of 100 nm. The material dispersions of silicon and silica have been taken into account in the simulation. We calculate the wavelength-dependence of the coupled power Tij from the ith TE mode at the input waveguide to the jth TE mode at the output waveguide. Here the excess loss for each mode is given by ELii=10log10Tii, while the intermode cross talk is calculated by CTij=10log10Tij(ji). To guarantee adequate simulation accuracy, we set the FDTD mesh sizes as dx=25  nm, dy=25  nm, and dz=22  nm. The theoretical transmission spectra of the MWB for TE0, TE1, TE2, and TE3 modes are shown in Figs. 6(a)–6(d), respectively. The simulation results at 1550 nm show that excess losses for TE0, TE1, TE2, and TE3 modes are only 0.010, 0.010, 0.018, and 0.022 dB, respectively, and intermode cross talks are all below 24  dB. From these diagrams, we can find that throughout a broad wavelength range of 1.16–1.66 μm, the losses of all modes are below 0.1 dB and the cross talks are all below 20  dB.

    3D FDTD simulation results of shape-optimized MWB in physical space. (a)–(d) The profiles of the Hz field component for (a) TE0, (b) TE1, (c) TE2, and (d) TE3 modes at 1550 nm wavelength, respectively, are shown. The insets show the input/output Hz field profiles.

    Figure 5.3D FDTD simulation results of shape-optimized MWB in physical space. (a)–(d) The profiles of the Hz field component for (a) TE0, (b) TE1, (c) TE2, and (d) TE3 modes at 1550 nm wavelength, respectively, are shown. The insets show the input/output Hz field profiles.

    Simulated transmission spectra of shape-optimized multimode bend waveguide. (a)–(d) The spectra for (a) TE0, (b) TE1, (c) TE2, and (d) TE3 modes, respectively, are shown. At 1550 nm, the excess losses for TE0, TE1, TE2, and TE3 modes are 0.010, 0.010, 0.018, and 0.022 dB, respectively, and intermode cross talks are all below −24 dB. Throughout a broad wavelength range of 1.16–1.66 μm, the losses of all modes are below 0.1 dB and the cross talks are all below −20 dB.

    Figure 6.Simulated transmission spectra of shape-optimized multimode bend waveguide. (a)–(d) The spectra for (a) TE0, (b) TE1, (c) TE2, and (d) TE3 modes, respectively, are shown. At 1550 nm, the excess losses for TE0, TE1, TE2, and TE3 modes are 0.010, 0.010, 0.018, and 0.022 dB, respectively, and intermode cross talks are all below 24  dB. Throughout a broad wavelength range of 1.16–1.66 μm, the losses of all modes are below 0.1 dB and the cross talks are all below 20  dB.

    To evaluate the permitted width fabrication deviation, we introduce a width deviation Δw of the MWB in physical space by moving each point on the boundary shape curve along the outer normal vector at a distance of Δw/2. The value of Δw varies from 100 to 140 nm with a 20 nm interval. The simulated transmission efficiencies of the shape-optimized MWB at 1550 nm wavelength vary with Δw, as shown in Fig. 7. When the required intermode cross talk is less than 20  dB, the theoretical permitted width fabrication deviation of our MWB is ±50  nm. This proves that our MWB has large fabrication tolerance.

    Simulated transmission efficiencies of shape-optimized MWB at 1550 nm wavelength vary with the width fabrication deviation Δw. (a)–(d) The efficiencies for (a) TE0, (b) TE1, (c) TE2, and (d) TE3 modes, respectively, are shown. For the intermode cross talk below −20 dB, the theoretical permitted width fabrication deviation of our MWB is ±50 nm.

    Figure 7.Simulated transmission efficiencies of shape-optimized MWB at 1550 nm wavelength vary with the width fabrication deviation Δw. (a)–(d) The efficiencies for (a) TE0, (b) TE1, (c) TE2, and (d) TE3 modes, respectively, are shown. For the intermode cross talk below 20  dB, the theoretical permitted width fabrication deviation of our MWB is ±50  nm.

    3. FABRICATION AND MEASUREMENT

    The proposed shape-optimized MWB was fabricated on an SOI wafer with 220 nm thick top silicon layer and 2 μm thick buried oxide (BOX) layer. The EBL and inductively coupled plasma (ICP) dry etching were applied to fabricate the photonics-integrated circuit (PIC). The micrographs of the fabricated PICs with and without MWB are shown in Figs. 8(a) and 8(b), respectively. A magnified view of MWB is shown in Fig. 8(c). The mode multiplexer and demultiplexer are designed for multiplexing and demultiplexing of the four TE modes (TE0, TE1, TE2, and TE3) [4].

    Microscopic view of the fabricated PICs and device. (a) The referenced PIC without MWB, including grating couplers, mode multiplexers, and demultiplexers; (b) PIC with MWB; (c) magnified microscopic view of the shape-optimized MWB.

    Figure 8.Microscopic view of the fabricated PICs and device. (a) The referenced PIC without MWB, including grating couplers, mode multiplexers, and demultiplexers; (b) PIC with MWB; (c) magnified microscopic view of the shape-optimized MWB.

    The mode multiplexer and demultiplexer are both formed by the asymmetry directional coupler (ADC) [31,32] between a narrow single-mode waveguide and a wide multimode bus waveguide. The ADC converts the TE0 mode of the narrow single-mode waveguide to one higher-order TE mode of the wide multimode bus waveguide, and vice versa. These ADCs are cascaded to realize MDM of the four TE modes. The bus waveguides with different widths are connected through adiabatic tapers, which are long enough to avoid any intermode coupling or radiation loss. Then the bus waveguide is expanded to the final width of 2 μm to access the two cascaded 90° MWBs. After the two MWBs, a mirrored circuit is used to demultiplex different waveguide modes back to TE0 mode in different output channels. Grating couplers are used to couple light in and out of the PIC from a single-mode fiber.

    A wideband amplified spontaneous emission (ASE) light source covering the C and L bands is used in our measurement. After the output fiber, a commercial 3 dB Y splitter is connected to split the output light into two paths. One is linked to the optical power meter for monitoring and tuning the alignment between the single-mode fibers and the grating couplers. The other is connected to the optical spectrum analyzer (OSA) to measure the transmission spectrum for every mode. To obtain the transmission efficiency of the MWB accurately, the losses of the vertical grating couplers, mode multiplexer, and demultiplexer should be deducted from the measured spectrum of the PIC with MWB. This can be easily done by measuring and subtracting the transmission spectra of the PIC with and without MWB.

    The measured transmission spectra for each mode of the shape optimized MWB are shown in Fig. 9. At central wavelength of 1550 nm, the excess losses for TE0, TE1, TE2, and TE3 modes are 0.232, 0.327, 0.477, and 0.546 dB, respectively, and intermode cross talks are all below 17  dB. Since the ASE light source used in the measurement has limited bandwidth of 1520–1620 nm and the 3 dB bandwidth of the grating coupler is only 50  nm, we could not measure the spectrum of the MWB in the whole wavelength range of 1200–1700 nm. Nevertheless, the transmission spectra for the four TE modes are very flat across the measured wavelength range, which demonstrates our shape-optimized MWB has wide operation bandwidth. The measured transmission efficiencies are slightly lower than the simulated results, which is possibly caused by the waveguide sidewall etching roughness and width random deviation in the fabrication.

    Measured transmission spectra of the shape-optimized MWB for (a) TE0, (b) TE1, (c) TE2, and (d) TE3 modes, respectively. At 1550 nm wavelength, the excess losses for TE0, TE1, TE2, and TE3 modes are 0.232, 0.327, 0.477, and 0.546 dB, respectively, and intermode cross talks are all below −17 dB.

    Figure 9.Measured transmission spectra of the shape-optimized MWB for (a) TE0, (b) TE1, (c) TE2, and (d) TE3 modes, respectively. At 1550 nm wavelength, the excess losses for TE0, TE1, TE2, and TE3 modes are 0.232, 0.327, 0.477, and 0.546 dB, respectively, and intermode cross talks are all below 17  dB.

    Table 1 compares the theorical performance of the reported MWBs on a silicon platform. Our MWB demonstrated in this paper has lowest excess losses, widest bandwidth, and sharpest bending among all the four modes channel MWBs. By slightly enlarging the bending radius and the waveguide width, our TO-optimized MWB can also be extended for TM modes or more TE modes. Furthermore, our MWB can be fabricated by one-step lithography and the dry-etching process with large fabrication tolerance.

    Ref.MethodMR (μm)EL (dB)CT (dB)BW (nm)
    [14]TO378.82.5//
    [15]SWGs3100.5−30100
    4200.5−26100
    [13]Euler curves4450.1−25100
    [16]PSO250.12−23100
    [17]SWG tapers445.80.4−20100
    [18]Corner bend270.18−36420
    10350.54−24420
    This workTO-optimized4170.1−20500

    Table 1. Comparison of Reported MWBs on Silicon Platform (Theoretical Performance)a

    4. CONCLUSION

    In this paper, we demonstrate a compact and broadband MWB supporting four TE modes simply by optimizing the waveguide boundary shape with TO. By designing and optimizing the MWB boundary shapes in virtual space of TO, the loss and intermode cross talk caused by the gradient index distribution are well suppressed. When transformed back to the physical space, the shape-optimized MWB can be formed. The theoretical simulation results show that this shape-optimized MWB exhibits low loss (<0.1  dB) and low cross talk (<20  dB) across an ultrabroad bandwidth from 1.16 to 1.66 μm, and the permitted waveguide width fabrication deviation is as large as ±50  nm. The proposed MWB was fabricated on a commercial SOI wafer. The measured excess losses for the four lowest-order TE modes are less than 0.6 dB, and the intermode cross talks are all below 17  dB. This shape-optimized MWB offers great potential to construct densely integrated and wideband MDM systems for on-chip ultrahigh capacity communications. Our proposed TO design method can be expanded to combine with more advanced optimization methods, such as genetic algorithm, particle swarm optimization, and topology optimization. Furthermore, in the future, this method will inspire more integrated TO devices without complicated fabrication.

    Acknowledgment

    Acknowledgment. We thank Dr. Cheng Zeng and the engineer Mr. Pan Li in the Center of Micro-Fabrication and Characterization (CMFC) of WNLO for the support in EBL and ICP dry etching.

    References

    [1] D. J. Richardson, J. M. Fini, L. E. Nelson. Space-division multiplexing in optical fibres. Nat. Photonics, 7, 354-362(2013).

    [2] G. F. Li, N. Bai, N. B. Zhao, C. Xia. Space-division multiplexing: the next frontier in optical communication. Adv. Opt. Photon., 6, 413-487(2014).

    [3] R. Ryf, S. Randel, A. H. Gnauck, C. Bolle, A. Sierra, S. Mumtaz, M. Esmaeelpour, E. C. Burrows, R. J. Essiambre, P. J. Winzer, D. W. Peckham, A. H. McCurdy, R. Lingle. Mode-division multiplexing over 96  km of few-mode fiber using coherent 6 × 6 MIMO processing. J. Lightwave Technol., 30, 521-531(2012).

    [4] D. Dai, J. Wang, Y. Shi. Silicon mode (de)multiplexer enabling high capacity photonic networks-on-chip with a single-wavelength-carrier light. Opt. Lett., 38, 1422-1424(2013).

    [5] G. Chen, Y. Yu, X. Zhang. Monolithically mode division multiplexing photonic integrated circuit for large-capacity optical interconnection. Opt. Lett., 41, 3543-3546(2016).

    [6] S. Wang, H. Wu, H. K. Tsang, D. Dai. Monolithically integrated reconfigurable add-drop multiplexer for mode-division-multiplexing systems. Opt. Lett., 41, 5298-5301(2016).

    [7] L. W. Luo, N. Ophir, C. P. Chen, L. H. Gabrielli, C. B. Poitras, K. Bergmen, M. Lipson. WDM-compatible mode-division multiplexing on a silicon chip. Nat. Commun., 5, 3069(2014).

    [8] S. Wang, X. Feng, S. Gao, Y. Shi, T. Dai, H. Yu, H. K. Tsang, D. Dai. On-chip reconfigurable optical add-drop multiplexer for hybrid wavelength/mode-division-multiplexing systems. Opt. Lett., 42, 2802-2805(2017).

    [9] D. X. Dai, J. Wang, S. T. Chen, S. P. Wang, S. L. He. Monolithically integrated 64-channel silicon hybrid demultiplexer enabling simultaneous wavelength- and mode-division-multiplexing. Laser Photon. Rev., 9, 339-344(2015).

    [10] J. Wang, S. L. He, D. X. Dai. On-chip silicon 8-channel hybrid (de)multiplexer enabling simultaneous mode-and polarization-division-multiplexing. Laser Photon. Rev., 8, L18-L22(2014).

    [11] C. L. Li, D. J. Liu, D. X. Dai. Multimode silicon photonics. Nanophotonics, 8, 227-247(2019).

    [12] D. X. Dai, J. E. Bowers. Silicon-based on-chip multiplexing technologies and devices for Peta-bit optical interconnects. Nanophotonics, 3, 283-311(2014).

    [13] X. Jiang, H. Wu, D. Dai. Low-loss and low-crosstalk multimode waveguide bend on silicon. Opt. Express, 26, 17680-17689(2018).

    [14] L. H. Gabrielli, D. Liu, S. G. Johnson, M. Lipson. On-chip transformation optics for multimode waveguide bends. Nat. Commun., 3, 1217(2012).

    [15] H. Wu, C. L. Li, L. J. Song, H. K. Tsang, J. E. Bowers, D. X. Dai. Ultra-sharp multimode waveguide bends with subwavelength gratings. Laser Photon. Rev., 13, 1800119(2019).

    [16] C. Sun, Y. Yu, G. Chen, X. Zhang. Ultra-compact bent multimode silicon waveguide with ultralow inter-mode crosstalk. Opt. Lett., 42, 3004-3007(2017).

    [17] H. N. Xu, Y. C. Shi. Ultra-sharp multi-mode waveguide bending assisted with metamaterial-based mode converters. Laser Photon. Rev., 12, 1700240(2018).

    [18] Y. Wang, D. Dai. Multimode silicon photonic waveguide corner-bend. Opt. Express, 28, 9062-9071(2020).

    [19] J. B. Pendry, D. Schurig, D. R. Smith. Controlling electromagnetic fields. Science, 312, 1780-1782(2006).

    [20] U. Leonhardt. Optical conformal mapping. Science, 312, 1777-1780(2006).

    [21] H. Chen, C. T. Chan, P. Sheng. Transformation optics and metamaterials. Nat. Mater., 9, 387-396(2010).

    [22] L. Xu, H. Chen. Conformal transformation optics. Nat. Photonics, 9, 15-23(2014).

    [23] L. H. Gabrielli, J. Cardenas, C. B. Poitras, M. Lipson. Silicon nanostructure cloak operating at optical frequencies. Nat. Photonics, 3, 461-463(2009).

    [24] J. Valentine, J. Li, T. Zentgraf, G. Bartal, X. Zhang. An optical cloak made of dielectrics. Nat. Mater., 8, 568-571(2009).

    [25] T. Zentgraf, J. Valentine, N. Tapia, J. Li, X. Zhang. An optical ‘Janus’ device for integrated photonics. Adv. Mater., 22, 2561-2564(2010).

    [26] S. Y. Li, Y. Y. Zhou, J. J. Dong, X. L. Zhang, E. Cassan, J. Hou, C. Y. Yang, S. P. Chen, D. S. Gao, H. Y. Chen. Universal multimode waveguide crossing based on transformation optics. Optica, 5, 1549-1556(2018).

    [27] Q. Wu, J. P. Turpin, D. H. Werner. Integrated photonic systems based on transformation optics enabled gradient index devices. Light: Sci. Appl., 1, e38(2012).

    [28] P. Markov, J. G. Valentine, S. M. Weiss. Fiber-to-chip coupler designed using an optical transformation. Opt. Express, 20, 14705-14713(2012).

    [29] Z. Chang, X. M. Zhou, J. Hu, G. K. Hu. Design method for quasi-isotropic transformation materials based on inverse Laplace’s equation with sliding boundaries. Opt. Express, 18, 6089-6096(2010).

    [30] Y. G. Ma, N. Wang, C. K. Ong. Application of inverse, strict conformal transformation to design waveguide devices. J. Opt. Soc. Am. A, 27, 968-972(2010).

    [31] T. H. Pan, S. Y. Tseng. Short and robust silicon mode (de)multiplexers using shortcuts to adiabaticity. Opt. Express, 23, 10405-10412(2015).

    [32] D. Guo, T. Chu. Silicon mode (de)multiplexers with parameters optimized using shortcuts to adiabaticity. Opt. Express, 25, 9160-9170(2017).

    Shuyi Li, Lifeng Cai, Dingshan Gao, Jianji Dong, Jin Hou, Chunyong Yang, Shaoping Chen, Xinliang Zhang, "Compact and broadband multimode waveguide bend by shape-optimizing with transformation optics," Photonics Res. 8, 1843 (2020)
    Download Citation