Abstract

We theoretically study the problem of Raman gain maximization in uniform silicon photonic-crystal waveguides supporting slow optical modes. For the first time, an exact solution to this problem is obtained within the framework of the undepleted-pump approximation. Specifically, we derive analytical expressions for the maximum signal gain, optimal input pump power, and optimal length of a silicon Raman amplifier and demonstrate that the ultimate gain is achieved when the pump beam propagates at its maximum speed. If the signal’s group velocity can be reduced by a factor of 10 compared to its value in a bulk silicon, it may result in ultrahigh gains exceeding 100 dB. We also optimize the device parameters of a silicon Raman amplifier in the regime of strong pump depletion and come up with general design guidelines that can be used in practice.

1. Introduction

The rapid progress in the field of silicon photonics realized during recent years draws the reality of all-optical signal processing on a nanoscale closer than ever [13]. Of paramount importance for the integration of optical functionalities on a silicon platform is the problem of efficient signal amplification through stimulated Raman scattering (SRS) [49]. Even though the Raman gain coefficient is relatively large for the silicon material, its use for signal amplification is hampered considerably by its large two-photon absorption (TPA) coefficient (~0.5 GW/cm2) in the wavelength region near 1550 nm [10] and the effects of free-carrier absorption (FCA) that come into play once free carriers are generated via TPA.

Several techniques have been successfully employed to mitigate the nonlinear losses and increase the Raman gain in silicon-on-insulator (SOI) waveguides. The most straightforward method lies in the reduction of free-carrier lifetime. It is typically implemented through a reverse-bias p-i-n diode, which uses silicon waveguide as an intrinsic region and removes carriers from the region near its mode center through the applied static electric field [1113]. This approach has been used to make silicon Raman amplifiers (SRAs) operating continuously. In the pulsed regime, in addition to gain enhancement, this technique allows operation of an SRA at a bit rate as high as 100 Gbit/s. Once free carriers are efficiently removed and they stop affecting optical propagation, TPA becomes the only nonlinear effect that limits performance of an SRA [14].

Signal gain can also be improved by successively optimizing both the length of the SOI waveguide and the input pump intensity [15, 16]. Existence of the unique length that provides the maximum gain for given material parameters is evident from the following consideration. If the waveguide is too short, signal does not have enough time to build up its energy through interaction with the pump. In a very long waveguide, on the other hand, the gain approaches zero, as the signal is constantly exposed to TPA and linear losses that decay with optical intensity slower than the SRS does. This means that the silicon waveguide should end at the point where the gain peaks. Inside an SRA of optimal length, the signal gain can still have a local minimum, in the event that FCA-induced energy dissipation predominates signal amplification near the input facet of the waveguide. If that is the case, one can increase the gain by reducing the intensity of the input pump beam. The existence of an optimal pump intensity for a fixed length of SRA can also be understood with a similar argument.

Another way of enhancing the amplified signal lies in engineering (tapering) the lateral cross-section of a SRA [1722]. The objective of tapering is to provide an optimal balance between the nonlinear losses and SRS, which is possible because they scale differently with the effective mode area. By developing a semianalytical method for calculating the optimal tapering profile of an SRA, we have shown that linear tapering of the amplifier is capable of increasing gain almost up to its ultimate limit that is reached with an irregular optimal tapering [17].

A somewhat more challenging, but quite efficient, gain-improvement technique relies on the slow-light enhancement of nonlinear optical phenomena in silicon photonic crystal (SPhC) waveguides [2329]. A reduction in the group velocities of the pump and signal beams inside a SPhC waveguide leads to an increase in beams intensities and allows them to stay longer inside the amplifier, thereby improving the efficiency of the SRS process. This speed reduction, unfortunately, also results in the enhancement of TPA and FCA, making the problem of gain maximization in slow-light SRAs more intricate compared to that in ordinary SRAs. The only attempt to analyze gain scaling with group velocity and find the optimal parameters for a SPhC-based SRA belongs to Krause et al. [30]. In this paper, we extend this study and present an exact solution to the problem of gain maximization in slow-light SRAs within the undepleted-pump approximation. We also perform a numerical optimization of various SRA parameters in the regime of strong pump depletion, in order to come up with the design guidelines that can be applied in practice.

2. Formulation of the Problem

Consider amplification of a continuous wave (CW) signal at the frequency 𝜔𝑠 by a CW pump at the frequency 𝜔𝑝 inside a SPhC waveguide, which is characterized by the slow-down factors 𝑆𝑠 and 𝑆𝑝 at frequencies 𝜔𝑠 and 𝜔𝑝, respectively. If the pump copropagates with the signal, then their “slowed-down” intensities 𝑄𝑝(𝑠)(𝑧)=𝑆𝑝(𝑠)𝐼𝑝(𝑠)(𝑧), where 𝐼𝑝(𝑠)(𝑧) are the ordinary intensities, satisfy the following system of coupled differential equations [23, 30]:1𝑆𝑝𝑄𝑝d𝑄𝑝d𝑧=𝛼𝑝𝛽𝑝𝑄𝑝𝛾𝑝𝑄𝑠𝜉𝑝𝑄2𝑝+4𝑄𝑝𝑄𝑠+𝜔𝑝/𝜔𝑠𝑄2𝑠,1𝑆𝑠𝑄𝑠d𝑄𝑠d𝑧=𝛼𝑠𝛽𝑠𝑄𝑠+𝛾𝑠𝑄𝑝𝜉𝑠𝑄2𝑝+4𝑄𝑝𝑄𝑠+𝜔𝑝/𝜔𝑠𝑄2𝑠.(1) In these equations, the terms characterized by the coefficients 𝛼𝑝(𝑠), 𝛽𝑝(𝑠), and 𝜉𝑝(𝑠) account for linear losses, TPA, and FCA of the pump and signal beams, respectively. The FCA coefficients are defined as 𝜉𝑝(𝑠)=𝜎0𝜔0/𝜔𝑝(𝑠)2𝜏𝑐𝛽𝑝2𝜔𝑝,(2) where 𝜎0=1.45×1021m2, 𝜔0=2𝜋𝑐/(1.55𝜇m), 𝜏𝑐 is the effective free-carrier lifetime, and is the reduced Planck's constant. The nonlinear parameters, 𝛾𝑝 and 𝛾𝑠, have been defined such that they describe both the cross-TPA and SRS effects. Assuming that the pump frequency exceeds the signal frequency exactly by the Raman shift 𝜔𝑝𝜔𝑠=ΩR(2𝜋)15.6 THz, we have𝛾𝑝=gR+2𝛽𝑝,𝛾𝑠=gR2𝛽𝑠𝜂,(3) where gR is the Raman gain coefficient and 𝜂=𝜔𝑝/𝜔𝑠.

Equations (1) are similar in form to the coupled-intensity equations obtained without the slow-light effects [3133], but they contain two additional parameters 𝑆𝑝 and 𝑆𝑠. Strictly speaking, these parameters are related to each other, since they both depend on geometric and material characteristics of the same photonic crystal. However, owing to a great number of such characteristics and the possibility of varying most of them in different ways within wide limits, we may think of the slow-down factors for pump and signal modes as independent of each other. These factors can be adjusted to take different value (from 1 to as large as 1000) by engineering the dispersion of the SPhC waveguide [29, 34]. We also assume that the input pump intensity 𝐼𝑝0 may take any value between zero and infinity because it can be varied by changing either the pump power or the SRA cross-section.

Hence, the problem of gain maximization in a slow-light SRA, for a given intensity of the input signal, consists of searching for the optimal values of the waveguide length 𝐿, input pump intensity 𝐼𝑝0, and slow-down factors 𝑆𝑠 and 𝑆𝑝. If some of these parameters are fixed by practical needs or restricted by fabrication capability, the optimization problem becomes simpler because it involves fewer than four parameters.

3. Exact Solution: Undepleted-Pump Approximation

The signal gain can be readily optimized analytically in the undepleted-pump regime, which usually holds for relatively short SRAs with gR10 GW/cm2 and for low input signal powers. We now consider this scenario in detail, for it is quite instructive and helps us to illuminate the peculiarities of the optimization problem in the general case.

In the approximation of an undepleted pump, valid as long as 𝑄𝑝𝑄𝑠 everywhere inside the SPhC waveguide, (1) reduces to1𝑆𝑝𝑄𝑝d𝑄𝑝d𝑧𝛼𝑝𝛽𝑝𝑄𝑝𝜉𝑝𝑄2𝑝,1(4a)𝑆𝑠𝑄𝑠d𝑄𝑠d𝑧𝛼𝑠+𝛾𝑠𝑄𝑝𝜉𝑠𝑄2𝑝.(4b)

Equation (4b) shows that the signal gain grows with increasing 𝑆𝑠, confirming that one can improve the performance of an SRA by slowing down the signal mode as much as possible [30]. The situation is less clear with the optimal slow-down factor for the pump. Indeed, when the input pump intensity is fixed, an increase in the time the pump spends inside the amplifier may result in either increase or decrease in the value of 𝑄𝑝, which, in turn, may enhance or reduce the local signal gain. Therefore, the three parameters to be optimized in the undepleted-pump regime are the amplifier length, the input pump intensity, and the slow-down factor for the pump field.

To find the values of the three optimal parameters, we use the implicit solution of (4a) and (4b). It can be written compactly with the shortened notations 𝑄0𝑄𝑝(0) and 𝑄𝐿𝑄𝑝(𝐿) in the form of the following coupled equations [35]: 𝛼𝑝+𝛽𝑝𝑄𝐿+𝜉𝑝𝑄2𝐿𝛼𝑝+𝛽𝑝𝑄0+𝜉𝑝𝑄20=𝑄2𝐿𝑄20exp2𝛼𝑝𝑆𝑝𝐿+𝛽𝑝𝑄0𝐿e𝐿,(5a)e=𝜇𝛽𝑝𝑄0ln𝜇𝐹𝐿+1𝜇𝐹𝐿1𝜇𝐹01𝜇𝐹0,+1(5b)where 𝜇=(14𝛼𝑝𝜉𝑝/𝛽2𝑝)1/2, 𝐹𝑗=1+2𝜉𝑝𝑄𝑗/𝛽𝑝, and 𝐿e is the effective length of the amplifier.

Using (4b), (5a), and (5b), it is easy to show that the signal gain 𝐺𝑠 is given by the relation𝑆𝑝𝑆𝑠ln𝐺𝑠𝛼=𝑠𝜂2𝛼𝑝𝑆𝑝𝐿+𝛾𝑠+𝜂2𝛽𝑝𝑄0𝐿e𝜂2𝑄ln0𝑄𝐿.(6)

We maximize 𝐺𝑠 with respect to the three optimization parameters, 𝐿, 𝐼𝑝0, and 𝑆𝑝, by requiring that the following three derivatives vanish: 𝜕𝐺𝑠𝐿,𝐼𝑝0,𝑆𝑝𝜕𝐿=0,(7a)𝜕𝐺𝑠𝐿,𝐼𝑝0,𝑆𝑝𝜕𝐼𝑝0=0,(7b)𝜕𝐺𝑠𝐿,𝐼𝑝0,𝑆𝑝𝜕𝑆𝑝=0.(7c)Upon employing (4b)–(6), these requirements lead to the following three relations: 𝛼𝑠𝛾𝑠𝑄𝐿+𝜉𝑠𝑄2𝐿𝑄=0,(8a)0+𝑄𝐿=𝛾𝑠𝜉𝑠,𝑆(8b)𝑝𝑆𝑠ln𝐺𝑠=𝛾𝑠𝜉𝑠𝑄0+𝑄𝐿𝛼𝑝+𝛽𝑝𝑄0+𝜉𝑝𝑄20𝑄0𝑄𝐿𝛼𝑠𝛾𝑠𝑄𝐿+𝜉𝑠𝑄2𝐿𝑆𝑝𝐿.(8c)

The combination of (8a)–(8c) leads to the equality ln𝐺𝑠=0, which is in conflict with (6). This implies that not more than two of the parameters 𝐿, 𝐼𝑝0, and 𝑆𝑝 can be optimized at the same time. Next, substitution of (8b) into (8c) results in ln𝐺𝑠𝛼=𝑠𝛾𝑠𝑄𝐿+𝜉𝑠𝑄2𝐿𝑆𝑠𝐿.(9) This equation is not satisfied in real SRAs since, according to (4b), it requires pump power to be constant along the amplifier. Consequently, the gain provided by an SRA of a given length peaks simultaneously with either the slow-down factor 𝑆𝑝 or the input pump intensity. The optimal parameters of an SRA depend on the optimization scenario. In what follows, we consider the three possible optimization schemes.

3.1. Optimal 𝐿 and 𝐼𝑝0 for a Fixed 𝑆𝑝

When the slow-down factor of the pump is fixed, (8a) and (8b) yield the critical values for the quantities 𝑄0 and 𝑄𝐿 of the form𝑄0=𝛾𝑠+𝛾2𝑠4𝛼𝑠𝜉𝑠2𝜉𝑠,𝑄𝐿=𝛾𝑠𝛾2𝑠4𝛼𝑠𝜉𝑠2𝜉𝑠.(10) It is evident from these expressions that the inequality 𝛾2𝑠>4𝛼𝑠𝜉𝑠 must be satisfied to allow for amplification of signal. This inequality sets the upper limit for the product 𝛼𝑠𝜏𝑐 [30]. For 𝛾2𝑠>4𝛼𝑠𝜉𝑠, one can see from (4b) that the derivative d𝑄𝑠/d𝑧 remains positive as long as 𝑄0>𝑄𝑝>𝑄𝐿. The expressions in (10) simply determine the pump intensities corresponding to the onset and the end of local net gain and thus enable signal amplification along the whole length of a SPhC waveguide.

With the initial and final intensities given by (10), the optimal length of the SRA can be explicitly found from (5a) while the optimal input intensity of the pump is simply given by 𝐼𝑝0=𝑄0/𝑆𝑝.

It is important to note that the optimal values of 𝑄0 and 𝑄𝐿 do not depend on the parameter 𝑆𝑝. Assuming that the signal exhibits a net gain, we conclude from (6) that the maximum signal gain corresponds to the lowest possible value of 𝑆𝑝, that is, to 𝑆𝑝=1. According to (4a), this value results in the slowest attenuation of the pump along the amplifier.

Hence, the optimal values of three parameters that maximize the signal gain for this optimization scheme are given by 𝐼𝑝0=𝑄0,𝑆𝑝=1,(11a)1𝐿=2𝛼𝑝𝑄ln20𝑄2𝐿𝛼𝑝+𝛽𝑝𝑄𝐿+𝜉𝑝𝑄2𝐿𝛼𝑝+𝛽𝑝𝑄0+𝜉𝑝𝑄20𝛽𝑝𝑄0𝐿e,(11b)where 𝑄0 and 𝑄𝐿 are given by (10).

3.2. Optimal 𝐿 and 𝑆𝑝 for a Fixed 𝐼𝑝0

If we keep the input pump power constant, then (8a) yields the optimal value for 𝑄𝐿 given in (10), where the minus before the radical is chosen to provide the maximal signal gain. With this value, the critical value of 𝑄0>𝑄𝐿 is found by simultaneously solving (5a), (5b), (6), and (8c). It should be noted that 𝑄0 obtained in this way does not generally coincide with that in (10).

Once the optimal quantities 𝑄0 and 𝑄𝐿 are known, the optimal slow-down factor 𝑆𝑝 is given by 𝑆𝑝=𝑄0/𝐼𝑝0 and the optimal length of the SRA is calculated from (5a).

On substitution of (8a) into (8c), the signal gain may be written asln𝐺𝑠=𝛾𝑠𝜉𝑠𝑄0+𝑄𝐿𝛼𝑝+𝛽𝑝𝑄0+𝜉𝑝𝑄20𝑄0𝑄𝐿𝑆1𝑠𝐼𝑝0.(12) This expression shows that the peak net gain grows in proportion to the input pump power. Since the maximum pump power is equal to 𝑄0 when 𝑆𝑝=1, we again arrive at the set of optimal parameters given in (11a) and (11b).

3.3. Optimal 𝐼𝑝0 and 𝑆𝑝 for a Fixed 𝐿

Consider finally the situation in which the length of an SRA is preset. As we have seen above, in this case (8b) and (8c) are inconsistent. Therefore, one may either find the optimal input pump power, 𝐼𝑝0(𝑆𝑝,𝐿), for given values of 𝑆𝑝 and 𝐿 by solving (5a), (5b), and (8b) or the optimal slow-down factor of the pump, 𝑆𝑝(𝐼𝑝0,𝐿), for given values of 𝐼𝑝0 and 𝐿 by solving (5a), (5b), and (8c). Of course, if after that we also optimize the length of the amplifier, we arrive at the results obtained for the previous two optimization schemes.

Since the values of 𝑄0 obtained with different optimization schemes are different, the optimal amplifier length and input pump power in (11a) and (11b) are also different. From physical considerations it is, however, clear that the optimal values of 𝑄0 and 𝑄𝐿 providing the ultimate signal gain are those given in (10).

In Figure 1, we plot the maximum signal gain as a function of free-carrier lifetime that can be achieved in the undepleted-pump regime with a “fast” signal mode (𝑆𝑠=1). The same figure shows the optimal SRA length and the optimal pump intensity required for the maximum gain. The material parameters are chosen to be 𝜆𝑠=1.55𝜇m, 𝜆𝑝=1.434𝜇m, 𝛼𝑠=𝛼𝑝=0.5–2 dB/cm, 𝛽𝑠=𝛽𝑝=0.5GW/cm2, and gR=20 GW/cm2. It is seen that the signal gains of >20 dB can be realized with pump powers below 1 GW/cm2 if 𝜏𝑐1 ns and linear losses do not exceed 1 dB/cm. These gains can be increased further by a factor of 𝑆𝑠, which can be more than an order of magnitude, provided the signal beam is sufficiently slowed down by dispersion engineering. Unfortunately, such strong amplification requires relatively long SPhC waveguides (longer than 8 cm in the above example).

Equations (6) and (11b) show that the length of a SRA can be reduced to a few micrometers without noticeably affecting the gain, by simultaneously increasing linear losses and reducing free-carrier absorption in such a manner that the product 𝛼𝑠𝜏𝑐 remains constant (so that the values of 𝑄0 and 𝑄𝐿 do not change). For instance, gains in excess of 20 dB (for 𝑆𝑠=1) can be achieved in a 1.5 mm long SRA with 𝜏𝑐=0.02 ns and 𝛼𝑠=50 dB/cm. The optimum length of a SRA also becomes shorter when the amplifier operates in the regime of considerable pump depletion, which we analyze in the next section.

4. Numerical Solution: General Case

As we saw in the previous section, in the undepleted-pump regime, the signal mode should be slowed down as much as possible, while the pump mode should travel with the highest speed. This is not the case when the intensity of input signal, 𝐼𝑠0, approaches the ratio 𝑄𝐿/(𝑆𝑠𝐺𝑠) and the pump’s depletion can no longer be ignored [36]. The examples below illustrate the peculiarities of SRA optimization in the regime of heavy pump depletion, in which signal gain is calculated numerically from (1).

Figure 2 shows the contour plots of the net signal gain in an 8 mm long SRA for 𝐼𝑝0=30 MW/cm2 and 𝐼𝑠0=0.1 MW/cm2; the depletion [37] of the pump beam at the amplifier output equals 76%. It is seen that the maximum gain of approximately 17.7 dB is achieved for the slow-down factors 𝑆𝑠7.8 and 𝑆𝑝3.4. More importantly, this gain value exceeds the maximum gain of 0.85 dB that occurs in the absence of the slow-light effects (𝑆𝑠=𝑆𝑝=1) by more than a factor of 20. It should also be noted that, owing to a relatively flat gain peak, the range of tolerable slow-down parameters is rather broad, a feature that facilitates the design of SRAs.

The gain corresponding to the optimal slow-down factors can be increased by optimizing the input pump power. Red curves in Figures 2 and 3 show how the peak gain and the optimal factors 𝑆𝑠 and 𝑆𝑝 vary with 𝐼𝑝0 for the same SRA of 8 mm length. One can see that the maximum peak gain of 31.7 dB is reached for 𝐼𝑝03.5 GW/cm2, 𝑆𝑠1.6, and 𝑆𝑝=1. Although the obtained values of 𝐼𝑝0 and 𝑆𝑝 coincide with those in the undepleted-pump regime, the maximum gain and the optimal amplifier length differ from the predictions in Figure 1 (114 dB and 11 cm, resp.). This is a consequence of strong pump depletion (about 84% in the present case).

A further enhancement of gain is generally possible by adjusting the length of the SPhC waveguide. Figure 3 shows this possibility by varying the amplifier length in the range of 2 to 12 mm. However, as the gain profiles corresponding to 2, 4, and 12 mm long SRAs suggest, 𝐿8 mm is the optimal length of the amplifier for this specific set of material parameters. For a typical effective mode area of 0.35 μm2, pump power of 10 W is required at the input of the optimized SRA. Since it is very difficult to get such high input power levels in a CW regime, the actual maximum gain will be lower than 31.7 dB and will be set by the available pump power.

Comparison of different curves in Figure 3 leads to several conclusions that are important for designing SRAs with optimum performance. First, the optimum slow-down factors for pump and signal can be reduced by increasing amplifier’s length, because SRS is more efficient in longer waveguides. Second, as long as the ultimate signal gain is not achieved, larger pump intensities can compensate for operating with faster pump and signal modes. This is possible, since the reduction in time that light spends inside the amplifier is canceled by its stronger interaction with silicon. Third, strong FCA occurring at high intensities reduces pump depletion and, as a result, requires the use of a slowed signal mode (see the blue dashed curve). This regime is to be avoided in practice because the same peak gain can be achieved with a smaller pump intensity.

As a concluding remark, it is worth noting that our results are valid not only for SPhC waveguides, but also for any other types of slow-light silicon waveguides, regardless of the relative group delay that occurs upon light propagation through them.

5. Conclusions

We have theoretically analyzed the problem of gain maximization in silicon Raman amplifiers operating with slow-light pump and signal modes. Using the undepleted-pump approximation, we derived expressions for the maximum signal gain, optimal input pump power, and optimal length of the amplifier with a fixed cross-section. We also showed that the signal gain is maximized when the pump beam moves at its maximum speed. The signal gain grows indefinitely and can exceed 100 dB, provided the signal beam travels more than 10 times slower than its speed in a bulk silicon. In the regime of substantial pump depletion, the parameters of an amplifier were optimized numerically and a relatively large CW gain of 31.7 dB was predicted for an 8 mm long waveguide pumped at an intensity of 3.5 GW/cm2. This result suggests that optimization of slow-light silicon photonic crystal waveguides is a promising approach to the realization of highly efficient Raman amplifiers built using the SOI platform.

Acknowledgments

This work was supported by the Australian Research Council, through its Discovery Grant scheme under grants nos. DP0877232 and DP110100713. The work of G. P. Agrawal is also supported by the US National Science Foundation award ECCS-0801772.