Active mode-locking of mid-infrared quantum cascade lasers with short gain recovery time

Yongrui Wang; Alexey Belyanin

doi:10.1364/OE.23.004173

1. Introduction

Generation of ultrashort pulses in mid-infrared (mid-IR) quantum cascade lasers (QCLs) is one of the last remaining fundamental challenges in QCL physics and laser physics in general. The gain recovery time of a typical QCL is of the order of 1 ps, which is much shorter than the cavity round-trip time of order 50 ps. This effectively prohibits passive mode locking with a saturable absorber of any kind. Indeed, any initial intensity fluctuation in a cavity will quickly get damped if the gain recovers to its small-signal value immediately after the passage of the pulse, leading to preferential amplification of the tails of the pulse.

The situation with active mode locking (AML) is not so obvious, although there is a widespread belief that a long gain recovery time is essential in this case as well. AML has been observed in terahertz QCLs where the gain recovery time is believed to be significantly longer than in the mid-IR QCLs; see e.g. [1, 2]. AML was also observed in two-section mid-IR QCLs of a ”super-diagonal” design in which the laser transition is diagonal in real space so that the upper-state lifetime can be as long as 40–50 ps [3, 4]. The latter papers also pointed out that active mode locking in monolithic Fabry-Perot lasers appears to be limited by the spatial hole burning (SHB), i.e. a transient population grating created by a standing wave pattern inside the cavity. Scattering of the laser field off the population grating leads to the proliferation of modes with uncorrelated phases, which destroys single-pulse operation and leads to a chaotic output [3–5]. To avoid the detrimental effects of the SHB, [6] proposed active modulation of an external ring cavity QCL where the field propagates in only one direction and the population grating is not formed. The active device modeled in [6] was assumed to have a vertical laser transition and a short gain recovery time. Simulations in [6] showed a stable and robust generation of mode-locked pulses, but the pulse duration was very long for modulation at the round-trip period. Reaching the output pulse duration of ∼ 5 ps required the gain modulation with ultrashort current pulses of 20–30 ps, which would be difficult to synthesize.

The studies in [6] revealed that short gain recovery does not prevent mode locking as long as the net gain window is opened for a short time as compared to the cavity round-trip time. To achieve this, the modulation of only a small part of the laser cavity seems to be crucial. Indeed if the gain responds to modulation instantaneously, modulation of the whole cavity will destroy an isolated pulse as it will experience both gain and loss during each roundtrip. According to our simulations below, in this case the steady-state output is either intensity modulated or simply dies out.

Opening a short gain window is straightforward in an external cavity. However, for many applications a compact monolithic design is preferable. Here we show that generation of mode-locked pulses of a few ps duration is feasible in a monolithic Fabry-Perot cavity even for QCLs with a short gain recovery time, as long as only a short section of the cavity is modulated at the cavity round-trip time. The cavity design is similar to the one employed in [3, 4]. However in contrast to those works we focus on a high-performance active region design with a vertical laser transition and a short gain recovery time. We compare the performance of lasers with long and short gain recovery times and find that standard QCLs with gain relaxation times of the order of 1 ps are in fact preferable for the AML. Indeed, in super-diagonal lasers electrons accumulate in the upper laser state and the modulation of injection is rather weak. The gain modulation could possibly be achieved by varying bias, which leads to variation of the oscillator strength. However, the DC bias modulation in the active region of a super-diagonal laser is reduced by a strong space charge effect and is accompanied by the modulation of the upper state lifetime and the shift of the intersubband transition energy. These effects can reduce the resulting gain modulation and are not easy to control and utilize. In vertical transition lasers the density of upper-state electrons is much smaller than the electron density in the injector. The dominant effect of the bias modulation is the modulation of injection current and gain.

We use a more realistic model of the QCL active region and transport which includes resonant tunneling injection, electron distribution over in-plane k-vectors, and space charge. We also perform time-domain and space-domain simulations of the propagating field, in which we can directly follow the pulse formation and the pulse instabilities in time. Note that all previous studies cited above used the frequency-domain modal approach and a two- or three-level laser model without current injection. We find the single-pulse operation to be quite robust to the variation of parameters such as the DC bias, bias modulation, and the length of the modulated section. At the same time the phase coherence and pulse duration are sensitive to even slight variations the modulation period: at the level of 0.3%. Moreover the optimal modulation period values for the shortest pulse and longest phase coherence time do not coincide.

2. The model of the active region and main equations

We use a four-subband model to describe the QCL active region. It includes the injector/extractor states (hereafter the ground state), upper laser state, and lower laser state, denoted as g, u and l, respectively. A schematic diagram of the model and the two-section cavity are shown in Fig. 1.

Fig. 1 (a) A schematic of the active region model. The current J (thick horizontal line) is due to resonant tunneling between aligned k_‖ states of injector g and the upper laser subband u. The distance δ is the separation between the centroids of electron states in g and u, and Δ is the detuning between the bottoms of these two subbands. The bias electric field is defined to be zero when g and u are aligned. Dashed arrows denote the non-radiative transitions, mostly due to LO phonons. (b) A schematic of a two-section cavity with an RF modulation at the cavity round-trip time applied to a shorter section.

Download Full Size | PDF

The system dynamics is described by coupled density matrix and Maxwell equations:

\begin{array}{l} \partial_{t} n_{g} = \frac{n_{u}}{T_{u g}} + \frac{n_{l}}{T_{\lg}} - J + D \frac{\partial^{2} n_{g}}{\partial z^{2}}, \\ \partial_{t} n_{u} = J - \frac{n_{u}}{T_{u l}} - \frac{n_{u}}{T_{u g}} - i \frac{d E}{ℏ} (ρ_{u l} - ρ_{u l}^{*}) + D \frac{\partial^{2} n_{u}}{\partial z^{2}}, \\ \partial_{t} n_{l} = \frac{n_{u}}{T_{u l}} - \frac{n_{l}}{T_{\lg}} + i \frac{d E}{ℏ} (ρ_{u l} - ρ_{u l}^{*}) + D \frac{\partial^{2} n_{l}}{\partial z^{2}}, \\ \partial_{t} ρ_{u l} = - (i ω + \frac{1}{T_{2}}) ρ_{u l} - i \frac{d E}{ℏ} (n_{u} - n_{l}), \\ \partial_{z}^{2} E - \frac{n^{2}}{c^{2}} \partial_{t}^{2} E = \frac{Γ d}{ε_{0} c^{2} L_{p}} \partial_{t}^{2} (ρ_{u l} + ρ_{u l}^{*}), \end{array}

where n_g, n_u and n_l are the sheet densities of the corresponding states, ρ_ul is the off-diagonal density matrix element for a laser transition, d is the dipole moment, E is the electric field of the laser mode, n is the refractive index, L_p is the thickness of one period, Γ is the overlap factor, ħω is the resonance energy of optical transition, and D is the diffusion coefficient, which we will take to be zero in simulations in order to consider the most unfavorable for AML case of a strong SHB. The injection current J is given by [7]

J = \frac{e Ω^{2} γ}{ℏ (Δ^{2} + γ^{2})} {θ (Δ) (n_{g} - n_{u} e^{- | Δ | / k_{B} T}) + θ (- Δ) (n_{g} e^{- | Δ | / k_{B} T} - n_{u})} .

Here Ω and γ are the coupling energy and the broadening of the states across the injection barrier, θ(x) is the Heaviside function, and Δ is the energy detuning from alignment. If we define the bias at alignment to be 0, then Δ = eE_extδ, where E_ext is the bias electric field, and δ is the separation between the centroids of states g and u.

Next, we generalize the ansatzes made in [5]:

\begin{array}{l} E (z, t) = \frac{1}{2} [E_{+} (z, t) e^{- i (ω t - k z)} + E_{+}^{*} (z, t) e^{i (ω t - k z)}] \\ + \frac{1}{2} [E_{-} (z, t) e^{- i (ω t + k z)} + E_{+}^{*} (z, t) e^{i (ω t + k z)}], \\ ρ_{u l} (z, t) = η_{+} e^{- i (ω t - k z)} + η_{-} e^{- i (ω t + k z)}, \\ n_{g} (z, t) = n_{g 0} + n_{g 2} e^{2 i k z} + n_{g 2}^{*} e^{- 2 i k z}, \\ n_{u} (z, t) = n_{u 0} + n_{u 2} e^{2 i k z} + n_{u 2}^{*} e^{- 2 i k z}, \\ n_{l} (z, t) = n_{l 0} + n_{l 2} e^{2 i k z} + n_{l 2}^{*} e^{- 2 i k z}, \end{array}

where k = nω/c. E₊ and E₋ are amplitudes of the fields traveling in right and left directions respectively. All populations contain the slowly varying average part and the grating part, denoted by subscripts 0 and 2. The injection current J is a linear combination of n_g and n_u when the space charge effect is not included. So, J(z, t) should have the same form as electron populations:

J (z, t) = J_{0} + J_{2} e^{2 i k z} + J_{2}^{*} e^{- 2 i k z} .

Although we include the space charge in our work, if we ignore the effect of the population grating on energy detuning Δ as the first order approximation, the conclusion above still holds.

Substituting these expressions into the density matrix equations, and making the slowly varying envelope approximation, we get the following equations for the envelope functions,

\partial_{t} n_{g 0} = \frac{n_{u 0}}{T_{u g}} + \frac{n_{l 0}}{T_{\lg}} - J_{0},

\partial_{t} n_{g 2} = \frac{n_{u 2}}{T_{u g}} + \frac{n_{l 2}}{T_{\lg}} - J_{2} - 4 k^{2} D n_{g 2},

\partial_{t} n_{u 0} = J_{0} - \frac{n_{u 0}}{T_{u l}} - \frac{n_{u 0}}{T_{u g}} + i \frac{d}{2 ℏ} [E_{+} η_{+}^{*} + E_{-} η_{-}^{*} - c . c .],

\partial_{t} n_{u 2} = J_{2} - \frac{n_{u 2}}{T_{u l}} - \frac{n_{u 2}}{T_{u g}} - 4 k^{2} D n_{u 2} + i \frac{d}{2 ℏ} [E_{+} η_{-}^{*} - E_{-}^{*} η_{+}],

\partial_{t} n_{l 0} = \frac{n_{u 0}}{T_{u l}} - \frac{n_{l 0}}{T_{\lg}} - i \frac{d}{2 ℏ} [E_{+} η_{+}^{*} + E_{-} η_{-}^{*} - c . c .],

\partial_{t} n_{l 2} = \frac{n_{u 2}}{T_{u l}} - \frac{n_{l 2}}{T_{u g}} - 4 k^{2} D n_{l 2} - i \frac{d}{2 ℏ} [E_{+} η_{-}^{*} - E_{-}^{*} η_{+}],

\partial_{t} η_{+} = - i \frac{d}{2 ℏ} [(n_{u 0} - n_{l 0}) E_{+} + (n_{u 2} - n_{l 2}) E_{-}] - \frac{η_{+}}{T_{2}},

\partial_{t} η_{-} = - i \frac{d}{2 ℏ} [(n_{u 0} - n_{l 0}) E_{-} + (n_{u 2}^{*} - n_{l 2}^{*}) E_{+}] - \frac{η_{-}}{T_{2}},

(\frac{n}{c} \partial_{t} + \partial_{z}) E_{+} = i \frac{Γ d ω}{n ε_{0} c L_{p}} η_{+} - l_{w} E_{+},

(\frac{n}{c} \partial_{t} - \partial_{z}) E_{-} = i \frac{Γ d ω}{n ε_{0} c L_{p}} η_{-} - l_{w} E_{-} .

Here the waveguide loss l_w is added in the equations for the fields. These are the partial differential equations we solve numerically. The derivatives are approximated by the 2nd order expression, e.g. a quantity f(t) is updated at each step as

f ((n + 1) Δ t) = f (n Δ t) + Δ t \dot{f} (n Δ t) + \frac{1}{2} {(Δ t)}^{2} \ddot{f} (n Δ t) .

In experiment the laser oscillations in a cavity grow starting from the spontaneous noise, i.e. quantum fluctuations of the field and polarization. To include them in our modeling in a rigorous way we would have to perform second quantization of the field and electrons which would make the problem intractable. To adequately simulate the spontaneous noise we add a complex random source term to the right-hand side of the polarization equations Eqs. (11) and (12) which has zero average and is delta-correlated in space and time within the accuracy of the space-time numerical grid. The maximum magnitude of the random source is βn_u0/T_rad where T_rad ≃ 10⁻⁷ s is the spontaneous radiative lifetime of the laser transition and β is the small geometric factor measuring the fraction of spontaneous emission coupled to the transverse laser mode. We also added a random delta-correlated electric field along the cavity at the initial moment of time to model initial fluctuations of the field. We performed numerous simulations with different values of the β parameter and initial noise. Figure 5 below shows one example of simulations for β = 10⁻⁴ and the initial random field of magnitude 10⁻³ V/cm. We found that the magnitude of fluctuations only affects the field dynamics during the initial transient stage of the laser field formation. The resulting pulses are not sensitive to the spontaneous noise level. We start observing noticeable changes to the shape of the pulses only when β approaches 1.

3. Base set of parameters

For QCLs with a short gain recovery time, the parameters used in our simulation are: (1) Resonant tunneling: Ω = 1.5 meV, γ = 6.58 meV, δ = 10 nm; (2) Lifetimes: T_ul = 1 ps, T_ug = 3 ps, T_lg = 0.1 ps, T₂ = 0.05 ps, D = 0; (3) Optical parameters: n = 3.3, l_w = 10 cm⁻¹, Γ = 1, d = e × 2 nm, λ₀ = 6.2 μm; (4) QCL structure parameters: L_p = 580Å, n_doping = 8.0 × 10¹⁰ cm⁻², T_e = 300 K, cavity length L_act = 3 mm, facet reflection coefficients R_l = 0.5 and R_r = 0.5. Those parameters are typical for mid-IR QCLs. The threshold bias field obtained is 8.31 kV/cm below the alignment bias. In this paper, we use the bias at alignment as a reference, defined to be 0, then the threshold bias field is V_th = −8.31 kV/cm. For a laser with a long gain recovery time T_ul = 50 ps, three parameters in the active region are different from the short gain recovery time design, namely Ω = 0.15 × 1.5 meV, T_ul = 50 ps, and T_ug = 100 ps. The corresponding threshold bias is V_th = −9.31 kV/cm, which is close to the threshold in the T_ul = 1 ps case.

We choose the modulated section near the left facet, so that the modulation needs to be applied at the cavity round-trip time. The modulated section length is denoted by l_mod, in units of L_act. The DC bias V_DC is applied to the remaining part of the cavity. The bias on the modulated section is written as

V_{Mod} = V_{Mod, DC} + V_{Mod, Amp} sin (2 π t / T_{\mod}) .

All biases are normalized to the threshold bias V_th. The modulation period T_mod and all times on the plots are normalized to the phase roundtrip time T_round = 2n_rL_act/c = 66 ps. In order to study the dependence of the QCL performance on modulation parameters, we choose the following base set of parameters: l_mod = 0.1, V_DC = 1.0, V_Mod,DC = 1.0, V_Mod,Amp = 0.5, and T_mod = 1.0. We will change only one parameter at a time. Note that the threshold bias V_th is negative, so V > 1 (V < 1) means below (above) threshold.

4. The effects of injection pumping and spatial hole burning

First, we analyze the effect of resonant tunneling. In previous works [3, 4], the current or gain were assumed to be uniform in each section of the cavity and sinusoidally modulated. This assumption looks reasonable at small modulation amplitudes but becomes problematic with deeper modulation. In particular a large modulation amplitude as compared to the DC pumping level in [4] would create large negative swings in the current and gain, which is unrealistic for QCLs. Here we adopt a more realistic point of view that the injection current to the upper laser state results from electron transport under the applied bias and cannot be directly controlled. We assume that the input parameter that we can control and modulate is an applied bias whereas the current is calculated from Eqs. (2) and (5)–(14) assuming the Fermi-Dirac distribution of electrons in the injector with a given temperature T_e. In this case a large modulation amplitude would move the injector out of resonance with the upper laser state, thus decreasing the gain but not flipping its sign and not reversing the current.

Figure 2 shows the spectrum of the bias, injection current, and gain (c) near the left facet assuming a sinusoidal modulation of the bias at a relatively large amplitude V_Mod,Amp = 0.9. The zero frequency part has been subtracted, since it only contains information about the average of the corresponding quantities. From this picture we can see that the response of the current to the sinusoidal modulation of the bias is noticeably nonlinear as it contains a number of higher harmonics. For a smaller modulation amplitude V_Mod,Amp = 0.5 only the second harmonic remains in the spectrum at a magnitude smaller than −10 dB; therefore a sinusoidal modulation of the current would be a reasonable approximation.

Fig. 2 The spectrum of the modulated bias (a), injection current (b), and gain (c) at the point adjacent to the left facet. All parameters are taken at base values except for a higher modulation amplitude V_Mod,Amp = 0.9. Although the modulation of the bias is at a single frequency, the injection current and gain contain higher harmonics.

Download Full Size | PDF

The effect of the spatial hole burning (SHB) described by the population and injection current gratings at the half-wavelength period in Eqs. (3) and (4) is significant, although in contrast to statements in previous studies, we do not find SHB in a monolithic Fabry-Perot cavity to be the main obstacle for the AML. The effect of SHB is more nuanced as illustrated in Figs. 3 and 4 which show the dependence of the pulse shape on the DC bias for lasers with short and long relaxation time. Here we keep the DC part of the bias in the modulation section V_mod,DC equal to the bias in the DC section V_DC, to minimize the number of free parameters. For a laser with T_ul = 1 ps, both the intensity and pulse duration increase when the DC bias is increased. Also, rapid oscillations within the pulse envelope are developed at high enough bias. Comparison with simulations carried out in the absence of the populations grating, n_(u,l,g)2 = 0 and J₂ = 0 (Fig. 4) show that the substructure is due to the SHB instability because it is not present in the absence of SHB. The same conclusion is true for the case of T_ul = 50 ps, where the multiple pulsations disappear when n_(u,l,g)2 = 0. At the same time, an overall increase in the pulse duration is due to an increase in the net gain window and resulting amplification of the tails of the pulse. We again notice that in a laser with a short gain recovery time the single pulse regime is more robust.

Fig. 3 The output field amplitude over one roundtrip time taken close to the end of the simulation time for different values of the DC bias V_DC in units of threshold bias for lasers with (a) T_ul = 1 ps and (b) T_ul = 50 ps. The DC biases in the modulated section and DC section are set to be equal.

Download Full Size | PDF

Fig. 4 Same as Fig. 3, but without population grating.

Download Full Size | PDF

5. The effect of the modulation period on the pulse duration and phase coherence

An apparent optimal choice of the modulation period is the phase round-trip time, i.e. normalized T_mod = 1 in Eq. (16). However, this turns out to be not the best choice for the long-time stability of the output. In Fig. 5 we plot the output field at the left facet for the whole simulation range of 10,000 roundtrips. As one can see, it takes a long time before the output intensity eventually stabilizes, although there is still a slight long-term variation in the magnitude of the pulses. Furthermore, although the amplitude of the output field is nearly periodic, the real (Fig. 5(b)) and imaginary (not shown) parts of the field are not. This means that the spectrum of the output field is not phase locked. The autocorrelation function for the real part of the field, defined as

R_{E} (τ) = \int_{t_{1}}^{t_{2}} Re [E (t)] Re [E (t + τ)] d t / \int_{t_{1}}^{t_{2}} (Re {[E (t)]}^{2}) d t,

experiences oscillations between 1 and 0 with a period of about 200 round-trips.

To understand the reason for these long-time variations, in Fig. 6 we show the output laser field from the left facet for the base set of parameters over the time interval of 5 periods, taken close to the end of the simulation time. Along with the laser field we show the injection current and net waveguide gain g − l_w (excluding mirror losses) near the left facet. The duration of intensity pulses is about 10% of the round-trip time, which is around 7 ps. Although the injection current and gain are shown for a point near the left facet, this point can represent the whole modulated section since the gain saturation effect is small at these intensities. From Fig. 6(a) we can see that the gain follows the injection current, while the peak of the pulse is delayed with respect to the gain maximum, namely the pulse is experiencing the tailing edge of the gain. This can be explained by the fact that the group velocity of the pulse is smaller than the phase velocity and therefore the gain modulation is slightly out of phase with the circulating pulse. This can be visualized in a movie ( Media 1, Figs. 6(d) and 6(e)) which shows the dynamics of the field intensity and gain over the two modulation periods. Therefore, the gain modulation is slightly out of resonance with the circulating pulse which causes long-period oscillations in the output.

Fig. 5 (a) The output field amplitude and (b) real part of the field at the left facet for the whole simulation range of 10,000 roundtrips and the base set of parameters. After the output is stable, there are still some small oscillations in the pulse amplitude. The real part of the field is not periodic in the modulation period, indicating that the phase is not locked.

Download Full Size | PDF

Fig. 6 (a) The amplitude of the output laser field on the left facet. (b) The injection current at the point adjacent to the left facet. (c) The difference between the gain and waveguide loss g − l_w at the point adjacent to the left facet. The gain follows the injection current almost instantaneously, due to the short gain recovery time. The peak of the pulse has a delay with respect to the maximum of the gain. (d, e) A frame in the movie ( Media 1, size: 3.50 MB) of the evolution of field intensity (d) and g − l_w (e) in the same laser.

Download Full Size | PDF

To eliminate the phase mismatch we slightly increase the modulation period to T_mod = 1.003 to compensate for the group delay of the pulses. As shown in Fig. 7, this leads to the stable output after the laser field build-up time of about 200 round-trips. Moreover, the real and imaginary parts of the field are also exactly periodic as can be verified by the autocorrelation function which has a constant value equal to 1. At the same time, the pulse duration is longer in this case as illustrated in Fig. 8 which shows the pulse shape as a function of the modulation period T_mod. When T_mod = 1, the pulses generated by a laser with a short gain recovery time have the shortest duration. When T_mod = 1.003, the intensity of the pulse increases significantly as it has the best overlap with the maximum of the gain in a modulated section. At the same time, the time window of the net gain is increased in this case, leading to a longer pulse duration.

Fig. 7 (a) The output field amplitude and (b) real part of the field at the left facet over 1,000 roundtrips in a laser with a short gain recovery time for the base set of parameters but a slightly longer modulation period T_mod = 1.003, which matches the group roundtrip time of the pulse. Both the field amplitude and the real part of the field become strictly periodic after about 200 roundtrips.

Download Full Size | PDF

Fig. 8 The output field amplitude over one roundtrip time taken close to the end of the simulation time for different values of the modulation period T_mod measured in units of phase roundtrip time for lasers with (a) T_ul = 1 ps and (b) T_ul = 50 ps and the base set of parameters.

Download Full Size | PDF

6. The output dependence on the modulation amplitude and length of the modulated section

Two other parameters that strongly affect the output pulses are the modulation amplitude and length of the modulated section. Figure 9 shows the effect of varying modulation amplitude. As expected, for both T_ul = 1 ps and T_ul = 50 ps cases, the peak intensity of the pulse increases with the modulation amplitude. At the same time, the pulse becomes more asymmetric in the T_ul = 1 ps case, which is the result of the phase mismatch between the modulation period and group round-trip time. For a laser with a long gain recovery time, the T_mod = 1 modulation period is closer to resonance with the group round trip time and the pulses are more symmetric. Also the gain saturation effect is stronger, so the tailing edge of pulse is not amplified.

Fig. 9 The output field amplitude over one roundtrip time taken close to the end of the simulation time for different values of the modulation amplitude V_Mod,Amp in Eq. (16), measured in units of threshold bias for lasers with (a) T_ul = 1 ps and (b) T_ul = 50 ps and the base set of parameters. No significant output is generated when V_Mod,Amp ≤ 0.1.

Download Full Size | PDF

With increasing length of the modulated section, the pulse experiences a stronger round-trip gain but at the same time the window of the net gain becomes wider which leads to the pulse broadening. Fig. 10 shows the dependence of the output on the length of a modulated section. In the T_ul = 1 ps case, when l_mod increases from 0.1 to 0.5, the intensity and duration of the pulse increase at the same time. For longer l_mod the pulse eventually broadens to the whole cavity and its peak intensity drops. However, the pulsation still exists even when the modulated section is 90% of the whole chip. In the T_ul = 50 ps case, one only gets a good pulse around l_mod = 0.3, and this pulse is shorter as compared to the T_ul = 1 ps design. However, with increasing l_mod the output quickly turns into multiple chaotic pulsations. It is indeed expected that a laser with a long population relaxation time comparable to the round-trip time is more prone to instabilities than the laser with a short recovery time in which the pulsations of the population inversion are strongly damped.

Fig. 10 The output field amplitude over one roundtrip time taken close to the end of the simulation time for different lengths of the modulated section l_mod measured in units of the total cavity length L_act for lasers with (a) T_ul = 1 ps and (b) T_ul = 50 ps and the base set of parameters.

Download Full Size | PDF

When the whole QCL cavity is modulated, no isolated pulses are obtained in both cases, and only the intensity is modulated. We checked that this conclusion remains true for the parameters typical for THz QCLs. The observation of isolated pulses in [1, 2] could be the result of the nonuniform injection of modulated pumping along the cavity. Indeed the THz waveguide forms a microstrip transmission line for the microwave modulation which forms a standing wave along the cavity.

7. Conclusion

In conclusion, we showed that standard high-performance mid-IR QCLs with a short gain recovery time can be actively modulated to generate mode-locked pulses of a few ps duration when the modulation is applied only to a short section of a monolithic Fabry-Perot cavity to create a net gain window for a circulating pulse. The performance of QCL with a short gain recovery time is more robust to varying parameters as compared to QCLs with a long gain recovery time. By fine tuning the modulation period to the group round-trip time of the pulse, one can phase lock the electric field of the pulse and use its spectrum as a frequency comb.

The peak power of the generated pulses grows with increasing modulation amplitude and DC bias, but is limited by the concomitant pulse broadening due to widening gain window, as is seen in Figs. 9 and 3. For a base set of parameters the peak power of the pulses is about 16 times higher than the CW output power for the DC bias equal to the time-averaged modulated bias of Eq. (16).

Acknowledgments

This work has been supported in part by NSF grants ECCS-1230517 and EEC-0540832, and by the AFOSR grant FA9550-14-1-0376. We are grateful to IQCLSW 2014 participants for numerous discussions and to anonymous reviewers for their comments that helped to improve the paper.

References and links

1. S. Barbieri, M. Ravaro, P. Gellie, G. Santarelli, C. Manquest, C. Sirtori, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Coherent sampling of active mode-locked terahertz quantum cascade lasers and frequency synthesis,” Nat. Photonics 5, 306–313 (2011). [CrossRef]

2. J. Freeman, J. Maysonnave, N. Jukam, P. Cavali, K. Maussang, H. Beere, D. Ritchie, J. Mangeney, S. Dhillon, and J. Tignon, “Direct intensity sampling of a modelocked terahertz quantum cascade laser,” Appl. Phys. Lett. 101, 181115 (2012). [CrossRef]

3. C. Y. Wang, L. Kuznetsova, V. M. Gkortsas, L. Diehl, F. X. Kärtner, M. A. Belkin, A. Belyanin, X. Li, D. Ham, H. Schneider, P. Grant, C. Y. Song, S. Haffouz, Z. R. Wasilewski, H. Liu, and F. Capasso, “Mode-locked pulses from mid-infrared quantum cascade lasers,” Opt. Express 17, 12929–12943 (2009). [CrossRef] [PubMed]

4. V.-M. Gkortsas, C. Wang, L. Kuznetsova, L. Diehl, A. Gordon, C. Jirauschek, M. A. Belkin, A. Belyanin, F. Capasso, and F. X. Kärtner, “Dynamics of actively mode-locked quantum cascade lasers,” Opt. Express 18, 13616–13630 (2010). [CrossRef] [PubMed]

5. A. Gordon, C. Y. Wang, L. Diehl, F. X. Kärtner, A. Belyanin, D. Bour, S. Corzine, G. Höfler, H. C. Liu, H. Schneider, T. Maier, M. Troccoli, J. Faist, and F. Capasso, “Multimode regimes in quantum cascade lasers: From coherent instabilities to spatial hole burning,” Phys. Rev. A 77, 053804 (2008). [CrossRef]

6. A. K. Wójcik, P. Malara, R. Blanchard, T. S. Mansuripur, F. Capasso, and A. Belyanin, “Generation of picosecond pulses and frequency combs in actively mode locked external ring cavity quantum cascade lasers,” Appl. Phys. Lett. 103, 231102 (2013). [CrossRef]

7. R. Terazzi, T. Gresch, A. Wittmann, and J. Faist, “Sequential resonant tunneling in quantum cascade lasers,” Phys. Rev. B 78, 155328 (2008). [CrossRef]

Active mode-locking of mid-infrared quantum cascade lasers with short gain recovery time

Abstract

1. Introduction

2. The model of the active region and main equations

3. Base set of parameters

4. The effects of injection pumping and spatial hole burning

5. The effect of the modulation period on the pulse duration and phase coherence

6. The output dependence on the modulation amplitude and length of the modulated section

7. Conclusion

Acknowledgments

References and links

Supplementary Material (1)

Cited By

Figures (10)

Equations (17)

Optics Express