1. Introduction
In photonic integrated circuits, energy is transferred in different structures, and even when it is transferred in the same structure, the structural parameters will change, making it necessary to design a connection structure to transfer energy from one module to another [
1,
2,
3,
4,
5]. In order to realize the design of connection structures, adiabatic devices based on the adiabatic mode evolution play an increasingly important role [
6,
7,
8]. In adiabatic devices, the designated mode at the input slowly evolves into the designated mode at the output. The output mode can be the same or different from the input mode, which is referred to as transfer efficiency when the modes are the same and conversion efficiency when the modes are different. No other unwanted modes were significantly excited during this evolution. The power transfer/conversion efficiency from the designated input mode to the designated output mode is very high.
One way to achieve adiabatic mode evolution transfer is to connect different functional units linearly. However, the resulting device structure has a long size [
9], which runs counter to the direction of optical integration towards higher integration. The analytical method can reduce the size of the adiabatic device. However, the analytical method design for adiabatic devices is only suitable for simple device structures, and the analytical method is not universal [
6,
7,
8,
9,
10,
11,
12,
13,
14]. The methods proposed in refs. [
6,
7,
8] are only suitable for the design of adiabatic couplers and are not suitable for the design of adiabatic tapered waveguides or adiabatic mode couplers in this study. Refs. [
10,
11,
12,
13,
14] are all analytical designs of adiabatic tapered waveguides, in which ref. [
10] describes mode propagation using a ray model, thus perfecting the analytical design rules for tapered waveguides proposed by Milton and Burns. In ref. [
14], an equivalent waveguide concept was proposed. This concept converts a straight waveguide into an equivalent conformal structure via conformal mapping and is applied to the design and analysis of a fully adiabatic tapered waveguide. These analytical methods are only suitable for the design of specific types of adiabatic devices. They are not suitable for the design of other types of adiabatic devices, which greatly limits the wide range of applications of these analytical methods. In the last three years, in order to design adiabatic devices with complicated three-dimensional (3D) geometries, numerical design methods [
15,
16,
17] have been proposed. Numerical methods can handle various types of adiabatic devices and can also be used to design device structures with geometric shapes, which is convenient for applications. The numerical design methods for the adiabatic device [
15,
16,
17] are calculated accurately. The methods are universal, but they require simulations for each section. Determining how many sections are divided, how many simulations are required, and selecting device length from simulation results increase the complexity of the design of the adiabatic device.
Therefore, in order to eliminate the need for simulations for each section, this study proposes a functional segmentation design method for the design of adiabatic taper waveguides and adiabatic mode converters. For the design of adiabatic taper waveguides, the device size of the quadratic function design is much smaller than that of the conventional linear-shape connection design (see
Section 2). For the design of adiabatic mode converters, the cubic function segmentation design method in this study is better than the analytical method in ref. [
9] (see
Section 3).
2. Quadratic Functional Segmentation Design of Adiabatic Taper Waveguides
Figure 1 shows two strip waveguides with different waveguide widths. The silicon core in the strip waveguide in
Figure 1a has a narrow waveguide width
WI, and the silicon core in the strip waveguide in
Figure 1b has a wide waveguide width
WF. The purpose is to design an adiabatic taper waveguide structure that connects the narrow and wide strip waveguides in
Figure 1 and efficiently transfers energy from one strip waveguide to the other with the shortest possible device length. The wavelength of the light beam is set to 1.55 μm.
A quadratic function is used to shorten the device length of the adiabatic taper waveguide. As shown in
Figure 2, assuming the
x-axis is horizontal and the
y-axis is vertical, the expression for the quadratic function is:
where A is a constant.
As shown in
Figure 2, the width is denoted by
WI =
yI at the left end and
WF =
yF at the right end. The absolute length is
LTOT =
xF −
xI. Similar to the numerical design methods for adiabatic devices that require segmentation, the functional segmentation design presented here also requires the entire structure to be divided into several distinct sections. A section is randomly intercepted in the quadratic function, and the coordinates of the initial and ending ends of the section in the upper part function are (
xi,
yi) and (
xi + 1,
yi + 1), respectively. Substituting the coordinates (
xi,
yi) and (
xi + 1,
yi + 1) into Equation (1) yields:
Since the width
Wi = 2
yi and the width
Wi + 1 = 2
yi + 1, substituting into Equations (2) and (3), respectively, yields:
Squaring both sides of Equations (4) and (5) yields:
Subtract Equation (6) from Equation (7) to obtain:
Considering
Wi =
WI and
Wi + 1 =
WF, we can obtain the expression for the coefficient A:
Substituting the expression (9) of the coefficient A into the Equation (8) gives the length of each section:
Equation (10) is used to calculate the length for each section in adiabatic taper waveguides. This equation is free from axial limitations and only depends on the structural parameters to be designed. For each section, the input and output widths Wi and Wi + 1 are connected linearly.
As shown in
Figure 1, the width of SiO
2 is
W0 = 10 μm, the thickness of SiO
2 is
h0 = 1.5 μm, the thickness of silicon is
h = 300 nm, and the width of silicon increases from
WI = 2 μm to
WF = 8 μm. The uniform segmentation is used to divide the width from
WI = 2 μm to
WF = 8 μm intervals (Δ
W = 0.1 μm) into
N = 60 sections. The higher the number of segments, the closer the shape is to the shape of the true quadratic function, but this also leads to higher computational demands. The widths of the input and output of the
i-th section are
Wi = [2 + (
i − 1) × 0.1] μm and
Wi + 1 = (2 +
i × 0.1) μm, respectively, where
i increases from 1 to
N. According to Equation (10), the length for each section can be obtained.
Each section is constructed according to the input width
Wi and output width
Wi + 1 of each section. The length of each section is obtained by Equation (10), and then all the sections are spliced to form a complete waveguide shape. Then, the power transfer efficiency curve of the entire structure can be obtained by the Finite-difference time-domain (FDTD) method or the Eigenmode Expansion (EME) method, as shown in
Figure 3. The figure shows the lengths corresponding to the power transfer efficiency, which can be used in practical applications.
As shown in
Figure 3, when 95% power transfer efficiency is required, the length required for the quadratic function segmentation design is 36 μm, while the length required for the linear-shape connection is 54 μm. When 99% power transfer efficiency is required, the function segmentation design requires a length of 78 μm, while the linear-shape connection requires a length of 109 μm. As a result, the quadratic function segmentation design is designed with a much smaller device size than the conventional linear-shape connection.
Ref. [
15] defines the adiabatic regime. When the power loss is greater than 10%, the adiabatic performance is poor. We can say that we are in an acceptable adiabatic regime when the power loss is less than or equal to 10%. Therefore, we choose a minimum power transfer efficiency of 90%. We can also choose a power transfer efficiency of 95% when we need a better adiabatic regime. Certainly, if someone wishes to pursue a better adiabatic regime, they can also choose a higher power transfer efficiency, such as 99% or even higher.
3. Cubic Functional Segmentation Design of Adiabatic Mode Converters
Figure 4 shows that the strip waveguide has a lower cladding of SiO
2 and an upper cladding of Si
3N
4, which makes the strip vertically asymmetric and allows mode conversion between different modes. The goal is to design an adiabatic mode converter that efficiently converts one mode at the input to another mode at the output with the shortest possible device length. The wavelength of the light beam is set to 1.55 μm.
As shown in
Figure 5, considering the horizontal direction as the
x-axis and the vertical direction as the
y-axis, the expression for the cubic function is:
where A and B are constants.
Input width
WI = 2
yI, output width
WF = 2
yF, and absolute length
; substitute coordinates (
xi,
yi) and (
xi + 1,
yi + 1) into Equation (11) to obtain:
where
yi + 1 >
yi > 0,
xi and
xi + 1 are calculated from Equation (11).
Since the input width is
Wi = 2
yi and the output width is
Wi + 1 = 2
yi + 1, substituting Equations (12) and (13) yields the following:
Equations (14) and (15) are opened on both sides of the equation to yield:
Subtract Equation (16) from Equation (17) to obtain:
Considering
Wi =
WI and
Wi + 1 =
WF, give the expression for the coefficient A:
Substituting the expression (19) of the coefficient A into Equation (18) yields the individual fragments:
Equation (20) is used to calculate the length for each section, which is independent of the coordinates and is only related to the structural parameters to be designed.
As shown in
Figure 6, strip waveguides of different widths
W are simulated using a time-domain simulation algorithm to obtain the effective refractive index for each mode corresponding to each width. The thickness of the silicon core is
h2 = 220 nm. From the figure, it can be determined that the width
W0 = 0.75 μm of the mode hybrid region appears. To achieve mode conversion, the input width
WI and output width
WO of the waveguide need to meet the following conditions:
WI <
W0 <
WF. After some simulation, the Eigenmode Expansion (EME) simulator from Lumerical is used to determine whether mode conversion really exists. Finally, it is determined that the input waveguide width is
WI = 0.69 μm and the output waveguide width is
WF = 0.83 μm.
This example uses non-uniform segmentation to divide the width from W = 0.69 μm to W = 0.73 μm at intervals ΔWNi = 0.002 μm into 20 sections; divide the width from W = 0.73 μm to W = 0.78 μm at intervals ΔWHi = 0.001 μm into 50 sections; divide the width from W = 0.78 μm to W = 0.83 μm at intervals ΔWNi = 0.002 μm into 25 sections. This section satisfies that the interval ΔWHi of the hybrid region is less than twice the interval ΔWNi of the non-mode hybrid region, that is, ΔWHi < 2ΔWNi.
Each section is constructed according to the input width
Wi and output width
Wi + 1 of each section. The length of each section is obtained by Equation (20), and then all the sections are spliced to form a complete waveguide shape. Then, the power conversion efficiency curve of the entire structure can be obtained by the Finite-difference time-domain (FDTD) method or the Eigenmode Expansion (EME) method, as shown in
Figure 7. The figure shows the lengths corresponding to the power conversion efficiency, which can be used in practical applications.
As shown in
Figure 7, when 90% power conversion efficiency is required, the length required for the cubic functional segmentation design is 37 μm. When 99% power transfer efficiency is required, the cubic functional segmentation design requires a length of 59 μm. In
Figure 7, the analytical method design in ref. [
9] is added as a comparison, which can better reflect the advantages of the cubic functional segmentation technique in this study. It can be seen from the figure that the method in ref. [
9] has obvious oscillation, indicating a large energy loss. Therefore, the effect of the cubic functional segmentation design is better than that of the design in ref. [
9]. The choice of total length in practical applications requires balancing the power conversion efficiency with the length requirement. Small device lengths can be chosen to achieve higher integration design goals when the conversion efficiency is not too high.
The variation in the effective refractive index with the width of the waveguide is one way in which we can determine the existence of a mode hybridization region. In this way, it is determined that there is mode conversion in some areas; for example, it is determined that there is mode conversion near
W0 = 0.75 μm in
Figure 6. To achieve mode conversion, the input width
WI and output width
WO of the waveguide need to meet the following conditions:
WI <
W0 <
WF. After some simulation, the Eigenmode Expansion (EME) simulator from Lumerical is used to determine whether mode conversion really exists. Finally, it is determined that when the input waveguide width is
WI = 0.69 μm, and the output waveguide width is
WF = 0.83 μm, the mode conversion can indeed be realized. The conversion efficiency is shown in
Figure 7. It is proven that there is indeed a mode conversion from width
WI = 0.69 μm to width
WF = 0.83 μm.
Through various works [
1,
2], it has been shown that quadratic function shapes generally give more efficient adiabatic taper waveguides compared to linear-shape designs. Ref. [
9] shows that cubic function shapes often give efficient adiabatic mode converters. Therefore, a quadratic function segmentation method is developed for the design of adiabatic taper waveguides and a cubic function segmentation method for the design of adiabatic mode converters. Certainly, the cubic function segmentation method can also be used in the design of adiabatic taper waveguides, but the results are poor, as shown in
Figure 3. It can be seen from
Figure 3 that the transfer efficiency of the adiabatic taper waveguide designed with cubic function segmentation is much worse than that of the adiabatic taper waveguide designed with quadratic function segmentation and even worse than that of the linear-shape connection.