Abstract
1. INTRODUCTION
Photonic integrated circuits (PICs) have drawn increasing attention over the past two decades. Their primary goal is to integrate complex manipulation of light (such as routing, filtering, coupling, interfering) onto a single chip [1–3]. Today, a PIC is usually designed for one specific application, so that it can be compact and power-efficient [4]. The design methodology for these chips is similar to that of application-specific integrated circuits (ASICs) in the electronic domain, and thus this kind of PIC is usually referred to as an application-specific PIC (ASPIC).
In contrast, another mainstream type of electronic circuit is the field-programmable gate array (FPGA). These circuits are generic in concept, and their functionality is programmed by configuring the on-chip connectivity of the logical building blocks. The photonic counterpart of FPGA, the programmable photonic integrated circuit (PPIC) [4–15], has been introduced recently based on the idea of run-time manipulation of light after a chip has been fabricated. Such reconfigurability is usually made available by controlling the active components (e.g., optical PSs [4]) with electrical/thermal signals. Due to its programmability, a PPIC is suitable for various applications such as fast prototyping of ASPICs [4], building optical neural networks (ONNs) [16], and processing quantum information [17,18].
A PPIC is composed of a mesh of tunable basic units (TBUs) [12], also called analog optical gates [4]. The most common implementation of a TBU is a Mach–Zehnder interferometer (MZI) circuit [4,12]. Considering the interconnections of TBUs, PPICs can generally be classified into two categories: (i) forward-only topologies [6–11,19–21], and (ii) loop-back (recirculating) topologies [4,12–15]. In a forward-only PPIC, light propagates in one direction (e.g., from left to right). It has been proven that with particular forward-only structures, a PPIC can realize any unitary transformation [10,19,22]. When fixed-length delay lines are introduced, it is also possible to implement finite impulse response (FIR) digital filters [11]. Feed-forward PPICs are commonly used to implement ONNs for AI computing. The first notable experimental realization of an ONN was published in 2017 [16]. Later works have considered
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
However, without loop-back connections, a forward-only PPIC cannot realize a ring resonator or an infinite impulse response (IIR) digital filter. Such shortcomings have motivated researchers to consider recirculating-based PPICs [15,23,24]. The most common recirculating configurations are triangular, square, or hexagonal close packing [25]. However, while these loop-back meshes offer the possibility of implementing more complex connectivities as well as FIR and IIR filters, the configuration of those functions is mostly done by manually assigning and configuring the optical gates in the mesh. Such an
To address these issues, a few published results [13,14,26] have proposed methods to perform this task automatically. The authors in Ref. [26] proposed to use optimization techniques to synthesize optical ring resonators and MZIs on a hexagonal-mesh PPIC. In Ref. [14], the authors proposed an auto-routing method based on graph theory for a hexagonal-mesh PPIC, and multiobjective routing is demonstrated by Ref. [13]. However, these methods can be dramatically improved to overcome the following key limitations: (i) their application range is restricted and many light-processing functions are not considered; and (ii) since many optical PSs need to be optimized in a PPIC, this high-dimensional optimization problem is not efficient with current methods that rely on nongradient methods (e.g., particle swarm optimization (PSO) in Ref. [26]) or gradient methods with numerical differentiation (e.g., Eqs. (4) and (5) in the supplementary material of Ref. [26]).
In this paper, we address these two main points by relying on scattering matrix theory, together with efficient calculation of analytical gradients. Specifically, we propose an efficient method that can realize configurations for many different light-processing functions on a square-mesh PPIC, without requiring
2. THEORY OF SCATTERING MATRICES
Following Refs. [4,9,15], we consider the TBU structure as shown in Fig. 1 throughout this paper. As shown in the top row of Fig. 1, we assume that two time-harmonic optical inputs are provided, respectively, at the two left ports . Then the outputs can be calculated based on the transfer matrix ,
Figure 1.Simplified schematic of a TBU. It is made up of two 50%:50% DCs on the left and right, and two optical PSs parameterized by
If we reverse the direction of light propagation, as shown in the bottom row of Fig. 1, then the vector on the left-hand side of Eq. (1) will be , while will be on the right. Combining the two propagation cases together, we have the scattering matrix relation,
For our filter synthesis application, the model in Eq. (2) is insufficient: we will never obtain a frequency-dependent response using this model, because does not rely on the light frequency . To remedy this, we modify the previous transfer matrix by taking the role of the TBU waveguides into consideration,
The circuit schematic of the recirculating PPIC waveguide mesh is shown in Fig. 2. In this paper, we ignore the TBUs in the right-most column. We also assume that the top and bottom connections (the yellow lines) are ideal connections, i.e., their transfer function is identity. These two assumptions are made for mathematical simplicity and the purpose of demonstration; note that our method is applicable without these assumptions.
Figure 2.Schematic of an
Next, we introduce naming conventions for the ports and propagation directions. As shown in Fig. 3, we adopt the following conventions for the ports in this PPIC: (i) the letters “A” and “B” are used to denote the ports on the left and right edge of a vertical TBU, respectively; (ii) the subscript is used to express that the port is on the th row and th column, where and . As shown in Fig. 3, for any port, the light can propagate in two directions. We define going into and out of the vertical TBU device as “I” (i.e., orange arrows) and “O” (i.e., purple arrows), respectively. One minor subtlety arises when applying this direction naming convention to the top line shown in Fig. 2, since this top line does not associate with any vertical device. In this case, we consider there to be virtual vertical TBUs above this top line, and then apply our direction naming convention. Similarly, we consider there to be virtual vertical TBUs beneath the bottom line in Fig. 2 for the purpose of notation consistency.
Figure 3.Naming conventions for port and direction. Capitalized “A” and “B” should be regarded as port names, and lowercases “a” and “b” are the complex magnitudes ahead of
Following Eq. (3) and applying the scattering matrix to the two propagating directions shown in the right figure in Fig. 3, we have
Recall that in writing Eq. (5), we apply the scattering matrix method to the two propagation directions of a vertical TBU. We can do the same thing for a horizontal TBU, which gives the following equation:
For a specific column index , if we vary the row index and in Eqs. (6) and (8), and next stack all the resulting equations in one column, we obtain a scattering matrix for the mapping: , where . Similarly, if we vary the row index in Eq. (10), we can write down the scattering matrix for the mapping: . Combining these two steps gives us the scattering matrix representing the mapping: . Mathematically, that is to say,
If we repeat Eq. (12) times for different column indices , then we can obtain the overall scattering matrix for the mapping from to ,
Thus far, we have obtained a relation between the input and the output. The ultimate scattering matrix is related to the individual (or ) matrix of a vertical (or horizontal) TBU device via Eqs. (18), (15), and (13), in sequence. Furthermore, the relations from and matrices to the individual PSs are also clear via Eqs. (11), (7), and (4). Thus, we have obtained an analytical expression of defined by all PSs . Although it is difficult to explicitly write down the expression for every entry in the matrix, we do know the sequential operations to construct it. Most importantly, all of the operations involved (e.g., matrix-vector multiplication) are differentiable, so that we can easily calculate or for any of any TBU device. As demonstrated later, this will form the basis for our synthesis method.
Without loss of generality, we assume that our desired forward light propagation is from left to right in the PPIC shown in Fig. 4. Then we can regard the forward input at the left and the backward input at the right both as given constant vectors. Based on Eq. (17), we can now express the forward output at the right and the backward output at the left as
Figure 4.Illustration of the forward input at the left and the backward input at the right. Note that the directions
Several points are worth noting. First, both and are of size . This provides us with some flexibility to synthesize multiple light-processing functions simultaneously. For instance, we can feed an input wave from the top port of the first vertical TBU, i.e., equal to 1 and all other entries of equal to 0. Then the outputs at the second and third entries of can be used to synthesize two different light-processing functions. Second, might not be zero in Eq. (19) even if is zero, because the information brought by can recirculate back. This is revealed by the term at the second line in Eq. (19). Third, recall we assume that the yellow lines in Fig. 2 are ideal connections, leading to the zero-one matrix in Eqs. (8) and (13). If the yellow lines are instead not ideal, we just need to revise the zero-one matrix, while our derivation (as well as the later synthesis method) still holds.
We make two additional remarks related to the yellow direct connections in the top and bottom rows. First, from the application perspective, the yellow direct connections in the top and bottom rows of the mesh introduce a peculiarity. These connections break the connection symmetry of the mesh, and in particular break the clockwise/counterclockwise degeneracy of the square waveguide mesh. Normally, when injecting light in a square waveguide mesh, light will either circulate in a clockwise or counterclockwise direction inside a unit cell, but these circulations are not coupled. This means that in the scattering matrix of a square waveguide mesh, at least half of the elements are zero. By adding the connections in the top and bottom rows, these clockwise/counterclockwise circulations can be coupled and more generic mesh functions can be defined.
Second, from the calculation perspective, introducing the yellow direct connections lets us only need to provide the forward input (i.e., scalars) if we assume the backward input . However, without these yellow direct connections, we would have to set input values for those floating ports in the top and bottom rows; otherwise, the conditions are insufficient to determine the circuit response. Analytical gradients can still be calculated in such a case, but our derivation will need substantial modification.
3. HORIZONTAL RELAXATION
In the previous section, we derive the scattering matrix for a square-mesh PPIC in a general form. In this section, we consider a simplified case under the assumption of horizontal relaxation: all horizontal TBUs are configured with and , and thus operate in the bar state [4]. This implies that in Fig. 1, the light propagates from Port to , to , or in reverse, but does not go from to . Namely, when passing a horizontal TBU, the light is confined in the upper or lower arm.
As a starting point, we consider a square-mesh PPIC under this horizontal relaxation. Its schematic is shown in Fig. 5. In this case, the matrix defined in Eq. (20) is of size . When is odd, its expression is
Figure 5.Schematic of a
A natural thought would be to extend the square PPIC under horizontal relaxation to an square PPIC. Fortunately, due to the assumptions of our horizontal relaxation, this is straightforward. Specifically, in Eqs. (21) and (22), we have a matrix with a size of corresponding to . For , we will have a matrix with a size of . Its first and last entries on the main diagonal will still be , as in Eqs. (21) and (22). The middle part of will be filled with different in a similar way to in Eqs. (21) and (22), where the th corresponds to the th row in the PPIC. To intuitively understand this, notice that the horizontal relaxation actually confines the horizontal propagation of signals in the same arm. Thus, the light propagating in the first row will never go to the second row, which means the transfer functions of two different rows are decoupled.
An important implication from this example is that even under this simplifying horizontal relaxation, the final transfer function shown in Eq. (23), though it has an analytical form, does not provide a direct solution–that is, it does not provide us with a direct analytical filter synthesis method. This motivates our optimization-based synthesis method proposed in the next section.
4. REALIZATION OF LIGHT-PROCESSING FUNCTIONS
In this section, we explain how we utilize our derivation to efficiently synthesize light-processing functions on an square-mesh PPIC. Assume that we want to attain light-processing functions represented by the complex transfer functions specifying the magnitude and phase responses in a range . We choose frequency points in this desired angular frequency range with incremental step equal to . Then we can define an error or cost function,
However, the difficulty lies in the fact that this optimization problem is extremely high-dimensional. For an square-mesh PPIC as shown in Fig. 2, it has PSs in total. Considering a fairly small PPIC, there are already 420 PSs to tune. To the best of our knowledge, such a high-dimensional optimization problem is inefficient to solve unless using a gradient descent method with analytical gradients. Specifically, nongradient methods take a long time to converge, and gradient descent methods based on numerical differentiation require many function evaluations to calculate the gradient once. Importantly, in our case, we do have the analytical derivative or for any and based on our previous derivations, because the operations that relate (or ) to the variable Cost are all differentiable. As a result, we can use gradient descent optimization to minimize Eq. (28) to perform the synthesis task. For details about how to calculate the gradient, please refer to Appendix B.
We note that in some applications, the desired light-processing functions only have requirements on the magnitude, but with no constraints on the phase. In such cases, we can choose to be real functions representing the desired magnitude response and revise the cost in Eq. (28) as
5. NUMERICAL RESULTS
In all our numerical experiments, we choose , , , and . We do not take dispersion effects into account (i.e., is considered to be constant and independent of , which means that ). Real waveguides do have dispersion, but this does not affect the method, as long as the dispersion of can be described by an analytically derivable function (e.g., with the help of ). Moreover, we emphasize that in high refractive index contrast platforms, usually depends on frequency, and the dispersion effect causes a narrow free spectral range (FSR) in the PPIC. Before moving on, we define a value for later simplicity,
Our algorithm is implemented in Python, and all our numerical experiments are performed on the same RedHat Linux server with 16 Intel Xeon E7-4850 CPUs working at 2.1 GHz. The initial guess required by the gradient descent optimization is randomly generated, consistent with our claim that our synthesis method does not require human design knowledge. However, we emphasize that in most of our examples, we are optimizing an interferometric system with many phase variables, and thus the cost function for most configurations will have many local peaks and valleys. The specific configuration coming out of the optimization algorithm will therefore depend strongly on the initial condition. Table 1 comprehensively lists the detailed information of all our experiments. In the following paragraphs, we comment on each case. Detailed Information for All Our Experiments All are performed on a In Cases 1, 2, and 3, the synthesized results have no interference; we use phase accumulation to depict how many TBUs the light path passes through (e.g., Input Port Output Port Target (s) Cost Results Run Time Phase Acc/FSR No. 1, routing Mag, phase Eq. ( Fig. 0.27 min No. 2, splitting Mag Eq. ( Fig. 1.09 min No. 3, splitting (c) Mag, phase Eq. ( Fig. 0.63 min No. 4, splitting (c) Mag, phase Eq. ( Fig. 4.48 min No. 5, filtering Mag Eq. ( Fig. 81.59 min No. 6, WDM Mag Eq. ( Fig. 108.25 min No. 7, WDM and filtering 1 at Mag Eq. ( Fig. 110.34 min
For Case 1, we consider routing the input light to an output port with minimum cost over the entire frequency band. Results are shown in Fig. 6. The synthesized path shown in Fig. 6(e) has gone through eight TBUs. Thus, according to Eq. (4), we know that the synthesized configuration has a phase accumulation corresponding to , or more specifically, that the output port has an dependence. This implies that we should witness a phase change of over a frequency range of , i.e., an interval with length 0.25 in the normalized frequency figure. This is indeed the case, as shown in Fig. 6(h). Also, if zooming in, Fig. 6(h) is exactly the same as Fig. 6(b). Since we have considered a loss term in our compact TBU model, the synthesized normalized power transmission shown in Fig. 6(g) cannot reach 0 dB. We see that the synthesized light path shown in Fig. 6 relies on the top line and passes through eight TBUs, and a quick calculation shows , consistent with Fig. 6(g). Last, but not least, Fig. 6(d) also demonstrates that only the first few rows have been adjusted by the optimization routine. This is as expected, since our input and output ports are both located at the top part of the mesh.
Figure 6.Case 1, routing. (a) and (b) show the target response
For Case 2, we consider equal power splitting to three output ports. Results are shown in Fig. 7. We note that due to reciprocity, combining three light inputs can also be readily solved. As shown in Fig. 7(d), the three light paths pass through 8, 10, and 17 TBUs, respectively, implying the three output responses should have phase accumulations corresponding to , , and . Namely, we will see a phase change of over a frequency range of , , and , respectively, corresponding to an interval with length 2/8, 2/10, and 2/17 in the normalized frequency figure. As shown in Fig. 7(g), these correctly reflect the 4, 5, and 8.5 periods in the interval of . One subtlety here is that when designing the target function , we consider the power loss due to and provide some margin in advance. Namely, the target magnitude chosen here is 0.5 on a linear scale [i.e., about in Fig. 7(a)], such that . Alternatively, choosing all three target magnitudes to be on a linear scale would be problematic. From a numerical perspective, the optimization routine would seek to push all three output magnitudes to , but since this is unattainable simultaneously due to the loss term , it could happen that the resulting three outputs would be unequal (e.g., ). Using a target magnitude that is attainable, as we do here, can prevent this issue. However, one side effect of a preprovided power loss margin is that it might encourage the light path to go through more TBUs. For instance, the zigzag light path with in this example is only one possible solution. It is obvious from Fig. 7(d) that this light path could propagate to the right bottom direction at the port with magnitude 0.57 in the middle, instead of going to the left bottom as it currently does.
Figure 7.Case 2, splitting. (a) Target equal three-way split magnitude response (normalized to the input); (b) heat map of all optimized PS values; (c) optimized configuration (
For Case 3, we consider coherent splitting. Namely, we want to split the input light to two output ports but now with identical phase. Results are shown in Fig. 8. As seen in Fig. 8(e), both light paths pass through 10 TBUs, implying that a frequency range of (i.e., an interval with length 0.2 in the normalized frequency figure) is required for a phase change of . This is also confirmed by Fig. 8(h). Moreover, we note that in a square mesh, without using the top or bottom line, the minimum number of TBUs required to propagate light from a port at left to a port at right is 10. Moreover, the optimization routine obtains a synthesized result that seems natural and readily understandable. Namely, we chose the output port row indices to be 3 and 7 in this case, while the input port row index was 1 (see Table 1). The resulting synthesized light path first goes from top left to the bottom right direction without any splitting, and then approximately stops at the middle between the output ports. Then it performs a 50%:50% power splitting and the resulting two light paths keep propagating without further splitting all the way to the output ports. This approach of first propagating to the middle followed by a 50%:50% splitting is a generic strategy to synthesize one-input to two-output coherent splitting and is automatically found by the optimization.
Figure 8.Case 3, coherent two-way splitting. (a) and (b), respectively, show the target equal magnitude split with equal phase response. (c) Heat map of all optimized PS values; (d) optimized configuration (
For Case 4, we consider a more complicated version of Case 3. Now, we attempt to do coherent splitting to four output ports. Results are shown in Fig. 9. Due to the structure of the square mesh, it is actually impossible to find four light paths all with the same length, meaning that the goal in this case is unachievable. Specifically, because the four output ports do not belong to the same clockwise/counterclockwise sub-mesh, there will be at least one light path difference. As shown in Fig. 9(e), it seems that the optimization attempts to utilize interference to approach this unattainable goal as closely as possible. From Fig. 9(g), we see that the red and blue curves both have an FSR equal to , while the green and cyan curves both have an FSR equal to . The difference of FSRs also indicates that we cannot achieve coherent splitting at an arbitrary frequency point, since these paths have a periodicity mismatch. This is also verified in Figs. 9(g) and 9(h). Our synthesized results do satisfy the given targets shown in Figs. 9(a) and 9(b): the optimization achieves coherent splitting in the normalized range , which corresponds to around a 25 GHz range in reality. However, we also notice that outside this range, the optimization cannot always achieve coherent splitting. An important note is that there are several rings in the synthesized configuration in this case and explains why we obtain a frequency-dependent response in Fig. 9(g). However, port magnitudes associated with some of the rings are smaller than 0.2, and thus are not drawn.
Figure 9.Case 4, coherent four-way splitting. (a) and (b), respectively, show the target equal four-way magnitude split and phase response. (c) Heat map of all optimized PS values; (d) optimized configuration (
For Case 5, we consider optical filtering. Results are shown in Fig. 10. As seen in Fig. 10(d), many rings have formed in the obtained configuration. We successfully achieve near 0 dB in the passband, and about in the stop band. The FSR is about , as depicted in Fig. 10(f).
Figure 10.Case 5, optical filtering. (a) Target magnitude response; (b) heat map of all optimized PS values; (c) optimized configuration (
As Case 6, we consider two-way wavelength division multiplexing (WDM), also called an optical interleaver, where the spectrum is separated into even and odd frequency channels over two outputs. From the results in Fig. 11(d), it is clear that many rings have formed in the optimized configuration. Moreover, Fig. 11(d) is plotted at the central frequency, and thus the other output port magnitudes are less than 0.3 and not drawn.
Figure 11.Case 6, WDM. (a) Target magnitude response; (b) heat map of all optimized PS values; (c) optimized configuration (
For Case 7, we consider synthesizing two light-processing functions (WDM and optical filtering) at the same time, given two in-phase inputs. Namely, we provide a complex input at and a complex input at . The two output ports for WDM are and , while that for filtering is . Results are shown in Fig. 12. Note that in Fig. 12(d), we see that some inner port magnitudes are larger than 1.0. This is possible because (i) the total input power is 2.0, and (ii) when a ring is formed, it can lead to the “intensity buildup” phenomenon [27,28] near resonance.
Figure 12.Case 7, simultaneously synthesizing two light-processing functions for two in-phase inputs. Figure caption is similar to that of Fig.
To better quantify the performance of our method, we implement two baseline methods for comparison: (i) differential evolution, a population-based gradient-free global optimization approach; and (ii) gradient descent optimization with numerical differentiation. Table 2 summarizes the run time of our method and the two baselines. We see that our proposed method achieves about computation time cost reduction compared with the implemented baseline methods. Run Time (in min) Comparison of Our Method with DE and ND DE is short for differential evolution, a population-based gradient-free global optimization approach. ND is short for gradient descent optimization with numerical differentiation. We stop DE/ND when the synthesized results attain similar cost values to our method or similar curve shapes in the magnitude or/and phase response figures. The “ Ours DE ND No. 1, routing 0.27 No. 2, splitting 1.09 No. 3, splitting (c) 0.63 No. 4, splitting (c) 4.48 No. 5, filtering 81.59 No. 6, WDM 108.25 No. 7, WDM and filtering 110.34
We emphasize that gradient descent optimization (with potentially nonconvex cost functions such as ours) is known to be only able to find local minima, and thus the specific configuration coming out of the optimization algorithm will depend strongly on the initial condition. To justify the practical utility of the proposed method, we also need to show that even with different initializations, the optimization routine can always yield a good result. Due to space limitations, we take Cases 1 and 5 as examples. We run our method with different initializations and plot the results in Figs. 13 and 14. These demonstrate the robustness of our method to random initialization. Note that for our applications, we do not necessarily need a global optimum, while a locally optimal configuration is already sufficient. Note that when the PPIC size further scales up, we would expect the optimization result to be more strongly impacted by the initialization, because more local optima might exist for a higher dimensional optimization problem.
Figure 13.Three different configurations are obtained under three random initializations; all satisfy the goal of routing in Case 1. Each row represents one synthesized configuration. Figure
Figure 14.Three different configurations are obtained under three random initializations; all satisfy the goal of filtering in Case 5. Each row represents one synthesized configuration. Figure
6. CONCLUSIONS AND FUTURE WORK
In this paper, we propose an efficient synthesis method that can be applied to realize configurations for a wide range of light-processing functions on a square-mesh PPIC. The key property that makes our method efficient is that we analytically derive the gradients of the mean squared error, or the log ratio, between target and realized circuit response with respect to the tunable PSs based on scattering matrix theory. Then, a gradient descent optimization can be carried out to synthesize the desired light-processing functions at time scales of minutes.
Another major source of nonideality comes from the thermal cross talk of heaters, which we have not considered in the main text. However, by using thermal eignemode decomposition [29], our proposed method remains applicable under a change of optimized variables and can account for thermal cross talk. A similar treatment can be adopted for other nonidealities of actuators in PPICs. Please refer to Appendix E for details.
Last, but not least, we emphasize that when nonidealities (e.g., process variation, dispersion effect, or beam-splitting error) are considered, a more complex variant of the proposed compact model might be needed. Please refer to Appendix F for details.
We also note that the matrix shown in Eq. (6) is not invertible when the two PSs have a difference (i.e., or vice versa). In this case, our approach will fail. This is consistent with our intuition: when the PSs have a difference, the vertical TBU is in the bar state, and knowing all port magnitudes related to “A” does not confer any knowledge on the port magnitudes related to “B” [see Fig. 3 and Eq. (6)]. In a real numerical implementation, this means that if the phase difference is close to , then the associated will be ill-conditioned, and our simulated frequency response at ports and the gradients might be inaccurate. Fortunately, in our optimized results, we have not encountered this issue. Observant readers might find that, for example, in Fig. 11(e), there exist a few vertical TBUs with a power coupling ratio reported as 0, implying that they are in bar states. However, this is because we only display up to integer percentage when drawing the figure for space reasons. To support the claim made in the main text, we have printed the condition number of defined in Eq. (18) as a way to examine numerical stability when running the program. But we warn of the need to pay attention to this case. In the future, this problem should readily be solved when we expand our approach to support any connections, since at that time, our scattering matrix will be set up based on graph theory, without requiring inverting .
Acknowledgment
Acknowledgment. Xiangfeng Chen and Wim Bogaerts received funding from the ERC.
APPENDIX A: JUSTIFICATION FOR THE COMPACT MODEL
Here we justify the validity of the compact TBU model given in Eq. (
APPENDIX B: EXAMPLE OF CALCULATING THE GRADIENT
According to Eq. (
As an example, if the we consider corresponds to the vertical TBU at Row 1 and Column 0, then only depends on this , so that we have
APPENDIX C: COLOR CELL ORDER IN HEAT MAP
We replot Figs.
Figure 15.Demonstration of how we plot the heat map, using Figs.
APPENDIX D: POWER COUPLING AND COMMON PS
Recall the compact model of a TBU as shown in Eq. (
APPENDIX E: CONSIDERING THERMAL CROSS TALK AND OTHER NONIDEALITIES
Recall that in the main text, we solve Eq. (
With the help of Eq. (
APPENDIX F: THE COMPACT MODEL UNDER NONIDEALITY
Our compact model, shown in Eqs. (
For any given instantiation of process variations (e.g., using random sampling), the compact model remains differentiable, and thus our proposed method can be used to generate a solution for that instance. Evaluation of performance degradation or yield (e.g., using the Monte Carlo method), becomes possible. Future research might consider the robust synthesis problem, building on the method proposed here.
References
[1] L. Chrostowski, M. Hochberg. Silicon Photonics Design: From Devices to Systems(2015).
[5] Z. Gao, X. Chen, Z. Zhang, U. Chakraborty, W. Bogaerts, D. S. Boning. Automatic realization of light processing functions for programmable photonics. IEEE Photonics Conference, 1-2(2022).
[6] D. A. B. Miller. Self-aligning universal beam coupler. Opt. Express, 21, 6360-6370(2013).
[7] D. A. Miller. Self-configuring universal linear optical component. Photon. Res., 1, 1-15(2013).
[8] C. Taballione, T. A. Wolterink, J. Lugani, A. Eckstein, B. A. Bell, R. Grootjans, I. Visscher, J. J. Renema, D. Geskus, C. G. Roeloffzen, I. A. Walmsley. 8 × 8 programmable quantum photonic processor based on silicon nitride waveguides. Frontiers in Optics, JTu3A.58(2018).
[23] J. Capmany, I. Gasulla, D. Pérez. The programmable processor. Nat. Photonics, 10, 6-8(2016).
[28] J. Heebner, R. Grover, T. Ibrahim. Optical Microresonator Theory(2008).
Set citation alerts for the article
Please enter your email address