arxiv: v3 [quant-ph] 25 Sep 2012

Transcription

1 Three-Level Laser Dynamics with the Atoms Pumpe by Electron Bombarment arxiv: v3 [quant-ph] 25 Sep 212 Fesseha Kassahun Department of Physics, Ais Ababa University, P. O. Box 33761, Ais Ababa, Ethiopia (Submitte on 25 Sep 212) Abstract We analyze the quantum properties of the light generate by a three-level laser with a close cavity an couple to a vacuum reservoir. The three-level atoms available in the cavity are pumpe from the bottom to the top level by means of electron bombarment an we carry out our analysis by putting the noise operators associate with the vacuum reservoir in normal orer. The maximum quarature squeezing of the light generate by the laser, operating far below threshol, is foun to be 5% below the coherent-state level. We have also establishe that the quarature squeezing of the output light is equal to that of the cavity light an has the same value in any frequency interval. This implies that the quarature squeezing of the laser light is an intrinsic property of the iniviual photons. Keywors: Stimulate emission, Photon statistics, Quarature squeezing, Noiseless vacuum reservoir 1. Introuction A three-level laser is a quantum optical system in which light is generate by three-level atoms insie a cavity usually couple to a vacuum reservoir. In one moel of a such laser, three-level atoms initially prepare in a coherent superposition of the top an bottom levels are injecte into a cavity an then remove after they have ecaye ue to spontaneous emission [1,2]. In another moel, the top an bottom levels of the three-level atoms injecte into a cavity are couple by coherent light [3,4]. The statistical an squeezing properties of the light generate by three-level lasers have been investigate by several authors [5-11]. It is foun that a three-level laser in either moel generates squeeze light uner certain conitions. It appears to be quite ifficult to prepare the atoms in a coherent superposition of the top an bottom levels before they are injecte into the laser cavity. In aition, it shoul certainly be har to fin out that the atoms have ecaye spontaneously before they are remove from the cavity. On the other han, the egree of squeezing of the light generate by the three-level laser, with the top an bottom levels couple by coherent light, is relatively large when the mean photon number is relatively small [4]. aress: fesseha.kassahun@aau.eu.et 1

2 Moreover, the quantum analysis of a three-level laser is usually carrie out by incluing the interaction of the atoms insie the cavity with the vacuum reservoir outsie the cavity. It may be possible to justify the feasibility of such interaction for a laser with an open cavity into which an from which atoms are injecte an remove. However, there cannot be any vali justification for leaving open the laser cavity in which the available atoms are pumpe to the top level by electron bombarment. Therefore, the aforementione interaction is not feasible for a laser in which the atoms available in a close cavity are pumpe to the top level by means of electron bombarment. We seek here to analyze the quantum properties of the light emitte by the three-level atoms available in a close cavity an pumpe to the top level at a constant rate. Thus taking into account the interaction of the three-level atoms with a resonant cavity light an the amping of the cavity light by a vacuum reservoir, we obtain the photon statistics of the cavity light an the quarature squeezing of the cavity (output) light. We also etermine the quarature squeezing of the cavity (output) light in any frequency interval. We carry out our calculation by putting the noise operators associate with the vacuum reservoir in normal orer an without consiering the interaction of the three-level atoms with the vacuum reservoir outsie the cavity. 2. Operator ynamics We consier here the case in whichn egenerate three-level atoms in cascae configuration are available in a close cavity. We enote the top, mile, an bottom levels of these atoms by a k, b k, an c k, respectively. In aition, we assume the cavity moe to be at resonance with the two transitions a k b k an b k c k, with irect transition between levels a k an c k to be ipole forbien. The interaction of one of the three-level atoms with the cavity moe can be escribe at resonance by the Hamiltonian [ ] Ĥ = ig (ˆσ a k + ˆσ k b )ˆb ˆb (ˆσ a k + ˆσk b ), (1) where an ˆσ k a = b k k a (2) ˆσ k b = c k k b (3) are lowering atomic operators, ˆb is the annihilation operator for the cavity moe, an g is the coupling constant between the atom an the cavity moe. We assume that the laser cavity is couple to a vacuum reservoir via a single-port mirror. In aition, we carry out our calculation by putting the noise operators associate with the vacuum reservoir in normal orer. Thus the noise operators will not have any effect on the ynamics of the cavity moe operators. We can therefore rop the noise operator an write the quantum Langevin equation for the operatorˆb as ˆb t = 2ˆb i[ˆb,ĥ], (4) whereis the cavity amping constant. Then with the ai of Eq. (1), we easily fin ˆb t = a 2ˆb g(ˆσ k + ˆσk b ). (5) Furthermore, applying the relation  = i [Â,Ĥ] (6) t 2

3 along with Eq. (1), one can reaily establish that t ˆσk a = g ˆηk bˆb g ˆη aˆb +g ˆb k ˆσ c k, (7) t ˆσk b = g ˆηk cˆb g ˆη bˆb g ˆb k ˆσ c k, (8) t ˆηk a = g ˆσ k a a k, (9) where t ˆηk b = g ˆσ k k b ˆb g ˆσ ˆb +g ˆb ˆσ a b k g ˆb ˆσ a, k (1) t ηk c = g σ kˆb g ˆb ˆσ k b b, (11) ˆσ k c = c k k a, (12) ˆη k a = a k k a, (13) ˆη k b = b k k b, (14) ˆη k c = c k k c. (15) We see that Eqs. (7)-(11) are nonlinear ifferential equations an hence it is not possible to fin exact time-epenent solutions of these equations. We inten to overcome this problem by applying the large-time approximation [12]. Then using this approximation scheme, we get from Eq. (5) the approximately vali relation 2g ˆb = (ˆσk a + ˆσk b ). (16) Eviently, this turns out to be an exact relation at steay state. Now introucing Eq. (16) into Eqs. (7)-(11), we get t ˆσk a = γ c ˆσ k a, (17) where t ˆσk b = γ c 2 ˆσk b +γ c ˆσ a k, (18) t ˆηk a = γ c ˆη a k, (19) t ˆηk b = γ c ˆη b k +γ c ˆη a, k (2) t ˆηk c = γ c ˆη b k, (21) γ c = 4g2. (22) We prefer to call the parameter efine by Eq. (22) the stimulate emission ecay constant. The three-level atoms available in the cavity are pumpe from the bottom to the top level by means of electron bombarment. The pumping process must surely affect the ynamics of ˆη k a an ˆη k c. If r a represents the rate at which a single atom is pumpe from the bottom to the top 3

4 level, then ˆη a k increases at the rate ofr a ˆη c k an ˆηk c ecreases at the same rate. In view of this, we rewrite Eqs. (19) an (21) as an t ˆηk a = γ c ˆη k a +r a ˆη k c (23) t ˆηk c = γ c ˆη b k r a ˆη c k. (24) We next sum Eqs. (17), (18), (2), (23), an (24) over then three-level atoms, so that in which t ˆm a = γ c ˆm a, (25) t ˆm b = γ c 2 ˆm b +γ c ˆm a, (26) t ˆN a = γ c ˆN a +r a ˆN c, (27) t ˆN b = γ c ˆN b +γ c ˆN a, (28) t ˆN c = γ c ˆN b r a ˆN c, (29) ˆm a = ˆm b = ˆN a = ˆN b = ˆN c = N ˆσ a k, (3) k=1 N ˆσ b k, (31) k=1 N ˆη a, k (32) k=1 N ˆη b k, (33) k=1 N ˆη c, k (34) k=1 with the operators ˆN a, ˆNb, an ˆN c representing the number of atoms in the top, mile, an bottom levels. In aition, employing the completeness relation ˆη a k + ˆη b k + ˆηk c = Î, (35) we easily arrive at ˆN a + ˆN b + ˆN c = N. (36) Furthermore, using the efinition given by (3) an setting for anyk ˆσ k a = b a, (37) we have ˆm a = N b a. (38) 4

5 Following the same proceure, one can also check that where Moreover, using the efinition ˆm b = N c b, (39) ˆm c = N c a, (4) ˆN a = N a a, (41) ˆN b = N b b, (42) ˆN c = N c c, (43) ˆm c = N ˆσ c. k (44) k=1 ˆm = ˆm a + ˆm b (45) an taking into account Eqs. (38)-(43), it can be reaily establishe that In the presence ofn three-level atoms, we rewrite Eq. (5) as ˆm ˆm = N( ˆN a + ˆN b ), (46) ˆmˆm = N( ˆN b + ˆN c ), (47) ˆm 2 = N ˆm c. (48) ˆb t = 2ˆb+λˆm, (49) in which λ is a constant whose value remains to be fixe. Applying the steay-state solution of Eq. (5), we get [ˆb,ˆb ] k = γ c (ˆηk c ˆη k a) (5) an on summing over all atoms, we have where [ˆb,ˆb ] = γ c ( ˆN c ˆN a ), (51) [ˆb,ˆb ] = N [ˆb,ˆb ] k (52) k=1 stans for the commutator ofˆb anˆb when the cavity moe is interacting with all the N threelevel atoms. On the other han, using the steay-state solution of Eq. (49), one can easily verify that ( ) 2λ 2 [ˆb,ˆb ] = N ( ˆN c ˆN a ). (53) Thus on account of Eqs. (51) an (53), we see that an in view of this result, Eq. (49) can be written as λ = ± g N (54) ˆb t = 2ˆb+ g ˆm. (55) N 5

6 3. Photon statistics We next seek to calculate the mean an variance of the photon number at steay state. However, we nee first to establish several important relations. Hence with the ai of (36), one can put Eq. (27) in the form t ˆN a = (γ c +r a ) ˆN a +r a (N ˆN b ) (56) an application of the large-time approximation scheme to Eq. (28) yiels ˆN b = ˆN a. (57) Thus on taking into account this result, we fin the steay-state solution of Eq. (56) to be We also notice that at steay state Eq. (27) leas to ˆN a = r an γ c +2r a. (58) ˆN c = γ c r a ˆN a. (59) We now procee to calculate the expectation value of the atomic operator ˆm c. We assume that the state vector of a three-level atom, put in the form ψ k = C a a k +C b b k +C c c k, (6) can be use to etermine the expectation value of an atomic operator forme by a pair of ientical energy levels or by two istinct energy levels between which transition with the emission of a photon is ipole forbien. One can thus reaily establish that an C a Ca = ˆη a, (61) C c C c = ˆη c, (62) ˆσ k c = C ac c. (63) In orer to have a mathematically manageable analysis, we take ˆσ c k to be real, so that an on subtractingc a C c from both sies of this equation, we have Hence for all possible values ofc a anc c, we see that an C a C c = C ac c (64) C a (Cc C c)+c c (C a Ca ) =. (65) C a = C a (66) C c = C c. (67) Now on account of (61) an (62) together with (66) an (67), Eq. (63) takes the form ˆσ c ˆη k = a k ˆηk c. (68) 6

7 Finally on summing overk from 1 up ton an taking into account Eqs. (32) an (34), we have ˆm c = ˆN a ˆN c. (69) Furthermore, aing Eqs. (25) an (26), we have t ˆm = 1 2 γ c ˆm γ c ˆm a. (7) In orer to inclue the effect of the pumping process, we rewrite this equation as ˆm t = 1 2 µˆm+ 1 2 µˆm a + ˆF m (t), (71) in which ˆF m (t) is a noise operator with vanishing mean anµis a parameter whose value remains to be etermine. Employing the relation ˆm ˆm ˆm = ˆm + ˆm ˆm, (72) t t t along with Eq. (71), we easily fin from which follows t ˆm ˆm = µ ˆm ˆm +µ ˆm aˆm a + ˆF m(t)ˆm + ˆm ˆFm (t), (73) t ˆN a + ˆN b = µ ˆN a + ˆN b +µ ˆN a + 1 N ˆF m (t)ˆm + 1 N ˆm ˆFm (t). (74) On the other han, on account of Eqs. (27) an (28), we see that This can be rewritten as t ˆN a + ˆN b = γ c ˆN b +r a ˆN c. (75) t ˆN a + ˆN b = (γ c +2r a ) ˆN a + ˆN b +(γ c +2r a ) ˆN a +r a N. (76) Hence comparison of Eqs. (74) an (76) shows that an We observe that Eq. (78) is equivalent to Upon casting Eq. (25) into the form one can also easily verify thatµhas the value given by (77) an µ = γ c +2r a (77) ˆF m (t)ˆm + ˆm ˆFm (t) = r a N 2. (78) ˆF m(t)ˆf m (t ) = r a N 2 δ(t t ). (79) t ˆm a = 1 2 µˆm a + ˆF a (t), (8) ˆF a (t)ˆf a (t ) = r a N 2 δ(t t ). (81) 7

8 With the atoms consiere to be initially in the bottom level, the expectation value of the solution of Eq. (8) happens to be ˆm a (t) =. (82) Then on account of this result, the expectation value of the solution of Eq. (71) turns out to be ˆm(t) =. (83) On the other han, the expectation value of the solution of Eq. (55) is expressible as ˆb(t) = ˆb() e t/2 + g t e t/2 e t /2 ˆm(t ) t. (84) N Now in view of Eq. (83) an the assumption that the cavity light is initially in a vacuum state, Eq. (84) goes over into ˆb(t) =. (85) We observe on the basis of Eqs. (55) an (85) thatˆb is a Gaussian variable with zero mean. The steay-state solution of Eq. (55) is given by 2g ˆb = ˆm. (86) N Now using Eq. (86) along with (46), the mean photon number of the cavity light is expressible as ˆb ˆb ( ) γ c = ˆN a + ˆN b. (87) It proves to be convenient to refer to the regime of laser operation with more atoms in the top level than in the bottom level as above threshol, the regime of laser operation with equal number of atoms in the top an bottom levels as threshol, an the regime of laser operation with less atoms in the top level than in the bottom level as below threshol. Thus accoring to Eq. (59) for the laser operating above threshol γ c < r a, for the laser operating at threshol γ c = r a, an for the laser operating below threshol γ c > r a. We note that for the laser operating well above threshol, Eq. (87) reuces to n = γ c N (88) an for the laser operating at threshol, we fin in which n = ˆb ˆb. n = 2γ c N, (89) 3 The variance of the photon number for the cavity light is expressible as ( n) 2 = ˆb ˆbˆb ˆb ˆb ˆb 2 (9) an using the fact thatˆb is a Gaussian variable with zero mean, we reaily get ( n) 2 = ˆb ˆb ˆbˆb + ˆb 2 ˆb 2. (91) Employing once more Eq. (86) an taking into account (47), we fin ˆbˆb = γ ( ) c ˆN b + ˆN c. (92) 8

9 In aition, taking into account (59), one can express Eq. (69) in the form γc ˆm c = r ˆN a. (93) a Hence with the ai of Eqs. (86), (48), an (93), we easily get ˆb 2 = γ c γc r ˆN a. (94) a Now on account of Eqs. (91), (92), an (94) along with (57) an (59), we arrive at ( n) 2 = 1 4 n2 (2η +η +2), (95) withη = γ c /r a. We then see that for the laser operating well above threshol ( n) 2 = 1 2 n2, (96) wherenis given by Eq. (88). An for the laser operating at threshol, we get in whichngiven by Eq. (89). ( n) 2 = 5 4 n2, (97) 4. Quarature squeezing In this section we wish to calculate the quarature squeezing of the cavity an output light in any frequency interval. 4.1 Global quarature squeezing The squeezing properties of the cavity light are escribe by two quarature operators efine by ˆb+ =ˆb +ˆb (98) an ˆb = i(ˆb ˆb). (99) It can be reaily establishe that It then follows that [ˆb,ˆb + ] = 2i γ c ( ˆN a ˆN c ). (1) b + b γ c The variance of the quarature operators is expressible as ˆN a ˆN c. (11) ( b ± ) 2 = ± [ˆb ±ˆb] 2 [ ˆb ± ˆb ] 2, (12) 9

10 so that on account of (85), we have ( b ± ) 2 = ˆb ˆb + ˆbˆb ± ˆb 2 ± ˆb 2. (13) Now employing (87), (92), an (94), we arrive at ( b + ) 2 = γ c ( N + ˆN a +2 ) r ˆN a a γc (14) an ( b ) 2 = γ ( ) c N + ˆN γc a 2 r ˆN a. (15) a We recall that the light generate by a two-level laser operating well above threshol is coherent, the quarature variance of which is given by [12] ( b ± ) 2 c = γ c N. (16) We calculate the quarature squeezing of the cavity (output) light relative to the quarature variance of the cavity (output) coherent light. We then efine the quarature squeezing of the cavity light by S = ( b ±) 2 c ( b ) 2 ( b ± ) 2. (17) c Hence employing (15), (16), an (58), one can put Eq. (17) in the form S = 2 η 1 η +2. (18) We observe that, unlike the mean photon number, the quarature squeezing oes not epen S η Fig. 1 A plot of Eq. (18) versusη. on the number of atoms. This implies that the quarature squeezing of the cavity light is inepenent of the number of photons. The plot in Fig. 1 shows that the maximum squeezing of the 1

11 cavity light is 5 % below the coherent-state level an this occurs when the three-level laser is operating below threshol atγ c = 4r a. On the other han, we efine the quarature squeezing of the output light by )2 S out = ( bout ± )2 c ( bout ( b out ± )2 c, (19) where ( b out ± )2 c is the quarature variance of the output coherent light. Since all calculations in this analysis are carrie out by putting the vacuum noise operators in normal orer, one can write ˆbout = ˆb. (11) We then easily see that an Now in view of Eqs. (19), (111), (112), an (17), we arrive at ( b out ± ) 2 c = ( b ± ) 2 c (111) ( b out )2 = ( b ) 2. (112) S out = S. (113) We observe that the quarature squeezing of the output light is equal to that of the cavity light. 4.2 Local quarature squeezing We finally seek to obtain the quarature squeezing of the cavity (output) light in a given frequency interval. To this en, we first etermine the spectrum of quarature fluctuations of the cavity light. We efine this spectrum for a single-moe light with central frequencyω by S ± (ω) = 1 ˆb± π Re (t),ˆb ± (t+τ) Upon integrating both sies of Eq. (114) overω, we get in which ss ei(ω ω )τ τ. (114) S ± (ω)ω = ( b ± ) 2, (115) ( b ± ) 2 = ˆb ± (t),ˆb ± (t) ss (116) is the quarature variance of the light moe at steay state. On the basis of the result given by Eq. (115), we assert thats ± (ω)ω is the steay-state quarature variance of the light moe in the interval betweenω anω +ω. In view of Eq. (85), we note that ˆb ± (t),ˆb ± (t+τ) = ˆb ± (t)ˆb ± (t+τ). (117) We now procee to etermine the two-time correlation function that appears in Eq. (117) for the cavity light. To this en, we realize that the solution of Eq. (55) can also be written as ˆb(t+τ) =ˆb(t)e τ/2 + g N e τ/2 τ e τ /2ˆm(t+τ )τ. (118) 11

12 On the other han, the solution of Eq. (71) is expressible as τ ( ) ˆm(t+τ) = ˆm(t)e µτ/2 +e µτ/2 e µτ /2 1 2 µˆm a(t+τ )+ ˆF m (t+τ ) τ. (119) Applying the large-time approximation scheme to Eq. (8), we get ˆm a (t+τ) = 2 µ ˆF a (t+τ), (12) so that on introucing this into Eq. (119), there follows τ ( ) ˆm(t+τ) = ˆm(t)e µτ/2 +e µτ/2 e µτ /2 ˆF a (t+τ )+ ˆF m (t+τ ) τ. (121) Now combination of Eqs. (118) an (121) yiels ˆb(t+τ)=ˆb(t)e τ/2 + 2gˆm(t) [ ] e µτ/2 e τ/2 + g τ e τ/2 N( µ) N τ e ( µ)τ /2 τ ( ) τ e µτ /2 ˆF a (t+τ )+ ˆF m (t+τ ). (122) Hence using this result, we reaily arrive at ˆb (t)ˆb(t+τ) = ˆb (t)ˆb(t) e τ/2 + 2g N( µ) ˆb (t)ˆm(t) [e µτ/2 e τ/2 ]. (123) Applying once more the large-time approximation, one gets from Eq. (55) With this substitute into Eq. (123), there emerges ( ˆb (t)ˆb(t+τ) = ˆb (t)ˆb(t) ˆm(t) = N 2g ˆb(t). (124) µ e µτ/2 µ µ e τ/2 ). (125) Following a similar proceure, one can also reaily establish that ( ˆb(t)ˆb (t+τ) = ˆb(t)ˆb (t) µ e µτ/2 µ ) µ e τ/2 (126) an ( ˆb(t)ˆb(t+τ) = ˆb 2 (t) µ e µτ/2 µ ) µ e τ/2. (127) Therefore, on account of Eqs. (125), (126), an (127), there follows ( ˆb ± (t),ˆb ± (t+τ) ss = ( b ± ) 2 µ e µτ/2 µ ) µ e τ/2. (128) Now on introucing Eq. (128) into Eq. (114) an carrying out the integration, we fin the spectrum of the minus quarature fluctuations for the cavity light to be [ ( ) S (ω) = ( b ) 2 µ/2π µ (ω ω ) 2 +[µ/2] 2 µ ( )] /2π µ (ω ω ) 2 +[/2] 2, (129) 12

13 in which( b ) 2 is given by (15). We realize that the variance of the minus quarature in the interval betweenω = λ anω = λ is expressible as λ ( b ) 2 ±λ = S (ω )ω, (13) in whichω = ω ω. Hence taking into account Eq. (129), we reaily get λ ( b ) 2 ±λ = z(λ)( b ) 2, (131) wherez(λ) is given by z(λ) = 2/π ( ) 2λ µ tan 1 µ 2µ/π µ tan 1 ( 2λ ). (132) On account of the result escribe by (131), the quarature variance of the output coherent light can be written in the same frequency interval as ( b ± ) 2 c±λ = z(λ)( b ±) 2 c. (133) We efine the quarature squeezing of the cavity light in theλ ± frequency interval by so that on account of (131), (133), an (17), there follows S ±λ = ( b ±) 2 c±λ ( b ) 2 ±λ ( b ± ) 2, (134) c±λ S ±λ = S. (135) This shows that the quarature squeezing of the cavity light in a given frequency interval is equal to that of the cavity light in the entire frequency interval. Finally, efining the quarature squeezing of the output light in the aforementione frequency interval by ± ) 2 c±λ ( bout ) 2 ±λ (136) an taking into account the fact that an we arrive at S±λ out = ( bout ( b out ± )2 c±λ ( b out ± ) 2 c±λ = z(λ)( bout ± ) c (137) ( b out ) 2 ±λ = z(λ)( bout ), (138) S out ±λ = Sout. (139) We see that the quarature squeezing of the output light in a certain frequency interval is the same as that of the output light in the entire frequency interval. 13

14 7. Conclusion We have consiere a three-level laser in which the three-level atoms available in a close cavity are pumpe from the bottom to the top level by means of electron bombarment. We have carrie out our analysis by putting the vacuum noise operators in normal orer an applying the large-time approximation scheme. The proceure of normal orering the noise operators reners the vacuum reservoir to be a noiseless physical entity. We maintain the stanpoint that the notion of a noiseless vacuum reservoir woul turn out to be compatible with observation. Base on the efinition given by Eq. (22), we infer that an atom in the top or mile level an insie a close cavity emits a photon ue to its interaction with the cavity light. We certainly ientify this process to be stimulate emission. We have seen that the light generate by the three-level laser operating uner the conition γ c > 1 4 r a is in a squeeze state, with the maximum quarature squeezing being 5% below the coherent-state level. This occurs when the three-level laser is operating below threshol at γ c = 4r a. We have also establishe that the quarature squeezing of the laser light has the same value in any frequency interval. On the basis of this result, we come to the conclusion that the quarature squeezing of the laser light is an intrinsic property of the iniviual photons. References [1] M.O. Scully, K. Wokiewicz, M.S. Zubairy, J. Bergou, N. Lu, J. Meyer ter Vehn, Phys. Rev. Lett. 6, 1832 (1988). [2] Fesseha Kassahun, Funamentals of Quantum Optics (Lulu Press Inc., North Carolina, 21). [3] N. Lu, F.X. Zhao, J. Bergou, Phys. Rev. A 39, 5189 (1989). [4] Eyob Alebachew an K. Fesseha, Opt. Commun. 265, 314 (26). [5] M.O. Scully an M.S. Zubairy, Opt. Commun. 66, 33 (1988). [6] N. Lu an S.Y. Zhu, Phys. Rev. A 41, 2865 (199). [7] N.A. Ansari, J. Gea-Banacloche, an M.S. Zubairy, Phys. Rev. A 41, 5179 (199). [8] C.A. Blockely an D.F. Walls, Phys. Rev. A 43, 549 (1991). [9] N.A. Ansari, Phys. Rev. A 48, 4686 (1993). [1] J. Anwar an M.S. Zubairy, Phys. Rev. A 49, 481 (1994). [11] K. Fesseha, Phys. Rev. A 63, (21). [12] Fesseha Kassahun, Opt. Commun. 284, 1357 (211). 14