Microscopic density matrix analysis on the linewidth enhancement factor (LEF) of both mid-infrared (mid-IR) and Terahertz (THz) quantum cascade lasers (QCLs) is reported, taking into account of the many body Coulomb interactions, coherence of resonant-tunneling transport and non-parabolicity. A non-zero LEF at the gain peak is obtained due to these combined microscopic effects. The results show that, for mid-IR QCLs, the many body Coulomb interaction and non-parabolicity contribute greatly to the non-zero LEF. In contrast, for THz QCLs, the many body Coulomb interactions and the resonant-tunneling effects greatly influence the LEF resulting in a non-zero value at the gain peak. This microscopic model not only partially explains the non-zero LEF of QCLs at the gain peak, which observed in the experiments for a while but cannot be explicitly explained, but also can be employed to improve the active region designs so as to reduce the LEF by optimizing the corresponding parameters.
© 2013 Optical Society of America
Quantum cascade lasers (QCLs) are unipolar laser sources relying on intersubband transitions in multiple quantum well systems. The light emission can be tuned across the mid-infrared (mid-IR, from 3 to 20 μm) and Terahertz (THz, from 1.2 to 5 THz, or 60 to 250 μm) ranges of the electromagnetic spectrum. As compact and coherent radiation sources, they have received considerable attentions since their first demonstrations [1, 2].
Linewidth enhancement factor (LEF, α) plays an important role in determining the optical emission linewidth of the laser systems  and the frequency responses . It is defined as the ratio of the carrier-induced variation of the real and imaginary parts of the susceptibility, and describes the amplitude-phase coupling . Owing to this coupling, a variation of the intensity of optical field will induce an excess perturbation of the phase of a laser mode, which influences the performance of lasers e.g. causes the laser linewidth increases well beyond the Schawlow-Townes limit. Typical value of LEF in interband semiconductor diode lasers is about 2-7 , resulting from the asymmetric differential gain spectrum caused by two bands associated with the laser transition with opposite curvature in the k-space. In contrast, due to the intersubband transition characteristics, QCLs are expected to have a narrow and symmetric gain spectrum, and hence, according to the Kramers-Kronig relation, resulting in a zero LEF at the peak gain wavelength predicted by Faist et al . However, the currently reported experiments [6–11] have demonstrated that the LEFs in both mid-IR and THz QCLs are not zero, although they are much smaller than those of diode lasers. For example, a LEF of up to −0.5 and 0.5 at the gain peak was reported at a lasing wavelength of 8.22 μm  and ~116 μm (~2.55 THz) , respectively. In addition, a strong dependence of the LEF on frequency detuning was observed, and a large value of up to −2 was reported for distributed feedback (DFB) mid-IR QCLs .
Although several experimental works have been carried out in studying the LEF of QCLs, few theoretical investigations on LEF has been reported  and the reasons behind the observed non-zero LEF at the gain peak has not been completely understood for both mid-IR and THz QCLs. In the previous work, the non-zero LEF was investigated in mid-IR QCLs by considering the refractive index change due to the device self-heating and the transitions not involved in the laser action based on the experiments and the macroscopic theoretical analysis . The results show that device self-heating is the dominate factor compared with the other factors. In addition to the device self-heating effect, other mechanisms such as many body Coulomb interactions, coherence of resonant tunneling (RT) and non-parabolicity can cause the non-zero LEF of QCLs at the gain peak, which cannot be considered in the macroscopic picture. Furthermore, a theoretical analysis on the LEF of THz QCLs is lacked. Since the active region structures of THz QCLs are different from those of mid-IR QCLs, they shall show different characteristic of non-zero LEF. To disclose more physical underlining mechanisms of non-zero LEF value in both mid-IR and THz QCLs, a microscopic model is required.
Recently, we have developed a microscopic density matrix model to examine the optical gain of THz QCLs demonstrating that the many body effects, which lead to renormalization of band structure and Rabi frequency, have significant modifications to the gain profile . In this paper, we extend this microscopic model to further investigate the role of many body Coulomb interactions, coherence of resonant tunneling and non-parabolicity on the LEF of both mid-IR and THz QCLs. The results show that, for mid-IR QCLs, both the many body Coulomb interactions and non-parabolicity contribute to the non-zero LEF among these three factors. Furthermore, the α-spectrum is blue-shifted as the injection current increases, which may explain the recent experimental result of α-value increases with the bias. In contrast, for THz QCLs, both the many body Coulomb interactions and the coherence of resonant tunneling greatly influence the LEF deviating from zero at the gain peak. Our results demonstrate the importance of the microscopic characteristics on the LEF of QCLs, and partially explain the discrepancy between the experimental results and the previous models.
2. LEF in mid-IR QCLs
2.1. The microscopic model
For mid-IR QCLs, in order to make a comparison with experiments, we use the active region with the three-well vertical design which has been used in the LEF measurement experiment [6, 13]. This structure is also very similar with the one used for another LEF experiment in . Figure 1 shows the conduction band diagram and magnitude squared envelope wave functions of this design in a “tight-binding” scheme. The coupling between the periods, achieved by tunneling (Ω51 is the coupling strength), is shown through the injector barriers. The energy states within one period are coupled through scattering processes, but at the injector barrier, the transport is modeled by tunneling. In order to conveniently treat the many body effects, we derive the dynamic equations of motion in the second quantized representation. The Hamiltonian of the system of mid-IR QCL in Fig. 1, which characterizes the electron-light coupling, the tunneling effects, free electrons and electron-electron Coulomb interactions, can be written as12, 14]. The parasitic coupling between levels 1 and 4 is neglected, which is a reasonably good approximation for mid-IR QCLs . The relaxation of electrons in level 3 into level 1 is characterized by an effective scattering rate.
Using the semiclassical laser theory and Maxwell equations, the gain G and carrier-induced refractive index change δn are given by
The linewidth enhancement factor α can be obtained by the ratio of the change in the real part of the refractive index change δn to the change in the gain G with respect to the carrier density N0 
2.2. Results and discussions
Figure 2 shows the α-spectrum at different biases calculated from the microscopic model including the many body Coulomb interactions, coherence of resonant-tunneling transport and non-parabolicity. The points shown in the figure indicate the positions of the gain peak. In contrast to the macroscopic calculations of zero LEFs the value of α-spectrum at gain peak cannot be neglected due to the asymmetry of gain spectrum caused by the interplay of many body interactions, non-parabolicity and tunneling effects. The obtained LEF is around −0.22. The sign of the α-value at gain peak reflects the asymmetry of gain spectrum relative to the lasing central frequency ωλ. The experimental α-value of this active region design at gain peak is around −0.5 . The absolute experimental value is 0.28 larger than our simulation result of –0.22. This discrepancy can be attributed to the refractive index change due to device self-heating effect demonstrated in , which needs a more complicated calculation on the thermal effects and is not taken into account in our simulations. Therefore, the microscopic contributions induced by many body interaction, non-parabolicity and coherence of resonant tunneling and refractive index change due to device self-heating are all important to determine the LEF of QCLs. Moreover, although the LEF at gain peak is smaller than 1, the LEF shows a large value as the lasing frequency is slightly away from the central frequency. This means that LEF strongly depends on the frequency detuning. Therefore, for the waveguide structure with grating, a laser will show a relatively large LEF at lasing frequency. In addition, for DFB-QCLs, the emission wavelengths will be shifted to longer wavelengths as the injection current increases , the value of LEF at the operation wavelength will rise as the external bias increases. Because the α-spectrum is blue-shifted as the external bias increases seen from Fig. 2, the LEF will increase as the injection current rises. This could explain the experimental observations in a QCL in , where the active region is very similar with the one used in this paper (see Fig. 1), that the LEF increases with the injection current. It is noted that the measured LEF in  does not include the contribution from the device self-heating effect due to self-mixing measurement technique used in this experiment.
Figure 3 shows the details of α-parameter at the gain peak position, the 0.18 μm blue-shift and red-shift positions relative to the peak under different biases computed from the microscopic model “many body + non-parabolicity” (considering both many body and non-parabolicity effects), the microscopic free-carrier model (considering both free carriers and non-parabolicity but neglecting the renormalization of band structure and Rabi frequency), the macroscopic density-matrix model with and without including the resonant-tunneling effect, respectively. As shown in Fig. 3(a), the many body Coulomb interaction, coherence of resonant tunneling and non-parabolicity all induce a finite value of LEF at the gain peak. Furthermore, by the comparison of these contributions to LEF, the many body Coulomb interaction and non-parabolicity parameter have more important effects on the LEF at gain peak. This is because the many body Coulomb interaction and non-parabolicity all tend to distort the shape of gain spectrum, as shown in Fig. 3(b), where the non-parabolicity can more greatly modify the gain spectrum as compared with Coulomb interaction.
The non-parabolicity in In0.53Ga0.47As quantum wells is higher than that in GaAs quantum wells. Therefore, the gain spectrum in mid-IR QCLs can be more greatly influenced by the non-parabolicity. If the non-parabolicity (proportional to the ratio of the effective mass at the upper laser level and the lower laser level) slightly increases, the LEF will increase according to Fig. 4(a).The increase is attributed to the influence of non-parabolicity on symmetry of gain spectrum, as shown in Fig. 4(b). Therefore, according to the above simulations, it is expected that, for similar structures, mid-IR QCL emitting at a shorter wavelength and with a higher non-parabolicity has a larger α-value.
3. LEF in THz QCLs
The present best temperature performance of THz QCLs is obtained by using resonant-phonon (RP) design  at an operation temperature of up to ~200 K. We consider this design in this paper. Figure 5 shows the conduction band diagram in the “tight-binding” scheme. The details of this structure can be found in . The Hamiltonian of this system can be expressed asEq. (4), we can obtain the LEF of THz QCLs.
Figure 6 shows the α-spectrum at different biases calculated from the microscopic model including the many body Coulomb interactions, coherence of resonant-tunneling transport and non-parabolicity. The points indicate the positions of the gain peak. Similarly, the value of α-spectrum of THz QCLs at gain peak cannot be neglected due to the asymmetry of gain spectrum caused by the interplay of many-body interaction, non-parabolicity and resonant-tunneling effects. Its value is around −0.7 when electric pumping does not exceed the designed bias. The sign of α-value at gain peak reflects the asymmetry of gain spectrum relative to the lasing central frequency ωλ. It is noted that, when the bias goes above the designed value and becomes 13.3 kV/cm, the LEF becomes −0.3. This is mainly attributed to the symmetric changes of gain spectrum with the extraction detuning due to tunneling effects i.e. the energy splitting due to the coupling between the lower laser level and the extraction level, as shown in Fig. 5 in . It has shown that the peak position and the lineshape of gain spectrum strongly follow the variations of extraction detuning. When the extraction detuning changes to a positive value from a negative one, the peak frequency is redshifted. Therefore, the absolute value of LEF is reduced. In addition, although the LEF at gain peak is smaller than 1, the LEF shows a large value as the lasing frequency is slightly away from the central frequency, which means that LEF of THz also strongly depends on the frequency detuning.
Figure 7 shows the details of α-parameter at frequencies of gain peak, 0.2 THz red-shift and 0.2 THz blue-shift relative to the peak position under different biases, computed from the microscopic model with “many body + non-parabolicity”, microscopic free-carrier model, the macroscopic models with and without resonant tunneling, respectively. According to the macroscopic models with and without resonant tunneling, the coherence of resonant tunneling contributes to a significant increase of the LEF at gain peak. Since tunneling exhibits an increasing broadening and modification for gain spectrum when injection level (level 1) and upper laser level (level 4), and lower laser (level 3) and extraction level (level 2) are in resonance simultaneously, the absolute α value rises with the increasing bias in the regime of negative injection and extraction detunings according to macroscopic model with resonant tunneling. Once the operation bias exceeds the designed bias, the sign of LEF is changed due to the variation of the symmetry of gain spectrum relative to the lasing central frequency ωλ. Furthermore, the non-parabolicity, in contrast to mid-IR QCLs, can only induce a slight influence on the LEF according to the comparisons between free-carrier model and macroscopic one with resonant tunneling, but the many body Coulomb interaction causes a large variation of LEF at gain peak with the comparison of the model “many-body + non-parabolicity” and free-carrier one due to its strong modifications to gain spectrum, as shown in Fig. 8. (The Coulomb interaction include the Hartree-Fock, dephasing and scattering contributions, more details can be seen in ). Overall, the non-zero LEF of THz QCLs is mainly due to the combined impacts from the Coulomb interaction and coherence of resonant-tunneling effects.
In conclusion, our study shows that the interplay of the many body interactions, coherence of resonant-tunneling effects and non-parabolicity for both mid-IR and THz QCL play an important role on the non-zero LEF at gain peak. A strong dependence of the LEF on the lasing frequency detuning is observed. The results show that, for mid-IR QCLs, the many body interactions and non-parabolicity all play an important key role in non-zero LEF at gain peak, where non-parabolicity has a more significant influence. In contrast, the many body Coulomb interactions and coherence of resonant tunneling have a significant impact on the LEF of THz QCLs, but the non-parabolicity only induces a little effect. Although the simulated LEF at gain peak is smaller than the measured value, our microscopic model can well explain some experimental observations e.g. the increase of LEF with injection current for DFB QCLs. The discrepancy between the experimental value and the proposed theoretical models can be attributed to the refractive index change due to device self-heating, which should also be considered in the future explorations.
According to the following Hershberger equation and the anticommutations relations of fermionic operator 12]. The influence of the subband dispersion, namely the nonparabolicity, is represented by using the effective mass of electrons . For subband j, we have . For our mid-IR structure, we estimate that [18, 19].
Similarly, one can obtain the following equations of motion for the polarizationand electron occupation for THz QCLs
Since levels (1′, 4) and (2, 3) are coherently coupled by the tunneling, the coherences corresponding to the levels (1′, 3), (2, 4) and (1′, 2) have a time-harmonic character due to the time-harmonic (3, 4) coherent coupling. Therefore, we try to look for solutions in the form of20].
We would like to acknowledge financial support from Nanyang Technological University (NTU) (M58040017), Defence Research and Technology Office, Singapore (TL-9009105606-01), and Ministry of Education, Singapore (MOE2011-T2-2-147). Support from the CNRS International-NTU-Thales Research Alliance (CINTRA) Laboratory, UMI 3288, Singapore 637553, is also acknowledged.
References and links
1. R. Köhler, A. Tredicucci, F. Beltram, H. E. Beere, E. H. Linfield, A. G. Davies, D. A. Ritchie, R. C. Iotti, and F. Rossi, “Terahertz semiconductor-heterostructure laser,” Nature 417(6885), 156–159 (2002). [CrossRef]
3. C. H. Henry, “Theory of the linewidth of semiconductor-lasers,” IEEE J. Quantum Electron. 18(2), 259–264 (1982). [CrossRef]
4. T. Chattopadhyay and P. Bhattacharyya, “Role of linewidth enhancement factor on the frequency response of the synchronized quantum cascade laser,” Opt. Commun. 309, 349–354 (2013). [CrossRef]
5. M. Osinski and J. Buus, “Linewidth broadening factor in semiconductor-lasers - an overview,” IEEE J. Quantum Electron. 23(1), 9–29 (1987). [CrossRef]
6. M. Lerttamrab, S. L. Chuang, C. Gmachl, D. L. Sivco, F. Capasso, and A. Y. Cho, “Linewidth enhancement factor of a type-I quantum-cascade laser,” J. Appl. Phys. 94(8), 5426–5428 (2003). [CrossRef]
7. R. P. Green, J. H. Xu, L. Mahler, A. Tredicucci, F. Beltram, G. Giuliani, H. E. Beere, and D. A. Ritchie, “Linewidth enhancement factor of terahertz quantum cascade lasers,” Appl. Phys. Lett. 92(7), 071106 (2008). [CrossRef]
8. 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]
9. J. von Staden, T. Gensty, W. Elsässer, 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]
10. N. Kumazaki, Y. Takagi, M. Ishihara, K. Kasahara, A. Sugiyama, N. Akikusa, and T. Edamura, “Detuning characteristics of the linewidth enhancement factor of a midinfrared quantum cascade laser,” Appl. Phys. Lett. 92(12), 121104 (2008). [CrossRef]
11. J. Kim, M. Lerttamrab, S. L. Chuang, C. Gmachl, D. L. Sivco, F. Capasso, and A. Y. Cho, “Theoretical and experimental study of optical gain and linewidth enhancement factor of type-I quantum-cascade lasers,” IEEE J. Quantum Electron. 40(12), 1663–1674 (2004). [CrossRef]
12. T. Liu, K. E. Lee, and Q. J. Wang, “Microscopic density matrix model for optical gain of terahertz quantum cascade lasers: Many-body, nonparabolicity, and resonant tunneling effects,” Phys. Rev. B 86(23), 235306 (2012). [CrossRef]
13. C. Gmachl, F. Capasso, J. Faist, A. L. Hutchinson, A. Tredicucci, D. L. Sivco, J. N. Baillargeon, S. N. G. Chu, and A. Y. Cho, “Continuous-wave and high-power pulsed operation of index-coupled distributed feedback quantum cascade laser at λ≈8.5 μm,” Appl. Phys. Lett. 72(12), 1430–1432 (1998). [CrossRef]
14. W. W. Chow, S. W. Koch, and M. S. I. I. I. Semiconductor-Laser Physics, (Springer-Verlag, Berlin, 1994).
15. S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]
16. C. Gmachl, F. Capasso, A. Tredicucci, D. L. Sivco, J. N. Baillargeon, A. L. Hutchinson, and A. Y. Cho, “High-power, continuous-wave, current-tunable, single-mode quantum-cascade distributed-feedback lasers at lambda ~5.2 and lambda ~7.95 μm,” Opt. Lett. 25(4), 230–232 (2000). [CrossRef]
17. S. Fathololoumi, E. Dupont, C. W. I. Chan, Z. R. Wasilewski, S. R. Laframboise, D. Ban, A. Mátyás, C. Jirauschek, Q. Hu, and H. C. Liu, “Terahertz quantum cascade lasers operating up to ~200 K with optimized oscillator strength and improved injection tunneling,” Opt. Express 20(4), 3866–3876 (2012). [CrossRef]
19. S. Panda, B. K. Panda, and S. Fung, “Effect of conduction band nonparabolicity on the dark current in a quantum well infrared detector,” J. Appl. Phys. 101(4), 043705 (2007). [CrossRef]
20. E. Dupont, S. Fathololoumi, and H. C. Liu, “Simplified density-matrix model applied to three-well terahertz quantum cascade lasers,” Phys. Rev. B 81(20), 205311 (2010). [CrossRef]