Photonics Research, Volume. 12, Issue 11, 2676(2024)

Integrated photonic modular arithmetic processor

Yuepeng Wu1, Hongxiang Guo1,2、*, Bowen Zhang1, Jifang Qiu1, Zhisheng Yang1, and Jian Wu1,3、*
Author Affiliations
  • 1School of Electronic Engineering, Beijing University of Posts and Telecommunications, Beijing 100876, China
  • 2e-mail: hxguo@bupt.edu.cn
  • 3e-mail: jianwu@bupt.edu.cn
  • show less

    Integrated photonic computing has emerged as a promising approach to overcome the limitations of electronic processors in the post-Moore era. However, present integrated photonic computing systems face challenges in achieving high-precision calculations, consequently limiting their potential applications, and their heavy reliance on analog-to-digital (AD) and digital-to-analog (DA) conversion interfaces undermines their performance. Here we propose an innovative photonic computing architecture featuring scalable calculation precision and, to our knowledge, a novel photonic conversion interface. By leveraging the residue number system (RNS) theory, the high-precision calculation is decomposed into multiple low-precision modular arithmetic operations executed through optical phase manipulation. Those operations directly interact with the digital system via our proposed optical digital-to-phase converter (ODPC) and phase-to-digital converter (OPDC). Through experimental demonstrations, we showcase a calculation precision of 9 bits and verify the feasibility of the ODPC/OPDC photonic interface. This approach paves the path towards liberating photonic computing from the constraints imposed by limited precision and AD/DA converters.

    1. INTRODUCTION

    Since the breakdown of Dennard scaling [1], the quest for faster and more energy-efficient computing units has been met with significant challenges. A prevailing consequence has been the limitation of current processors in clock frequencies and utilization [2]. Stemming from inherent physical problems faced by CMOS transistor technology, it has resulted in computing performance being a major bottleneck for numerous potential applications [3]. Efforts have therefore been attempted to enhance the conventional transistor-based computing paradigm with alternative platforms. Among these, promising advances in the energy efficiency of integrated optoelectronic devices and their compatibility with CMOS [47] have led to an increased interest in integrated photonic computing, which allows operations at unprecedented bandwidths [8].

    Most published works on integrated photonic computing employ intensity as the “calculation quantity”, as illustrated in Fig. 1(a), encompassing both the digital and the analog types. The former, conducting bit-wise operations with optical logic gates, enables superior performance to their electric counterparts [9,10]. However, optical intensity signals suffer losses from both the material absorption and the intrinsic property of non-reciprocal Boolean operation [11], making it challenging to be cascaded to form a photonic arithmetic unit with a high bit-width [12] using photonic digital computing architectures.

    (a) General structure of emerging photonic computing. The digital-system-compatible and efficient calculation is performed by exploiting converters and the intrinsic alignment between a controllable physical phenomenon and a specific arithmetic method. (b) Overview of residue number system. (c) Schematic diagram of the phase-based photonic modular arithmetic processor (PPMAP) architecture. (d) Mapping between operands and phase when taking modulus mi=12 as an example, along with a demonstration of basic modular arithmetic performed based on this mapping.

    Figure 1.(a) General structure of emerging photonic computing. The digital-system-compatible and efficient calculation is performed by exploiting converters and the intrinsic alignment between a controllable physical phenomenon and a specific arithmetic method. (b) Overview of residue number system. (c) Schematic diagram of the phase-based photonic modular arithmetic processor (PPMAP) architecture. (d) Mapping between operands and phase when taking modulus mi=12 as an example, along with a demonstration of basic modular arithmetic performed based on this mapping.

    Photonic analog computing, on the other hand, utilizes the manipulation of optical analog quantities to efficiently execute algebraic operations. Substantial advancements in intensity-based photonic analog computing have been made in conducting operations, such as matrix multiplication [1315], convolution operations [1618], and neural network model inference [1921]. Despite the successes of analog photonic computing in several domains, these architectures are generally optimized for low-precision operations and applications that can tolerate high noise levels. This optimization limits the ability to fully leverage the flexibility and efficiency of analog signals to meet the broader demands of wider application domains. Consequently, the development of an analog-based photonic arithmetic processor architecture with higher and scalable precision becomes imperative.

    However, utilizing analog optical signals for high-precision computations faces two primary challenges. First, optical signals themselves inherently exhibit lower precision. The analog intensity signal is susceptible to noise, which limits its calculation precision due to a relatively poor signal-to-noise ratio (SNR) [22]. The calculation precision of published works (usually 4–5 bits [23,24]) falls short of meeting general requirements, with no efficient method available to improve their calculation precision as compared to digital computation systems [25,26]. Second, the mandatory use of analog-to-digital converters (ADCs) and digital-to-analog converters (DACs) incurs additional costs and confinements. A photonic analog computing system operating at ultra-high bandwidth necessitates high-precision ADCs and DACs matching its speed. This requirement poses a significant challenge as there is a well-known tradeoff between bandwidth and resolution for the electric AD/DA converters [27,28], while also posing significant energy consumption overhead [29,30] to the integrated photonic computing system. These “precision challenge” and “converter challenge” have hindered the practical implementation of the intensity-based methods and have necessitated the pursuit of alternative solutions for integrated photonic computing.

    Phase is another controllable physical quantity in integrated photonics apart from the intensity and has not yet been fully explored in the field of photonic computing, which can be controlled more efficiently by an integrated phase modulator (PM) [31,32] and exhibiting intrinsic periodicity. That feature of the optical phase aligns naturally with the modular arithmetic in which numbers are folded by a specific modulus [33,34]. For modulus arithmetic, leveraging a non-traditional numerical representation format, i.e., residue number system (RNS), allows for the decomposition of high-precision computations into multiple parallel low-precision operations, as depicted in Fig. 1(b). The correspondence between the modular arithmetic and the optical phase opens up an avenue for the efficient extension of calculation precision in integrated photonic computing systems [35]. Meanwhile, representing information in the phase domain offers the opportunity to alleviate the dependence on electric AD/DA converters through the design of optical counterparts. In this paper, we propose an integrated phase-based photonic modular arithmetic processor (PPMAP) architecture aimed at addressing the precision and converter challenges prevalent in photonic computing. Utilizing the unique characteristics of phase, its components, PPMAP units, enable accurate execution of basic modular arithmetic operations and support multi-operand operations, which is non-trivial for specialized accelerators [36]. The precision of this architecture can be effectively elevated using the RNS theory. In a proof-of-concept experiment, the three-unit PPMAP achieved 9-bit calculation precision, which is unprecedented for intensity-based photonic computing architecture. Moreover, based on the multi-operand characteristics of PPMAP architecture and our previous work on a phase-shifted optical quantizer [37,38], we demonstrate the accurate loading and extraction of information into and from the phase domain with the seamlessly integrated photonic conversion interfaces, which makes an advancement in eliminating the constraints of AD/DA converters. Further simulations on a 4-unit 15-bit PPMAP indicate comparable reliability to electronic processors under reachable SNR conditions.

    2. PPMAP ARCHITECTURE

    The basic structure of the PPMAP unit is depicted in Fig. 1(c). A beam splitter (BS) evenly splits the power of an incident monochromatic laser into two beams, which propagate through two separate paths. Each path is equipped with a set of cascaded phase modulators (PMs). Finally, the phase difference between the signals in the two paths is extracted by a phase-to-digital converter (PDC) linked at the end. Since both beams originate from the same laser source, in the absence of a voltage bias on the modulator, they exhibit a stable static phase difference denoted as φ0. When a voltage is applied to one of the modulators, it introduces a dynamic phase difference denoted as Δφ. In the case of simultaneously applying voltage to the PMs, the overall dynamic phase difference is the sum of their individual contributions.

    To perform modular arithmetic operations under modulus m, m discrete phase points from the analog quantity Δφ[0,2π) are designated to represent the m elements in the set {0,1,2,,m1}. These points must include the phase of 0 as the identity element and be uniformly spaced with an interval of 2π/m to satisfy the closure property, which ensures that the result of the addition operation remains within the same set [33]. Consequently, the integers in modulo m are mapped to the discrete dynamic phase differences ΔφZ={0,2π/m,4π/m,,2π2π/m}. This mapping allows for modular addition or subtraction operations by loading operands into the corresponding phase modulators on the same or different paths, as shown in Fig. 1(d), which illustrates an example of m=12. To perform a modular addition operation with operands “10” and “7”, followed by subtraction with “3”, we regulate two PMs on the same path and one PM on the other path to induce phase changes of 5π/3, 7π/6, and π/2. The closure property allows us to derive the results once we identify the specific discrete phase points from ΔφZ. The PDC then determines the resulting phase Δφ=|7π/3|2π=π/3, corresponding to the calculation result of “|10+73|12=2”. This efficient implementation of modular addition also enables modular multiplication, using the same structure but with the inclusion of an index-sum multiplication method [39] that performs modular multiplication via a simple modular addition and a table lookup operation (see Appendix G).

    Additionally, by leveraging the theory of the residue number system (RNS), which represents the integer operand, denoted as X, by its remainders of a set of pairwise coprime moduli {m1,,mN}, i.e., the RNS format, XRNS=X%{m1,,mN}={|X|m1,,|X|mN},the original operation with high bit-width operands can be broken into multiple parallel modular arithmetic operations with much lower bit-width operands and different moduli, as long as the result is within the range M=i=1Nmi. By orchestrating multiple PPMAP units with different moduli according to the RNS theory, as shown in Fig. 1(c), we construct the PPMAP architecture with a bit-width of log2i=1Nmi, thereby realizing calculation with scalable precision.

    3. RESULTS

    A. High-Precision Integrated Photonic Calculation

    Figures 2(a) and 2(b), respectively, illustrate the experimental setup and the integrated photonic chip utilized for the proof-of-concept experiments. The latter was fabricated on an SOI platform and corresponds to the fundamental structure depicted in Fig. 1(c). A continuous wave (CW) laser is coupled into the integrated waveguide through a grating coupler and is evenly separated via a subsequent 1×2 multimode interferometer (MMI). Those two laser beams then traverse two adjacent waveguides where four PMs load operands into the produced phase difference. The magnitude of the phase difference induced by a specific voltage value applied to the PM can be precisely controlled by the method described in Appendix B. Then, the two laser beams with a phase difference carrying the calculation result of those operands enter into a 2×K MMI to identify and quantify phase information. Figure 2(c) displays a micrograph of the MMI-structured optical PDC (OPDC) with five output channels (i.e., K=5), where the two input laser beams take place with interference based on the self-imaging effect [40]. That makes the optical power at each output channel capture different profiles of the phase information Δφ. The relationship between the phase changes induced by varied control voltages applied to a modulator, PM0 [as depicted in Fig. 2(d)] and the optical power of the five channels is shown in Fig. 2(e). By performing threshold decision on the optical power signal of a channel using comparators, the 360-deg phase range is divided into multiple phase intervals. This division occurs due to a change in the binary codewords (0s and 1s) generated by that channel. By performing one decision on each channel’s optical power, the phase range can be divided into up to 2K intervals. The specific phase interval in which the current phase lies can be determined according to the resulting Gray code digital output [41], as shown in Fig. 2(c) for the case of m=10. To discern the phase-encoded information accurately, we introduce an additional static phase bias φb and pre-optimize it along with the decision thresholds for any modulus m2K (see Appendix C).

    Phase manipulation on the integrated chip. (a) Schematic diagram of the experimental setup. PCB, printed circuit board; MVS, multi-channel voltage source; MPM, multi-channel optical power meter. (b) Micrograph of the integrated photonic chip and (c) MMI-structured OPDC, with five output channels. (d) The relationship between the voltage applied to PM0 and the magnitude of consequential phase change, and (e) the corresponding normalized optical power of the output signals. The zero phase reference point is chosen as φb, while the optical power signals from multiple channels were normalized relative to optimized thresholds, ensuring that all thresholds are aligned at 0.5 on the graph. (f) The multi-channel optical power signals along with their corresponding threshold decision results (red or green for zero or one, respectively) at each time slot, and (g) the phase signal reconstructed from the optical power signals.

    Figure 2.Phase manipulation on the integrated chip. (a) Schematic diagram of the experimental setup. PCB, printed circuit board; MVS, multi-channel voltage source; MPM, multi-channel optical power meter. (b) Micrograph of the integrated photonic chip and (c) MMI-structured OPDC, with five output channels. (d) The relationship between the voltage applied to PM0 and the magnitude of consequential phase change, and (e) the corresponding normalized optical power of the output signals. The zero phase reference point is chosen as φb, while the optical power signals from multiple channels were normalized relative to optimized thresholds, ensuring that all thresholds are aligned at 0.5 on the graph. (f) The multi-channel optical power signals along with their corresponding threshold decision results (red or green for zero or one, respectively) at each time slot, and (g) the phase signal reconstructed from the optical power signals.

    Figures 2(d)–2(g) illustrate the experimental process of the loading and extracting phase information for the case of m=10. The voltages corresponding to different discrete phases in φb+ΔφZ are applied to the modulator in varying time slots. By recovering the phase value from the multiple-channel optical power signals in Fig. 2(f), it can be observed that the phase is controlled to be in the middle of the expected phase intervals, as shown in Fig. 2(g). Furthermore, based on the digital results of the threshold decision to the optical power of multiple channels, which was simulated on the host computer in this experiment, we were able to accurately obtain the integer value corresponding to each specific discrete phase for every time slot.

    The methods of phase information loading and extraction described above can be applied to various moduli, affording the flexibility to instantiate the PPMAP units corresponding to different moduli using a single integrated photonic device.

    To demonstrate the superior calculation precision and process of our PPMAP architecture, we conducted experimental demonstrations of a 9-bit three-operand addition operation using moduli {10, 9, 7}. During the experiment, the same physical PPMAP unit depicted in Fig. 2 performed all the calculations involving different moduli, i.e., implementing three PPMAP units utilized one device via time division multiplexing. Seven sets of independent operand samples were selected and each underwent the operation of A+B+C and the bit-width of the binary operands is log2(10×9×7)=9. After converting into RNS format, the original multi-operand addition operation transformed into three corresponding modular arithmetic operations with varying moduli of {10, 9, 7}.

    For each specific modulus, the operations on different sample sets perform consecutively at varying time slots. For implementing the operation of a single set in the phase domain, we simultaneously apply the voltages corresponding to the converted operands to multiple PMs. As shown in Fig. 3(b), we apply a constant voltage to PM0 to assign it the responsibility of providing the bias phase φb. PM1, PM2, and PM3 are also assigned to accurately map operands A, B, and C to their respective phase change amounts. The mapping method for m=9 is illustrated in Fig. 3(f) as an example. The applied voltage signals and their consequent calculation results in the phase domain are illustrated in Figs. 3(c) and 3(d). For a specific time slot, the phase values fall within the phase interval corresponding to the exact computed results, thus achieving the desired calculations in the phase domain. Figure 3(e) displays the integer signals generated by the OPDC, which infers the precise extraction of the RNS-format calculation results from the phase domain.

    Demonstration of PPMAP for algebraic operation (a) Schematic of the PPMAP architecture in the demonstration experiment, which contains three units. (b) The operand or function assigned to each modulator. (c)–(e) The parallel photonics-accelerated modular arithmetic operations in moduli {10, 9, 7}. In different time slots, applying (c) voltage signals corresponding to RNS-formatted operands to the modulators to perform operations on different sample sets, resulting in (d) output phase signals and (e) quantized results obtained from the OPDC. (f) The method of mapping operands to corresponding phases: since φb has been provided by PM0, we directly control the voltage of PM2 to induce a phase change of ΔφX=X·2π/m for operand X and control the voltage of PM1 and PM3 to generate a phase change equal to 2π−ΔφX when mapping X to phase ΔφX.

    Figure 3.Demonstration of PPMAP for algebraic operation (a) Schematic of the PPMAP architecture in the demonstration experiment, which contains three units. (b) The operand or function assigned to each modulator. (c)–(e) The parallel photonics-accelerated modular arithmetic operations in moduli {10, 9, 7}. In different time slots, applying (c) voltage signals corresponding to RNS-formatted operands to the modulators to perform operations on different sample sets, resulting in (d) output phase signals and (e) quantized results obtained from the OPDC. (f) The method of mapping operands to corresponding phases: since φb has been provided by PM0, we directly control the voltage of PM2 to induce a phase change of ΔφX=X·2π/m for operand X and control the voltage of PM1 and PM3 to generate a phase change equal to 2πΔφX when mapping X to phase ΔφX.

    When the result is obtained in RNS format, a non-positional and multiple-radix number system, the PPMAP successfully completes a high-precision arithmetic operation. Although the outcomes are represented differently from the binary format, which is a positional and single-radix number system, they are both valid results of the same operation [34]. To clearly demonstrate that the PPMAP has performed the correct calculation, we convert the RNS format into conventional notations on an electronic processor. In Fig. 3(g), we illustrate the process of the Chinese Remainder Theorem (CRT) algorithm, a basic method for the format conversion. Here, M=Πi=1Nmi, Mi=M/mi, and |·|mi1 represents the modular multiplicative inverse, ensuring that ||Mi|mi1·Mi|mi1. The RNS results from the photonic processor and the binary results from the electronic one are consistent, except for the result in set5. This discrepancy occurs because the computation in set5 produces a result that exceeds the valid range of [0, 630], leading to a numerical overflow that warped 835 to |835|630=205. It should be emphasized that format conversion is not obligatory for performing actual computing tasks. Recent advances in computer architecture have confirmed that all computer components efficiently function within the RNS format [42], and computing architectures built on the RNS format can perform a variety of applications, such as matrix multiplication, fast Fourier transform, and combinatorial optimization, without frequent format conversion [43].

    While our initial proof-of-concept experiment employed thermo-optic modulators on an SOI platform at a moderate speed, it is crucial to highlight the adaptability of the proposed PPMAP architecture across diverse platforms. For instance, when implemented using electro-optic modulators on a commercially available thin-film lithium niobate platform, the PPMAP architecture seamlessly achieves high-precision calculation with bandwidths extending into the range of several tens of GHz [38].

    B. Optical Digital-to-Phase Converters

    The thorough resolution of the “converter challenge” necessitates not only the OPDC but also dedicated electronic-photonic DAC modules that effectively convert digital operands into the analog phase domain. As a DAC essentially performs a weighted summation with fixed weight on a digital operand, denoted as X, and concurrently transfers the result into analog quantity Ψ, X=(xL1xL2x0)2ΨX(i=0L12ixi),where xi=0 or 1, and L represents the bit-width of X, we can leverage the just-demonstrated multi-operand feature of PPMAP architecture to implement the optical digital-to-phase converter (ODPC).

    An ODPC comprises a collection of PMs positioned along the paths of the PPMAP unit. It receives the digital signal of one operand from the parallel bus. In the basic design of the ODPC, it is equipped with L PMs that are controlled by each binary digit xi in the X. Each PM in the ODPC is engineered to exhibit a specific phase response to the voltage amplitude of the digital signal. In particular, when the phase response of the ith PM is designed as Δφi=2π|2i|m/m, the OPDC enables the direct transformation of the L synchronously applied digital signal into the resulting summation in the phase domain: ΔφX=i=0L1xiΔφi=i=0L1xi|2i|m2πm=2πmi=0L1|2ixi|m=2πm|X|m.

    When the digital operand is in RNS format, satisfying X<m, the OPDC produces the phase value Δφ=2πX/m, performing the role of a DAC in the phase domain. Furthermore, when the digital operand X exceeds the modulus m, the ODPC performs the digital-to-analog conversion and the format conversion simultaneously.

    The experimental demonstration of a 4-bit ODPC along with an OPDC, as shown in Fig. 4, validated the successful operation of this photonic conversion interface. During this experiment, digital operands representing 4-bit binary codes ranging from 0 to 15 were sequentially fed into the corresponding phase modulators PM0, PM1, PM2, and PM3. With m=7, all modulators generated phase changes of 0 or 2π|2i|7/7 based on the values of the digital signals (0 or 1) at each bit position i, while PM0 provided an additional bias phase φb simultaneously. The experimental results show that the OPDC successfully performed the modulo-7 remainder calculation on the 4-bit binary operands and accurately transferred the results to analog phase values.

    Feasibility verification of the photonic conversion interface. This experiment demonstrates a PPMAP unit with m=7 consisting of one OPDC and one 4-bit ODPC. Four PMs are controlled by four parallel digital signals and convert the 4-bit operands into phases. In this proof-of-concept experiment, the digital signals are scaled corresponding to the bit position to produce the expected phase change.

    Figure 4.Feasibility verification of the photonic conversion interface. This experiment demonstrates a PPMAP unit with m=7 consisting of one OPDC and one 4-bit ODPC. Four PMs are controlled by four parallel digital signals and convert the 4-bit operands into phases. In this proof-of-concept experiment, the digital signals are scaled corresponding to the bit position to produce the expected phase change.

    In the generalized ODPC design, low-bit-width DACs are utilized to achieve a larger overall bit-width (see Appendix D). As integrated intensity modulators, like microring resonators, exhibit poor linearity, achieving mapping of operands with a given precision into transmissivities necessitates the DAC with higher precision [44]. In contrast, the favorable linearity and the characteristics of phase modulation enable the decomposition of high-precision DAC tasks into G parallel low-bit-width tasks, which can reduce power consumption and delays, from O(eL) to O(G·eL/G) and O(eL/G) [30], endowing the PPMAP for more efficient processing capabilities.

    4. DISCUSSION

    A. Precision and Reliability

    The PPMAP architecture offers superior calculation precision compared to intensity-based photonic computing methods. While intensity-based methods heavily rely on the system’s SNR for calculation precision, PPMAP leverages the RNS theory to distribute the calculation precision across multiple physical implementations. This decomposition substantially reduces the dependence on SNR and enables a relatively high calculation precision. Section 3.A has demonstrated the superior calculation precision of PPMAP (9 bits) compared to intensity-based methods. Now, the focus is on further analyzing the impact of noise on both approaches by conducting simulations to compare their performance in a given 15-bit multiplication scenario.

    Analysis on precision and reliability. (a) The relative calculation error of PPMAP and intensity-based optical computing methods in simulating a 15-bit multiplication task under varying SNR levels. (b) The relative frequency of errors obtained from simulation results for PPMAP, along with its probability model.

    Figure 5.Analysis on precision and reliability. (a) The relative calculation error of PPMAP and intensity-based optical computing methods in simulating a 15-bit multiplication task under varying SNR levels. (b) The relative frequency of errors obtained from simulation results for PPMAP, along with its probability model.

    To further quantify the reliability of the PPMAP architecture, a probability model [Eq. (4)] was developed (detailed in Appendix E): P=1i=1Nj=1Ki12(1+erf(3sin(πKi(j12))22·10SNR20)).

    As shown in Fig. 5(b), the predicted probabilities of system errors at different SNR levels aligned well with the simulated relative frequency of errors. With the addition of computational error correction mechanisms, such as triple modular redundancy (TMR) [47], the error probability P can be reduced to impressively low levels. For example, with an SNR of approximately 34 dB and employing TMR, the system can achieve an error probability as low as 3P22P33×1020 [48], indicating a mean time between failure (MTBF) [49] of around 10 years when continuously operating at a frequency of 100 GHz. This level of reliability is comparable with commercial electronic processors [50], which makes PPMAP a promising candidate for performing high-precision and reliable computation tasks.

    In high-speed scenarios, the PPMAP system demands high-quality modulators and detectors to maintain a sufficiently high SNR at elevated bandwidths. This contributes to reducing the overall error probability and supporting more modulus values m, further enhancing calculation precision.

    B. Computational Efficiency and Scalability

    As analyzed in the previous section, the RNS-based PPMAP architecture overcomes the precision limitations of current photonic computing systems. However, it is essential to analyze whether this high precision is achieved at the cost of significantly reduced computational efficiency and scalability. To address this concern, we conducted a qualitative comparison between the PPMAP architecture and other recent optical computing architectures. We analyzed the computational volume (i.e., the equivalent number of basic operations the system executes per symbol period) and system overhead [including the number of required transmitters (Tx) and receivers (Rx), which respectively consist of modulators/detectors, DACs/ADCs, etc.] as functions of the scale factor D (or D1D2).

    As shown in Table 1, the computational volume and system overhead for optical computing architectures, including PPMAP, increase linearly with the system scale. Compared to other architectures, PPMAP reduces computational volume by up to approximately 50%, while the system overhead does not exceed N times that of the architecture with the minimal overhead. This indicates that, regardless of system scale, PPMAP achieves a linear improvement in precision relative to N with no more than a 2N decrease in computational efficiency. This ensures that the PPMAP architecture maintains computational efficiency comparable to other photonic computing methods while also providing scalability consistent with other integrated solutions.

    Therefore, for general-purpose computing, PPMAP can serve as a high-speed general-purpose processor with high computational efficiency, executing basic arithmetic operations to accelerate various tasks. For application-specific tasks, PPMAP’s scalability allows for tailored designs to further enhance efficiency.

    C. Optoelectronic Signal Converters

    Despite the use of moderate bandwidth in the experimental demonstration, the implementation of ODPC/OPDC addresses the performance challenges encountered by electronic signal converters in high-speed, high-precision optical computing applications. This is attributed to the unique characteristics of the optical system employed by ODPC/OPDC, which decompose a single high-precision ADC or DAC task into multiple parallel low-precision tasks, such as codeword decision or digital modulation.

    This decomposition is pivotal as it breaks the precision-speed tradeoff constraint inherent to electronic signal converters, enabling optoelectronic converters to process high-precision signals at low-precision conversion speeds. Moreover, with recent advancements in the performance of integrated photonic devices [51,52], the current analog bandwidth limitations [29] can be alleviated by incorporating optoelectronic devices. The proposed OPDC/ODPC architectures adopt an integrated design compatible with existing digital systems. By leveraging advanced optoelectronic packaging and fabrication techniques [53], these specially designed optoelectronic signal converters can outperform conventional off-the-shelf electronic signal converters.

    5. CONCLUSION

    We have proposed PPMAP, an integrated photonic computing architecture that performs modular arithmetic based on the optical phase and RNS. The architecture is capable of addition, subtraction, and multiplication and supports multi-operand operations. PPMAP can achieve relatively high precision, and a 9-bit calculation was experimentally demonstrated. Further simulation extends its calculation precision to 15 bits and indicates that its reliability can approach the level of commercial processors under the SNR condition of 34 dB when fortified with error correction mechanisms. Additionally, this architecture harmoniously integrates our proposed ODPC and OPDC, proposing a solution to the long-standing constraints instituted by electronic AD/DA interfaces in photonic computing, and an experimental demonstration confirms the feasibility of this photonic conversion interface. Future research will explore the integration of higher-performance integrated photonic devices into the functional modules of PPMAP to improve the system’s calculation precision in high-speed scenarios.

    APPENDIX A: EXPERIMENTAL METHOD

    We fabricated the photonic integrated chip on a commercially available 180-nm SOI platform via the multi-project wafer (MPW) process. Within the 500-Å-thick titanium nitride layer, we constructed thermo-optic modulators with a length of 0.7 mm, measuring a resistance of 1.5 kΩ and consuming power of 10  mW. We conducted testing on our photonic integrated chip using a single-mode fiber vertical coupling platform [54], optical spectrum analyzer (Yokogawa AQ6370B), and broadband light source (Golight SLED) prior to the optical and electrical packaging process. Measurements indicated that the overall on-chip losses for the five output ports were 8.7, 8.9, 9.2, 8.5, and 9.1 dB, centered around the wavelength of 1535 nm. After the optical packaging process, a pair of single-mode waveguide arrays facilitated the coupling of the laser light into and out of the optical chip through on-chip grating couplers with a total loss of 6 dB.

    We employed a tunable semiconductor laser (Santec TSL 570) to provide the incident CW laser at the wavelength of 1535 nm, and a multi-port optical power meter (Santec MPM 210) to measure the multi-channel optical power signals. After the electrical packaging of the integrated photonic chip, metal wires connect the on-chip electrodes to the FPC (flexible printed circuit) adapter board. We employed a multi-channel programmable voltage source (nicslab XDAC-40MUB-R4G8) to control the applied voltages on the electrodes. The experimental equipment and the host computer were connected via GPIB interface and Ethernet cables, which facilitated the transmission of control signals and the acquisition of experimental data. In-house software was employed to perform real-time monitoring of the optical power signals from each port, which allowed one to evaluate system stability and monitor the polarization state of the incident light for controlling the coupling efficiency.

    APPENDIX B: PHASE RECOVERY FROM MULTI-CHANNEL OPTICAL POWER SIGNALS

    To illustrate the phase difference between two paths, φΔφ+φ0, from multi-channel optical power signals I, it is necessary to have knowledge of the corresponding multi-channel optical power states at different values of φ beforehand. That requires utilizing voltage as an intermediate quantity and obtaining two mappings: (1) mapping between the applied voltage to PM and resulting I, and (2) mapping between the voltage and its produced φ.

    An experimental procedure was conducted to access the relationship between I and φ. The voltage applied to PM0 was varied from 0 V to 10 V with a voltage increment of 0.01 V. Then the first mapping is obtained by recording the multi-channel optical power values at different voltage levels. The optical power for the k-th channel and produced phase variance at voltage index i are denoted as Ik,i and Δφi. Moreover, the establishment of the second mapping is based on the prior knowledge of the photonic device. Theoretically, the variation of output optical power in 2×K MMI follows a cosine function with respect to the phase [37]: Ik,i=PchipK·(1+cosφik)=PchipK·(1+cos(Δφi+φ0k)),φ0k={φ0k2·2πKfor  even    kφ0+k12·2πKfor  odd    k,where Pchip is the power of the input laser coupled to the chip. The multiple-channel signals are analyzed individually. First, The signal from the k-th channel, denoted as Ik, is scaled to [1,1] to conform to the input range of the arccosine function. As the arccosine function is not monotonically changing with respect to the phase Δφ, the relationship of φik is estimated based on the absolute value of the variation of the arccosine function: I^k=2×Ikmax(Ik)1,φik={φikarccosI^k,ii=0φi1k+|φi+1kφik|i0.

    As that accumulation-based approach introduces errors that accumulate over time, a further process is employed to enhance the precision. This process utilizes the ratio of normalized intensity signals with adjacent voltages, denoted as Rk,i=I^k,i+1/I^k,i: φik=arctan(cos(δik)Rk,isin(δik))+Nπ,φ=φ0+Δφi¯=φ0+1Kk=1Kφikφ0k,where δik=φi+1kφik is the finite difference of phase with respect to voltage, which along with the integer value N is determined by the previous estimation obtained through Eq. (B2). This process provides finer estimates within the π interval based on the coarse-grained results provided in the previous step, and averages the results obtained from each channel to obtain the final voltage-phase relationship.

    With the obtained mapping between I and φ, the phase at a given moment can be determined from the measured optical power. Expressly, the current φ is specified as the phase value corresponding to the previously recorded optical power with the minimum Euclidean norm compared to the currently measured optical power. It is worth noting that there are waveform spikes observed in the obtained phase signals at slot transitions. These spikes are due to synchronization issues related to the multi-channel voltage sources and multi-channel optical detectors in our experimental setup, as well as the inherent folding of phase signals when exceeding 2π. However, these spikes do not affect the accurate operation of the PPMAP, given the high-precision signal delay calibration and the use of synchronous clocks.

    APPENDIX C: OPDC THRESHOLD OPTIMIZATION

    Given the multi-channel optical power signals I(φ)=[I1(φ),I2(φ),,IK(φ)], which were obtained from the experimental procedure in Appendix B, and the thresholds of comparators for the K channels t=[t1,,tK], the intersection of tk with its corresponding IK(φ) determines two boundary points in the 360-deg phase range. The sorted boundaries of all channels are denoted as pexpt=pexpt(I,t), which represents the actual phase boundary points. And in the ideal scenario, the phase boundaries pideal should be shifted by half-interval 2π2m to the discrete phase points φ0+ΔφZ, where the total static phase φ0 is adjusted by applying an extra bias phase φb. Thus the ideal phase boundaries can be described as pideali=φ0+φb+i·2πm+πm.

    To accurately extract the phase information and reduce the probability of incorrect quantization operation, an optimization procedure was performed to find the optimal values of φb and t that minimize the discrepancy between pideal(φb) and pexpt(I,t). This optimization aims to enhance the robustness of OPDCs by placing the most error-prone phase point in a less vulnerable position. The optimization problem can be formulated as follows: [φb,t]I=argminφb,tmaxi|pideali(φb)pexpti(I,t)|.

    The optimization objective is to minimize the loss pideal(φb)pexpt(I,t). A heuristic optimization method, the differential evolution (DE) algorithm [55], is employed. The DE solver provides multiple trial values for t at each iteration. A search algorithm was performed for each trial threshold to determine the minimized loss under a specific search range of φb with a given precision. Then those minimum values were passed back to the DE solver for further exploration in the next iteration.

    In the practical experiments, the DE solver had a population size of 20 and was run for 40 iterations. When modulus m=2K, the intervals tend to be evenly divided. When modulus m<2K, there is a need to merge the remaining phase intervals to unify them into uniformly divided intervals of m. This merging process involves mapping multiple different codewords to the same integer. Figures 6(a) and 6(b) display the optimized results for m=9 and m=7. It is worth noting that for m=7, it is possible to utilize information from five optical power signals to quantify the phase information. However, signals from only channels 1, 3, 4, and 5 were employed to reduce the length of the codewords in our experimental setup. The experimental results demonstrate that this approach still yields satisfactory outcomes.

    The relationship between multi-channel optical power values and normalized optical power values with respect to phase change, as well as the optimized thresholds for m=9 (a) and 7 (b). The normalization method involved scaling the original optical power values for each channel such that the optimized thresholds were uniformly scaled to 0.5.

    Figure 6.The relationship between multi-channel optical power values and normalized optical power values with respect to phase change, as well as the optimized thresholds for m=9 (a) and 7 (b). The normalization method involved scaling the original optical power values for each channel such that the optimized thresholds were uniformly scaled to 0.5.

    APPENDIX D: GENERALIZED DESIGN OF ODPC

    The incorporation of high-linearity phase modulators, e.g., lithium niobate modulators [56], enhances the capabilities of the ODPC by allowing it to be combined with low-bit-width DACs, thus achieving a larger overall bit-width in the limitation of the given number of PMs. To accomplish this, the L-bit digital operand is partitioned into G groups corresponding to G DACs, where the i-th DAC has a bit-width of Wi (with Σi=1GWi=L). The phase response of the i-th PM to the least significant bit (LSB) voltage of the i-th DAC is set to Δφi=2π|2wi|m/m, where wi=Σj=1i1Wj represents the bit position of the LSB of the i-th DAC in the original L-bit operand. By mapping the digital operand X using the DACs and the corresponding PMs, the resulting ΔφX is as shown in Eq. (3): ΔφX=i=1Gj=0Wi12jxjΔφi=2πmi=1Gj=0Wi12wi|2jxj|m=2πmj=0L1|2jxj|m=2πm|X|m.

    The basic design illustrated in Fig. 4 is the specific case in which Wi is limited to one. The PPMAP unit, which incorporates the generalized ODPCs, is depicted in Fig. 7(a).

    (a) Structure diagram of a PPMAP unit that adopts the generalized design of OPDC. (b) Example of the grouping strategy of input digital signals, which is based on the bit weight |2i|m.

    Figure 7.(a) Structure diagram of a PPMAP unit that adopts the generalized design of OPDC. (b) Example of the grouping strategy of input digital signals, which is based on the bit weight |2i|m.

    For a given electro-optic coefficient of the electro-optic material in the modulator and a given modulus, m, the objective of the grouping strategy is to minimize the magnitude of the modulation phase. According to Eq. (D1), the LSB phase response of the i-th modulator for the DAC should be 2π|2wi|m/m. To achieve this, we need to choose a grouping method for the L-bit digital signal such that |2wi|m is as close to zero as possible. Being close to zero implies being in proximity to either 0 or m because a digital value |2i|m close to m can be achieved by modulating the phase of the digital value m|2i|m on the opposite waveguide. Figure 7(b) illustrates examples of bit grouping for moduli of 15, 19, 21, 30, 31, and 33, where the maximum modulation phase of any PM, (2Wi1)×min(|2i|m,m|2i|m), is designed to around 2π. We have observed that the bit weights |2i|m can be close to zero by employing an appropriate grouping strategy.

    APPENDIX E: ANALYSIS OF CALCULATION PRECISION AND RELIABILITY

    Calculation precision analysis. In Fig. 5(a), the relative error (RE) is defined as the average absolute value to error over the dynamic range. The RE of a calculation with n-bit precision is defined as the n-bit relative quantization error of the calculation result: REn=minaRaδ+a|x|·1δdxδ·(2n1)=12n+242(n+2).

    And the theoretical RE of the intensity-based method under Gaussian noise with uniform intensity is given by REI=|x|f(x)dxAmax=20x·1σ2πex22σ2dxAmax=22π(σAmax),where f(x) represents the probability density function of the normal distribution with standard deviation σ, and Amax is the maximum signal amplitude. When Eqs. (E1) and (E2) are combined, we can get the relationship between the calculation precision expressed in the form of the bit length, n, and the relative standard deviation σrel: n=log2(σAmax)[log2(22π)+2]log2σrel1.67.

    The ratio of σ and Amax, i.e., σrel, is determined by its relationship with SNR and peak-to-average power ratio (PAPR), i.e., SNR=20log10(Amax/σ)PAPR. Consequently, the relationships of RE between SNR and precision n, is inferred: log10REI=log10(22π)PAPR20SNR20=(n+2)·log102.

    The logarithm of the RE decreases linearly with SNR and n. Assuming the calculation results follow the uniform distribution, it has a PAPR of 4.77 dB (i.e., Amax=3Arms), and Eq. (E4) can be simplified to SNR6.02·n+5.31.

    This equation shows that as n increases, the required SNR increases linearly with a slope of 6.02, corresponding to the numerical simulation shown in Fig. 5(a). That linear term indicates that as the calculation precision increases, the required SNR will reach an unattainable level and imply an exponential increase in energy consumption [22].

    It is worth noting that the precision metrics in other referenced papers are not consistently defined. For example, in Ref. [25], precision is estimated by dividing SNR by 6.02, resulting in a relatively rough estimation. Discrepancies exist in Refs. [23,24], where the former estimates 5-bit precision based on a σrel of 0.008, and the latter estimates 6-bit precision based on a σrel of 0.02, resulting in an apparent contradiction. Therefore, we use Eqs. (E3) and (E5) to convert other indicators into a unified form of relative quantization error for quantitative analysis. It is also worth to note that some works claim a precision of up to 9 bits [44,46], but this refers to weight control precision rather than calculation precision, so they were not included in the comparison presented in Table 1.

    Reliability analysis. For the PPMAP architecture, the probability of the occurrence of errors, denoted as P, is equal to the probability that an error occurs in any of the N PPMAP units. The probability of the i-th PPMAP unit not encountering errors, denoted as Piunit, is determined by accurately achieving threshold decisions for the output optical power on all the Ki channels. The success rate of the j-th channel is represented as Pi,jchan: P=1i=1NPiunit=1i=1Nj=1KiPi,jchan.

    In the simulation, the chosen modulus, i.e., 7, 11, 19, 23, is denoted as m. In this case, an index-sum modular multiplication operation is implemented via a modular addition operation with moduli m1 [57], equal to twice the designed Ki in each PPMAP unit. That makes for all discrete phase points in ΔφZ, the relative distance, |dj|/Amax, between their corresponding ideal output optical power and decision threshold in all Ki channels given by [41] {|dj|Amax|jZ,1jKi}={12sin(πKi(j12))|jZ,1jKi}.

    Considering the relative variance noise Amax/σ=3Arms/σ=3·10SNR20, the success rate for the j-th channel in the i-th unit is Pi,jchan=Φ(32sin((j12)πKi)·10SNR20),where Φ(x)  =  (1+erf(x/2))/2, i.e., the cumulative distribution function of standard Gaussian distribution N(0,1). Then considering the commutative law of the continued multiplication in Eq. (E6), Πj=1KiPi,jchan=Πj=1KiPi,jchan, the error model described in Eq. (4) can be derived.

    The deviation between the simulated relative frequency and the probability model around 31 dB in Fig. 5 is only a few samples out of the 40,000 experienced errors near an SNR of 31 dB, which does not satisfy the prerequisites of the law of large numbers. Apart from this specific range, the relative frequency obtained in our simulations aligns well with our probability model.

    To further enhance the reliability, optimizations to the codeword decision scheme can be implemented, along with the incorporation of computational error correction codes, such as exploring the use of redundant RNS (RRNS) error-correcting codes that can improve reliability by adding extra redundant PPMAP units. These methods can reduce or tolerate errors from individual PPMAP units, allowing the PPMAP system to support a greater number of units and thereby improve calculation precision.

    APPENDIX F: EXTENSIBILITY OF OPDC

    In this section, we introduce the extended design of the OPDC, focusing on two key aspects: first, the extensibility of operating speed to accommodate high-speed scenarios by demonstrating the quantization on a lithium niobate platform; second, the extensibility of the number of ports beyond the five-port MMI configuration used in our experiments.

    Extensibility in speed. To verify the operational speed of OPDC, we conducted experiments using a thin-film lithium niobate (LN) photonic integrated chip. Figure 8(a) illustrates the schematic representation of the PPMAP structure on the LN platform, which, in essence, mirrors that depicted in Fig. 2(b) but includes multiple (M) lithium niobate phase modulators (LN-PMs) that simultaneously occupy the two waveguides. We demonstrated the loading of high-precision, RNS-formatted operands using moduli of 7 and 10. In this experiment, the LN-PPMAP utilized was configured with M=1. The chip used in the experiment is consistent with those in Refs. [38,58]. We attained a half-wave voltage of 2.8 V at the central wavelength of 1529 nm and an overall bandwidth of 5 GHz after electrical packaging. A voltage source provided a static bias for PPMAP, while an AWG generated a 1 Gbaud PAM4 signal that was amplified and fed into the LN-PM as operands. The AWG output amplitudes were adjusted to deliver 2.40 V and 1.68 V to the modulators, ensuring that the PAM signal’s output codewords matched the PPMAP’s phase levels.

    (a) Schematic diagram of PPMAP chip on lithium niobate platform. (b) Experimental setup for the performance demonstration experiment. AWG, arbitrary waveform generator; OCS, oscilloscope; AMP, voltage amplifier; BPF, band-pass filter. (c) Waveforms recovered in the phase domain, displaying the computational results achieved at 1 GHz under moduli 7 and 10. (d) For m=10, the output intensity waveforms of OPDC.

    Figure 8.(a) Schematic diagram of PPMAP chip on lithium niobate platform. (b) Experimental setup for the performance demonstration experiment. AWG, arbitrary waveform generator; OCS, oscilloscope; AMP, voltage amplifier; BPF, band-pass filter. (c) Waveforms recovered in the phase domain, displaying the computational results achieved at 1 GHz under moduli 7 and 10. (d) For m=10, the output intensity waveforms of OPDC.

    (a) Structural rendering, (b) micrograph, and (c) the relationship between optical power and phase change of each channel, along with the phase intervals (taking m=18 as an example), for the nine-channel OPDC.

    Figure 9.(a) Structural rendering, (b) micrograph, and (c) the relationship between optical power and phase change of each channel, along with the phase intervals (taking m=18 as an example), for the nine-channel OPDC.

    Besides horizontal expansion by adding the number of output channels, OPDC can be vertically expanded with a cascaded structure [59], which possesses relatively low wavelength sensitivity. On the other hand, as the number of ports increases, manufacturing and design errors can lead to a reduced working wavelength range for the OPDC, a problem that is more pronounced on the lithium niobate platform [60]. To mitigate this effect, a compact MMI with sub-wavelength gratings can be employed [61], which has proven effective in counteracting the impacts of etching and substrate non-uniformity and is a focus of our ongoing parallel research.

    APPENDIX G: INDEX-SUM MULTIPLICATION

    In this section, we briefly introduce the process of index-sum multiplication. Since the overall multiplication is equivalent to independent modular multiplication within each modulus for numbers in RNS format, we only discuss the multiplication within a single modulus here. Index-sum multiplication necessitates the modulus m to be representable as pN, where p is a prime number and N is an integer. For illustration, let us consider the modulus m=19. In this case, any integer A0 can be expressed as A=|gα|m.

    Here, α is termed the index, with α[0,m2]. This expression demonstrates a crucial correspondence. Through the base number g, each A is mapped to a value of α in the “index field”. The number g is known as the primitive root of modulus m [33]. Based on the properties of primitive roots, each integer A corresponds uniquely and exclusively to a specific α, and each α is also uniquely mapped to a different A.

    Index-sum multiplication requires a precomputed lookup table (LUT) of length m to record the correspondence between A and α. Since RNS decomposes the modulus into very small numbers, the LUT’s size is minimal, often less than 10 bytes. Let us consider a more specific example where we have chosen the prime root of m=19 as g=2.

    Taking |8×9|19 (equal to 15) as an example, when implementing |8×9|19, we look up Table 2 to obtain the mapping of A and α: 8=|23|19,9=|28|19,that is, A=8α=3 and A=9α=8. Subsequently, by adding the indices of the two operands modulo m1, |3+8|191=11,we obtain the multiplication result in the index field. Benefited from the uniqueness of the mapping between A and α, we can subsequently map the result back to the integer field (inverse lookup) to get the final results: α=11A=15. It is worth noting that in Table 2, there is no corresponding index for A=0. This is because, in index-sum multiplication, any number multiplied by zero is pre-determined as zero without further executing other steps.

    Mapping Table of GF(19) with a Primitive Root of 2

    A0123456789101112131415161718
    α01132161463817121557114109

    In summary, the index-sum multiplication process involves mapping all operands into the index field, performing addition operations on indices instead of multiplication, and then mapping the result back from the index field.

    APPENDIX H: EVALUATION ON POWER CONSUMPTION

    To demonstrate the advantages and potential of PPMAP, we estimated its power consumption and compared it with electronic processors. As depicted in Fig. 10, which illustrates the relationship between transistor power consumption and operating frequency for a 16-bit multiplier under different process technologies [62], it is evident that while increasing the power supply voltage can boost processor performance, it comes at the expense of exponentially higher power consumption. This results in a practical limitation to the frequency achievable in electronic processors, typically constrained to the GHz range.

    Performance comparison of electronic processor under various processes and PPMAP. The black dots represent performance data for different VDD values under various processes, derived from the scaling method and data presented in Refs. [62,64]. The solid line corresponds to the fitting and extrapolation of the former. The power consumption model of PPMAP assumes a fixed capacitance, with parameters obtained from Refs. [7,63], while PPMAP* provides an estimate under the assumption of a one-order-of-magnitude improvement in the energy efficiency of optoelectronic devices.

    Figure 10.Performance comparison of electronic processor under various processes and PPMAP. The black dots represent performance data for different VDD values under various processes, derived from the scaling method and data presented in Refs. [62,64]. The solid line corresponds to the fitting and extrapolation of the former. The power consumption model of PPMAP assumes a fixed capacitance, with parameters obtained from Refs. [7,63], while PPMAP* provides an estimate under the assumption of a one-order-of-magnitude improvement in the energy efficiency of optoelectronic devices.

    In contrast, optoelectronic systems, including technologies like PPMAP, can provide high bandwidth at relatively low power consumption, enabling operating frequencies of tens of GHz. We present a simplified model illustrating the efficiency of a 15-bit multiplication operation using PPMAP, consistent with Section 4. The parameters in this simplified model are derived from published works, encompassing the power consumption of modulators [7], lasers, and detectors [63]. The overhead of OPDC and ODPC is included in the power consumption model by converting the costs into the equivalent number of optoelectronic components required. Figure 10 demonstrates the potential of PPMAP to surpass the frequency limitations of CMOS processors with lower power consumption, especially at high frequencies.

    It is important to note that Fig. 10 does not account for interconnect wires for electronic processors, which constitutes a significant portion of energy consumption and delays in advanced processes [4,65]. Considering the benefits brought by optical interconnects, the performance of PPMAP is expected to be even more pronounced. Furthermore, PPMAP supports direct multi-operand operations efficiently, a capability that is challenging to achieve in digital circuits. A more detailed and comprehensive energy consumption analysis is part of our future research.

    [1] J. Hennessy, D. Patterson. A new golden age for computer architecture: domain-specific hardware/software co-design, enhanced security, open instruction sets, and agile chip development. ACM/IEEE 45th Annual International Symposium on Computer Architecture (ISCA), 27-29(2018).

    [2] H. Esmaeilzadeh, E. Blem, R. St. Amant. Dark silicon and the end of multicore scaling. Proceedings of the 38th Annual International Symposium on Computer Architecture, 365-376(2011).

    [27] R. J. Van de Plassche. CMOS Integrated Analog-to-Digital and Digital-to-Analog Converters, 742(2013).

    [33] W. K. Nicholson. Introduction to Abstract Algebra(2012).

    [34] A. R. Omondi, A. B. Premkumar. Residue Number Systems: Theory and Implementation, 2(2007).

    [35] Y. Wu, H. Guo, B. Zhang. Calculating with phase opens up the high-precision and high-reconfigurability integrated photonic computing. Optical Fiber Communication Conference, M3J–4(2023).

    [36] H. Park, S. Kim. Hardware accelerator systems for artificial intelligence and machine learning. Advances in Computers, 122, 51-95(2021).

    [42] S. Srikanth, P. G. Rabbat, E. R. Hein. Memory system design for ultra low power, computationally error resilient processor microarchitectures. IEEE International Symposium on High Performance Computer Architecture (HPCA), 696-709(2018).

    [47] P. H. Hochschild, P. Turner, J. C. Mogul. Cores that don’t count. Proceedings of the Workshop on Hot Topics in Operating Systems, 9-16(2021).

    [49] J. Lienig, H. Bruemmer, J. Lienig. Reliability analysis. Fundamentals of Electronic Systems Design, 45-73(2017).

    [50] A. Sari, M. Psarakis. Scrubbing-based SEU mitigation approach for systems-on-programmable-chips. International Conference on Field-Programmable Technology, 1-8(2011).

    [61] L. Xu, Y. Wang, D. Patel. Ultra-broadband and ultra-compact optical 90 hybrid based on 2 × 4 MMI coupler with subwavelength gratings on silicon-on-insulator. Optical Fiber Communication Conference, M3I–7(2018).

    [65] A. Khursheed, K. Khare. Nano Interconnects: Device Physics, Modeling and Simulation(2021).

    Tools

    Get Citation

    Copy Citation Text

    Yuepeng Wu, Hongxiang Guo, Bowen Zhang, Jifang Qiu, Zhisheng Yang, Jian Wu, "Integrated photonic modular arithmetic processor," Photonics Res. 12, 2676 (2024)

    Download Citation

    EndNote(RIS)BibTexPlain Text
    Save article for my favorites
    Paper Information

    Category: Integrated Optics

    Received: Apr. 22, 2024

    Accepted: Sep. 2, 2024

    Published Online: Nov. 1, 2024

    The Author Email: Hongxiang Guo (hxguo@bupt.edu.cn), Jian Wu (jianwu@bupt.edu.cn)

    DOI:10.1364/PRJ.527762

    Topics