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

Small-signal analysis of bidirectional operating characteristics in a Raman ring laser with external optical injections

Open Access Open Access

Abstract

Employing the small-signal analysis, we study the bidirectional operating characteristics in a Raman ring laser. Bidirectional operation is found not stable in a silicon ring because of the existence of two-photon absorption. Using polar crystals which have different forward and backward Raman gain or introducing external optical injections helps to establish stable bidirectional operation in a Raman ring laser.

©2010 Optical Society of America

1. Introduction

The gyroscope, which is used to detect the rotation, is one of the most potential applications of the ring laser. Soon after the first Raman ring laser on silicon-on-insulator was demonstrated [1], some groups study its gyro characteristics [2]. However, to make a gyro work, one of the preconditions is the coexistence of two laser beams travel in opposite directions in the ring cavity [3]. In this paper, we employed the small-signal analysis to study the bidirectional operating characteristics in a Raman ring laser. The material in discussion is not only silicon but also polar crystals, e.g. GaP, which provides different forward and backward Raman gains.

2. Modeling

In our modeling, we have four counter propagating waves inside the resonator: two quasi-TM polarizing (dominant vertical component of electric field) pump beams (“p”) and two quasi-TE polarizing (dominant horizontal component) Stokes waves (“s”). The starting point in the development of the ring laser model is the time differential intensity equations for each clockwise (CW, “−”) or counterclockwise (CCW, “+”) propagating wave [4]:

I˙p,±=(vg/L)ηp(Ip,±extIp,±)Ip,±/τpcv{vg[(λs/λp)gR(κfIs,±+κbIs,)+β(Ip,±+2Ip,+2Is,++2Is,)+φ¯λp2Neff]}Ip,±I˙s,±=(vg/L)ηs(Is,±extIs,±)Is,±/τscv+{vg[gR(κfIp,±+κbIp,)β(Is,±+2Is,+2Ip,++2Ip,)φ¯λs2Neff]}Is,±

In the equations, τp,scvis the life time of a pump or Stokes photon in the ring resonator, whose wavelengths in vacuum areλp,srespectively. It is assumed that both the pump and Stokes waves have the same group velocityvgin the ring cavity with a length L. The stimulated Raman scattering (SRS) with the gain factorgR, two-photon absorption with the absorption factor β, and the free carrier absorption with the effective factorφ¯are taken into concern in this model. The forward and backward SRS coefficientsκf,bcan be set different because the forward and backward SRS efficiencies in some polar crystals [5], e.g. GaP [6], can be different though they are all the same in most materials, e.g. silicon [2,4]. The effective free carrier number can be expressed as follows

Neff=βτeff2ωp(Ip,++Ip,+Is,++Is,)2
whereωpis the energy of a photon andτeffis the effective carrier life time [7,8]. The external optical injectionsIp,±extandIs,±extare coupled into the ring though a coupler, whose cross coupling coefficients for pump and Stokes lights areηpandηs, respectively. In general, ηs is designed as small as possible to reduce the loss of the Stokes radiation, butηpis large to enhance the pump injection. The backscattering is weak and has little effect to the intensity equations, so it is omitted in the modeling. By defining
k=λs/λp,G=vggR,B=vgβ,T=L/vg,F=φ¯λp2Bτeff/(2ωp),τp,scp=ηp.s/T,1/τp,s=1/τp,scv+1/τp,scp
the equations are simplified as

I˙p,±=Ip,±ext/τpcpIp,±/τp+[kG(κfIs,±+κbIs,)B(Ip,±+2Ip,+2Is,++2Is,)F(Ip,++Ip,+Is,++Is,)2]Ip,±I˙s,±=Is,±ext/τscpIs,±/τs+[G(κfIp,±+κbIp,)B(Is,±+2Is,+2Ip,++2Ip,)k2F(Ip,++Ip,+Is,++Is,)2]Is,±

Without Stokes light emission, the pump intensities areIp,±(0). And considering the near threshold situation, both the Stokes waves are so weak that small-signal analysis works. So Is,+/Ip,±(0)1, Is,/Ip,±(0)1, and the pump intensity can be approximately written in the form

Ip,±=Ip,±(0)[1+a±+(Is,+/Ip,±(0))+a±(Is,/Ip,±(0))+ο(Is,+/Ip,±(0))]Ip,±(0)+a±+Is,++a±Is,

By substituting (4) and (5) into the following equations,

I˙p,±=a±+I˙s,++a±I˙s,
omitting the higher order terms, and then by identifying the coefficients of Is,±on both sides, it can be deduced that
(1/τp1/τs+GW++FU+)a+++2VIp,+(0)a++(κfkG+2V)Ip,+(0)=0(1/τp1/τs+GW++FU)a++2VIp,(0)a+++(κbkG+2V)Ip,(0)=0(1/τp1/τs+GW+FU+)a++2VIp,+(0)a+(κbkG+2V)Ip,+(0)=0(1/τp1/τs+GW+FU)a+2VIp,(0)a++(κfkG+2V)Ip,(0)=0
where

U±=(Ip,+(0)+Ip,(0))[(1k2)(Ip,+(0)+Ip,(0))+2Ip,±(0)]W±=κfIp,±(0)+κbIp,(0),V=B+F(Ip,+(0)+Ip,(0))

AndIp,±(0)satisfy the following equations

[1/τp+B(2Ip,(0)+Ip,±(0))+F(Ip,+(0)+Ip,(0))2]Ip,±(0)=Ip,±ext/τpcpa±+Is,+ext/τscpa±Is,ext/τscp

The coefficients can be obtained by solving the linear Eqs. (7). With the help of these coefficients, we substitute the pump intensities in the time differential equations of Stokes light, and rewrite the equations in the follow forms with second-order (third-order on the point of electric field) terms, such as the self- and cross-saturation effect.

I˙s,±=Is,±ext/τscp+{ρ±+BIs,±+[G(κfa±++κba+)θ+]Is,++[G(κfa±+κba)θ]Is,}Is,±
where the zero-order gain/loss
ρ±=1/τs+G(κfIp,±(0)+κbIp,(0))2B(Ip,+(0)+Ip,(0))k2F(Ip,+(0)+Ip,(0))2
and the saturation factor

θ±=2B(1+a+±+a±)+2k2F[(a+±+a±+1)Ip,+(0)+(a±+a+±+1)Ip,(0)]

3. Discussions

Consider the simplest situation, the injection is symmetric, that isIp,±ext=IpextandIs,±ext=Isext. An obvious consequence is that both pump intensities Ip,±(0)are the sameIp(0) when Stokes lights do not emit. So the equations for these counter running Stokes waves become

I˙s,±=(ραslfIs,±αcrsIs,)Is,±+Isext/τscp
where the self- and cross-saturation factors are defined as
αslf=G[(κf+κb)aB+(κf2+κb2)aG+(κfκb)2aΔ]+θBαcrs=G[(κf+κb)aB+2κfκbaG(κfκb)2aΔ]+θ
and

ρ=1/τs+(κf+κb)GIp(0)4BIp(0)4k2F(Ip(0))2θ=2[12aB(κf+κb)aG](B+2k2FIp(0))aB=2VIp(0)/(1/τp+ρ+6VIp(0)),aG=kGIp(0)/(1/τp+ρ+6VIp(0)),aΔ=2VkG(Ip(0))2/[(1/τp+ρ+6VIp(0))(1/τp+ρ+2VIp(0))]

For the first case, suppose there is no injection at the wavelength of the Stokes light, that isIsext=0. At the stationary solutionIs,±(s) which satisfy the stationary condition

(ραslfIs,±(s)αcrsIs,(s))Is,±(s)=0
we consider small deviations of the intensities about them

Is,±=Is,±(s)+ε±

Substituting them into (13), we find

ddt(ε+ε)=(αslfIs,+(s)αcrsIs,+(s)αcrsIs,(s)αslfIs,(s))(ε+ε)Θ(ε+ε)

The stationary solution is stable only if both the eigenvalues of the stability matrixΘ are negative. For a second order matrix, this requires its trace is negative and its determinant is positive. This leads to

αslf2αcrs2>0,αslf>0

It is easy to prove that bothαslfandαcrsare positive in a normally operating laser, and

αslfαcrs=G(κfκb)2(aG+2aΔ)B=(κfκb)2kG2Ip(0)1/τp+ρ+2VIp(0)B

This reveals that if the backscattering between the two Stokes beams is weak, such as in a laser gyroscope, the ring laser with materials of which SRS has no directional preference, e.g. silicon, does not support stable bidirectional operating state. Becauseαslfis always smaller thanαcrsin this case. This can be explained in a positive feedback way. When the ring laser works at the bidirectional operating stationary point, as the forward and backward Raman gains are all the same, the only difference between the gain/loss factors of these two oppositely running Stokes beams originates from the difference between self- and cross-TPA. Once intensities of these two Stokes beams depart a little from the stationary point, the departure will be enlarged because the cross-TPA is twice as large as self-TPA. This positive feedback makes the bidirectional stationary operating point unstable.

But for those polar crystals whose Raman scattering has directional preference, e.g. longitude phonon scattering in GaP only occurs in the backward direction, stable bidirectional operating can be realized if

(κfκb)2kG2Ip(0)B(1/τp+ρ+2VIp(0))>0

Combining with the lasing condition: the zero order net gain ρ>0, we have

0<(κf+κb)G[1/(τsIp(0))+4B+4k2FIp(0)]<(κfκb)2kG2/B[1/(τpIp(0))+2B+4FIp(0)]

So a difference between the forward and backward gains is needed to realize stable bidirectional operation. The left diagram in Fig. 1 gives an example of the dependence of the operating regime to the forward and backward scattering gains. The dashed line is the threshold. Only in the region above it where ρ>0, laser light emits. This region is divided into three areas by the solid curves, which are obtained by setting the two expressions in (22) to be equal. In the middle area where the stability condition, the right part of (22), is not satisfied, the ring laser works in unidirectional operating regime. Only in the two side areas where there is a distinct difference between the forward and backward gains, stable bidirectional operation can be obtained. To justify this conclusion, we simulate A~D four situations by the side of the threshold line. The numerical simulation is based on the discretisation of (1) by simply replacing the time derivatives with finite differences. The simulation results agree with the distribution of these work points. In the cases A and B, though the laser cavity is symmetric pumped, there is always one suppressed Stokes mode, CW or CCW, depending on the initial condition. It should be noticed that case A is the situation that the material have the same forward and backward scattering gain. But in the case C and D, both CW and CCW Stokes modes can coexist and be stable in the ring cavity after sufficiently long simulating time.

 figure: Fig. 1

Fig. 1 Theory and numerical simulation results illustrating the occurrence of the two different ring laser regimes (unidirectional and bidirectional). The left diagram is the small-signal analysis results from (22). Figures A, B, C, and D are, respectively, the simulation results corresponding to the points A, B, C, and D in the left diagram. The numerical simulation is in the case of symmetric pumped. Values of other parameters, common to these diagrams, are: λs = 1550nm, λp = 1646nm, φ¯ = 6×10−10, τeff = 0.1ns, β = 0.6cm2/GW, vgτp = 2.35cm, vgτs = 8cm, Isext = 0, Ipext/vgτpcp = 0.15GW/cm3, and the corresponding Ip(0) = 0.1648GW/cm2.

Download Full Size | PDF

As discussed above, when there is no Stokes light injected, those materials of which SRS has no directional preference do not support stable bidirectional emission. Now we add such injections. By settingκf=κb=1, we have

ρ=1/τs+2GIp(0)4BIp(0)4k2F(Ip(0))2,αcrs=αslf+B=2G(aB+aG)+θ

Consider the stationary solutionIs,+(s)=Is,(s)=Is(s), which satisfy the stationary condition

0=(2αslf+B)(Is(s))2+ρIs(s)+Isext/τscp

Appling the above vibrations analysis, we have the stability matrix

Θ=(A+αslfIs(s)αcrsIs(s)αcrsIs(s)A+αslfIs(s)),whereA=Isext/(Is(s)τscp)

Again, the stationary solution is stable only if the trace is negative and the determinant is positive. Because all these variables are positive, it only requires that

Isext/τscp>B(Is(s))2

Combining with the stationary condition, we finally have

Is(s)>ρ/(2αslf),Isext>Bρ2τscp/(4αslf2)

Even in a laser with the same forward and backward Raman gain and no backscattering, bidirectional operating is possible if the injection is strong enough. The conclusion is illustrated in Fig. 2 . As show in the left figure, when the extent injections are not strong enough, the competition between these two opposite propagating Stokes modes induces asymmetric intensities of them. Once the injections are strong enough, as shown in the right, both CW and CCW modes have the same intensity in the end.

 figure: Fig. 2

Fig. 2 Simulation results ring laser regimes with injections at Stokes wavelength: Isext/vgτscp = 2kW/cm3 (left) and 3kW/cm3 (right), κfgR/β=κbgR/β = 5 and values of other parameters are the same with those in Fig. 1.

Download Full Size | PDF

All the discussion above is based on the small-signal analysis, which requires the Raman ring laser work in the near threshold region. Once the intensity of the Stokes emission is no longer much smaller than that of the pump, all of these conclusions are invalid. In such region, further study is still needed.

4. Conclusions

Employing the small-signal analysis, we studied the bidirectional operating characteristics in a Raman ring laser. If without injections, bidirectional operating only occurs in the ring made up of polar crystals, and it requires that the difference between forward and backward stimulated Raman gains is large enough. Further, with strong enough external injections, bidirectional operating can also be obtained in a ring laser made up by general crystal, e.g. silicon, which has the same forward and backward gains.

References and links

1. H. Rong, Y. H. Kuo, S. Xu, A. Liu, R. Jones, M. Paniccia, O. Cohen, and O. Raday, “Monolithic integrated Raman silicon laser,” Opt. Express 14(15), 6705–6712 (2006). [CrossRef]   [PubMed]  

2. F. De Leonardis and V. M. N. Passaro, “Modeling and Performance of a Guided-Wave Optical Angular-Velocity Sensor Based on Raman Effect in SOI,” J. Lightwave Technol. 25(9), 2352–2366 (2007). [CrossRef]  

3. W. W. Chow, J. Gea-Banacloche, L. M. Pedrotti, V. E. Sanders, W. Schleich, and M. O. Scully, “The ring laser gyro,” Rev. Mod. Phys. 57(1), 61–104 (1985). [CrossRef]  

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

5. R. Loudon, “The Raman effect in crystals,” Adv. Phys. 50(7), 813–864 (2001). [CrossRef]  

6. T. Saito, K. Suto, J. Nishizawa, and M. Kawasaki, “Spontaneous Raman scattering in [100], [110], and [11-2] directional GaP waveguides,” J. Appl. Phys. 90(4), 1831–1835 (2001). [CrossRef]  

7. H. Rong, A. Liu, R. Nicolaescu, M. Paniccia, O. Cohen, and D. Hak, “Raman gain and nonlinear optical absorption measurements in a low-loss silicon waveguide,” Appl. Phys. Lett. 85(12), 2196–2198 (2004). [CrossRef]  

8. F. De Leonardis, V. Dimastrodonato, and V. M. N. Passaro, “Modelling of a DBR laser based on Raman effect in a silicon-on-insulator rib waveguide,” Semicond. Sci. Technol. 23(6), 064008 (2008). [CrossRef]  

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 (2)

Fig. 1
Fig. 1 Theory and numerical simulation results illustrating the occurrence of the two different ring laser regimes (unidirectional and bidirectional). The left diagram is the small-signal analysis results from (22). Figures A, B, C, and D are, respectively, the simulation results corresponding to the points A, B, C, and D in the left diagram. The numerical simulation is in the case of symmetric pumped. Values of other parameters, common to these diagrams, are: λ s = 1550nm, λ p = 1646nm, φ ¯ = 6×10−10, τ e f f = 0.1ns, β = 0.6cm2/GW, v g τ p = 2.35cm, v g τ s = 8cm, I s e x t = 0, I p e x t / v g τ p c p = 0.15GW/cm3, and the corresponding I p ( 0 ) = 0.1648GW/cm2.
Fig. 2
Fig. 2 Simulation results ring laser regimes with injections at Stokes wavelength: I s e x t / v g τ s c p = 2kW/cm3 (left) and 3kW/cm3 (right), κ f g R / β = κ b g R / β = 5 and values of other parameters are the same with those in Fig. 1.

Equations (27)

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

I ˙ p , ± = ( v g / L ) η p ( I p , ± e x t I p , ± ) I p , ± / τ p c v { v g [ ( λ s / λ p ) g R ( κ f I s , ± + κ b I s , ) + β ( I p , ± + 2 I p , + 2 I s , + + 2 I s , ) + φ ¯ λ p 2 N e f f ] } I p , ± I ˙ s , ± = ( v g / L ) η s ( I s , ± e x t I s , ± ) I s , ± / τ s c v + { v g [ g R ( κ f I p , ± + κ b I p , ) β ( I s , ± + 2 I s , + 2 I p , + + 2 I p , ) φ ¯ λ s 2 N e f f ] } I s , ±
N e f f = β τ e f f 2 ω p ( I p , + + I p , + I s , + + I s , ) 2
k = λ s / λ p , G = v g g R , B = v g β , T = L / v g , F = φ ¯ λ p 2 B τ e f f / ( 2 ω p ) , τ p , s c p = η p . s / T , 1 / τ p , s = 1 / τ p , s c v + 1 / τ p , s c p
I ˙ p , ± = I p , ± e x t / τ p c p I p , ± / τ p + [ k G ( κ f I s , ± + κ b I s , ) B ( I p , ± + 2 I p , + 2 I s , + + 2 I s , ) F ( I p , + + I p , + I s , + + I s , ) 2 ] I p , ± I ˙ s , ± = I s , ± e x t / τ s c p I s , ± / τ s + [ G ( κ f I p , ± + κ b I p , ) B ( I s , ± + 2 I s , + 2 I p , + + 2 I p , ) k 2 F ( I p , + + I p , + I s , + + I s , ) 2 ] I s , ±
I p , ± = I p , ± ( 0 ) [ 1 + a ± + ( I s , + / I p , ± ( 0 ) ) + a ± ( I s , / I p , ± ( 0 ) ) + ο ( I s , + / I p , ± ( 0 ) ) ] I p , ± ( 0 ) + a ± + I s , + + a ± I s ,
I ˙ p , ± = a ± + I ˙ s , + + a ± I ˙ s ,
( 1 / τ p 1 / τ s + G W + + F U + ) a + + + 2 V I p , + ( 0 ) a + + ( κ f k G + 2 V ) I p , + ( 0 ) = 0 ( 1 / τ p 1 / τ s + G W + + F U ) a + + 2 V I p , ( 0 ) a + + + ( κ b k G + 2 V ) I p , ( 0 ) = 0 ( 1 / τ p 1 / τ s + G W + F U + ) a + + 2 V I p , + ( 0 ) a + ( κ b k G + 2 V ) I p , + ( 0 ) = 0 ( 1 / τ p 1 / τ s + G W + F U ) a + 2 V I p , ( 0 ) a + + ( κ f k G + 2 V ) I p , ( 0 ) = 0
U ± = ( I p , + ( 0 ) + I p , ( 0 ) ) [ ( 1 k 2 ) ( I p , + ( 0 ) + I p , ( 0 ) ) + 2 I p , ± ( 0 ) ] W ± = κ f I p , ± ( 0 ) + κ b I p , ( 0 ) , V = B + F ( I p , + ( 0 ) + I p , ( 0 ) )
[ 1 / τ p + B ( 2 I p , ( 0 ) + I p , ± ( 0 ) ) + F ( I p , + ( 0 ) + I p , ( 0 ) ) 2 ] I p , ± ( 0 ) = I p , ± e x t / τ p c p a ± + I s , + e x t / τ s c p a ± I s , e x t / τ s c p
I ˙ s , ± = I s , ± e x t / τ s c p + { ρ ± + B I s , ± + [ G ( κ f a ± + + κ b a + ) θ + ] I s , + + [ G ( κ f a ± + κ b a ) θ ] I s , } I s , ±
ρ ± = 1 / τ s + G ( κ f I p , ± ( 0 ) + κ b I p , ( 0 ) ) 2 B ( I p , + ( 0 ) + I p , ( 0 ) ) k 2 F ( I p , + ( 0 ) + I p , ( 0 ) ) 2
θ ± = 2 B ( 1 + a + ± + a ± ) + 2 k 2 F [ ( a + ± + a ± + 1 ) I p , + ( 0 ) + ( a ± + a + ± + 1 ) I p , ( 0 ) ]
I ˙ s , ± = ( ρ α s l f I s , ± α c r s I s , ) I s , ± + I s e x t / τ s c p
α s l f = G [ ( κ f + κ b ) a B + ( κ f 2 + κ b 2 ) a G + ( κ f κ b ) 2 a Δ ] + θ B α c r s = G [ ( κ f + κ b ) a B + 2 κ f κ b a G ( κ f κ b ) 2 a Δ ] + θ
ρ = 1 / τ s + ( κ f + κ b ) G I p ( 0 ) 4 B I p ( 0 ) 4 k 2 F ( I p ( 0 ) ) 2 θ = 2 [ 1 2 a B ( κ f + κ b ) a G ] ( B + 2 k 2 F I p ( 0 ) ) a B = 2 V I p ( 0 ) / ( 1 / τ p + ρ + 6 V I p ( 0 ) ) , a G = k G I p ( 0 ) / ( 1 / τ p + ρ + 6 V I p ( 0 ) ) , a Δ = 2 V k G ( I p ( 0 ) ) 2 / [ ( 1 / τ p + ρ + 6 V I p ( 0 ) ) ( 1 / τ p + ρ + 2 V I p ( 0 ) ) ]
( ρ α s l f I s , ± ( s ) α c r s I s , ( s ) ) I s , ± ( s ) = 0
I s , ± = I s , ± ( s ) + ε ±
d d t ( ε + ε ) = ( α s l f I s , + ( s ) α c r s I s , + ( s ) α c r s I s , ( s ) α s l f I s , ( s ) ) ( ε + ε ) Θ ( ε + ε )
α s l f 2 α c r s 2 > 0 , α s l f > 0
α s l f α c r s = G ( κ f κ b ) 2 ( a G + 2 a Δ ) B = ( κ f κ b ) 2 k G 2 I p ( 0 ) 1 / τ p + ρ + 2 V I p ( 0 ) B
( κ f κ b ) 2 k G 2 I p ( 0 ) B ( 1 / τ p + ρ + 2 V I p ( 0 ) ) > 0
0 < ( κ f + κ b ) G [ 1 / ( τ s I p ( 0 ) ) + 4 B + 4 k 2 F I p ( 0 ) ] < ( κ f κ b ) 2 k G 2 / B [ 1 / ( τ p I p ( 0 ) ) + 2 B + 4 F I p ( 0 ) ]
ρ = 1 / τ s + 2 G I p ( 0 ) 4 B I p ( 0 ) 4 k 2 F ( I p ( 0 ) ) 2 , α c r s = α s l f + B = 2 G ( a B + a G ) + θ
0 = ( 2 α s l f + B ) ( I s ( s ) ) 2 + ρ I s ( s ) + I s e x t / τ s c p
Θ = ( A + α s l f I s ( s ) α c r s I s ( s ) α c r s I s ( s ) A + α s l f I s ( s ) ) , where A = I s e x t / ( I s ( s ) τ s c p )
I s e x t / τ s c p > B ( I s ( s ) ) 2
I s ( s ) > ρ / ( 2 α s l f ) , I s e x t > B ρ 2 τ s c p / ( 4 α s l f 2 )
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.