Chapter 5Active Mode LockingFor simplicity, we assume, that the laser operates in the transverse fundamen-tal modes and, therefore, we only have to treat the longitudinal modes of thelaser similar to a simple plane parallel Fabry-Perot resonator (Figure: 5.1)We consider one polarization of the field only,however,as we will say laterfor some mode-locked laser polarization dynamics will become important.The task of mode-locking is to get as many of the longitudinal modeslasing in a phase synchronous fashion, such that the superposition of allmodes represents a pulse with a spatial extent much shorter than the cavity.The pulse will then propagate at the group velocity corresponding to thecenterfrequency of thepulse.XXX3Cavity Length, LFigure 5.1:Fabry-Perot resonator173
Chapter 5 Active Mode Locking For simplicity, we assume, that the laser operates in the transverse fundamental modes and, therefore, we only have to treat the longitudinal modes of the laser similar to a simple plane parallel Fabry-Perot resonator (Figure: 5.1). We consider one polarization of the field only, however, as we will say later for some mode-locked laser polarization dynamics will become important. The task of mode-locking is to get as many of the longitudinal modes lasing in a phase synchronous fashion, such that the superposition of all modes represents a pulse with a spatial extent much shorter than the cavity. The pulse will then propagate at the group velocity corresponding to the center frequency of the pulse. Figure 5.1: Fabry-Perot resonator 173
174CHAPTER5.ACTIVEMODELOCKING5.1The Master Equation of Mode Locking5Lets consider for the moment the cold cavity (i.e.there is only a simplelinear medium in the cavity no lasing). The most general solution for theintracavity field is a superpositon of left- and rightward running wavesZEnej(nt+Kn2)E(lef)(z,t) = Re ^ >(5.1)n=0andej(Ont-KnE(right)(z,t)= Re(5.2)The possible values for the wavenumbers are Kn = n/L, resulting from theboundary conditions on metallic mirrors or periodicity after one roundtrip inthe cavity. If the mirrors are perfectly reflecting, the leftward and rightwardmoving waves Eqs.(5.1) and (5.2) contain the same information and it issufficient to treat only one of them.Usually one of the cavity mirrors isnot perfectly refecting in order to couple out light, however, this can beconsidered a perturbation to the ideal mode structure.We consider the modes in Eq.(5.2) as a continuum and replace the sumby an integralE(K)e(2(K)t-K2)dKE(right)(z,t) (5.3)Be2TwithE(Km)= Em2L(5.4)Eq.(5.3) is similar to the pulse propagation discussed in chapter 2 and de-scribes the pulse propagation in the resonator. However, here it is ratheran initialvalueproblem,ratherthan aboundary valueproblem.Note,thewavenumbers of the modes are fixed, not the frequencies. To emphasize thiseven more, we introduce a new time variable T = t and a local time framet'=t-z/ug.o, instead of the propagation distance z, where ug.o is the groupvelocity at the central wave number Kno of the pulseaw(5.5)Ug.0ak
174 CHAPTER 5. ACTIVE MODE LOCKING 5.1 The Master Equation of Mode Locking Lets consider for the moment the cold cavity (i.e. there is only a simple linear medium in the cavity no lasing). The most general solution for the intracavity field is a superpositon of left- and rightward running waves E(lef t) (z, t) = Re (X∞ n=0 Eˆnej(Ωnt+Knz) ) , (5.1) and E(right) (z, t) = Re (X∞ n=0 Eˆnej(Ωnt−Knz) ) . (5.2) The possible values for the wavenumbers are Kn = nπ/L, resulting from the boundary conditions on metallic mirrors or periodicity after one roundtrip in the cavity. If the mirrors are perfectly reflecting, the leftward and rightward moving waves Eqs.(5.1) and (5.2) contain the same information and it is sufficient to treat only one of them. Usually one of the cavity mirrors is not perfectly reflecting in order to couple out light, however, this can be considered a perturbation to the ideal mode structure. We consider the modes in Eq.(5.2) as a continuum and replace the sum by an integral E(right) (z, t) = 1 2π Re ½Z ∞ K=0 Eˆ(K)ej(Ω(K)t−Kz) dK¾ (5.3) with Eˆ(Km) = Eˆm2L. (5.4) Eq.(5.3) is similar to the pulse propagation discussed in chapter 2 and describes the pulse propagation in the resonator. However, here it is rather an initial value problem, rather than a boundary value problem. Note, the wavenumbers of the modes are fixed, not the frequencies. To emphasize this even more, we introduce a new time variable T = t and a local time frame t 0 = t − z/υg,0, instead of the propagation distance z, where υg,0 is the group velocity at the central wave number Kn0 of the pulse υg,0 = ∂ω ∂k ¯ ¯ ¯ ¯ k=0 = µ∂k ∂ω¶−1 ¯ ¯ ¯ ¯ ¯ ω=0 . (5.5)
1755.1.THEMASTEREQUATIONOFMODELOCKINGFor introduction of a slowly varying envelope, we shift the frequency andwavenumber by the center frequency wo = no and center wave numberko = Knok=K-Kno,(5.6)(5.7)w(k) = 2(Kno +k)-no,E(k) = E(Kno +k),(5.8)The temporal evolution of the pulse is than determined by1.E(k)e;(a(k)t-k2) dkj(wot-koz)E(right)(z,t) =Re(5.9)2元Analogous to chapter 2, we define a slowly varying field envelope, that isalready normalized to the total power flow in the beamAeff1E(k)e;(k)t-k2) dk.A(z,t) =(5.10)2Z02元With the retarded time t' and time T, we obtain analogous to Eq. (2.184).AefiE(k)ei((k)-ug,ok)T+kug.ot dk.A(T,t) =(5.11)2Z02元whichcanbewrittenasQA(T,t)A(T,t),TR(5.12)i>OTOtI(GDD)n=22with thedispersion coefficients per resonator round-trip Tr=Ua.02Lan-1ug(k)Dn :(5.13)nlug,tOkn-The dispersion coefficients (5.13) look somewhat suspicious, however, it isnot difficult to show, that they are equivalent to derivatives of the roundtripphase r(2) = n(2)2L in the resonator at the center frequency1 angm(2)Dn=(5.14)n!02nl0=wo
5.1. THE MASTER EQUATION OF MODE LOCKING 175 For introduction of a slowly varying envelope, we shift the frequency and wavenumber by the center frequency ω0 = Ωn0 and center wave number k0 = Kn0 k = K − Kn0 , (5.6) ω(k) = Ω(Kn0 + k) − Ωn0 , (5.7) Eˆ(k) = Eˆ(Kn0 + k), (5.8) The temporal evolution of the pulse is than determined by E(right) (z, t) = 1 2π Re (Z ∞ −Kn0→−∞ Eˆ(k)ej(ω(k)t−kz) dk) ej(ω0t−k0z) . (5.9) Analogous to chapter 2, we define a slowly varying field envelope, that is already normalized to the total power flow in the beam A(z, t) = rAef f 2Z0 1 2π Z ∞ −∞ Eˆ(k)ej(ω(k)t−kz) dk. (5.10) With the retarded time t 0 and time T, we obtain analogous to Eq. (2.184). A(T,t0 ) = rAef f 2Z0 1 2π Z ∞ −∞ Eˆ(k)ej((ω(k)−υg,0k)T +kυg,0t0 dk. (5.11) which can be written as TR ∂A(T,t0 ) ∂T ¯ ¯ ¯ ¯ (GDD) = j X∞ n=2 Dn µ −j ∂n ∂t0 ¶n A(T, t0 ), (5.12) with the dispersion coefficients per resonator round-trip TR = 2L υg,0 Dn = 2L n!υn+1 g,0 ∂n−1υg(k) ∂kn−1 ¯ ¯ ¯ ¯ ¯ k=0 . (5.13) The dispersion coefficients (5.13) look somewhat suspicious, however, it is not difficult to show, that they are equivalent to derivatives of the roundtrip phase φR(Ω) = Ω c n(Ω)2L in the resonator at the center frequency Dn = − 1 n! ∂nφ(n) R (Ω) ∂Ωn ¯ ¯ ¯ ¯ ¯ Ω=ω0 , (5.14)
176CHAPTER5.ACTIVEMODELOCKINGSofar, only the lossless resonator is treated. The gain and loss can be mod-elled by adding a term likeA(T,t)= -IA(T,t)(5.15)TRaTI(loss)where l is the amplitude loss per round-trip. In an analogous manner we canwrite for thegainA(T,t)A(T,t),(5.16)TR--1OTat/asnwhere g(T) is the gain and and D, is the curvature of the gain at the maxi-mumoftheLorentzian lineshape.g(T)Dg=(5.17)22Dg is the gain dispersion. g(T) is an average gain, which can be computedThe dis-from the rate equation valid for each unit cell in the resonator.tributed gain obeys the equation[A(z, t)2Og(z,t)g - 9o(5.18)otELTLAeff,T the upper state lifetimewhere E is the saturation energy E=and o the gain cross section. For typical solid-state lasers, the intracavitypulse energy is much smaller than the saturation energy. Therefore, the gainchanges within one roundtrip are small. Furthermore, we assume that thegain saturates spatially homogeneous, g(z,t') = g(t'). Then, the equation forthe average gain g(T) can be found by averageing (5.18) over one round-tripandweobtainW(T)0g(T)g - 90(5.19)'ELTRaTTLwhere W(T) is the intracavity pulse energy at time t = TTR/2IA(T,t)’dt ~A(T,t)Pdt(5.20)W(T) 2
176 CHAPTER 5. ACTIVE MODE LOCKING Sofar, only the lossless resonator is treated. The gain and loss can be modelled by adding a term like TR ∂A(T, t0 ) ∂T ¯ ¯ ¯ ¯ (loss) = −lA(T, t0 ) (5.15) where l is the amplitude loss per round-trip. In an analogous manner we can write for the gain TR ∂A(T,t0 ) ∂T ¯ ¯ ¯ ¯ (gain) = µ g(T) + Dg ∂2 ∂t02 ¶ A(T,t0 ), (5.16) where g(T) is the gain and and Dg is the curvature of the gain at the maximum of the Lorentzian lineshape. Dg = g(T) Ω2 g (5.17) Dg is the gain dispersion. g(T) is an average gain, which can be computed from the rate equation valid for each unit cell in the resonator. The distributed gain obeys the equation ∂g(z, t) ∂t = −g − g0 τ L − g |A(z, t)| 2 EL , (5.18) where EL is the saturation energy EL = hνL 2∗σL Aef f , τ L the upper state lifetime and σL the gain cross section. For typical solid-state lasers, the intracavity pulse energy is much smaller than the saturation energy. Therefore, the gain changes within one roundtrip are small. Furthermore, we assume that the gain saturates spatially homogeneous, g(z, t0 ) = g(t 0 ). Then, the equation for the average gain g(T) can be found by averageing (5.18) over one round-trip and we obtain ∂g(T) ∂T = −g − g0 τ L − g W(T) ELTR , (5.19) where W(T) is the intracavity pulse energy at time t = T W(T) = Z TR/2 t0=−TR/2 |A(T,t0 )| 2 dt0 ≈ Z ∞ −∞ |A(T,t0 )| 2 dt0 . (5.20)
1775.2.ACTIVEMODELOCKINGBYLOSSMODULATIONHighOutputCouplerReflectorAOMGainFigure 5.2:Actively modelocked laser with an amplitude modulator(Acousto-Optic-Modulator).Taking all effects into account, the linear ones: loss, dispersion, gain andgain dispersion, as well as the nonlinear ones like saturable absorption andself-phase modulation, weend up with the master equation of modelockingaA(T,t)ar-IA(T,t)+ZA(T,t)TROTotn=2102(5.21)A(T,t)+ g(T22 0t/2q(T,t)A(T,t) - js|A(T,t)?A(T,t)To keep notation simple, we replace t' by t again. This equation was firstderived by Haus [4] under the assumption of small changes in pulse shapeper round-trip and per element passed within one round-trip.5.2Active Mode Locking by Loss Modula-tionActive mode locking was first investigated in 1970 by Kuizenga and Siegmanusing a gaussian pulse analyses, which we want to delegate to the exercises[3]. Later in 1975 Haus [4] introduced the master equation approach (5.21),We follow the approach of Haus, because it also shows the stability of thesolution.We introduce a loss modulator into the cavity,for example an acousto-optic modulator, which periodically varias the intracavity loss according to
5.2. ACTIVE MODE LOCKING BY LOSS MODULATION 177 Figure 5.2: Actively modelocked laser with an amplitude modulator (Acousto-Optic-Modulator). Taking all effects into account, the linear ones: loss, dispersion, gain and gain dispersion, as well as the nonlinear ones like saturable absorption and self-phase modulation, we end up with the master equation of modelocking TR ∂A(T,t0 ) ∂T = −lA(T,t0 ) + j X∞ n=2 Dn µ j ∂n ∂t ¶n A(T, t0 ) + g(T) µ 1 + 1 Ω2 g ∂2 ∂t02 ¶ A(T,t0 ) (5.21) − q(T, t0 )A(T,t0 ) − jδ|A(T, t0 )| 2 A(T,t0 ). To keep notation simple, we replace t 0 by t again. This equation was first derived by Haus [4] under the assumption of small changes in pulse shape per round-trip and per element passed within one round-trip. 5.2 Active Mode Locking by Loss Modulation Active mode locking was first investigated in 1970 by Kuizenga and Siegman using a gaussian pulse analyses, which we want to delegate to the exercises [3]. Later in 1975 Haus [4] introduced the master equation approach (5.21). We follow the approach of Haus, because it also shows the stability of the solution. We introduce a loss modulator into the cavity, for example an acoustooptic modulator, which periodically varias the intracavity loss according to