1. Introduction
Silicon-on-insulator (SOI), a mature and promising platform for silicon photonics, has been commonly employed to realize compact, high-performance, and high-yield photonic integrated circuits (PICs) [
1,
2,
3,
4,
5]. To construct such on-chip PICs, power splitters are fundamental and indispensable components, which play a key role in separating the input light power into different output ports [
6,
7,
8,
9,
10]. For the power splitter, different waveguide structures have been reported, e.g., typical Y-junctions [
6], adiabatic tapers [
7], directional couplers (DCs) [
8], photonic crystals (PCs) [
9], and multimode interference (MMI) couplers [
10]. With overall consideration of the device size, performance, and fabrication tolerance, MMI couplers are the best choice since Y-branches require a large device size, while DCs and PCs have limited working bandwidths and tight fabrication tolerances, respectively. Therefore, MMI-based power splitters are vital components applied in on-chip Mach–Zehnder interferometer (MZI) modulators [
11,
12], optical switches [
13,
14], and other on-chip optical circuits requiring a light power-splitting function [
15,
16]. Within these application cases, a smaller device size and better device performance have become the new requirements for silicon-based power splitters.
To further shorten the device size of MMI-based power splitters (e.g., typical MMI length of ~5.2 μm under a waveguide width of 2.5 μm [
17,
18]), we should enlarge the effective index difference of the transmitted interference modes within the MMI region, which needs to change the refractive index distribution of the conventional MMI structure. Fortunately, a subwavelength grating (SWG) structure whose grating pitch is obviously shorter than the grating Bragg period behaves as a homogenous medium without a reflection and diffraction effect for the input light (e.g., λ = 1.55 μm) [
19,
20,
21,
22]. Furthermore, the effective refractive index of SWG structure can be easily changed by varying its grating duty cycle, offering a new degree of freedom for the device design [
19,
20,
21,
22]. Therefore, by adding such an SWG structure into the conventional MMI waveguide, the refractive index distribution of the whole structure becomes no longer uniform in the waveguide region, contributing to a reduction in MMI length [
19]. Under this condition, various device schemes have been proposed. For instance, through embedding a row of SWGs in the central region of the MMI, the required MMI length was clearly reduced to 3.8 μm (waveguide width 2.2 μm) for power splitting [
23]. Meanwhile, a dual polarization operation could also be achieved with the help of the tunability of the effective refractive index for an SWG structure [
23]. To further shrink the MMI length, a tapered SWG structure was embedded in the central region of the MMI along the propagation direction, and the MMI length was reduced to only 1.92 μm (waveguide width 2.0 μm) with an insertion loss (IL) of 0.39 dB, which is quite beneficial for on-chip compact integration [
24]. In addition to these schemes, an SWG-based MMI power splitter was shown to be a useful structure, where the MMI region was formed by SWGs with identical or different orientations for the grating component [
25]. The corresponding MMI lengths could even be reduced to 3.2 μm (waveguide width 2.8 μm) with a large bandwidth of 400 nm (IL < 1 dB) [
25]. However, we should note that the IL and working bandwidth of the power splitter still need to be efficiently reduced and increased, respectively, since the fundamental power splitter will be greatly required and heavily used in current and future on-chip large-scale PICs for splitting light power [
26,
27,
28].
In this paper, we propose a silicon-based MMI power splitter, where the key MMI region is embedded by an inverse-tapered SWG in the center and two rows of uniform SWG slots on both sides along the propagation direction. The embedded central SWG is employed to change the whole refractive index distribution of the multimode silicon waveguide, which can help enlarge the effective index difference between excited modes and further reduce the MMI length. The embedded bilateral SWG slots are used to increase the coupling efficiency with the output waveguides via mode matching, corresponding to a reduction in power splitting loss. Through these two main techniques applied in the MMI structure for the power splitter, the key MMI length could be reduced to 3.2 μm (waveguide width 2.5 μm), while the obtained IL was only 0.08 dB at a wavelength of 1.55 μm, along with a quite low reflection loss (RL) of <−35 dB. Note that the device working bandwidth could be over 550 nm (covering the whole optical communication band) even under a strict criterion of IL < 0.6 dB. In addition, the fabrication tolerances of the key structural parameters were also analyzed. We hope that the proposed device can find important applications in the field of on-chip PICs.
2. Device Structure and Principle
Figure 1 illustrates the three-dimensional schematic of our proposed SWG-assisted MMI power splitter, which consisted of an inverse-tapered SWG and two rows of uniform SWG slots embedded in the MMI center and on both sides, respectively. For a conventional MMI power splitter, as the input fundamental TE
0 mode is injected into the MMI region from the central port, the new TE
0 and TE
2 mode are excited dominantly within the MMI region, thus leading to a clear mode interference along the propagation direction. Once the accumulated phase difference between these two excited modes (TE
0 and TE
2 mode) is equal to π, a double image of the input TE
0 mode with the same phase can be achieved after mode interference [
17]. Finally, using two waveguides to separate and output such a double image, we can realize the function of power splitting. To efficiently shorten the MMI length and reduce the power splitting loss, we embedded an inverse-tapered SWG in the MMI center, where the width of the SWG was linearly tapered from
wg1 (=200 nm) to
wgn in a period number of
n, and we embedded two rows of uniform SWG slots (period number N) on both sides of the MMI region, whereby their central positions along the propagation direction were aligned with those of the output waveguides. These embedded SWG structures had the same grating pitch and duty cycle of Λ = 200 nm and
a/Λ = 0.5, respectively. The gap length between the left side of the inverse-tapered SWG (uniform SWG slot) and the left side of the MMI region was
LS (
LE), shown in the inset of
Figure 1. Moreover, the length and width of MMI region were
LM and
WM (=2.5 μm), respectively. The widths of the central input waveguide and bilateral output waveguides were tapered from
w1 to
w2 in a length of
LI and from
w3 to
w4 in a length of
LO, respectively, where the input and output waveguide widths were the same (
w1 =
w4 = 500 nm). In addition, we also introduced two right-angled cutting structures on both sides of the input port of MMI region symmetrically to further enhance the device performance, where the corresponding length and width of the cutting structures were
LT and (
WM −
w2)/2, respectively. The whole device was analyzed and simulated on the basis of a commercial SOI wafer with a 220 nm thick top silicon layer and a 2 μm thick buried oxide layer (SiO
2), while the upper cladding was also SiO
2 with a thickness of 2 μm.
Figure 2a shows the calculated effective indices
neff of guided modes under different waveguide widths, where the TE polarization modes were analyzed. From this figure, we can easily find the effective index of every supported mode for the waveguide width. For example, the silicon waveguide could support seven modes (from TE
0 to TE
6) when the waveguide width was 2.5 μm. For the multimode waveguide transmission, only even modes (TE
0, TE
2, TE
4, TE
6) could be adequately excited for the central injection from the input single-mode waveguide [
17]. Meanwhile, the TE
0 and TE
2 modes of the multimode waveguide would dominate among the excited modes according to our calculations. On the basis of these two modes, we could roughly determine the required multimode waveguide length for the power splitting, where more precise length should be optimized using a numerical method. By adding SWG structures into the conventional MMI waveguide, the effective index difference of dominant excited modes (TE
0 and TE
2 mode) within the MMI region was increased, leading to a reduction in beat length
Lπ [
17,
18,
19].
where
nTE0,
nTE2 are the mode effective indices of the excited TE
0 and TE
2 modes, and
λ is the working wavelength. If we set the MMI length
LM equal to the beat length
Lπ, we can obtain two TE
0 modes (double image of the input TE
0 mode) with the same phase and same power, corresponding to the power splitting. The SWG structure can be treated as a homogeneous medium, and its equivalent refractive index
nSWG can be estimated as follows [
19]:
where
nc and
ncl denote the refractive index of the core (silicon) and cladding (silica) of the grating, respectively. Using Equation (2), we can roughly obtain the equivalent refractive index of the SWG structure.
Figure 2b shows the calculated effective indices
neff of guided modes for the multimode waveguide (
WM = 2.5 μm) embedded with the central SWG structure as a function of the SWG width
wg. Note that the effective index difference between dominant excited TE
0 and TE
2 modes is really increased with the increase in embedded SWG width
wg. This is the reason that the required MMI length can be reduced compared with the conventional MMI structure. In the analyses below, we study the embedded SWG structures in detail to achieve larger bandwidth and lower loss compared with previous reports [
23,
24,
25].
3. Results and Discussion
For the device performance analyses, the three-dimensional finite-difference time-domain (3D-FDTD) method was employed to help design and optimize the device parameters [
29,
30].
Figure 3 shows the device transmission of several typical silicon-based MMI power splitters as a function of the MMI length
LM under the same MMI width (
WM = 2.5 μm). As shown in
Figure 3, three types of MMI power splitter were considered, conventional MMI structure, uniform SWG embedded MMI structure, and inverse-tapered SWG embedded MMI structure. One can clearly find that the variation of device transmission became small within the calculation range of
LM when the MMI region was embedded with a uniform or inverse-tapered SWG structure, corresponding to an increased MMI length tolerance. Moreover, we can also observe that the optimum MMI length was obviously decreased from 5.0 μm to 3.6 μm (3.3 μm) due to the embedded uniform (inverse-tapered) SWG structure in the center of the MMI region along the propagation direction, contributing to on-chip compact integration. If we further added linearly tapered structures connecting the input/output waveguides with the MMI region, the device transmission could be obviously increased, while the variation range of transmission could be reduced, as shown in
Figure 3. By comparison, we can find that the inverse-tapered SWG embedded MMI structure with input and output tapers obtained the best power splitting function together with largest tolerance for MMI length
LM, where the device transmission was about −0.25 dB at
LM = 3.3 μm.
To perform a more detailed analysis of the inverse-tapered SWG structure, the embedded SWG parameters were defined as follows:
where
k is the width increment of the inverse-tapered SWG along the
y-direction;
k was chosen as 50 nm during structural optimization. More details about the choice of
k can be found in the tolerance analyses described later. Therefore,
LM,
LS, and
n are closely related.
Figure 4a shows the detailed device transmission with MMI length
LM, where the optimal half period number
m (
n =
m/2) of the SWG structure under different
LM is also plotted. From
Figure 4, we can find that the device transmission was higher than −0.5 dB within the whole range from
LM = 2.6 to 4.0 μm, and the optimal period number almost increased with the increase in
LM. Note that the device transmission reached a high value with relatively small fluctuation when
LM varied from 3.2 μm to 3.5 μm, marked by the oval in
Figure 4a. By consideration of the device size and required grating number, we set
LM to 3.2 μm for the subsequent analysis, where the optimal half period number
m was 18 (
n = 9).
Figure 4b shows the device transmission with gap length
LS between the left side of the inverse-tapered SWG and the left side of the MMI region, where
LM was equal to 3.2 μm. According to Equation (2), the period number of the embedded SWG structure decreased as the gap length
LS increased within the calculation range shown in
Figure 4b. Meanwhile, the obtained transmission curve presented a nearly flat top as
LS changed from 1.1 to 1.4 μm, and we chose
LS = 1.4 μm, corresponding to the highest transmission shown in
Figure 4b.
According to above analyses, the best device transmission of the MMI power splitter embedded with the inverse-tapered SWG was limited to about −0.25 dB at the wavelength of 1.55 μm. To further enhance the device transmission or reduce the power splitting loss, we embedded two rows of uniform SWG on both sides of the MMI region, where the period number was N. Such embedded uniform SWG structures were employed to match the separated modes (double image of the input TE
0 mode) within the MMI region with the output waveguide modes. Therefore, we set the central positions of embedded uniform SWGs along the propagation direction aligned with those of the output waveguides. Thus, the period number N and relative position
LE of such a uniform SWG was also determined.
Figure 5 illustrates the device transmission as a function of
LS, where the optimal period number N and relative position
LE of the embedded uniform SWG are plotted for every gap length
LS. The grating parameters of the added uniform SWG were the same as those of the inverse-tapered SWG (grating pitch Λ = 200 nm, duty cycle
a/Λ = 0.5), and the etching width of the uniform SWG in the
y-direction was set as 100 nm, corresponding to the square etching slots. As shown in
Figure 5, the optimal value of
LS was reduced compared with no embedded uniform SWG structure shown in
Figure 5b. The reason is that the embedded uniform SWG could reduce the effective index on both sides of the MMI waveguide, and such a structure also required a further reduction in the effective index on the central side, leading to a reduction in
LS or an increase in period number
n of the inverse-tapered SWG. Therefore,
LS decreased from 1.4 to 0.9 μm, corresponding to the optimal values of N = 5 and
LE = 1.45 μm. For the device transmission of our proposed power splitter, its value was obviously reduced to −0.08 dB at a wavelength of 1.55 μm upon embedding the extra uniform SWG structures on both sides of the MMI region. Such a low transmission loss of the power splitter would be very beneficial to construct on-chip photonic devices (e.g., MZI modulators [
11,
12] and optical switches [
13,
14]) and large-scale PICs. During the calculation process shown in
Figure 5, the choice of optimal values of N and
LE for every gap length
LS was also important. Here, we used
LS = 0.9 μm as an example, and the results are plotted in
Figure 6, where the period number of the embedded uniform SWG was set as N = 4, 5, and 6. Note that the obtained device transmissions revealed some fluctuations as the relative position
LE increased from 0.85 to 2.15 μm, and the largest transmission loss was almost lower than 0.5 dB, as shown in
Figure 6. These transmission fluctuations may have resulted from the interaction of the separated modes with the uniform SWG structure embedded on both sides of the multimode waveguide, where different interact positions led to different transmission losses due to the wave features of input light. In addition, we also introduced two right-angled cutting structures on both sides of the input MMI side in a symmetric manner (
LT = 0.25 μm) to further reduce the device reflection loss (RL < −35 dB) and stabilize the light evolution through the proposed device. Therefore, through these structural designs and optimizations applied in the MMI region for power splitting, the device transmission loss was reduced to only 0.08 dB, while the required MMI length was 3.2 μm, representing an improvement compared to the conventional MMI power splitter without an embedded SWG structure.
Next, we conduct wavelength spectrum analyses for the proposed device, where IL is used to denote the device transmission loss.
Figure 7 shows the wavelength dependence of IL for the proposed device and conventional MMI power splitter, where the inset shows the device schematics, and key parameters are also marked. For a better comparison, we chose quite a large wavelength range calculated from 1.15 to 1.95 μm, and the MMI width was set as the same (
WM = 2.5 μm) for both devices, where the material dispersions of silicon and silica were also considered [
31]. From
Figure 7, we can clearly find that the allowable working bandwidth of the conventional MMI power splitter was very large (1360–1770 nm) when keeping IL < 0.6 dB. This is the main reason for the MMI structure’s superiority over the DC or PC structure for power splitting [
8,
9,
10]. By comparison, the allowable working bandwidth of our proposed power splitter could be extended from 1240 nm to 1800 nm covering the whole optical communication band if IL < 0.6 dB was also satisfied. Therefore, the obtained working bandwidth (560 nm) of the power splitter based on our proposed structure could be even higher than that of the commonly used broadband MMI power splitter (410 nm), revealing the ultra-broadband feature of our proposed device.
Table 1 compares our proposed power splitter with other MMI power splitters embedded with an SWG structure reported recently, where the MMI dimension, IL, RL, and allowable working bandwidth were all considered. It can be noted that the proposed power splitter had obvious advantages of ultra-broadband, low IL, and low RL, while the required MMI dimension was comparable with other reports. Therefore, the present device can be employed as an efficient and broadband power splitting component applied in on-chip PICs.
For the device fabrication, we only needed one-step lithography and etching processes on a commercial SOI wafer with 220 nm-thick top silicon layer, since the proposed device had a uniform etching depth (220 nm). Moreover, the required minimum linewidth was 100 nm, which can be easily achieved by current E-beam lithography [
1,
32]. Within our proposed device, the most important structures were the inverse-tapered SWG embedded in the MMI center and the two rows of uniform SWG embedded on both sides of the MMI region. Here, we mainly considered the device performance affected by the lateral shift of the embedded SWG structure in the
y-direction with respect to the MMI waveguide, and the results are shown in
Figure 8, where Δ
wS and Δ
wE represent the lateral shift (
y-direction) of the inverse-tapered SWG and uniform SWG from their designated positions mentioned above. Moreover, we introduce a new index called the splitting ratio (SR) to characterize the ratio of light power received at output ports due to structural parameter deviations,
where
POutput1 and
POutput2 are the receiving power at the two output ports illustrated in
Figure 1. Considering the structural symmetry of our proposed power splitter, performance variation due to lateral shift of Δ
wS was only calculated on one side (positive
y-direction, Δ
wS > 0), as shown in
Figure 8a,b, and the results due to a shift on the other side (negative
y-direction, Δ
wS < 0) were the same. To keep SR higher than 0.9, Δ
wS should be controlled within the interval [−30 nm, 30 nm]. Moreover, for the two rows of uniform SWGs embedded on both sides of the MMI region, we considered one row shifting from its optimal position due to fabrication imperfections, while the other row’s effect on device performance was considered identical due to structural symmetry. As shown in
Figure 8c,d, the available variation of Δ
wE was within the interval [−140 nm, 200 nm] when keeping SR > 0.9. By comparison, the central SWG structure had a tighter tolerance range with regard to device fabrication because such a structure had a strong influence on the mode splitting in the MMI region. Thus, the variation range of SR for the central inverse-tapered SWG shift was larger than that for the bilateral uniform SWG shift. Under such conditions, we further analyzed the effect of the central inverse-tapered SWG shift on the power ratio between the two output ports, as shown in
Figure 8e. Note that we could achieve a change in power ratio from nearly 30:70 to 70:30 by only shifting the central inverse-tapered SWG structure relative to the silicon waveguide along the
y-direction. Such a characteristic would introduce new applications for the power splitter. We also analyzed the effect of the width increment
k of the inverse-tapered SWG on the device performance, as shown in
Figure 9a. As shown in this figure, the highest IL was lower than 0.3 dB within the whole calculation range from
k = 20 nm to
k = 100 nm, and the optimum value of
k was located at the position of
k = 40 nm or
k = 50 nm. This is why we chose
k as 50 nm in the previous analyses.
Figure 9b shows the effect of grating size variation on device performance (IL). Here, we considered the grating width variation along both
x- and
y-directions. For example, a grating size of 100% corresponded to the optimum widths along the
x- and
y-directions, while a grating size of 120% corresponded to a width increment of 20% relative to the optimum widths along the
x- and
y-directions. According to the results, the available grating width variation range should be controlled within the range of 53% to 135% relative to the optimum values by keeping IL < 0.3 dB, where 53% indicates a grating width decrement of 47% relative to the optimum width. Such a relatively large width variation range is beneficial for device fabrication. Therefore, these obtained tolerance ranges of the key embedded SWG structures in the MMI region need to be guaranteed during practical device fabrication [
32].
On the basis of the abovementioned device design and optimization processes, the final device parameters were as follows:
w1 =
w4 = 0.5 μm,
w2 = 0.94 μm,
w3 = 0.6 μm,
LI = 1.38 μm,
LO = 2.0 μm,
WM = 2.5 μm,
LM = 3.2 μm,
LS = 0.9 μm,
LE = 1.45 μm,
LT = 0.25 μm,
m = 23, N = 5,
wg1 = 0.2 μm, and waveguide thickness = 220 nm.
Figure 10 plots the electric field evolution through the designed power splitter, where the working wavelength was 1.55 μm. From
Figure 10, we can clearly observe that the input fundamental TE mode could be evenly separated into the two output ports, and the electric field evolution was quite stable without any fluctuations or ripples. In comparison to the conventional MMI structure, no clear periodic interference pattern could be observed. The reason is that the embedded SWG region had a lower refractive index compared with the silicon waveguide, and the inverse-tapered shape of the SWG structure broke the period interference behavior in the conventional MMI structure. Furthermore, we can find that its evolution pattern was similar to that of a typical Y splitter. Thus, we cannot find a standard period interference pattern in the SWG assisted MMI region, but its working principle is still based on the MMI effect since the power splitter based on the Y junction nearly cannot be realized in a length of only 3.2 μm. Both a small branch angle and a large conversion length are required for the Y junction-based power splitter. For the proposed power splitter, its total device length is ~6.5 μm when both considering the input and output tapers, while the MMI length is only 3.2 μm. Using such a device, we can efficiently realize the power splitting function with an ultrabroad bandwidth, low insertion loss, and low reflection loss in a compact size, which can be used as the fundamental component for constructing other photonic devices and would be very promising for building on-chip large-scale PICs [
26,
27,
28].