2.1. Working Principle and Structure
The “checkerboard” imager shown in
Figure 1 is a classic structure of the photonics integrated interference imaging system [
8]. The lens array of the “checkerboard” imager is set in a
square grid arrangement, and three layers of PIC arrays are used. The first layer of 2D PIC after the lens array is used to acquire and split light from the target. The second layer of 3D PIC is used to transmit and match the light of each pair of lenses, bands, and sub-field. The third layer of 2D PIC is used for interferometry, where a balanced four quadrature detector is the waveguide chip with a phase measurement function, acquiring the mutual coherence factor
of the target.
The schematic diagram of the balanced four quadrature detector with a wide bandwidth phase measurement function proposed is shown in
Figure 2a. This device is mainly composed of a 90° phase shifter and 3 dB couplers labeled as A, B, C, and D. The schematic diagram of the 3 dB coupler is shown in
Figure 2b, and
Figure 2c is a schematic diagram of the connection section shown in the block diagram in
Figure 2a. The role of the 90° phase shifter is to introduce a 90° phase difference in the optical path.
The amplitudes of the two input beams injected from input port in1 and input port in2 are
and
, respectively, and the phases are
and
, respectively. The input optical fields of the two input ports of the device can be expressed as follows:
The transmission matrix of the 3 dB coupler can be expressed as [
25,
26]:
In Equation (3),
z is the coupling length,
K is the coupling coefficient, and the expression for the coupling coefficient
K is [
27]:
In Equation (4), , , , , is the effective refractive index of the waveguide, is the waveguide radius, and is the waveguide spacing.
The two outputs of coupler A can be represented by the following matrix:
In Equation (5), is the output of the upper branch of coupler A, and is the output of the lower branch of coupler A.
Similarly, we can obtain the two outputs of coupler B:
In Equation (6), is the output of the upper branch of coupler B, and is the output of the lower branch of coupler B.
From Equations (5) and (6), we can obtain the inputs of coupler C and coupler D. The input matrices for coupler C and coupler D are given by Equation (7) and Equation (8), respectively.
In Equation (7), is the input of the upper branch of coupler C, and is the input of the lower branch of coupler C. In Equation (8), is the input of the upper branch of coupler D, and is the input of the lower branch of coupler D.
Therefore, we can obtain the output matrices of coupler C and coupler D as Equation (9) and Equation (10), respectively:
In Equation (9), is the output of the upper branch of coupler C, and is the output of the lower branch of coupler C. In Equation (10), is the output of the upper branch of coupler D, and is the output of the lower branch of coupler D.
Based on Equations (9) and (10), the light intensity values of the output beams of the four output ports of the device can be expressed as follows:
For monochromatic light, if the 3 dB coupler has a splitting ratio of 0.5, from Equations (11) to (14), the relationship between the phase difference
and the light intensities
,
,
,
at the four output ports can be expressed as:
For broadband light, the dispersion effect of the device causes a deviation between the calculated phase difference using Equation (15) and the true value. To minimize phase measurement error while increasing bandwidth, it is necessary to carefully select the design parameters of the device.
2.2. Selection of Materials
At present, there are various types of materials used in the production of optical waveguides, including silicon [
28], silicon dioxide [
29,
30], lithium niobate [
16], triple-five compounds such as indium phosphide [
18], and polymers [
31,
32]. Each material offers distinct advantages and disadvantages. The choice of material depends on the specific functional requirements of the device to achieve optimal performance.
Table 1 summarizes the characteristics of these materials. From
Table 1, silica material has several advantages, such as minimal output phase deviation, a well-established fabrication process, low dispersion and lower cost [
33,
34]. Additionally, it is suitable for a broad bandwidth range [
35]. Therefore, this paper selects silica-based optical waveguides for a series of simulations and optimizations.
Figure 3 illustrates the cross-section of the optical waveguide made from silicon dioxide. The chosen three-dimensional structure is a widely used buried optical waveguide in waveguide-type optical devices. It consists of a silicon substrate with a silicon dioxide cladding layer (refractive index 1.445) and a core layer made of germanium-doped silicon dioxide. The refractive index of the core layer, which can be adjusted by varying the germanium doping concentration, is set to 1.456 to match the mode field with a 1550 nm single-mode fiber. This results in a relative refractive index difference of approximately 0.75% between the core and cladding layers [
36].
2.3. Design of Waveguide Cross-Section Size
To achieve single-mode transmission, it is essential to carefully design the waveguide’s cross-sectional size, as single-mode transmission helps prevent mode dispersion and minimizes signal distortion. Based on the previous section, this paper adopts a center wavelength of 1550 nm, with silicon dioxide as the material. The refractive index of the cladding layer is 1.445, while the core layer has a refractive index of 1.456. Additionally, to reduce waveguide birefringence and transmission loss, this paper designs the waveguide with a square cross-section, where the width-to-height ratio is 1:1 [
37]. Simulations are conducted to assess single-mode transmission for wavelengths of 1400 nm, 1550 nm, and 1700 nm. The results, shown in
Figure 4, plot the waveguide width (in μm) on the horizontal axis and the effective refractive index on the vertical axis. The simulation results indicate that when the wavelength is 1400 nm, the first-order mode appears in the waveguide when the waveguide width reaches approximately 6 um. When the wavelength is 1550 nm, the first-order mode appears in the waveguide when the width reaches approximately 6.5 um. When the wavelength is 1700 nm, the first-order mode appears in the waveguide when the width reaches approximately 7 um. Therefore, to ensure single mode transmission in the wavelength range of 1400 nm to 1700 nm, the waveguide width should be approximately 6 µm or less. The simulations indicate that to ensure single-mode transmission within the wavelength range of 1400 nm to 1700 nm, the waveguide width should be approximately 6 μm or less. Since smaller waveguide dimensions increase fabrication complexity, a slightly larger cross-section is acceptable while still meeting single-mode requirements. Based on the simulation results and current fabrication capabilities, this paper selects a waveguide cross-section size of w = h = 6 μm. The schematic of the silicon dioxide rectangular optical waveguide model is depicted in
Figure 5. The cross-sectional view of the silicon dioxide rectangular waveguide model is shown in
Figure 5a, where the
x-axis represents the width of the waveguide and the
y-axis represents the height of the waveguide. The transverse mode field distribution of the straight waveguide is shown in
Figure 5b, and the light field is well confined within the rectangular waveguide and the rectangular optical waveguide has
= 1.450865.
2.4. Design of 3 dB Coupler
The 3 dB coupler is a crucial component for forming a waveguide chip with a phase measurement function. Its schematic structure is illustrated in
Figure 2b, and includes the input port region, the S-bend region, the coupling region, and the output port region. The beam splitting ratio is a vital parameter of the 3 dB coupler. The expression for the beam splitting ratio SR is [
24]:
In Equation (16), is the coupling length, and is the coupling coefficient. From Equations (4) and (16), it can be observed that the parameters related to the beam splitting ratio are the wavelength , waveguide radius , waveguide spacing , coupling length, core refractive index , and cladding refractive index .
Assuming that the change of the input wavelength mainly leads to the change of the phase difference introduced by the delay line, but dose not only affects the beam splitting ratio of the 3 dB coupler. Based on Equations (4)–(10), the simulation results show the variation of phase deviation with input phase difference and input wavelength, as shown in
Figure 6. When the beam splitting ratio of coupler A changes in the range of 0.1–0.9, the beam splitting ratio of coupler B, C, and D is set to 0.5, the input phase difference varies from −π to π, and the phase deviation varies in the range from −0.6570 rad to 0.5635 rad, as shown in
Figure 6a. That is, the phase deviation is within 37.6°. When the beam splitting ratio of coupler B changes in the range of 0.1–0.9, the beam splitting ratio of coupler A, C, and D is set to 0.5, the input phase difference varies in the range of −π to π, and the phase deviation varies in the range from −0.5731 rad to 0.6329 rad, as shown in
Figure 6b. That is, the phase deviation is within 36.3°. When the beam splitting ratio of coupler C changes in the range of 0.1–0.9, the beam splitting ratio of coupler A, B, and D is set to 0.5, the input phase difference varies in the range of −π to π, and the phase deviation varies in the range from −0.3718 rad to 0.3493 rad, as shown in
Figure 6c. That is, the phase deviation is within 21.3°. When the beam splitting ratio of coupler D changes in the range of 0.1–0.9, the beam splitting ratio of coupler A, B, and C is set to 0.5, the input phase difference varies in the range of −π to π, and the phase deviation varies in the range from −0.3297 rad to 0.3145 rad, as shown in
Figure 6d. That is, the phase deviation is within 18.9°. From the data above, it can be seen that the variation in beam splitting ratio of coupler A and coupler B has a more significant impact on the output of the chip than coupler C and coupler D. In addition, it also means that when the beam splitting ratio of each coupler changes in the range of 0.1–0.9, the phase measurement deviation basically meets the accuracy requirements of 0.18π (32.4°). This provides a rough range for our next step in the parameter design of the 3 dB coupler.
The structural design of the 3 dB coupler mainly includes the parameter design of the input and output port region, the S-bend region, and the coupling region.
2.4.1. Design of the Waveguide Spacing in the Input and Output Port Region
For the 3 dB coupler, energy coupling occurs when the core layer spacing of the two input waveguides is too close. The same problem also exists in the core layer of the two output waveguides. To minimize non-essential energy coupling, it is essential to design a suitable waveguide spacing that avoids unintended coupling, thus reducing device loss and size. For a rectangular silica waveguide with a width and height of 6 μm, the relationship between the coupling length and the waveguide spacing is shown in
Figure 7. In this figure, the horizontal axis represents the waveguide spacing, and the vertical axis represents the coupling length, with the unit in μm. From
Figure 7, when the waveguide spacing is 20 μm, the waveguide coupling length reaches 4 m, and the coupling effect can be considered negligible. Therefore, this paper just needs to make sure that the spacing between the waveguide in the input and output port regions of the 3 dB coupler is designed to be greater than 20 μm.
2.4.2. Design of S-Bend Region Parameters
The S-bend region connects the input or output port regions with the coupling region. If the parameters of the S-bend region are not optimized, it can lead to significant transmission losses. Therefore, careful design and optimization of the S-bend parameters are crucial.
In this study, the cosine-type bending waveguide is employed due to its lower transmission loss compared to other designs [
38]. The cosine-type bending waveguide structure is illustrated in
Figure 8. The cosine-type bending waveguide has two important design parameters, namely the length L projected in the x-direction and the height h projected in the y-direction.
To balance device size reduction and prevent optical signal interference between adjacent paths, the value of h is usually set to 50 μm. When the center wavelength of the incident signal is λ = 1550 nm and the waveguide height is 50 μm, the relationship between the output power of the cosine-type bending waveguide and the bending waveguide length L is shown in
Figure 9. The horizontal axis of the figure is the bending waveguide length L in μm. The vertical axis is the normalized output power of the waveguide. As observed, increasing the bending waveguide length L leads to a gradual increase in output power and a reduction in loss. However, once the bending waveguide length reaches 1000 μm or more, the loss levels off. From the curve in
Figure 9, it can be seen that the normalized output power reaches its maximum value at 1200 μm, which means that the insertion loss is minimized at this point. Therefore, in this paper, the length of the bending waveguide is set to 1200 μm.
When the bending waveguide length is set to 1200 μm,
Figure 10 illustrates the output power of the cosine-type bending waveguide as a function of the input wavelength. The horizontal axis represents the input wavelength in μm, while the vertical axis represents the normalized output power of the waveguide. The figure shows that the output power of the waveguide remains relatively constant within the input wavelength range of 1.4 μm–1.7 μm.
2.4.3. Design of Coupling Region Parameters
The structure of the coupling region, which is the core unit of the 3 dB coupler, is shown in
Figure 11. The image shows two design parameters related to the coupling region, namely the waveguide spacing
and the coupling length
.
This paper first determines the waveguide spacing, which is mainly limited by the waveguide coupling. For the determined waveguide structure, the minimum coupling length of the 3 dB coupler increases exponentially as the waveguide spacing grows, and usually the value of the waveguide spacing is 6 μm~8 μm. Considering the goal of minimizing the size of the device, the waveguide spacing of the 3 dB coupler is designed as 6 μm in this paper. Then the value of the coupling length can be obtained. Based on the refractive index of the waveguide material and the waveguide cross-sectional dimensions, the coupling between the rectangular waveguides is simulated for an incident light with a center wavelength of λ = 1550 nm. The simulation structure is shown in
Figure 11, and the results are shown in
Figure 12. The left diagram in
Figure 12 depicts the coupling transmission of the mode field, revealing that the mode field is gradually coupled to the right waveguide during the transmission process of the left waveguide. The mode coupling is accompanied by the optical energy also produces alternating coupling between the two waveguides as shown in the right diagram in
Figure 12. It can be observed that when the coupling length reaches approximately 2600 μm (the intersection of the blue line and the green line in the right diagram of
Figure 12), the output values of the two waveguides are equal, which means that the coupled waveguide can satisfy a splitting ratio of 1:1.
Considering the influence of the bending waveguide on the coupling length, this paper simulates the transmission of optical signals with a wavelength of 1550 nm in the 3 dB coupler. Simulate the 3 dB coupler with the previously determined design parameters.
Figure 13 shows the transmission results at an incident wavelength of 1550 nm. The output values of the two output ports are inconsistent (the values of the blue line and green line in the right diagram of
Figure 13 when
takes the maximum value), and have a significant difference which does not meet the design requirement of a 1:1 splitting ratio. This is because the equivalent coupling length generated by the bending waveguide in the S-bend region was not taken into consideration during the initial design process. Therefore, it is necessary to account for this factor to achieve the desired performance of the coupler.
To further optimize the design parameters of the 3 dB coupler, a parameter scanning analysis of the coupling length is conducted. According to the previously determined design parameters of the bending waveguide, the variation of the two outputs (output1, output2) power of the 3 dB coupler with the coupling length is shown in
Figure 14. The horizontal axis represents the coupling length in μm, while the vertical axis represents the normalized output power of the waveguide. From
Figure 14, the straight waveguide length of the 3 dB coupler should be set to 2340 μm (the x-coordinate value of the intersection points between the black line and the red line in
Figure 14).
Figure 15 presents the final transmission simulation results for the 3 dB coupler at a wavelength of 1550 nm with coupling length of 2340 μm. As can be seen from the right diagram of
Figure 15, there is almost no difference between the output of the two output ports (the values of the blue line and green line in the right diagram of
Figure 15 when
takes the maximum value), meeting the design requirement of a 1:1 beam splitting ratio and achieving perfect beam splitting. As shown in
Figure 15, when the input power is 1 unit, the single-branch output power of the 3 dB coupler is approximately 0.48. Therefore, the sum of the two outputs of the 3 dB coupler is 0.96, the insertion loss of the 3 dB coupler is about 0.18 dB.
Based on the above determination of coupling length, numerical simulation analysis shows that when the input wavelength varies in the range of 1.4 μm–1.7 μm, the output of the 3 dB coupler (output1 and output2) varies between 0.14 and 0.83, as shown in
Figure 16. The horizontal axis represents the wavelength in μm, while the vertical axis represents the normalized output power of the waveguide. The splitting ratio of the 3 dB coupler numerically equals the normalized output power at the two output ports. In other words, the splitting ratio of the 3 dB coupler varies between 0.14 and 0.83, which satisfies the requirement of beam splitting ratio of the 3 dB coupler for the design of the broadband phase measurement chip, with a measurement accuracy better than 0.18π.
2.5. Design of Bending Waveguide Parameter
In addition to the bending waveguide in the S-bend region mentioned above, the connection section between the four 3 dB couplers also introduces the bending waveguide, as shown in
Figure 2c. This bending waveguide may cause significant losses, so the design of the semicircular arc bending waveguide used at the turn of the device must be optimized. For the simulation of the semicircular arc bending waveguide, when the incident wavelength is 1550 nm, the length of the bending waveguide is 1500 μm and the values of the waveguide width and height are both 6 μm, the effect of different curvature radius on the output of the waveguide is shown in
Figure 17. The horizontal axis of the figure is the waveguide curvature radius in μm, while the vertical axis is the output power of the waveguide. With the increase of the radius of curvature, the output power of the rectangular waveguide gradually increases. That is, the bending loss gradually decreases. From
Figure 17, when the radius of the waveguide reaches more than 5000 μm, the loss tends to stabilize. The design of the bending waveguide focuses on reducing the insertion loss of the device and the device length. For the above reasons, this paper sets the curvature radius of the semicircular arc waveguide to 5000 μm. According to
Figure 17, when the input power is 1 unit, the output power of the bending waveguide with a length of 1500 μm and a curvature radius of 5000 mm is approximately 0.98. From this, the loss of the bending waveguide can be calculated to be about 0.58 dB/cm.
The variation in output power of the semicircular arc bending waveguide across an input wavelength range of 1.4 μm to 1.7 μm as shown in
Figure 18. The horizontal axis of the figure is the wavelength in μm. The vertical axis is the output power of the waveguide. The figure shows that the output power of the waveguide remains relatively stable within the input wavelength range of 1.4 μm–1.7 μm.
2.6. Design of Crossed Waveguide Parameter
The 3 dB couplers labeled A and B, as well as the 3 dB couplers labeled C and D, involve a physical crossing of the waveguides during beam transmission. The structure of this configuration is illustrated in
Figure 2c.
Figure 19 shows the schematic structure of the crossed waveguide, with the crossing angle denoted as θ. The crossed waveguide introduces additional loss and crosstalk, therefore, the effect of the crossed waveguide on the performance of the device must be considered. To reduce the loss, the appropriate crossing angle must be selected. According to theory, the larger the crossing angle (right-angle crossing is optimal), the smaller the loss will be. However, the increase of the crossing angle will inevitably increase the size of the device, so it is necessary to carefully choose the appropriate crossing angle.
The loss of the crossed waveguide is simulated, with the layout of the crossed waveguides shown in
Figure 20. In this setup, two optical waveguides intersect: one is vertical, and the other crosses it at an angle θ.
Light with a wavelength of 1550 nm is incident from the vertical waveguide. As the crossing angle is varied, the change in optical power in the vertical waveguide is calculated, as depicted in
Figure 21. In the figure, the horizontal axis represents the crossing angle in degrees, and the vertical axis represents the output power of the vertical waveguide. The figure shows that the output power of the vertical waveguide gradually increases with the increase of the crossing angle, and the output power changes slowly when the crossing angle reaches more than 30°. Therefore, the crossing angle selected in this paper is 30°. As shown in
Figure 21, when the input power is 1 unit and the crossing angle is 30°, the output power of the vertical waveguide is approximately 0.98, from which the cross loss of the crossed waveguide can be calculated to be about 0.12 dB.
When the crossing angle of the waveguide is set to 30°, the curve of the output power of the vertical waveguide changes with the input wavelength 1.4 μm–1.7 μm is shown in
Figure 22. The horizontal axis of the figure is the input wavelength in μm and the vertical axis is the output power of the vertical waveguide. The figure shows that the output power of the waveguide remains relatively stable within the range of input wavelength from 1.4 μm to 1.7 μm.
2.7. Design of 90° Phase Shifter and Connection Part Parameter
From the structure shown in
Figure 2a, the design scheme of the device contains a 90° phase shifter. In this paper, the method of extending the length of the waveguide is used to introduce the phase shifter into the optical path, as shown in
Figure 23. In the figure,
represents the length difference between the two waveguides.
In an optical waveguide, the optical phase is expressed as:
where
is the effective refractive index of the waveguide, and
L is the waveguide length. Therefore, this paper can introduce a phase shift of 90° between the two waveguides by increasing the length of the waveguide:
Using the above Equation (17) this paper can get:
From the simulation, shown in
Figure 5b, the effective refractive index
of the waveguide is about 1.451 for a waveguide with a cross-section size of
, and the value of
is about 267 nm when the incident wavelength
= 1550 nm, that is, the extended waveguide length is about 267 nm.
A schematic diagram of the connection section is shown in
Figure 2c. By setting the length of the waveguide in the input and output areas of the 3 dB coupler and the distance between the upper and lower couplers of the two waveguides to 1000 μm, the length of the crossed waveguide can be calculated to be approximately 3863.70 μm, with a transverse length of about 3732.05 μm. To maintain the same optical path of the connection part, the other two connecting waveguides, except for the cross waveguide, are composed of four identical circular arcs, including the branch that introduces the 90° phase difference, as shown in
Figure 24.
Figure 24 indicates the radius R of each arc and the central angle β subtended by each arc. The length of the cross waveguide is 3863.70 μm, so the length of each arc is 965.93 μm. The bending radius of the semi-circular waveguide is 5000 μm, so the radian corresponding to each circular arc is converted to an angle of 11.0687°. When the 90° phase difference is introduced into the optical path, the length of the extended waveguide is about 267 nm, and the length of each arc is 966.18 μm. The bending radius of the semi-circular waveguide is 5000 μm, so the radian corresponding to each circular arc is converted to an angle of 11.0716°.
2.8. Design Results
The design results of the above parameters are shown in
Table 2.
As can be seen from
Figure 2a, the length of the chip is determined by the length of the two 3 dB couplers and the transverse length of the intermediate connected waveguide. The width of the chip is primarily determined by the width and spacing of the two 3 dB couplers and the bending radius of the semi-circular waveguide. After conducting calculations based on the data in
Table 2, it is concluded that the width of the chip is about 1230 μm, and the length of the chip is about 17212.05 μm. Considering the bending radius of the semi-circular waveguide, the size of the device is estimated to be around 17.5 mm × 1.6 mm.
Based on the previous description, we can know that the insertion loss of the 3 dB coupler is about 0.18 dB, the loss of the bending waveguide is about 0.58 dB/cm, and the cross loss of the crossed waveguide is about 0.12 dB. For a silicon dioxide straight waveguide with a refractive index difference of 0.75%, we take the transmission loss to be 0.03 dB/cm [
39,
40]. Combining the design results of the connection part, we can estimate the loss of the device at a wavelength of 1500 nm. When each of the two input ports of the structure shown in
Figure 2a is supplied with a power of 1 unit, we can determine that the output power at each output port is 0.42. Thus, the total input power of the device (
) is 2, and the total output power
is 1.68. It can be inferred that the total insertion loss of the device is approximately 0.8 dB.