Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Maximization of net optical gain in silicon-waveguide Raman amplifiers

Open Access Open Access

Abstract

We present a novel method for maximizing signal gain in continuously pumped silicon-waveguide Raman amplifiers made with silicon-on-insulator technology. Our method allows for pump-power depletion during Raman amplification and makes use of a variational technique. Its use leads to a system of four coupled nonlinear differential equations, whose numerical solution provides the optimal axial profile of the effective mode area along the waveguide length that maximizes the output signal power for a given amplifier length and a preset input (or output) cross-section area. In practice, the optimum profile can be realized by varying the cross-section area of a silicon waveguide along its length by tapering its width appropriately.

©2009 Optical Society of America

1. Introduction

Owing to the perceived benefits of active silicon-photonics devices, there is a relentless interest in improving their performance since the first experimental demonstration of optical amplification via stimulated Raman scattering (SRS) in silicon-on-insulator (SOI) waveguides [1–3]. In addition to linear losses, which is the only inherent impairment in fiber optics [4, 5], silicon is plagued with two-photon absorption (TPA) and free-carrier absorption (FCA), two processes that dissipates energy of optical waves propagating through a silicon waveguide [6, 7]. To reduce the impact of TPA-induced FCA on Raman amplification, several theoretical and experimental studies have been carried out during the last five years [8–17]. Among others, several recently developed analytical and semi-analytical results, which provide a clear physical insight into the free-carrier effects in SOI waveguides deserve specific attention [18–20]. Particularly, an analytical solution describing propagation of single pulse through the SOI waveguides has been obtained in [18]. An approximate expressions for pump and signal intensities evolution along the single-pass silicon amplifier have been derived in [19]. In addition, the problem of continuous-wave (CW) Raman amplification of a signal pulse have been addressed semi-analytically using variational techniques in [20]. The results of these studies can be used to improve the performance of SOI-based Raman amplifiers. In this paper, we present another semi-analytical method that allows one to minimize the detrimental effects of FCA during the Raman amplification process.

The first and the most straightforward way to tackle the FCA problem involves decreasing the photogenerated carrier density inside the modal area by either decreasing the carrier lifetime [8, 15] or sweeping the carriers out of the waveguide by utilizing an external electric field [12, 13, 15, 21]. The other method of increasing the gain of a silicon Raman amplifier consists of adjusting the cross-section area of the SOI waveguide as well as its length appropriately [10]. Clearly, the optimal choice of these parameters is always possible because, for each fixed waveguide length, there is an optimal cross-section since zero cross-section corresponds to zero output signal power and an infinite cross-section implies absence of amplification. Similarly, for each fixed cross-section area, there is an optimal waveguide length since an infinite length corresponds to zero output signal power and a zero length implies absence of amplification. Hence, it should be possible to devise an algorithm that finds the optimal cross-section area and optimal length of the SOI waveguide that maximize the signal gain for given input pump and signal powers.

Recently, waveguide tapering was suggested as an effective way to increase the efficiency of SOI-based Raman lasers [9–11]. The principle behind tapering consists of balancing of FCA and SRS processes inside the waveguide by varying its lateral dimensions. Numerical calculations have shown that width tapering allows one to significantly increase the output signal power in silicon Raman amplifiers, and reduce the lasing threshold of silicon Raman lasers, compared with the untapered situation [10,11]. Under the undepleted pump approximation, the optimal tapering profile for CW silicon Raman amplifiers can be found analytically [9]. In fact, within this approximation, the optimal waveguide cross-section area as well as the pump and signal intensities have simple exponential forms along the amplifier length. However, even though the undepleted pump approximation can be justified for some Raman lasers, it is inapplicable in the case of Raman amplifiers, which ideally should be designed for total pump depletion and maximum signal amplification.

To the best of our knowledge, only a standard generic algorithm has been used up to now to find approximately the optimum axial profile of the effective mode area (EMA), related to the cross-section area of a SOI-based Raman laser [10]. This method cannot generate the exact optimal EMA axial profile along the waveguide length because it does not provide a systematic way to converge toward it. In this paper, we develop a novel numerical scheme, based on the variational principle, which allows one to find the exact EMA axial profile for a Raman amplifier that maximizes the output signal power. We reduce the optimization problem to a boundary-value problem for the pump and signal powers and two auxiliary functions introduced judiciously. Using our method for different lengths and input EMAs of a SOI waveguide, one can find the unique amplifier dimensions (cross-section area and length) that provide the maximum output signal power.

2. Variational approach and resulting equations

We focus on a silicon Raman amplifier that is pumped continuously at a frequency ωp to amplify a CW signal at the frequency (ωs. Since optical phases do not play any important role during CW Raman amplification, we only need to consider the evolution of pump and signal powers (Pp and Ps) along the SOI waveguide. The two equations governing the SRS process in this situation are well known [6,7] and can be written as

dPpdz=αpPpβpPp2AeffζpsPpPsAeffσpτc(ρpPp2+ρsPs2+ρpsPpPs)PpAeff2,
dPsdz=αsPsβsPs2AeffζspPsPpAeffσsτc(ρpPp2+ρsPs2+ρpsPpPs)PsAeff2,

where α p(s) is the linear loss coefficient, β p(s) is the TPA coefficient, A eff is the EMA of the propagating mode, βps is the cross-TPA coefficient, τc is the effective carrier lifetime, and ζps and ζsp are the coefficients responsible for both cross-TPA and SRS

ζps=2βps+4gRγR2ΩRΩps(ΩR2Ωps2)2+4γR2Ωps2,
ζsp=ωsωp(2βps4gRγR2ΩRΩps(ΩR2Ωps2)2+4γR2Ωps2).

In these equations, gR = 76 cm/GW is the Raman gain coefficient, ΩR = 15.6 THz is the Stokes shift, γr = 105 GHz is the Raman-gain bandwidth, and Ωps = ωp - ωs. Finally, the FCA effects are included in Eq. (1) through σ p(s) = σ 0(λ p(s)/λ 0)2 with σ 0 = 1.45 × 10-21 m2 and λ 0 = 1550 nm. Also, ρ p(s) = β p(s)(2h̄ω p(s)) and ρps = 2βps/(h̄ωp).

It follows from Eq. (1b) that, to maximize the net gain of signal at the end of the waveguide, we should maximize the functional

G(Aeff)=0L1PsdPsdzdz=0Lg(z,Aeff)dz,

where L is the amplifier length and the local gain coefficient g(z,A eff) is given by

g(z,Aeff)=αsβsPs(z)Aeff(z)ζpsPp(z)Aeff(z)σsτcρpPp2(z)+ρsPs2(z)+ρpsPp(z)Ps(z)Aeff2(z).

In contrast to the case of undepleted pump approximation, such a maximization does not require the local g(z,A eff) to be maximum at each point z along the amplifier. This behavior results from a sophisticated differential relation between functions Pp(z) and Ps(z) which is clearly absent in the case of undepleted pump approximation.

We can solve the maximization problem using the standard method of the calculus of variations [22]. Let us assume that A eff(z) is the EMA axial profile for which G is stationary. Equating the variation δG to zero with respect to three variables, Pp, Ps, and A eff, we obtain the relation

0L(𝒜δPp+𝓑δPs)dz=0L𝒞δAeffdz,

where δPp, δPs, and δA eff represent small axial variations in the three quantities and we have introduced

𝒜=1Aeff(ζsp+σsτc2ρpPp+ρpsPsAeff),
𝓑=1Aeff(βs+σsτc2ρsPs+ρpsPpAeff),
𝒞=βsPs+ζspPpAeff2+2σsτcρpPp2+ρsPs2+ρpsPpPsAeff3.

To find the connection between δPp, δPs, and δA eff, we vary Pp and Ps in Eqs. (1) by small amount around fixed values and obtain

d(δPpdz)=apδPpbsδPscpδAeff,
d(δPsdz)=asδPsbpδPpcsδAeff,

where

ap=αp+2βpPpAeff+ζpsPsAeff+σpτc3ρpPp2+ρsPs2+2ρpsPpPsAeff2,
bp=𝒜Ps,cp=𝒞Pp.

The coefficients as, bs, and cs are obtained from ap, bp, and cp by interchanging subscripts ps and assuming that ρsp = ρps.

Since the variations δPp and δPs are nonlinearly related to the variation δA eff, we cannot solve Eq. (3) directly. However, it is still possible to express the optimal EMA axial profile, A eff(z), in terms of the pump and signal powers. To do this, we introduce two auxiliary functions, φ(z) and ψ(z), whose functional form is yet to be determined. However, apart from the existence of the derivative of these functions, Eqs. (1) do not impose any additional constraints on them. In particular, we can choose their end-point values (i.e., boundary values at the end of the amplifier) freely. Later, we make a specific choice in such a way that subsequent algebraic calculations are simplified.

Multiplying Eq. (4a) by φ and Eq. (4b) by ψ and integrating the resulting equations with respect to z from 0 to L, we obtain the following two equations:

0Lφd(δPp)dzdz=φδPp|0L0LφδPpdz=0Lφ(apδPp+bsδPscpδAeff)dz,
0Lψd(δPs)dzdz=ψδPs|0L0LψδPsdz=0Lψ(asδPs+bpδPpcsδAeff)dz,

where φ′ = /dz and ψ′ = dψ/dz.

We assume that the boundary conditions at the left end of the amplifier for the pump and the signal are fixed so that their variations are equal to zero, i.e., δPp(0) = δPs(0) = 0. Since we are looking for a variational solution that maximizes the signal power at the waveguide output (i.e., at the right end of the amplifier), δPs(L) = 0 as well. Because of these constraints, we do not need to set the boundary values for the auxiliary function ψ. However, we can simplify the algebra considerably by forcing φ(z) to be equal to zero at the right end of the waveguide, i.e., by choosing

φ(L)=0.

With this choice, the preceding two equations reduce to

0L[(φapφ)δPp+φbsδPs]dz=0LφcpδAeffdz,
0L[ψbpδPp+(ψasψ)δPs]dz=0LψcpδAeffdz.

Adding these equations termwise, we finally obtain the relation

0L[(φap+ψbpφ)δPp+(ψas+φbsψ)δPs]dz=0L(φcp+ψcs)δAeffdz.

A comparison of Eqs. (3) and (6) shows that the optimal EMA of the SOI waveguide can be found from a simple relation

φcp+ψcs=𝒞,

if the two auxiliary functions satisfy the differential equations

φ=apφ+bpψ𝒜,
ψ=asψ+bsψ𝓑.
 figure: Fig. 1.

Fig. 1. Left panel: Evolution of the pump and signal powers along SOI waveguides of three different length (L = 8, 10, 12 mm) with optimized EMA profiles. Right panel: Optimized EMA axial profiles for five waveguide lengths ranging from 8 to 12 mm. In all cases, P p0 = 1 W, P s0 = 0.01 W, and A 0 = 1μm2. See the text for other device parameters.

Download Full Size | PDF

Equation (7) can be solved easily to provide the following relation:

Aeff(z)=2τc[σs(1ψPs)σpφPp](ρpPp2+ρsPs2+ρpsPpPs)φPp(βpPp+ζpsPs)(1ψPs)(βsPs+ζspPp).

Irrespective of the input powers and waveguide parameters, this solution corresponds to the maximum of the functional G and thus maximizes net gain of the signal during Raman amplification. This can be proved by the following arguments. On the one hand, the output signal power, Ps(L), equals to zero for all functions A eff(z) such that A eff(z 0) = 0 for z 0 < L; by increasing the waveguide waist at z = z 0 from zero to some finite value, one can increase the output signal power. On the other hand, if for z > 0 the waveguide cross-section is so large that the nonlinear effects are negligible, Ps(L) can be increased by narrowing the waveguide and enhancing the Raman amplification. Since the above analysis has revealed only one optimum A eff(z), we conclude that it should give a maximum value of G (otherwise there should be at least two other optimal functional forms of A eff(z) that correspond to the maxima of G).

Hence, our variational analysis shows that, to find the optimal EMA axial profile of the SOI waveguide (which is related to the optimum tapering of the waveguide width or cross-section area), we should solve the system of four coupled nonlinear differential equations, given in Eqs. (1) and (8), assuming that A eff is given by (9). The 4 boundary conditions needed for solving these 4 equations come from the condition (5) and the three initial conditions: Pp(0) = P p0, Ps(0) = Ps0, and A eff(z = 0) = A 0. The last condition can also be used at z = L if the objective is to design a Raman amplifier with a preset EMA at the output end of the amplifier.

3. Numerical examples

We illustrate the versatility of our method by considering a silicon Raman amplifier with the following simulation parameters: λp = 1550 nm, λs = 1686 nm, αp = αs = 1 dB/cm, βp = βs = βps = 0.5 cm/GW, and τc = 1 ns. Evolution of the pump and signal powers inside optimized amplifiers of three different lengths is shown in the left panel of Fig. 1. The right panel in the same figure shows the optimum EMA axial profiles for five amplifiers designed with A eff = 1 μm2 at z = 0. The incident pump and signals powers were assumed to be P p0 = 1 W and Ps0 =0.01 W. As mentioned earlier, the signal power is maximized at the output end of the amplifier rather than at each point along the waveguide. The longer the amplifier, the larger the output modal area and the smaller the tilt of the EMA axial profile needed for optimization. This is because for longer amplifiers, Raman energy conversion should take place further from the input end to minimize losses. In fact, to provide such a delay in energy conversion, in some cases the modal area should even increase locally at a certain position of the waveguide. This feature is apparent in the right panel of Fig. 1 for the longest Raman amplifier for which the initially decreasing modal area increases slightly between 4 and 6 mm. Note also that the output signal power exceeds 0.4 W in all four cases, corresponding to a net gain by more than a factor of 400 in spite of the TPA and FCA. This is possible because our theory allows for significant pump depletion.

 figure: Fig. 2.

Fig. 2. Left panel: Solution of Eqs. (1) and (8) for a 3-mm-long waveguide with input values P p0 = 1 W, P s0 = 0.5 W, and A 0 = 1 μm2. Solid curves show the pump and signal powers; dashed and dotted curves represent the two auxiliary functions. Right panel: Optimized EMA profile (solid curve, right scale) and the pump and signal intensities (dashed curves, left scale) corresponding to the pump and signal powers shown on the left panel.

Download Full Size | PDF

The extent of pump depletion is apparent in the left panel of Fig. 2 where we show the optimized solution in the case of a 3 mm-long amplifier by plotting variations along the amplifier length of the pump and signal powers together with those of the two auxiliary functions. Now the incident powers, P p0 = 1 W and P s0 = 0.5 W, are so high that the EMA of the silicon waveguide should even increase by about 30% between 0.6 and 1.4 mm (see the right panel in Fig. 2). In this case, the pump and signal intensities (dashed curves in the right panel of Fig. 2) decrease for a while, because it reduces the pump depletion (flattens the pump power profile) and thus conserves pump’s power for subsequent energy conversion closer to the waveguide output. The main point to note is that, even in this highly saturated regime, pump is able to transfer about 85% of its power to the signal, in spite of losses resulting from TPA and FCA.

It is interesting to compare Raman amplification provided by optimized silicon waveguides with that occurring in the case of untapered as well as linearly tapered waveguides. Such a comparison is shown in Fig. 3. Solid curves represent the dependance of output signal power on the output EMA for four linearly tapered waveguides. The dotted line corresponds to the waveguide with a constant EMA of 1 μm2. The crosses represent the output signal powers realized with optimum tapering in the four cases displayed in Figs. 1 and 2. As one may expect, silicon waveguides with a constant EMA all along their length perform the worst in some cases. For example, the output signal power is only 35% of that realized with the optimum EMA profile for an untapered 8-mm-long amplifier. A smaller enhancement of the output signal power (around 10%) occurs for a 10-mm-long waveguide, and nearly no enhancement is found for 12-mm- and 3-mm waveguides. What is more important from a practical perspective is our observation that, for each optimized amplifier there exist a linearly tapered waveguide which gives nearly the same signal amplification (its parameters are determined by the maxima of the solid curves). Thus, for example, instead of reproducing the modal area profile shown in the right panel of Fig. 2, one can reach the same signal gain by fabricating a linearly-tapered amplifier with the output EMA approximately equal to 1.3 μm2.

 figure: Fig. 3.

Fig. 3. Output signal power as a function of output EMA for linearly tapered waveguides of four different lengths. The incident powers correspond to those used in Figs. 1 and 2, and A 0 = 1 μm2. The central dotted line represents a waveguide with constant EMA. The crosses show signal power and modal area at the end of the optimized waveguides in Figs. 1 and 2.

Download Full Size | PDF

An important note concerning the applicability of the proposed method should be made. Clearly, the boundary-value problem that optimizes the amplifier “knows” nothing about the physical meaning of functions involved. As a result, physically meaningless solutions can sometimes arise in the case of very long amplifiers. These solutions lead to divergences in the EMA axial profile at some critical points zi, where the denominator in Eq. (9) vanishes, and to negative values of A eff between two successive critical points. Further work needs to be done to improve the stability of the numerical algorithm.

4. Conclusion

In conclusion, this paper has focused on the optimization of silicon Raman amplifiers in the saturation regime in which a substantial amount of pump power is transferred to the signal. We presented a novel numerical method for calculating the optimum EMA axial profile that provides the largest signal gain for a given value of the EMA and pump and signal powers at the input end of the amplifier. The implementation of the method requires numerical solution of a system of nonlinear differential equations, which contains twice as many equations as the original system describing the Raman amplification in a constant cross-section waveguide. We have solved these equations in several cases of practical interest and have found the optimum EMA profile that can be realized in practice by tapering the width of an SOI waveguide appropriately. We have also shown that in some cases the same maximal signal gain can be achieved by a simple linear tapering of the SOI waveguide. Our method is also applicable to the optimization of silicon Raman lasers with its suitable extension.

Acknowledgment

This work was funded by the Australian Research Council through its Discovery Grant scheme under grant DP0877232. The work of GPA is also supported by the NSF award ECCS-0801772.

References and links

1. R. Claps, D. Dimitropoulos, and B. Jalali, “Stimulated Raman scattering in silicon waveguides,” Electron. Lett. 38, 1352 (2002). [CrossRef]  

2. R. Claps, D. Dimitropoulos, V. Raghunathan, Y. Han, and B. Jalali, “Observation of stimulated Raman amplification in silicon waveguides,” Opt. Express 11, 1731 (2003). [CrossRef]   [PubMed]  

3. B. Jalali, V. Raghunathan, D. Dimitropoulos, and O. Boyraz, “Raman-based silicon photonics,” IEEE J. Sel. Top. Quantum Electron. 12, 412 (2006). [CrossRef]  

4. C. Headley and G. P. Agrawal, Eds., Raman Amplification in Fiber-Optic Communication Systems, (Academic Press, San Diego, 2005).

5. G. P. Agrawal, Nonlinear Fiber Optics, (Academic, Boston, 2007).

6. Q. Lin, O. J. Painter, and G. P. Agrawal, “Nonlinear optical phenomena in silicon waveguides: Modeling and applications,” Opt. Express 15, 16604 (2007). [CrossRef]   [PubMed]  

7. X. Chen, N. C. Panoiu, and R. M. Osgood, “Theory of Raman-mediated pulsed amplification in silicon-wire waveguides,” IEEE J. Quantum Electron. 42, 160 (2006). [CrossRef]  

8. M. Krause, H. Renner, and E. Brinkmeyer, “Analysis of Raman lasing characteristics in silicon-on-insulator waveguides,” Opt. Express 12, 5703–5710 (2004). [CrossRef]   [PubMed]  

9. H. Renner, M. Krause, and E. Brinkmeyer, “Maximal gain and optimal taper design for Raman amplifiers in silicon-on-insulator waveguides,” in Integrated Photonics Research and Applications Topical Meetings (IPRA 2005), paper JWA3.

10. M. Krause, H. Renner, and E. Brinkmeyer, “Efficiency increase of silicon-on-insulator Raman lasers by reduction of free-carrier absorption in tapered waveguides,” in Conference on Lasers and Electro-Optics (CLEO 2005), paper CThB1.

11. M. Krause, H. Renner, and E. Brinkmeyer, “Efficient Raman lasing in tapered silicon waveguides,” Spectroscopy 21(1), 26–32 (2006).

12. D. Dimitropoulos, R. Jhaveri, R. Claps, J. C. S. Woo, and B. Jalali, “Lifetime of photogenerated carriers in silicon-on-insulator rib waveguides,” Appl. Phys. Lett. 86, 071115 (2005). [CrossRef]  

13. S. Fathpour, K. K. Tsia, and B. Jalali, “Energy harvesting in silicon Raman amplifiers,” Appl. Phys. Lett. 89, 061109 (2006). [CrossRef]  

14. A. Liu, H. Rong, M. Paniccia, O. Cohen, and D. Hak, “Net optical gain in a low loss silicon-on-insulator waveguide by stimulated Raman scattering,” Opt. Express 12, 4261 (2004). [CrossRef]   [PubMed]  

15. R. Claps, V. Raghunathan, D. Dimitropoulos, and B. Jalali, “Influence of nonlinear absorption on Raman amplification in Silicon waveguides,” Opt. Express 12, 2774 (2004). [CrossRef]   [PubMed]  

16. T. K. Liang and H. K. Tsang, “Role of free carriers from two-photon absorption in Raman amplification in silicon-on-insulator waveguides,” Appl. Phys. Lett. 84, 2745 (2004). [CrossRef]  

17. T. K. Liang and H. K. Tsang, “Nonlinear absorption and Raman scattering in silicon-on-insulator optical waveguides,” IEEE J. Sel. Top. Quantum Electron. 10, 1149 (2004). [CrossRef]  

18. I. D. Rukhlenko, M. Premaratne, C. Dissanayake, and G. P. Agrawal, “Nonlinear pulse evolution in silicon waveguides: An approximate analytic approach,” IEEE J. Lightwave Technol. 28, in press (2009).

19. I. D. Rukhlenko, M. Premaratne, C. Dissanayake, and G. P. Agrawal, “Continuous-wave Raman amplification in silicon waveguides: beyond the undepleted pump approximation,” Opt. Lett. 34, 536 (2009). [CrossRef]   [PubMed]  

20. S. Roy, S. K. Bhadra, and G. P. Agrawal, “Raman amplification of optical pulses in silicon waveguides: effects of finite gain bandwidth, pulse width, and chirp,” J. Opt. Soc. Am. B 26, 17 (2009). [CrossRef]  

21. B. Jalali, O. Boyraz, V. Raghunathan, D. Dimitropoulos, and P. Koonath, “Silicon Raman amplifiers, lasers and their applications,” in Active and Passive Optical Components for WDM Communications V, A. K. Dutta, Y. Ohishi, N. K. Dutta, and J. Moerk, eds., Proc. SPIE6014, 21–26 (2005).

22. C. Fox, An Introduction to the Calculus of Variations (Dover, New York, 1987).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (3)

Fig. 1.
Fig. 1. Left panel: Evolution of the pump and signal powers along SOI waveguides of three different length (L = 8, 10, 12 mm) with optimized EMA profiles. Right panel: Optimized EMA axial profiles for five waveguide lengths ranging from 8 to 12 mm. In all cases, P p0 = 1 W, P s0 = 0.01 W, and A 0 = 1μm2. See the text for other device parameters.
Fig. 2.
Fig. 2. Left panel: Solution of Eqs. (1) and (8) for a 3-mm-long waveguide with input values P p0 = 1 W, P s0 = 0.5 W, and A 0 = 1 μm2. Solid curves show the pump and signal powers; dashed and dotted curves represent the two auxiliary functions. Right panel: Optimized EMA profile (solid curve, right scale) and the pump and signal intensities (dashed curves, left scale) corresponding to the pump and signal powers shown on the left panel.
Fig. 3.
Fig. 3. Output signal power as a function of output EMA for linearly tapered waveguides of four different lengths. The incident powers correspond to those used in Figs. 1 and 2, and A 0 = 1 μm2. The central dotted line represents a waveguide with constant EMA. The crosses show signal power and modal area at the end of the optimized waveguides in Figs. 1 and 2.

Equations (24)

Equations on this page are rendered with MathJax. Learn more.

d P p d z = α p P p β p P p 2 A eff ζ ps P p P s A eff σ p τ c ( ρ p P p 2 + ρ s P s 2 + ρ ps P p P s ) P p A eff 2 ,
d P s d z = α s P s β s P s 2 A eff ζ sp P s P p A eff σ s τ c ( ρ p P p 2 + ρ s P s 2 + ρ ps P p P s ) P s A eff 2 ,
ζ ps = 2 β ps + 4 g R γ R 2 Ω R Ω ps ( Ω R 2 Ω ps 2 ) 2 + 4 γ R 2 Ω ps 2 ,
ζ sp = ω s ω p ( 2 β ps 4 g R γ R 2 Ω R Ω ps ( Ω R 2 Ω ps 2 ) 2 + 4 γ R 2 Ω ps 2 ) .
G ( A eff ) = 0 L 1 P s d P s d z d z = 0 L g ( z , A eff ) d z ,
g ( z , A eff ) = α s β s P s ( z ) A eff ( z ) ζ ps P p ( z ) A eff ( z ) σ s τ c ρ p P p 2 ( z ) + ρ s P s 2 ( z ) + ρ ps P p ( z ) P s ( z ) A eff 2 ( z ) .
0 L ( 𝒜 δ P p + 𝓑 δ P s ) d z = 0 L 𝒞 δ A eff d z ,
𝒜 = 1 A eff ( ζ sp + σ s τ c 2 ρ p P p + ρ ps P s A eff ) ,
𝓑 = 1 A eff ( β s + σ s τ c 2 ρ s P s + ρ ps P p A eff ) ,
𝒞 = β s P s + ζ sp P p A eff 2 + 2 σ s τ c ρ p P p 2 + ρ s P s 2 + ρ ps P p P s A eff 3 .
d ( δ P p d z ) = a p δ P p b s δ P s c p δ A eff ,
d ( δ P s d z ) = a s δ P s b p δ P p c s δ A eff ,
a p = α p + 2 β p P p A eff + ζ ps P s A eff + σ p τ c 3 ρ p P p 2 + ρ s P s 2 + 2 ρ ps P p P s A eff 2 ,
b p = 𝒜 P s , c p = 𝒞 P p .
0 L φ d ( δ P p ) d z d z = φδ P p | 0 L 0 L φ δ P p d z = 0 L φ ( a p δ P p + b s δ P s c p δ A eff ) d z ,
0 L ψ d ( δ P s ) d z d z = ψδ P s | 0 L 0 L ψ δ P s d z = 0 L ψ ( a s δ P s + b p δ P p c s δ A eff ) d z ,
φ ( L ) = 0 .
0 L [ ( φ a p φ ) δ P p + φ b s δ P s ] d z = 0 L φ c p δ A eff d z ,
0 L [ ψ b p δ P p + ( ψ a s ψ ) δ P s ] d z = 0 L ψ c p δ A eff d z .
0 L [ ( φ a p + ψ b p φ ) δ P p + ( ψ a s + φ b s ψ ) δ P s ] d z = 0 L ( φ c p + ψ c s ) δ A eff d z .
φ c p + ψ c s = 𝒞 ,
φ = a p φ + b p ψ 𝒜 ,
ψ = a s ψ + b s ψ 𝓑 .
A eff ( z ) = 2 τ c [ σ s ( 1 ψ P s ) σ p φ P p ] ( ρ p P p 2 + ρ s P s 2 + ρ ps P p P s ) φ P p ( β p P p + ζ ps P s ) ( 1 ψ P s ) ( β s P s + ζ sp P p ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.