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

Theoretical investigation of injection-locked high modulation bandwidth quantum cascade lasers

Open Access Open Access

Abstract

In this study, we report for the first time to our knowledge theoretical investigation of modulation responses of injection-locked mid-infrared quantum cascade lasers (QCLs) at wavelengths of 4.6 μm and 9 μm, respectively. It is shown through a three-level rate equations model that the direct intensity modulation of QCLs gives the maximum modulation bandwidths of ~7 GHz at 4.6 μm and ~20 GHz at 9 μm. By applying the injection locking scheme, we find that the modulation bandwidths of up to ~30 GHz and ~70 GHz can be achieved for QCLs at 4.6 μm and 9 μm, respectively, with an injection ratio of 5 dB. The result also shows that an ultrawide modulation bandwidth of more than 200 GHz is possible with a 10 dB injection ratio for QCLs at 9 μm. An important characteristic of injection-locked QCLs is the nonexistence of unstable locking region in the locking map, in contrast to their diode laser counterparts. We attribute this to the ultra-short upper laser state lifetimes of QCLs.

©2012 Optical Society of America

1. Introduction

Semiconductor laser sources with high modulation bandwidths are always desirable for high speed data transmission systems. However, the modulation bandwidth of direct modulated semiconductor laser is largely limited by relaxation resonance frequency determined by the interactions of the carriers and photons in the laser cavity. Optical injection locking scheme has been shown theoretically [1] and experimentally [2] as an effective approach to enhance the relaxation resonance frequency of semiconductor lasers, making it a useful method to increase the modulation bandwidth of lasers. With this technique, vertical-cavity surface-emitting lasers (VCSELs) with a 80 GHz intrinsic 3-dB bandwidth [3], quantum-dot (QD) DFB lasers with a 16.3 GHz bandwidth [4], and distributed reflector (DR) lasers with wirelike active regions [5] with a 15 GHz bandwidth have been demonstrated recently.

Quantum cascade lasers (QCLs) are the most promising semiconductor-based mid-infrared and terahertz sources. As they have intrinsic ultrafast intersubband transitions and versatile-designed wavelengths, they are very attractive for high-speed free-space optical (FSO) communication systems [6, 7] (two atmospheric transmission windows are in the regimes (3-5 μm) and (8-14 μm), respectively.) and Light Detection and Ranging (LIDAR) applications. Recent advancement of high-performance continuous-wave room-temperature operated QCLs [8, 9] has further brought interests of applying QCLs for FSO applications [10]. With direct modulation of 8.1 μm QCLs, enhanced link stability has been demonstrated in FSO communication system in fog weather due to the less Mie scattering at relatively longer wavelengths, showing that QCL is an ideal source for FSO communication in atmospheric conditions with restricted visibility [10]. Earlier theoretical studies show that the modulation bandwidth of QCLs can reach 100 GHz or even Terahertz regimes [1113]. However, this is misleading as the value of photon lifetime was inaccurately estimated for practical QCLs. Also in [14], the rate equations for QCLs were not correctly established, as the cascade characteristics of QCLs were not properly considered. Direct modulation of QCLs actually leads to a much lower modulation bandwidth ~10 GHz for some mid-infrared QCLs [15] and ~13 GHz for typical terahertz QCLs [16]. Further increase of the modulation bandwidth of QCLs through other approaches is highly desired.

In the paper, we theoretically and systematically investigate the injection-locking of QCLs, for the first time to our knowledge, at emitting wavelengths of 4.6 μm [17] and 9 μm [18] corresponding to the two atmospheric transmission windows, respectively. We first examine the direct intensity modulation characteristics of QCLs, using a three-level rate equations model. The results show that no resonant frequency appears in the frequency modulation response, and the maximum modulation bandwidths of round 7 GHz for QCLs at 4.6 μm and 20 GHz for QCLs at 9 μm are obtained. The theoretical analysis is in good agreement with the experimental observations [15]. In addition, we apply the optical injection locking scheme to increase the frequency modulation bandwidth of QCL. Rate equation analysis shows that the obtained modulation bandwidth of injection locked QCLs is about three times higher than that of the directly modulated QCLs for both lasers at 4.6 μm and 9 μm under a 5 dB injection ratio, with modulation bandwidths up to ~30 GHz and ~70 GHz, respectively. With a 10 dB injection ratio, a modulation bandwidth of over 200 GHz can be achieved with the injection locking scheme. A unique feature of injection locked QCLs is that there is no unstable locking range, as opposed to other semiconductor lasers [19]. We attribute this effect to the ultra-short lifetime of the upper laser state.

2. Theory

2.1 Direct intensity modulation of QCLs

In this study, we use a three-level rate equation model to describe the dynamic behavior of injection-locked QCLs. It is noted that neglecting the carrier number in the lower laser state is invalid in QCLs, especially for room temperature operation, also to accurately describe the dynamic behavior of QCLs, a self-consistent scattering model [20] is required, which however is too complicated for the present study. We denote the instantaneous carrier numbers in the lower and upper laser state by N2 and N3, respectively, and the photon number by P. Here the cavity is assumed to have only one longitudinal mode. We noticed that analysis of direct intensity modulation of QCLs with rate equation model has also been investigated in [14]; however, the gain coefficients for carrier rate equations are miscalculated due to the neglect of the cascade characteristics of QCLs. Therefore, the obtained modulation bandwidth in [14] is much larger than that demonstrated in experiments [15]. This conclusion is also supported by the analysis in Ref [21]. The resulted rate equations for QCLs should be:

dN3dt=JN3τ3GNp(N3N2)P
dN2dt=N3τ3+GNP(N3N2)PN2τ2
dPdt=G(N3N2)PPτP
where J denotes the current injected into the active region divided by electronic charge e, G is the optical gain coefficient of the entire active region, Np is the number of stages, τ2 and τ3 represent the lifetime of lower and upper laser levels, respectively, and τP is the lifetime of photon in the cavity, expressed as τP = νgwm) where αw, αm, and νg are the waveguide loss, mirror loss, and group velocity, respectively. The mirror loss can be calculated using αm = -ln(R1R2)/(2L), where R1 and R2 are the power reflectivity of the facets 1 and 2, respectively. The group velocity νg is given by νg = c/neff, where c and neff are the speed of light in vacuum and modal effective refractive index, respectively. In our analysis, G/NP instead of G as in [14] is correctly used in Eqs. (1) and (2) to represent the optical gain coefficient of one single active region period in QCLs. Denote G0 = G/NP, we have G0 = Γνgσ32/V, where Γ is the optical mode confinement factor for a single period, and V stands for the volume for a period, given by WLLP, where W and L are the width and length of the active region, and LP denotes the period length, respectively. The stimulated emission cross section σ32 is,
σ32=4πe2z322ε0nneffλ0(2γ32)
where e is the electronic charge, z32 is the dipole matrix element between the upper and the lower lasing levels, ε0 is the vacuum dielectric constant, λ0 is the free-space emission wavelength, and (2γ32) stands for the full width at half maximum (FWHM) of the optical transition spectrum.

2.1.1 Steady state analysis

We set the time derivatives of Eqs. (1)-(3) to zero, and denote the steady state values of J, N2, N3, and P as J0, N20, N30 and P0, then the steady state rate equations can be expressed as:

J0N30τ3G0(N30N20)P0=0
N30τ3+G0(N30N20)P0N20τ2=0
G(N30N20)P0P0τP=0

The steady state population of lower lasing level, N20 is given as

N20=J0τ2
with Eqs. (4) and (5), P0 can be obtained in the following equation,
J0(1η)N30N20τ3G0(N30N20)P0=0
where η=τ2/ τ3. Also, we note that at threshold,
(N30N20)=1NpG0τp
where Jth is the threshold current (s−1), given by Jth = 1/3τP (1-η). Substituting Eq. (9) into Eq. (8), one gets the steady state photon number P0,
P0=1G0τ3(J0Jth1)
with Eqs. (7) and (9), the upper laser level population N30 can be written as,

N30=1NPG0τP+J0τ2

2.1.2 Small-signal modulation analysis

The small-signal analysis is necessary in order to investigate the modulation response of QCLs. Assuming small variations ∆N2, ∆N3, ∆P and ∆J around the steady state values, and substituting N2=N20+∆N2, N3=N30+∆N3, P=P0=+∆P and J=J0+∆J into Eqs. (1)(3), we get the rate equations for the small deviations as follows:

dΔN3dt=ΔJΔN3τ3G0(ΔN3ΔN2)P0G0(N30N20)ΔP
dΔN2dt=ΔN3τ3+G0(ΔN3ΔN2)P0+G0(N30N20)ΔPΔN2τ2
dΔPdt=G(ΔN3ΔN2)P0+G(N30N20)ΔPΔPτP

Taking the Laplace transform of Eqs. (12)(14), and placing them into a matrix form, we have the following transformed matrix,

[s+f11f12f13f21s+f22f23f31f32s+f33][ΔN3ΔN2ΔP]=[ΔJ00]
where the matrix terms are expressed as

f11=1τ3+G0P0f21=(1τ3+G0P0)f31=GP0f12=G0P0f22=1τ2+G0P0f32=GP0f13=1τPNPf23=1τPNPf33=0

The normalized modulation response is

H(s)=ΔPΔJ(1η)τPNP=sA+1Bs3+Cs2+Ds+1
Where

A=(1η1)1τ3B=τPτ2NPGP0C=B(2G0P0+1τ3(1+1η))D=B(1τ2G0P0+1τ3τ2+2GP01τPNP)

2.2 Injection locking of quantum cascade lasers (QCLs)

2.2.1 Steady state analysis

The injection locking scheme involves two lasers: the slave laser (SL) is optically locked by the master laser (ML) [2]. The differential equation describing electric field inside the slave laser under injection locking was first proposed by Lang and Kobayashi, given as

ddtESL(t){jω(N)+12[NPG0(N)1τp]}ESL(t)=fdEML(t)
where fd is the coupling rate between the master laser and the slave laser, approximately expressed as: fd = νg(1-R)/(2LR1/2) where R stands for the power reflectivity of the injected cavity facet. We denote the electromagnetic fields in the slave laser and the master laser as,
ESL(t)=E0(t)ej(ω0t+ϕ0(t))
EML(t)=E1(t)ej(ω1t+ϕ1(t))
where E1(t) is taken as a constant, ϕ0(t) and ϕ1(t) are the phases of the two electromagnetic fields, (ϕ1(t) is usually set as zero for computing convenience); ω0 and ω1 are the angular frequencies of the slave laser and the master laser, respectively. Substituting Eq. (20) into (19), the differential equation can be split into field magnitude and phase rate equations.

ddtE0(t)=12NPG0ΔN(t)E0(t)+fdE1cosϕ(t)
ddtϕ(t)=12αNPG0ΔN(t)fdE1E0sinϕ(t)Δωinj

In the above equations, ∆N is the carrier number change due to light injection from the master laser, given by N-Nth. Different from rate equations for diode lasers, the carrier number in the above equations is the carrier number difference between the upper and the lower laser levels, i.e. N = N3-N2. From Eqs. (12) and (13), and taking into account the carrier number in the lower laser state, the differential equations governing the carrier dynamics are given as,

dN(t)dt=J(N+N2)2τ32E02G0N+N2τ2
dN2(t)dt=(N+N2)1τ3+E02G0NN2τ2
where E0(t), ϕ0(t), N(t), and N2(t) represent the slave laser’s field magnitude, phase, carrier number difference between the transition states, and carrier number in the lower laser state, respectively. The field magnitude E0(t) has been normalized, so that |E0(t)|2=P(t), where P(t) is regarded as the total photon number for a single longitudinal mode inside the cavity. ϕ(t) is defined as the phase difference between the master and the slave lasers, i.e. ϕ(t)=ϕSL(t)−ϕML(t), G0, NP, α, J, τ2, and τ3 are the slave laser’s gain coefficient for one stage, the number of period, linewidth enhancement factor, the injection current, the lower laser state lifetime, and the upper laser state lifetime, respectively. While fd, E1 and ∆ωinj represent the coupling rate, injected electrical field magnitude, and frequency detuning, respectively. The latter is expressed as ∆ωinj = ω1-ω0. For free-running laser, N = Nth, N2 = N2th, E = Efr, with Eq. (24), the carrier number in the lower laser state for free-running slave laser can be obtained,
N2th=(Efr2γP+Nthτ3)τ2τ3τ3τ2
where we define γP as,

γP=G0Nth

Remembering that cavity photon decay rate γP0 is given by GNth, so that γP = γP0/NP. With Eq. (23) and the expression of γP, the threshold current, defined as I/e (where I is current, and e is the electronic charge), and the free-running field magnitude of slave laser are found to be,

Jth=Nth1τ3τ2
Efr2=JNth1τ3τ2γPτ3τ3τ2

Under the injection-locking region, the carrier number difference N in the slave laser should decrease, due to the enhanced stimulated mission in the gain medium. The steady state values for E, ∆N, ϕ, and N2 under injection locking are denoted as E0, ∆N0, ϕ0 and N20 respectively. The corresponding expressions of ∆N0, ϕ0 and N20 are given below,

ΔN0=2fdE1NpG0E0cosϕ0
ϕ0=sin1(ΔωinjE0fdE11+α2)tan1α
N20=[E02(γP+G0ΔN0)+ΔN0τ3+Nthτ3]τ2τ3τ3τ2

Applying Eqs. (29)-(31) into Eqs. (23), and setting the time derivative to zero, one gets a cubic equation for the steady state value E0. Defining two new parameters,

T1=2τ2τ3τ3τ2
T2=2τ3+T1τ3
the equation is shown as,

E03(2G0NPγP+G0NPγPT1)E022G0fdE1cosϕ0(2+T1)E0(G0JNPNthG0NPT2)2fdE1cosϕ0T2=0

The steady state value of E0 can be solved by various numerical methods, given that the phase difference ϕ0 varies approximately from -π/2 to cot−1α in the injection locking range, and ∆N0 should be negative.

2.2.2 Small-signal modulation analysis

Assuming a small variances ∆E, ∆ϕ, ∆N and ∆N2 around the corresponding steady state value, and substituting E=E0+∆E, ϕ=ϕ0+ϕ, N=N0+∆N, N2=N20+∆N2 into Eqs. (21)(24), one gets the rate equations for the small signal modulation:

ddtΔE=12NPG0ΔN0ΔEfdE1sinϕ0Δϕ+12NpG0E0ΔN
ddtΔϕ=fdE1E02sinϕ0ΔEfdE1E0cosϕ0Δϕ+α2NPG0ΔN
ddtΔN=ΔJ(4E0γP+4E0G0ΔN0)ΔE(2τ3+2E02G0)ΔN+(1τ22τ3)ΔN2
ddtΔN2=(2E0γP+2G0ΔN0E0)ΔE+(1τ3+E02G0)ΔN+(1τ31τ2)ΔN2

Converting the rate equations into the matrix form and taking Laplace transform, we have the following matrix:

[Q11+sQ12Q13Q14Q21Q22+sQ23Q24Q31Q32Q33+sQ34Q41Q42Q43Q44+s][ΔEΔϕΔNΔN2]=[00ΔJ0]

In which the matrix elements are

Q11=12G0ΔN0NPQ12=fdE1sinϕ0Q13=12G0E0NPQ14=0Q21=fdE1E02sinϕ0Q22=fdE1E0cosϕ0Q23=α2G0NPQ24=0Q31=4E0γP+4E0G0ΔN0Q32=0Q33=2τ3+2E02G0Q34=2τ31τ2Q41=(2E0γP+2G0ΔN0E0)Q42=0Q43=(1τ3+E02G0)Q44=(1τ21τ3)

The direct modulation frequency response is given by,

H(s)=Gs2+Es+Fs4+As3+Bs2+Cs+D
where

A=Q11+Q22+Q33+Q44B=Q33Q44+Q11Q33+Q11Q44+Q22Q33+Q22Q44+Q11Q22Q21Q12Q31Q13Q34Q43C=Q11Q33Q44+Q22Q33Q44+Q33Q11Q22+Q44Q11Q22Q33Q21Q12Q44Q21Q12Q13Q31Q22Q13Q31Q44+Q12Q31Q23+Q13Q34Q41Q22Q43Q34Q11Q34Q43D=Q11Q22Q33Q44Q12Q21Q33Q44+Q12Q21Q43Q34Q13Q31Q22Q44+Q12Q31Q23Q44Q12Q23Q34Q41+Q13Q22Q34Q41Q11Q22Q34Q43E=(Q12Q23Q13Q22Q13Q44)/Q13F=(Q12Q23Q44Q13Q22Q44)/Q13G=Q13

3. Results and discussions

3.1 Direct modulations

First, we examine the direct modulation bandwidths of QCLs at different wavelengths in the mid-IR regime. Table 1 shows the parameters used in the simulations, of the state-of-the-art mid-infrared QCLs emitting at around 4.6 μm [17] and 9 μm [18], which correspond to the two atmospheric windows in the mid-infrared spectrum range. A simple way to analyze the modulation behavior of QCLs can be done by analyzing the zeros and poles of the normalized frequency response H(s). Tables 2 and 3 list the corresponding zeros and poles for the QCLs at 4.6 μm and 9 μm under different injection currents. Three poles can be obtained from the denominator of the frequency response, which are expressed as p1, p2 and p3 respectively. One zero will be obtained from the numerator. Also shown is the calculated 3-dB bandwidth of the modulation response. We also notice that the phenomena that no resonance frequency appears and the calculated direct modulation bandwidth are in good agreement with the experimental observations [15], if the experimental parameters [15] are used in the calculations.

Tables Icon

Table 1. Characteristic parameters of the state-of-the-art QCLs at 4.6 μm [17] and 9 μm [18] for direct modulation

Tables Icon

Table 2. Poles, zeros, and f3dB for QCLs at 4.6 μm

Tables Icon

Table 3. Poles, zeros, and f3dB for QCLs at 9 μm

Figures 1(a) and (b) show the normalized frequency responses of QCLs at 4.6 μm and 9 μm, for I = 1.5Ith, 2Ith, 3Ith, and 4Ith, respectively, corresponding to the injection currents from just above threshold to roll-over, respectively. Obtaining the poles and zeros of the frequency response, we can transfer it into the following form,

 figure: Fig. 1

Fig. 1 Normalized modulation response of mid-infrared QCLs. (a) QCLs emitting at 4.6 μm with parameters shown in Table.1 (b) QCLs emitting at 9 μm with parameters shown in Table.1. The 3-dB bandwidth is indicated by the black dashed line.

Download Full Size | PDF

H(s)=(s|z1|+1)(s|p1|+1)(s|p2|+1)(s|p3|+1)

This is the normal form of the Bode plot. Because the first pole is three to four magnitudes smaller than the other poles and zeros, it will mainly determine the modulation bandwidth, so H(s)≈1/(s/|p1|+1) in the Bode plot. For cubic equations, we have p1p2+p2p3+p1p3=D/B, and p1p2p3=−1/B, leading to 1/p1+1/p2+1/p3=−1/D. For p1 is much smaller than p2 and p3, we get |1/p1|≈|D|. Substituting |1/p1| with |D| in the expression of H(s), and keeping in mind that D is positive in the whole locking region; one can approximate the frequency response as,

H(s)=1sD+1
where D is given by (τP + 2τ2) + τP/(G0P0τ3). The definition of 3-dB bandwidth is defined as |H(s)| at half of its zero value, i.e. |H(s)|f3dB| = 1/2 in the case of normalized frequency response. An approximated expression for 3-dB bandwidth is given as,

f3dB32πD=32π[(τP+2τ2)+τpG0P0τ3]

Combining the D and f3dB, we can easily see that a larger value of τp leads to a decreased modulation bandwidth. This can also be seen from the simulated frequency responses of QCLs at 4.6 μm and 9 μm. Due to larger optical loss, e.g. free carrier absorption (proportional to λ2) and intersubband absorption in the waveguide, the photon lifetime is much shorter for longer emission wavelengths, leading to an increased modulation bandwidth. However, a larger τp means a decreased optical loss, leading to a higher optical output. Thus, the tradeoff between the modulation bandwidth and the output optical power has to be taken into account when designing QCLs for high speed modulation applications. Another feature of QCLs is the absence of resonance peak as normally shown in conventional diode lasers. The physics behind it lies in the ultrafast upper state lifetime compared with the photon lifetime, making QCLs an overdamped system, showing no resonance peak in the frequency modulation response.

3.2 Injection locking modulations

To further enhance the modulation bandwidth of QCLs in the mid-IR regime, we can employ injection locking scheme for QCLs, the theoretical analysis of which is shown in section 2.2. All the additional parameters used in injection locking calculation are listed in Table 4 , for QCLs emitting at 4.6 μm and 9 μm, respectively. The rest of the parameters which are common to direct modulated QCLs are shown in Tables 1. Though the linewidth enhance factor of QCLs is expected to be zero [22], the experiments showed different values (0.02 ± 0.2 in [23],-0.44 to 2.29 in [24], and −1.8 to −1.7 in [25]). Without loss of generality, we set the linewidth enhancement factor as unity in the calculations.

Tables Icon

Table 4. Characteristic parameters of the state-of-the-art QCL at 4.6 μm [17] and 9 μm [18] for injection locking modulation

Under injection locking scheme, four equations are involved in the calculations. This results in a four-order denominator and a two-order numerator in the frequency response. The typical locking maps for both of QCLs at 4.6 μm and 9μm are illustrated in Fig. 2 .(a) and (b), respectively, with J = 4Jth. In our simulations, the boundaries of phase in the injection locking range are approximately cot−1α to -π/2, from negative to positive detuning edge. This is derived from rate equations for the field magnitude E and phase difference ϕ, where the noise terms and spontaneous terms are neglected. The calculated locking maps of QCLs are similar to those of diode lasers, where the locking range increases linearly with the increase of the amplitude of the injected optical field, as demonstrated in Ref [26].

 figure: Fig. 2

Fig. 2 Locking maps for QCLs emitting at (a) 4.6 μm and (b) 9 μm. The parameters used in the simulations are shown in Tables 1 and 4.

Download Full Size | PDF

To examine the frequency response of the injection locked QCLs, we calculate the response curves of QCLs at 9 μm with an injection ratio of 5 dB across the locking range and a frequency spacing of 1 GHz. The corresponding waterfall plot is shown in Fig. 3 . Table 5 lists the accompanying poles, zeros and 3 dB bandwidth of the responses.

 figure: Fig. 3

Fig. 3 Normalized frequency response curves versus frequency detuning at a fixed injection ratio R = 5 dB, for the 9 μm QCL.

Download Full Size | PDF

Tables Icon

Table 5. Poles, Zeros and 3-dB bandwidth f3dB for the QCL in Fig. 3.

The frequency response exhibits resonance-like behaviors close to both edges of the locking range (e.g. the green and the black curves, respectively, in Fig. 3). However, the reasons behind these two are different. A further explanation can be sought by examining the poles of the frequency response. As illustrated in Table 5, when the frequency comes close to the positive frequency detuning edge (black curve in Fig. 3), two complex conjugate poles appear. In conventional semiconductor lasers, the imaginary part of the complex conjugate poles gives the resonance frequency of the response, while the real part defines the damping term, so that a peak appears in the frequency response. This is also the case for QCLs. At the negative detuning edge (black curve in Fig. 3), because all the poles are real (and negative), the peak in the frequency response has to be analyzed by the Bode Plots. According to the Bode Plots theory, zeros will increase the magnitude of the frequency response from its critical frequency and beyond 10 dB per decade, while poles decrease the magnitude at the same rate. For QCLs, the first zero z1 is smaller than any of the poles, making the peak in the frequency response apparent, even though all the poles are real.

For conventional semiconductor lasers, increasing slave laser’s bias current is an effective way to enhance the bandwidth of the modulation system [12]. Here we calculate the frequency response at a fixed injection ratio of 5 dB, with the injection current varying 1.5x, 2x, 3x, and 4xJth, respectively, and a frequency detuning 0.1 GHz away from the negative detuning edge, as shown in Fig. 4 . The inset is the corresponding pole/zero diagrams. Tables 6 and 7 list the associated poles and zeros of the modulation responses. The increased bandwidth can be attributed to the increased values of poles and zeros of the frequency response, and the much smaller z1 than any of the poles. Similar effects can also be seen by increasing the injection ratio, for which a much broader modulation bandwidth can be achieved. Our calculation shows that for injection ratio R = 10 dB and frequency detuning of 0.1 GHz away from the edge, over 200 GHz modulation bandwidth can be obtained. However, there is a tradeoff between the injection ratio and the slave laser’s injection current, even though both of which can effectively increase the modulation bandwidth. This is because the increased injection current leads to an increased slave laser output power, so that in order to achieve a high injection ratio, a master laser with a higher output power is required, which in some cases is limited in real applications. We noticed that by direct modulating the injection current using an RF source, Pierre et al demonstrated that the cavity resonance frequency of terahertz QCLs can be injection-locked [27]. However, their scheme is different to ours in that the injection signal in our case is the optical signal from the master laser.

 figure: Fig. 4

Fig. 4 Normalized modulation response of mid-infrared QCLs under injection locking scheme. (a) QCLs emitting at 4.6 μm. (b) QCLs emitting at 9 μm. The inset show the corresponding poles (x) and zeros (), with the arrow indicating direction of increasing the injection current of the slave laser. The 3 dB bandwidth is indicated by the dashed line.

Download Full Size | PDF

Tables Icon

Table 6. Poles, zeros, and f3dB for the QCL in Fig. 4 (a)

Tables Icon

Table 7. Poles, zeros, and f3dB for Fig. 4 (b)

It is also noted that there is no unstable locking range in QCLs based on the above rate equation analysis, compared to other types of semiconductor lasers. Mathematically, the stable locking condition is satisfied when the eigenvalues of the frequency response function’s denominators are located at the left half of the complex s-plane. In the entire locking region of QCLs, we find that the roots of the denominator are always negative, making the whole locking region a stable locking system. One of the experimental evidence is that the injected light induced pulsation effect [19], which is a unique sign of the unstable locking condition in diode lasers [28] caused by the reduced damping and the strong correlation between the phase and intensity of photons, has not been observed in QCLs. The reasons are as follows: first, there is no relaxation oscillation in QCLs which has been experimentally verified, owing to the picosecond carrier lifetime of the laser states [15]. It is equivalent to say that the damping for the relaxation oscillation is so high such that the carrier numbers and photons come to a steady state in a much shorter time than those in conventional diode lasers. Second, using the phasor diagram and small signal analysis (here we express the small change of N, P, and ϕ as δN, δP and δϕ), Henry pointed out [28] that the increase of δP causes an increased δϕ in relaxation oscillations. On the other hand, at the negative detuning edge, the increase of δϕ causes an enhanced mismatch of E0 and E1, leading to a larger cavity intensity change, i.e. larger δP. This forms a positive feedback between δϕ and δP, which is enhanced as the relaxation oscillations increase. However, there is no relaxation oscillation in QCLs, thus the interaction between δϕ and δP is reduced. Therefore the large damping effect and the reduced interactions between ϕ and P, make the unstable locking range in QCLs disappear, as shown in Fig. 2.

4. Conclusions

Theoretical investigation of injection-locked QCLs at wavelengths of 4.6 μm and 9 μm are carried out in the paper. By using a three-level rate equations model, we find that the maximum modulation bandwidths of round 7 GHz for QCLs at 4.6 μm and 20 GHz for QCLs at 9 μm were obtained. It is shown that by applying the injection locking scheme, the modulation bandwidths can be increased by around three times under a 5 dB injection ratio, compared to the direct modulation scheme. The frequency modulation responses were analyzed using the Bode diagram, based on which it shows that increasing the injection ratio and the injected current of the slave laser are effective ways to enhance the modulation bandwidth of the master QCLs. With a 10 dB injection ratio, more than 200 GHz modulation bandwidth can be obtained for QCLs at 9 μm. Unlike conventional semiconductor lasers, no unstable locking ranges appear in the locking map, which we attribute to the ultra-short lifetime of the upper laser state of QCLs, due to the nature of intersubband transitions.

References and links

1. A. Murakami, K. Kawashima, and K. Atsuki, “Cavity resonance shift and bandwidth enhancement in semiconductor lasers with strong light injection,” IEEE J. Quantum Electron. 39(10), 1196–1204 (2003). [CrossRef]  

2. E. K. Lau, H.-K. Sung, and M. C. Wu, “Frequency response enhancement of optical injection-locked lasers,” IEEE J. Quantum Electron. 44(1), 90–99 (2008). [CrossRef]  

3. E. K. Lau, X. Zhao, H.-K. Sung, D. Parekh, C. Chang-Hasnain, and M. C. Wu, “Strong optical injection-locked semiconductor lasers demonstrating > 100-GHz resonance frequencies and 80-GHz intrinsic bandwidths,” Opt. Express 16(9), 6609–6618 (2008). [CrossRef]   [PubMed]  

4. N. B. Terry, N. A. Naderi, M. Pochet, A. J. Moscho, L. F. Lester, and V. Kovanis, “Bandwidth enhancement of injection-locked 1.3 μm quantum-dot DFB laser,” Electron. Lett. 44(15), 904–905 (2008). [CrossRef]  

5. S. H. Lee, D. Parekh, T. Shindo, W. J. Yang, P. Guo, D. Takahashi, N. Nishiyama, C. J. Chang-Hasnain, and S. Arai, “Bandwidth enhancement of injection-locked distributed reflector lasers with wirelike active regions,” Opt. Express 18(16), 16370–16378 (2010). [CrossRef]   [PubMed]  

6. R. Martini, R. Paiella, C. Gmachl, F. Capasso, E. A. Whittaker, H. C. Liu, H. Y. Hwang, D. L. Sivco, J. N. Baillargeon, and A. Y. Cho, “High-speed digital data transmission using mid-infrared quantum cascade lasers,” Electron. Lett. 37(21), 1290–1291 (2001). [CrossRef]  

7. R. Martini, C. Bethea, F. Capasso, C. Gmachl, R. Paiella, E. A. Whittaker, H. Y. Hwang, D. L. Sivco, J. N. Baillargeon, and A. Y. Cho, “Free-space optical transmission of multimedia satellite data streams using mid-infrared quantum cascade lasers,” Electron. Lett. 38(4), 181–183 (2002). [CrossRef]  

8. A. Lyakh, C. Pflügl, L. Diehl, Q. J. Wang, F. Capasso, X. J. Wang, J. Y. Fan, T. Tanban-Ek, R. Maulini, A. Tsekoun, R. Go, and C. K. N. Patel, “1.6 W high wall plug efficiency, continuous-wave room temperature quantum cascade laser emitting at 4.6μm,” Appl. Phys. Lett. 92, 111110 (2008).

9. Y. Bai, S. R. Darvish, S. Slivken, W. Zhang, A. Evans, J. Nguyen, and M. Razeghi, “Room temperature continuous wave operation of quantum cascade lasers with watt-level optical power,” Appl. Phys. Lett. 92(10), 101105 (2008). [CrossRef]  

10. P. Corrigan, R. Martini, E. A. Whittaker, and C. Bethea, “Quantum cascade lasers and the Kruse model in free space optical communication,” Opt. Express 17(6), 4355–4359 (2009). [CrossRef]   [PubMed]  

11. C. Y. L. Cheung, P. S. Spencer, and K. A. Shore, “Modulation bandwidth optimization for unipolar intersubband semiconductor lasers,” IEE Proc.: Optoelectron. 144, 44–47 (1997). [CrossRef]  

12. C. Y. Cheung and K. A. Shore, “Self-consistent analysis of dc modulation response of unipolar semiconductor lasers,” J. Mod. Opt. 45(6), 1219–1229 (1998). [CrossRef]  

13. N. Mustafa, L. Pesquera, C. Y. L. Cheung, and K. A. Shore, “Terahertz bandwidth prediction for amplitude modulation response of unipolar intersubband semiconductor lasers,” IEEE Photon. Technol. Lett. 11(5), 527–529 (1999). [CrossRef]  

14. M. K. Haldar, “A simplified analysis of direct intensity modulation of quantum cascade lasers,” IEEE J. Quantum Electron. 41(11), 1349–1355 (2005). [CrossRef]  

15. R. Paiella, R. Martini, F. Capasso, C. Gmachl, H. Y. Hwang, D. L. Sivco, J. N. Baillargeon, A. Y. Cho, E. A. Whittaker, and H. C. Liu, “High-frequency modulation without the relaxation oscillation resonance in quantum cascade lasers,” Appl. Phys. Lett. 79(16), 2526–2528 (2001). [CrossRef]  

16. S. Barbieri, W. Maineult, S. S. Dhillon, C. Sirtori, J. Alton, N. Breuil, H. E. Beere, and D. A. Ritchie, “13 GHz direct modulation of terahertz quantum cascade lasers,” Appl. Phys. Lett. 91(14), 143510 (2007). [CrossRef]  

17. A. Lyakh, R. Maulini, A. Tsekoun, R. Go, C. Pflügl, L. Diehl, Q. J. Wang, F. Capasso, C. Kumar, and N. Patel, “3 W continuous-wave room temperature single-facet emission from quantum cascade lasers based on nonresonant extraction design approach,” Appl. Phys. Lett. 95(14), 141113 (2009).

18. Q. J. Wang, C. Pflügl, L. Diehl, F. Capasso, T. Edamura, S. Furuta, M. Yamanishi, and H. Kan, “High performance quantum cascade lasers based on three-phonon-resonance design,” Appl. Phys. Lett. 94(1), 011103 (2009). [CrossRef]  

19. F. Mogensen, H. Olesen, and G. Jacobsen, “Locking conditions and stability properties for a semiconductor laser with external light injection,” IEEE J. Quantum Electron. 21(7), 784–793 (1985). [CrossRef]  

20. D. Indjin, P. Harrison, R. W. Kelsall, and Z. Ikonić, “Self-consistent scattering theory of transport and output characteristics of quantum cascade lasers,” J. Appl. Phys. 91(11), 9019–9026 (2002). [CrossRef]  

21. F. Rana and R. J. Ram, “Current noise and photon noise in quantum cascade lasers,” Phys. Rev. B 65(12), 125313 (2002). [CrossRef]  

22. J. Faist, F. Capasso, D. L. Sivco, C. Sirtori, A. L. Hutchinson, and A. Y. Cho, “Quantum cascade laser,” Science 264(5158), 553–556 (1994). [CrossRef]   [PubMed]  

23. T. Aellen, R. Maulini, R. Terazzi, N. Hoyler, M. Giovannini, J. Faist, S. Blaser, and L. Hvozdara, “Direct measurement of the linewidth enhancement factor by optical heterodyning of an amplitude-modulated quantum cascade laser,” Appl. Phys. Lett. 89(9), 091121 (2006). [CrossRef]  

24. J. von Staden, T. Gensty, W. Elsäßer, G. Giuliani, and C. Mann, “Measurements of the α factor of a distributed-feedback quantum cascade laser by an optical feedback self-mixing technique,” Opt. Lett. 31(17), 2574–2576 (2006). [CrossRef]   [PubMed]  

25. M. Ishihara, T. Morimoto, S. Furuta, K. Kasahara, N. Akikusa, K. Fujita, and T. Edamura, “Linewidth enhancement factor of quantum cascade lasers with single phonon resonance-continuum depopulation structure on Peltier cooler,” Electron. Lett. 45(23), 1168–1169 (2009). [CrossRef]  

26. M. S. Taubman, T. L. Myers, B. D. Cannon, and R. M. Williams, “Stabilization, injection and control of quantum cascade lasers, and their application to chemical sensing in the infrared,” Spectrochim. Acta A Mol. Biomol. Spectrosc. 60(14), 3457–3468 (2004). [CrossRef]   [PubMed]  

27. P. Gellie, S. Barbieri, J. F. Lampin, P. Filloux, C. Manquest, C. Sirtori, I. Sagnes, S. P. Khanna, E. H. Linfield, A. G. Davies, H. Beere, and D. Ritchie, “Injection-locking of terahertz quantum cascade lasers up to 35 GHz using RF amplitude modulation,” Opt. Express 18(20), 20799–20816 (2010). [CrossRef]  

28. C. H. Henry, N. A. Olsson, and N. K. Dutta, “Locking range and stability of injection locked 1.54 μm InGaAsP semiconductor laser,” IEEE J. Quantum Electron. 21(8), 1152–1156 (1985). [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 (4)

Fig. 1
Fig. 1 Normalized modulation response of mid-infrared QCLs. (a) QCLs emitting at 4.6 μm with parameters shown in Table.1 (b) QCLs emitting at 9 μm with parameters shown in Table.1. The 3-dB bandwidth is indicated by the black dashed line.
Fig. 2
Fig. 2 Locking maps for QCLs emitting at (a) 4.6 μm and (b) 9 μm. The parameters used in the simulations are shown in Tables 1 and 4.
Fig. 3
Fig. 3 Normalized frequency response curves versus frequency detuning at a fixed injection ratio R = 5 dB, for the 9 μm QCL.
Fig. 4
Fig. 4 Normalized modulation response of mid-infrared QCLs under injection locking scheme. (a) QCLs emitting at 4.6 μm. (b) QCLs emitting at 9 μm. The inset show the corresponding poles (x) and zeros (), with the arrow indicating direction of increasing the injection current of the slave laser. The 3 dB bandwidth is indicated by the dashed line.

Tables (7)

Tables Icon

Table 1 Characteristic parameters of the state-of-the-art QCLs at 4.6 μm [17] and 9 μm [18] for direct modulation

Tables Icon

Table 2 Poles, zeros, and f3dB for QCLs at 4.6 μm

Tables Icon

Table 3 Poles, zeros, and f3dB for QCLs at 9 μm

Tables Icon

Table 4 Characteristic parameters of the state-of-the-art QCL at 4.6 μm [17] and 9 μm [18] for injection locking modulation

Tables Icon

Table 5 Poles, Zeros and 3-dB bandwidth f3dB for the QCL in Fig. 3.

Tables Icon

Table 6 Poles, zeros, and f3dB for the QCL in Fig. 4 (a)

Tables Icon

Table 7 Poles, zeros, and f3dB for Fig. 4 (b)

Equations (47)

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

d N 3 dt =J N 3 τ 3 G N p ( N 3 N 2 )P
d N 2 dt = N 3 τ 3 + G N P ( N 3 N 2 )P N 2 τ 2
dP dt =G( N 3 N 2 )P P τ P
σ 32 = 4π e 2 z 32 2 ε 0 n neff λ 0 ( 2 γ 32 )
J 0 N 30 τ 3 G 0 ( N 3 0 N 2 0 ) P 0 =0
N 30 τ 3 + G 0 ( N 3 0 N 2 0 ) P 0 N 20 τ 2 =0
G( N 30 N 20 ) P 0 P 0 τ P =0
N 20 = J 0 τ 2
J 0 ( 1η ) N 30 N 20 τ 3 G 0 ( N 30 N 20 ) P 0 =0
( N 30 N 20 )= 1 N p G 0 τ p
P 0 = 1 G 0 τ 3 ( J 0 J th 1 )
N 30 = 1 N P G 0 τ P + J 0 τ 2
dΔ N 3 dt =ΔJ Δ N 3 τ 3 G 0 ( Δ N 3 Δ N 2 ) P 0 G 0 ( N 30 N 20 )ΔP
dΔ N 2 dt = Δ N 3 τ 3 + G 0 ( Δ N 3 Δ N 2 ) P 0 + G 0 ( N 30 N 20 )ΔP Δ N 2 τ 2
dΔP dt =G( Δ N 3 Δ N 2 ) P 0 +G( N 30 N 20 )ΔP ΔP τ P
[ s+ f 11 f 12 f 13 f 21 s+ f 22 f 23 f 31 f 32 s+ f 33 ][ Δ N 3 Δ N 2 ΔP ]=[ ΔJ 0 0 ]
f 11 = 1 τ 3 + G 0 P 0 f 21 =( 1 τ 3 + G 0 P 0 ) f 31 =G P 0 f 12 = G 0 P 0 f 22 = 1 τ 2 + G 0 P 0 f 32 =G P 0 f 13 = 1 τ P N P f 23 = 1 τ P N P f 33 =0
H( s )= ΔP ΔJ( 1η ) τ P N P = s A +1 B s 3 +C s 2 +Ds+1
A=( 1 η 1 ) 1 τ 3 B= τ P τ 2 N P G P 0 C=B( 2 G 0 P 0 + 1 τ 3 ( 1+ 1 η ) )D=B( 1 τ 2 G 0 P 0 + 1 τ 3 τ 2 +2G P 0 1 τ P N P )
d dt E SL (t){ jω(N)+ 1 2 [ N P G 0 ( N ) 1 τ p ] } E SL ( t )= f d E ML ( t )
E SL ( t )= E 0 ( t ) e j( ω 0 t+ ϕ 0 ( t ) )
E ML ( t )= E 1 ( t ) e j( ω 1 t+ ϕ 1 ( t ) )
d dt E 0 ( t )= 1 2 N P G 0 ΔN( t ) E 0 ( t )+ f d E 1 cosϕ( t )
d dt ϕ( t )= 1 2 α N P G 0 ΔN( t ) f d E 1 E 0 sinϕ( t )Δ ω inj
dN( t ) dt =J( N+ N 2 ) 2 τ 3 2 E 0 2 G 0 N+ N 2 τ 2
d N 2 ( t ) dt =( N+ N 2 ) 1 τ 3 + E 0 2 G 0 N N 2 τ 2
N 2th =( E fr 2 γ P + N th τ 3 ) τ 2 τ 3 τ 3 τ 2
γ P = G 0 N th
J th = N th 1 τ 3 τ 2
E fr 2 = J N th 1 τ 3 τ 2 γ P τ 3 τ 3 τ 2
Δ N 0 = 2 f d E 1 N p G 0 E 0 cos ϕ 0
ϕ 0 = sin 1 ( Δ ω inj E 0 f d E 1 1+ α 2 ) tan 1 α
N 20 =[ E 0 2 ( γ P + G 0 Δ N 0 )+ Δ N 0 τ 3 + N th τ 3 ] τ 2 τ 3 τ 3 τ 2
T 1 = 2 τ 2 τ 3 τ 3 τ 2
T 2 = 2 τ 3 + T 1 τ 3
E 0 3 ( 2 G 0 N P γ P + G 0 N P γ P T 1 ) E 0 2 2 G 0 f d E 1 cos ϕ 0 ( 2+ T 1 ) E 0 ( G 0 J N P N th G 0 N P T 2 ) 2 f d E 1 cos ϕ 0 T 2 =0
d dt ΔE= 1 2 N P G 0 Δ N 0 ΔE f d E 1 sin ϕ 0 Δϕ+ 1 2 N p G 0 E 0 ΔN
d dt Δϕ= f d E 1 E 0 2 sin ϕ 0 ΔE f d E 1 E 0 cos ϕ 0 Δϕ+ α 2 N P G 0 ΔN
d dt ΔN=ΔJ( 4 E 0 γ P +4 E 0 G 0 Δ N 0 )ΔE( 2 τ 3 +2 E 0 2 G 0 )ΔN+( 1 τ 2 2 τ 3 )Δ N 2
d dt Δ N 2 =( 2 E 0 γ P +2 G 0 Δ N 0 E 0 )ΔE+( 1 τ 3 + E 0 2 G 0 )ΔN+( 1 τ 3 1 τ 2 )Δ N 2
[ Q 11 +s Q 12 Q 13 Q 14 Q 21 Q 22 +s Q 23 Q 24 Q 31 Q 32 Q 33 +s Q 34 Q 41 Q 42 Q 43 Q 44 +s ][ ΔE Δϕ ΔN Δ N 2 ]=[ 0 0 ΔJ 0 ]
Q 11 = 1 2 G 0 Δ N 0 N P Q 12 = f d E 1 sin ϕ 0 Q 13 = 1 2 G 0 E 0 N P Q 14 =0 Q 21 = f d E 1 E 0 2 sin ϕ 0 Q 22 = f d E 1 E 0 cos ϕ 0 Q 23 = α 2 G 0 N P Q 24 =0 Q 31 =4 E 0 γ P +4 E 0 G 0 Δ N 0 Q 32 =0 Q 33 = 2 τ 3 +2 E 0 2 G 0 Q 34 = 2 τ 3 1 τ 2 Q 41 =( 2 E 0 γ P +2 G 0 Δ N 0 E 0 ) Q 42 =0 Q 43 =( 1 τ 3 + E 0 2 G 0 ) Q 44 =( 1 τ 2 1 τ 3 )
H( s )=G s 2 +Es+F s 4 +A s 3 +B s 2 +Cs+D
A= Q 11 + Q 22 + Q 33 + Q 44 B= Q 33 Q 44 + Q 11 Q 33 + Q 11 Q 44 + Q 22 Q 33 + Q 22 Q 44 + Q 11 Q 22 Q 21 Q 12 Q 31 Q 13 Q 34 Q 43 C= Q 11 Q 33 Q 44 + Q 22 Q 33 Q 44 + Q 33 Q 11 Q 22 + Q 44 Q 11 Q 22 Q 33 Q 21 Q 12 Q 44 Q 21 Q 12 Q 13 Q 31 Q 22 Q 13 Q 31 Q 44 + Q 12 Q 31 Q 23 + Q 13 Q 34 Q 41 Q 22 Q 43 Q 34 Q 11 Q 34 Q 43 D= Q 11 Q 22 Q 33 Q 44 Q 12 Q 21 Q 33 Q 44 + Q 12 Q 21 Q 43 Q 34 Q 13 Q 31 Q 22 Q 44 + Q 12 Q 31 Q 23 Q 44 Q 12 Q 23 Q 34 Q 41 + Q 13 Q 22 Q 34 Q 41 Q 11 Q 22 Q 34 Q 43 E= ( Q 12 Q 23 Q 13 Q 22 Q 13 Q 44 ) /Q 13 F= ( Q 12 Q 23 Q 44 Q 13 Q 22 Q 44 ) / Q 13 G= Q 13
H( s )= ( s | z 1 | +1 ) ( s | p 1 | +1 )( s | p 2 | +1 )( s | p 3 | +1 )
H( s )= 1 sD+1
f 3dB 3 2πD = 3 2π[ ( τ P +2 τ 2 )+ τ p G 0 P 0 τ 3 ]
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.