{'fraction_non_alphanumeric': 0.02898000506200962, 'fraction_numerical': 0.015945330296127564, 'mean_word_length': 4.465848175687359, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Small-Signal Analysis of All-Si Microring Resonator Photodiode

All-silicon microring resonator photodiodes are attractive for silicon photonics integrated circuits due to their compactness, wavelength division multiplexing ability, and the absence of germanium growth. To analyze and evaluate the performance of the microring photodiode, we derived closed-form expression of the response transfer function with both electrical and optical behavior included, using a small-signal analysis. The thermo-optic nonlinearity resulting from optical loss and ohmic heating was simulated and considered in the model. The predicted response achieved close agreement with the experiment results, which provides an intuitive understanding of device performance. We analytically investigated the responsivity–bandwidth product and demonstrated that the performance is superior when the detuning frequency is zero.


Introduction
The silicon photonics platform has been considered the most attractive solution for high-speed network and artificial intelligence applications due to its low cost and high miniaturization [1,2]. Because the cutoff wavelength of bulk Si is only around 1.1 µm, the photodetection at telecommunication wavelengths is usually achieved by the heteroepitaxial growth of germanium (Ge) or heterogeneous integration of III-V materials [3][4][5][6]. Both methods require added material costs and extra process steps. All-Si photodetection based on defect state absorption has been demonstrated to achieve high quantum efficiency and high bandwidth [7,8]. However, the method needs special steps, such as defects implantation, which are not compatible with the standard process in foundries. Si photodetectors (PDs) based on two-photon absorption (TPA) and photon-assisted tunnelling (PAT) can absorb sub-bandgap wavelengths [9]. Applying a high bias voltage across the PN junction, the spatial distance between the conduction and valence band will be shortened due to the strong electric field. The electrons in the valence band can possibly tunnel into the conduction band [10]. However, the mechanisms are weak and PDs suffer from low responsivity. Microring resonator (MRR) PD benefits from optical resonant enhancement effect is a promising solution to improve responsivity [11]. Intel has reported a 112 Gb/s (56 GBaud PAM4) all-Si MRR PD with a responsivity of 0.23 A/W [12].
Thus far, there is not much dynamic analysis work to quantify the performance of the MRR PD. Some steady-state and dynamical discussions on MRR modulators show that the modulation rate and efficiency are limited by the resonator quality factor (Q) [13][14][15][16][17][18]. For the MRR PDs, Q value limits the rate at which the optical power can be injected in the ring and absorbed, which determines the modulation speed. At the same time, the resonant enhancement effect, and thus the responsivity, are related to the Q value. Bandwidth improvement always comes at the expense of a reduced responsivity and it is important to weigh the tradeoff in design [19]. The performance of MRR PD also depends on the input wavelength offset away from the resonance wavelength. Optical power losses and electrical power consumption within the small volume usually lead to heat accumulation as well as thermally-induced nonlinearities, which make the model more complicated because thermal power distribution in MRR can dynamically change the resonance wavelength [20][21][22][23]. For high-power applications, the thermal nonlinearity of the MRR PD cannot be ignored. To completely understand the performance of MRR PD, a fully dynamical model considering thermal nonlinear effects is required.
In this paper, we provide a small-signal analysis of all-Si MRR PD in terms of frequency response and responsivity. The perturbation treatment is applied to derive the expression of the response based on the MRR energy coupling model. We take thermal nonlinear effects into consideration by simulating the temperature variance and calculating the temperature induced resonance shift. The results predicted from the model match well with the measured data for our high-performance Ge-free MRR PD with the simple Si PN junction. The responsivity-bandwidth, which is an important parameter for PD [24,25], is also analyzed versus the detuning frequency. The verified model helps us understand the MRR PD and is appealing to future all-Si PD optimization.

Analytic Model
In this section, we aim to derive the small-signal response of the MRR PD. As shown in Figure 1, the transfer function of MRR PD can be divided into optical response H O (ω m ) = P r P in and electrical response H E (ω m ) = I L Is , where P r and P in are optical power in cavity and input optical power and I L and I s are load current and source current, respectively. P r is related to I s by P r = I s R , where R is the responsivity of the MRR at critical coupling at resonant wavelength. The input optical signal consists of a small sinusoidal signal of cos(ω m t) and an optical carrier of cos(ω in t), where ω m and ω in are modulation angular frequency and optical carrier frequency. Under stimulus, the RF part of the output signal oscillates at ω m .
Electronics 2022, 11, x FOR PEER REVIEW 2 of 12 width improvement always comes at the expense of a reduced responsivity and it is important to weigh the tradeoff in design [19]. The performance of MRR PD also depends on the input wavelength offset away from the resonance wavelength. Optical power losses and electrical power consumption within the small volume usually lead to heat accumulation as well as thermally-induced nonlinearities, which make the model more complicated because thermal power distribution in MRR can dynamically change the resonance wavelength [20][21][22][23]. For high-power applications, the thermal nonlinearity of the MRR PD cannot be ignored. To completely understand the performance of MRR PD, a fully dynamical model considering thermal nonlinear effects is required.
In this paper, we provide a small-signal analysis of all-Si MRR PD in terms of frequency response and responsivity. The perturbation treatment is applied to derive the expression of the response based on the MRR energy coupling model. We take thermal nonlinear effects into consideration by simulating the temperature variance and calculating the temperature induced resonance shift. The results predicted from the model match well with the measured data for our high-performance Ge-free MRR PD with the simple Si PN junction. The responsivity-bandwidth, which is an important parameter for PD [24,25], is also analyzed versus the detuning frequency. The verified model helps us understand the MRR PD and is appealing to future all-Si PD optimization.

Analytic Model
In this section, we aim to derive the small-signal response of the MRR PD. As shown in Figure 1, the transfer function of MRR PD can be divided into optical response  In order to obtain the optical response of the PD, we started with the rate equations as [26][27][28]: where ( ) is the energy amplitude travelling inside the resonator, is the angular resonant frequency of the ring, ( ) is the input amplitude of waveguide mode and ( ) is the modulated amplitude. is the amplitude decay rate which consists of a decay rate due to the attenuation in the resonator and the ring to bus waveguide coupling . The Q factor and photon lifetime can be obtained by = = . The energy coupling coefficient is related to the decay rate and field coupling coefficient by In order to obtain the optical response of the PD, we started with the rate equations as [26][27][28]: where a(t) is the energy amplitude travelling inside the resonator, ω r is the angular resonant frequency of the ring, E(t) is the input amplitude of waveguide mode and δE(t) is the modulated amplitude. γ r is the amplitude decay rate which consists of a decay rate due to the attenuation in the resonator γ l and the ring to bus waveguide coupling γ e . The Q factor and photon lifetime τ p can be obtained by Q = 2γ r ω r = 1 ω r τ p . The energy coupling coefficient µ is related to the decay rate γ e and field coupling coefficient k by µ 2 = 2γ e = k 2 T rt , where T rt is the round-trip transit time. For the steady state, the input field is written as E(t) = E 0 e jω in t and the energy amplitude is a(t) = a 0 e jω in t , where a 0 is: Electronics 2022, 11, 183 3 of 12 ∆ω = ω in − ω r = 2π∆ f is the angular frequency detuning of the input carrier frequency from the resonance frequency.
With the sinusoidal modulated signal, the energy amplitude will experience a small amount of change δa(t). To solve the equation, the δa(t) is assumed to be the product of a slowly varying envelope and an exponential item: δa(t) = δae (jω r −γ r )t . The modulated signal can be written with upper and lower sidebands as: δE(t) = 0.5 × δE 0 e j(ω in +ω m )t + e j(ω in −ω m )t . The output signal consists of beat signals generated between the carrier and sidebands and between two sidebands. Because the modulation depth is small, the high order harmonics can be neglected. The δa(t) can be written as: The output signal depends on the absorption coefficient due to PAT and optical power in the resonator. The transfer function can be written as: The transfer function can be divided as the DC part: . By substituting Equations (2) and (3), the DC and RF part of H O (ω m ) can be written as: where ε = j∆ω + γ r . The equivalent circuit of PD is shown in Figure 1. C j and R s represent series capacitance and resistance from the Si PN junction. C ox and R sub represent the capacitance from oxide layer and resistance from the high-resistivity substrate. C P , R P , and L P are parasitic components. All the parameters can be extracted from the S11 fitting and the electrical response can be modeled with equivalent resistance R eq and capacitance C eq as: The small-signal transfer function of the MRR PD is the product of H E (ω m ), H O (ω m ), and R. The MRR PD is a fourth-order system with a complex-conjugate pole pair, two poles from the RLC circuit and a zero. The pole pair and zero depend on the Q value and laser detuning frequency, which determine the speed of the PD. The product of H O-DC and R is the responsivity of the MRR PD.
The responsivity-bandwidth product, which is defined as the product of the 3 dB bandwidth and the responsivity, is an important figure of merit to characterize the performance of the PD. Figure 2 shows the relationship between the responsivity-bandwidth product versus the detuning frequency and decay rate. During the calculation, the responsivity was normalized as 1 A/W when γ r = 2π × 10G at zero detuning. The detailed normalized responsivity and bandwidth at zero detuning are shown as star markers in the Figure 2. It can be observed that the responsivity-bandwidth product decreases with the laser detuning frequency for a fixed decay rate. Therefore, the PD can achieve the best performance when the optical carrier frequency is set as the same as the resonance frequency. When ∆ω = 0, the 3 dB bandwidth for MRR equals to γ r 2π if we neglect the RC time limitation. It should be noted that critical coupled resonator is favorable for the MRR PD responsivity only if the PAT effect is much higher than the loss or is proportional to the power decay rate in cavity. The experimental relationship between frequency response and detuning frequency will be discussed in the following section. The bandwidth can be improved with a higher decay Electronics 2022, 11, 183 4 of 12 rate at a cost of lower responsivity. To obtain the high output signal, the MRR PD usually needs to work at a high optical power density. The heating power will accumulate in the small volume, which will increase the temperature of the structure and shift the resonance frequency [20]. To analyze the PD based on Equations (4) and (5), the thermal-induced shift needs to be calculated. laser detuning frequency for a fixed decay rate. Therefore, the PD can achieve the best performance when the optical carrier frequency is set as the same as the resonance frequency. When ∆ = 0, the 3 dB bandwidth for MRR equals to if we neglect the RC time limitation. It should be noted that critical coupled resonator is favorable for the MRR PD responsivity only if the PAT effect is much higher than the loss or is proportional to the power decay rate in cavity. The experimental relationship between frequency response and detuning frequency will be discussed in the following section. The bandwidth can be improved with a higher decay rate at a cost of lower responsivity. To obtain the high output signal, the MRR PD usually needs to work at a high optical power density. The heating power will accumulate in the small volume, which will increase the temperature of the structure and shift the resonance frequency [20]. To analyze the PD based on Equations (4) and (5), the thermal-induced shift needs to be calculated. The optical energy in the MRR can be dissipated due to mode mismatch, propagation loss, bend loss and free-carrier absorption (FCA) [29]. For the doped Si MRR, the FCA usually dominates the loss [30,31]. The heating power resulted from the optical loss at low bias voltage can be described as: where is the loss rate due to FCA and is the round-trip field attenuation factor, which can be obtained from the intrinsic loss induced decay rate = . Another dissipated power source is the ohmic heating, which can be determined by [32]: , represents the DC electrical power delivered to the PD and , represents the RF power delivered to the load.
, is approximately equal to the product of the current and applied voltage. Because the voltage is mainly applied at the PN junction in the waveguide, the heat will accumulate in this area, where it largely overlaps with the mode area. Due to limited thermal conductivity, the heat will change the temperature of the structure and consequently change the refractive index as well as the resonance wavelength. The change in resonant wavelength as a function of the waveguide temperature variation can be expressed as [33,34]: The optical energy in the MRR can be dissipated due to mode mismatch, propagation loss, bend loss and free-carrier absorption (FCA) [29]. For the doped Si MRR, the FCA usually dominates the loss [30,31]. The heating power resulted from the optical loss at low bias voltage can be described as: where α FCA is the loss rate due to FCA and a rt is the round-trip field attenuation factor, which can be obtained from the intrinsic loss induced decay rate a rt = e −γ l T rt . Another dissipated power source is the ohmic heating, which can be determined by [32]: P DC,in represents the DC electrical power delivered to the PD and P RF,out represents the RF power delivered to the load. P DC,in is approximately equal to the product of the current and applied voltage. Because the voltage is mainly applied at the PN junction in the waveguide, the heat will accumulate in this area, where it largely overlaps with the mode area. Due to limited thermal conductivity, the heat will change the temperature of the structure and consequently change the refractive index as well as the resonance wavelength. The change in resonant wavelength as a function of the waveguide temperature variation can be expressed as [33,34]: where ∈ and n g are the thermal expansion coefficient and group index of the MRR. dn dT is the index thermo-optic coefficient of silicon, which is~1.94 × 10 −4 K −1 at 1310 nm.
To analyze the small-signal response, including the nonlinear thermal optical effect of a MRR PD, we first need to obtain the static optical power in the MRR by calculating the resonator optical enhancement. The heating power and generated temperature variation can then be used to calculate the resonator frequency shift as well as the detuning frequency. Finally, Equations (5) and (6) will tell us the frequency response.

Experimental Verification
To verify the small-signal model in the previous section, we designed an all-Si MRR PD at Advanced Micro Foundry (AMF). The cross section of the MRR PD and a microscope image of the fabricated device are shown in Figure 3a,b. The investigated PD has an equivalent radius of 12 µm and was implemented on a 220 nm SOI wafer with a remaining slab height of 90 nm. A simple PN junction was formed within the 500 nm wide waveguide.
is the index thermo-optic coefficient of silicon, which is ~1.94 × 10 K at 1310 nm.
To analyze the small-signal response, including the nonlinear thermal optical effect of a MRR PD, we first need to obtain the static optical power in the MRR by calculating the resonator optical enhancement. The heating power and generated temperature variation can then be used to calculate the resonator frequency shift as well as the detuning frequency. Finally, Equations (5) and (6) will tell us the frequency response.

Experimental Verification
To verify the small-signal model in the previous section, we designed an all-Si MRR PD at Advanced Micro Foundry (AMF). The cross section of the MRR PD and a microscope image of the fabricated device are shown in Figure 3a,b. The investigated PD has an equivalent radius of 12 μm and was implemented on a 220 nm SOI wafer with a remaining slab height of 90 nm. A simple PN junction was formed within the 500 nm wide waveguide.  Figure 4a shows the measured transmission spectrum of the MRR around 1310 nm with low optical input power (−4.5 dBm) when the bias voltage is −6 V. The extinction ratio of the resonator is ~10 dB and Q factor is ~7000. The field coupling coefficient and the round-trip field attenuation factor can be extracted as 0.36 and 0.96, respectively, by fitting the measured spectrum using the equations as shown in [28]. The attenuation mainly comes from the high impurity concentration in the MRR waveguide area and PAT effect.  Figure 4a shows the measured transmission spectrum of the MRR around 1310 nm with low optical input power (−4.5 dBm) when the bias voltage is −6 V. The extinction ratio of the resonator is~10 dB and Q factor is~7000. The field coupling coefficient k and the round-trip field attenuation factor a rt can be extracted as 0.36 and 0.96, respectively, by fitting the measured spectrum using the equations as shown in [28]. The attenuation mainly comes from the high impurity concentration in the MRR waveguide area and PAT effect.  The measured photo and dark current are demonstrated in Figure 4b. The current was measured at its highest point (~1311.6 nm) when the optical input power is 6.5 dBm. The current increases with the reverse bias due to the enhanced PAT effect and the avalanche gain. At −6 V, the MRR PD shows a dark current of 0.5 μA and a photo current of 0.41 mA. The measured photo and dark current are demonstrated in Figure 4b. The current was measured at its highest point (~1311.6 nm) when the optical input power is 6.5 dBm. The current increases with the reverse bias due to the enhanced PAT effect and the avalanche gain. At −6 V, the MRR PD shows a dark current of 0.5 µA and a photo current of 0.41 mA.

Thermal-Induced Nonlinearity
For the MRR PD, the CW transmission spectrum can be significantly modified with thermal-induced nonlinearity. The optical transmission spectrum with high optical input power (6.5 dBm) was measured at −6 V and is shown in Figure 5. We observe the resonance wavelength redshifts until the wavelength scans up to~1311.6 nm. The nonlinear refractive index change can be attributed to the increase of temperature resulting from the ohmic heating and optical absorption. Beyond this wavelength, the heating power disappears, and the resonance shifts back to the cold cavity value. To verify the above analysis, Lumerical HEAT was used to simulate temperature variation when the bus waveguide power is 6.5 dBm and bias voltage is −6 V. The input power consists of the ohmic heating and optical loss. Ohmic heating is proportional to the product , which will accumulate around the junction. Total optical attenuation can be divided into absorption due to PAT that generates photocurrent and loss mainly due to FCA that generates heat. The power loss coefficient, which is about 6.5% per round trip, can be extracted from the transmission spectrum at low bias when the PAT effect is negligible. The detailed simulation parameters are shown in the Table 1. The simulated temperature variation at 1311.45 nm is shown in Figure 6. The temperature at the bottom of the simulation region (in the Si substrate) was set to the room temperature of 300 K. It can be observed that there is a strong localization of temperature variation near the waveguide area. The average temperature rises by 12.1 K and the corresponding resonance shift is 0.78 nm. We observe close agreement with the calculated value based on the transmission spectrum, which corroborates the thermal analysis. The temperature variation and induced resonant shift at other wavelengths were also simulated and are shown in Table 1. By linear fitting, the thermal impedance of the structure In order to calculate the small-signal response at different input wavelengths, the nonlinear resonance shift and detuning frequency are required. The easiest way to obtain these is based on the transmission power given by the spectrum. Because the nonlinear loss from TPA is small, we can assume that the loss coefficient does not change significantly with the optical power [22,35]. The calculated linear spectrums at different input wavelengths are shown as dotted lines in Figure 5. The difference between input and resonance wavelength can be calculated as the transmission power T ap and is known using the following formula [28]: where τ is the field transmission coefficient (τ = √ 1 − k 2 ) and ∅ rt is the round-trip phase, which is related to the detuning frequency. The calculated resonances at different input wavelengths are shown in Table 1.
To verify the above analysis, Lumerical HEAT was used to simulate temperature variation when the bus waveguide power is 6.5 dBm and bias voltage is −6 V. The input power consists of the ohmic heating and optical loss. Ohmic heating is proportional to the product I total V bias , which will accumulate around the junction. Total optical attenuation can be divided into absorption due to PAT that generates photocurrent and loss mainly due to FCA that generates heat. The power loss coefficient, which is about 6.5% per round trip, can be extracted from the transmission spectrum at low bias when the PAT effect is negligible. The detailed simulation parameters are shown in the Table 1. The simulated temperature variation at 1311.45 nm is shown in Figure 6. The temperature at the bottom of the simulation region (in the Si substrate) was set to the room temperature of 300 K. It can be observed that there is a strong localization of temperature variation near the waveguide area. The average temperature rises by 12.1 K and the corresponding resonance shift is 0.78 nm. We observe close agreement with the calculated value based on the transmission spectrum, which corroborates the thermal analysis. The temperature variation and induced resonant shift at other wavelengths were also simulated and are shown in Table 1. By linear fitting, the thermal impedance of the structure is 3.45 K/mW. The thermal impedance was also measured when the input optical power is −4.5 dBm and the thermal nonlinear effect can be neglected. As shown in Figure 7, the measured resonance shift versus heat power was 215 pm/mW and the corresponding thermal impedance was 3.37 K/mW, which is close to the simulated value. The small discrepancy may come from the small piece of contact metal on the rim that was not included in the simulation. The relatively high thermal impedance can be attributed to the poor thermal conductivity of SiO2 and thick BOX. The thermal impedance was also measured when the input optical power is −4.5 dBm and the thermal nonlinear effect can be neglected. As shown in Figure 7, the measured resonance shift versus heat power was 215 pm/mW and the corresponding thermal impedance was 3.37 K/mW, which is close to the simulated value. The small discrepancy may come from the small piece of contact metal on the rim that was not included in the simulation. The relatively high thermal impedance can be attributed to the poor thermal conductivity of SiO 2 and thick BOX.
The nonlinear coefficient n TO , due to the thermo-optic effect, can be defined as , where ∆n is the nonlinear index change and P r is the optical power in the MRR. A e f f is the effective mode area, which was simulated in Lumerical MODE as 0.15 µm 2 . The calculated nonlinear coefficient is 0.15 cm 2 /GW, which is comparable to the results in [20]. The value is more than 3 orders of magnitude higher compared to the Kerr effect in Si, which can be attributed to the photocurrent generated in the MRR PD [36].

Frequency Response
In the small-signal equivalent circuit model as shown in Figure 1, the series capacitance C j , oxide layer capacitance C ox , and parasitic capacitance C P are 5.7 fF, 5.8 fF, and 36.8 fF, respectively. The series resistance R s , substrate resistance R sub , and parasitic resistance R P are 91 Ω, 1967 Ω, and 22 Ω, respectively. Parasitic inductance L P is 40 pH. All the components are extracted from S parameter fitting. The load resistance R L is 50 Ω. Based on these values, the RLC-limited 3 dB bandwidth of the MRR PD is~66 GHz.
(a) (b) Figure 6. Simulated temperature variation when the input wavelength is 1311.45 nm in (a) X-Y plane and (b) Y-Z plane.
The thermal impedance was also measured when the input optical power is −4.5 dBm and the thermal nonlinear effect can be neglected. As shown in Figure 7, the measured resonance shift versus heat power was 215 pm/mW and the corresponding thermal impedance was 3.37 K/mW, which is close to the simulated value. The small discrepancy may come from the small piece of contact metal on the rim that was not included in the simulation. The relatively high thermal impedance can be attributed to the poor thermal conductivity of SiO2 and thick BOX. is the effective mode area, which was simulated in Lumerical MODE as 0.15 μm 2 . The calculated nonlinear coefficient is 0.15 cm 2 /GW, which is comparable to the results in [20]. The value is more than 3 orders of magnitude higher compared to the Kerr effect in Si, which can be attributed to the photocurrent generated in the MRR PD [36]. Substituting the laser detuning obtained in the last subsection and other MRR parameters in Equation (5), the frequency responses of the MRR PD are calculated and shown as the solid lines in Figure 8. The S21 were measured for input wavelength from 1310.65 nm to 1311.45 nm at −6 V with a commercial 65 GHz MZM and a vector network analyzer (VNA). The measurement results are also shown in Figure 8 with dotted lines, and they show excellent agreement with our calculated results. The data is normalized to the DC value at 1311.45 nm. The small ripples in the measured response are due to the multiple reflections of the signal at the GSG probe and the high-speed electrical lines. It can be observed that peaking occurs when the detuning frequency is large. This can be explained by the fact that the absolute value of the zero in the transfer function becomes larger than the pole pair. Although the bandwidth can be improved with detuning frequency, the output signal at high frequency is weaker and the responsivity-bandwidth product decreases as shown in Figure 2. It can be concluded that the best operating wavelength for MRR PD is the resonance wavelength to obtain the maximum output power at each frequency. When the detuning frequency is zero, the 3 dB bandwidth of the MRR PD can be simplified as The MRR PD eye diagrams were measured and are shown in Appendix A. Our PDs can support 80 Gb/s nonreturn-to-zero (NRZ) signal and 100 Gb/s PAM4 signal transmission. The summary of the all-Si PDs and APDs is shown in Appendix B. quency, the output signal at high frequency is weaker and the responsivity-bandwidth product decreases as shown in Figure 2. It can be concluded that the best operating wavelength for MRR PD is the resonance wavelength to obtain the maximum output power at each frequency. When the detuning frequency is zero, the 3 dB bandwidth of the MRR PD can be simplified as ( ) ( ) .

Discussion and Conclusions
For the first time, we modeled the small-signal response of MRR PDs that can predict the frequency response as well as the responsivity. We derived the thermal-and frequency-dependent transfer function based on the rate equation and perturbation theory.

Discussion and Conclusions
For the first time, we modeled the small-signal response of MRR PDs that can predict the frequency response as well as the responsivity. We derived the thermal-and frequencydependent transfer function based on the rate equation and perturbation theory. Both electrical and optical behavior of the PD are considered in the expression, which suggests that a MRR PD is a fourth-order system with one zero, a pair of complex-conjugate poles and two poles from the RLC circuit. The thermo-optic nonlinear effect was investigated, and the simulated thermal impedance has been verified with measurement. The model has been verified experimentally to be accurate and allows for the determination of an O/E response of a MRR PD with high precision. The responsivity-bandwidth product is adopted as a figure of merit to demonstrate the performance of PD is superior when the detuning frequency is zero.
The derived closed-form transfer function allows us to calculate the small-signal behaviors and optimize the MRR PD design for different aims. The resonant wavelength is verified to be the optimal operating wavelength. The responsivity of PD can be further improved if the PAT effect can be enlarged.
The wavelength-selective photodetection makes the MRR PD promising in the WDM Rx application. There is no need for an additional optical filter and a small resonator size can enable a small footprint for the Rx. All-Si photodetection without Ge growth can also reduce material costs and simplify the process complexity. In addition, if the MRR PD can work at a higher bias voltage, such as an APD, the Rx sensitivity can be improved. Further work is required to check its performance as an APD, and a more complex model can be established to analyze the response for an MRR APD.

Appendix A
The NRZ and PAM4 eye diagrams of the MRR PD were measured with the pseudorandom bit sequence 9 (PRBS9) signals, which were generated with a 92 GSa/s arbitrary waveform generator (AWG). The optical signal from an O-band tunable laser was modulated by a 65 GHz MZM. The signal was then amplified with a praseodymium-doped fiber amplifier (PDFA). The RF power from the all-Si MRR PD was detected after a bias tee by a 60 GHz digital communication analyzer oscilloscope (DCA). The eye diagrams were measured at resonant wavelength at −6.4 V without an equalizer and have been demonstrated in [35]. As shown in Figure A1, the MRR PD can support the NRZ signal with a data rate up to 80 Gb/s and a PAM4 signal with a data rate up to 100 Gb/s. The eye opening at PAM4 modulation is mainly limited by the responsivity of PD.

Appendix B
To better understand the performance of our MRR PDs, the comparison of different all-Si PDs/APDs is demonstrated in Table A1. Compared to other waveguide PDs, our MRR PDs achieved comparable bandwidth and high bit rate transmission with a compact size.