Effects of resonant tunneling and dynamics of coherent interaction on intrinsic linewidth of quantum cascade lasers

doi:10.1364/OE.20.017145

« Show journal navigation

Effects of resonant tunneling and dynamics of coherent interaction on intrinsic linewidth of quantum cascade lasers

Tao Liu, Kenneth E. Lee, and Qi Jie Wang »View Author Affiliations

Optics Express, Vol. 20, Issue 15, pp. 17145-17159 (2012)
http://dx.doi.org/10.1364/OE.20.017145

View Full Text Article

Acrobat PDF (1357 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Quantum Langevin equations
3. Laser intrinsic linewidths
4. Results and discussions
5. Conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (24)
Cited By
Figures (9)
Metrics

Abstract

A theoretical model for calculation of the intrinsic linewidth of QCLs is built on the basis of the quantum Langevin approach. It differs from the traditional rate equation model in that the resonant tunneling and the dynamics of coherent interaction can be considered. Results show that the coupling strength and the dephasing rate associated with resonant tunneling strongly affect the linewidth of THz QCLs in the incoherent resonant-tunneling transport regime but only induce little influence in the coherent regime. The dynamics of coherent interaction and resonant-tunneling transport show insignificant effects on the linewidth calculation of mid-infrared QCLs due to strong coupling in resonant tunneling. We also demonstrate that by properly designing the active regions of QCLs, one can reduce the intrinsic linewidth according to our model.

1. Introduction

Quantum cascade lasers (QCLs) have been a subject of great interest since their invention in 1994 [1

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]

]. Since then, the emission wavelength can be engineered across the mid-infrared (mid-IR) (~3-24 μm) and Terahertz (THz) (~1.2-5 THz, or ~60-250 μm) regions. A thorough understanding of the linewidth of QCLs is increasingly important, as it is related to many practical applications e.g. trace-gas absorption spectroscopy and optical free-space data communication [2

T. Liu and Q. J. Wang, “Fundamental frequency noise and linewidth broadening caused by intrinsic temperature fluctuations in quantum cascade lasers,” Phys. Rev. B 84(12), 125322 (2011). [CrossRef]

–7

L. Tombez, S. Schilt, J. Di Francesco, P. Thomann, and D. Hofstetter, “Temperature dependence of the frequency noise in a mid-IR DFB quantum cascade laser from cryogenic to room temperature,” Opt. Express 20(7), 6851–6859 (2012). [CrossRef] [PubMed]

]. Recently, an ultra-narrow intrinsic linewidth of ~510 Hz has been experimentally observed in a distributed feedback mid-IR QCL at I₀/I_th = 1.54 (I₀ is the operation current and I_th is the laser threshold current) [8

S. Bartalini, S. Borri, P. Cancio, A. Castrillo, I. Galli, G. Giusfredi, D. Mazzotti, L. Gianfrani, and P. De Natale, “Observing the intrinsic linewidth of a quantum-cascade laser: beyond the Schawlow-Townes limit,” Phys. Rev. Lett. 104(8), 083904 (2010). [CrossRef] [PubMed]

]. The experimental results are fitted by theoretical calculations derived from the rate equation model. Although the rate equation model can explain the reason behind the narrow intrinsic linewidth well [9

M. Yamanishi, T. Edamura, K. Fujita, N. Akikusa, and H. Kan, “Theory of the intrinsic linewidth of quantum-cascade lasers: hidden reason for the narrow linewidth and line-broadening by thermal photons,” IEEE J. Quantum Electron. 44(1), 12–29 (2008). [CrossRef]

], it cannot adequately describe the electron transport characteristics of QCLs, and therefore cannot accurately predict the actual linewidth value.

In the rate equation model, localization of wavefunctions due to the dephasing scattering is overlooked which can lead to unrealistic results and limit the usefulness of the model. As reported in Refs 10

S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]

–12

H. Callebaut and Q. Hu, “Importance of coherence for electron transport in terahertz quantum cascade lasers,” J. Appl. Phys. 98(10), 104505 (2005). [CrossRef]

, the effects of resonant tunneling and dynamics of coherent interaction (influenced by dephasing) are important when describing the transport of QCLs. The coherent (or incoherent) injection coupling plays an important role in the electrical injection process to populate the upper radiative level, which affects the electron populations associated with stimulated emission processes and hence the ratio of spontaneous emission coupled into the lasing mode to the net stimulated emission. As a result, the resonant tunneling and the dynamics of coherent interaction may finally affect the intrinsic linewidth of QCLs. In addition, the rate equation model uses the operating pumping current as an input parameter, thus knowledge on linewidth reduction through changing the active regions of QCLs, e.g. the coupling strength and the doping level, is restricted before obtaining experimental data.

The rate equation model also assumes that the polarization follows the other dynamic variables adiabatically, and therefore the dynamics of coherent interaction associated with laser transitions are neglected. But the polarization dynamics may help us interpret some physical phenomena such as electron memory effect induced intrinsic linewidth reduction [13

L. Davidovich, “Sub-Poissonian processes in quantum optics,” Rev. Mod. Phys. 68(1), 127–173 (1996). [CrossRef]

,14

G. P. Agrawal and C. M. Bowden, “Concept of linewidth enhancement factor in semiconductor lasers: its usefulness and limitations,” IEEE Photon. Technol. Lett. 5(6), 640–642 (1993). [CrossRef]

]. Therefore, in some circumstances a careful consideration of coherent interaction dynamics associated with the laser transition is essential.

Since the resonant-tunneling transport and dynamics of coherent interaction cannot be properly included in the classical rate equation model, one has to refer to the more fundamental quantum mechanical model. In this paper, we report a new model for calculation of intrinsic linewidth of QCLs based on the quantum Langevin equations. It differs from the results derived from the classical rate equation model in that the resonant-tunneling effects and dynamics of coherent interaction are considered. The results show that the coupling strength and the dephasing rate associated with resonant tunneling strongly affect the linewidth of QCLs in the incoherent resonant-tunneling transport regime but only induce little influence in the coherent regime. Finally, the intrinsic linewidth reduction through optimization of the doping density and lifetime of different energy levels is presented.

2. Quantum Langevin equations

We consider a three-level QCLs system, as shown in Fig. 1 .

Δ_{0} / 2

denotes the interaction of the resonant tunneling,

ℏ g

is the electron-light interaction. The electrons in level 1′ are injected into the upper laser level 3 by resonant-tunneling transport. Then, they relax from level 3 to level 2 by emitting photons. The Hamiltonian of this three-level QCLs system in the rotating-wave approximation can be obtained from Ref. 15

M. O. Scully and M. S. Zubairy, Quantum optics (Cambridge University Press, 1997).

by adding the resonant-tunneling term (the fourth item on the right side of the equation):

H = ℏ ω_{λ} c^{†} c + \sum_{j} ε_{j} b_{j}^{†} b_{j} + ℏ (g^{*} a^{†} b_{2}^{†} b_{3} + g b_{3}^{†} b_{2} a) - (Δ_{0} / 2) (b_{3}^{†} b_{1^{'}} + b_{1^{'}}^{†} b_{3})

(1)

where

ω_{λ}

is the single-mode field lasing angular frequency,

ε_{j}

is the energy of level j,

c^{†}

(

c

) and

b_{j}^{†}

(

b_{j}

) denote the creation (annihilation) operator of the laser field and electrons in the subband j, respectively. The coupling constant g is given by

g = \frac{e μ}{ℏ} \sqrt{\frac{ℏ ω_{λ} Γ}{2 ε_{0} n^{2} V_{m}}}

(2)

where

μ

is the radiative dipole matrix element,

Γ

is the confinement factor of the optical mode that overlaps with the entire active region,

ε_{0}

is the vacuum permittivity,

n

is the refractive index,

V_{m}

is the volume of one single gain stage in the active region.

Fig. 1 Schematic diagram showing the conduction-band of a three-level QCL structure and the square modulus of the wave functions. Levels 1′, 2 and 3 are the injector level, the lower laser level and the upper laser level, respectively. The inset showing the simplified energy levels with resonant tunneling and optical response.

Δ_{0} / 2

denotes the interaction of the resonant tunneling,

ℏ g

is the electron-light interaction.

The operators corresponding to the observables of our interest are the single-mode intracavity slowly varying laser field operator

a = c^{†} c \exp (i ω_{λ} t)

, electron polarization operator

σ_{23} = b_{2}^{†} b_{3} \exp (i ω_{λ} t)

and electron population operator

ρ_{j} = b_{j}^{†} b_{j}

in level j. By including the noise operators

f_{i j}

and

f_{j}

, the quantum Langevin equations which rule the dynamics of the laser are derived according to the Hershberger’s equation [15

M. O. Scully and M. S. Zubairy, Quantum optics (Cambridge University Press, 1997).

]

\frac{d a}{d t} = - \frac{κ}{2} a - i g^{*} σ_{23} + f_{a} (t)

(3a)

\frac{d σ_{23}}{d t} = - γ_{23} σ_{23} + i g (ρ_{3} - ρ_{2}) a + f_{23} (t)

(3b)

\frac{d σ_{3 1^{'}}}{d t} = - γ_{3 1^{'}} σ_{3 1^{'}} - i \frac{Δ_{3 1^{'}}}{ℏ} σ_{3 1^{'}} - \frac{i Δ_{0}}{2 ℏ} (ρ_{1^{'}} - ρ_{3}) + f_{3 1^{'}} (t)

(3c)

\frac{d ρ_{3}}{d t} = - γ_{s p} ρ_{3} - γ_{3} ρ_{3} + i (g^{*} a^{†} σ_{23} - g σ_{23}^{†} a) - \frac{i Δ_{0}}{2 ℏ} (σ_{3 1^{'}}^{†} - σ_{3 1^{'}}) + f_{3} (t)

(3d)

\frac{d ρ_{2}}{d t} = γ_{s p} ρ_{3} + (1 / τ_{32}) ρ_{3} - γ_{2} ρ_{2} + i (g σ_{23}^{†} a - g^{*} a^{†} σ_{23}) + f_{2} (t)

(3e)

Where

Δ_{3 1^{'}}

is the energy detuning between level 3 and 1′. γ_sp is the rate of spontaneous emission coupling into all the field modes except for the laser mode,

κ

is the field decay rate,

γ_{j}

is the total relaxation rate of level j,

1 / τ_{32}

is the relaxation rate from level 3 to 2,

γ_{i j} = γ_{i} / 2 + γ_{j} / 2 + T_{2}^{*}

is the polarization decay rate (dephasing rate) with a phenomenological pure dephasing rate

T_{2}^{*}

due to interface roughness and impurity scattering.

In the above equations, we neglect the polarization

σ_{2 1^{'}}

driven by higher quantum coherence effects. This assumption is, in general, valid for mid-IR designs due to a large value of

Δ_{32}

and also reasonably valid for THz QCLs due to thicker injector barriers [10

S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]

The Langevin noise operators are fully defined by their first- and second-order moments. For the field Langevin forces

f_{a} (t)

one gets [15

M. O. Scully and M. S. Zubairy, Quantum optics (Cambridge University Press, 1997).

,16

W. H. Louisell, Quantum statistical properties of radiation (Wiley, 1973).

]

〈 f_{a} (t) 〉 = 0

(4a)

〈 f_{a}^{†} (t) f_{a} (t^{'}) 〉 = κ n_{t h} δ (t - t^{'})

(4b)

〈 f_{a} (t) f_{a}^{†} (t^{'}) 〉 = κ (n_{t h} + 1) δ (t - t^{'})

(4c)

〈 f_{a}^{†} (t) f_{a}^{†} (t^{'}) 〉 = 〈 f_{a} (t) f_{a} (t^{'}) 〉 = 0

(4d)

where “+1” is the contribution from vacuum fluctuations and has no classical equivalent.

n_{t h}

is the average number of thermal photons in the laser cavity at temperature T_e (approximate electron temperature at the quasi-thermal equilibrium) and can be written as

n_{t h} = {[\exp (ℏ ω_{λ} / k_{B} T_{e}) - 1]}^{- 1}

(5)

For the electron population and polarization noise operators, the correlation functions can be derived according to the generalized Einstein relations [15

M. O. Scully and M. S. Zubairy, Quantum optics (Cambridge University Press, 1997).

,16

W. H. Louisell, Quantum statistical properties of radiation (Wiley, 1973).

], as explained in the Appendix A.

3. Laser intrinsic linewidths

3.1. Equivalent c-number Langevin equations

To solve the present problem, we have to convert the above operator equations into c-number equations [15

M. O. Scully and M. S. Zubairy, Quantum optics (Cambridge University Press, 1997).

]. For this we have to choose certain particular ordering for field and electron density operators, because the c-number variables commute with each other while the operators do not. Here we choose the normal ordering of field and electron density operators to be

a^{†}

σ_{23}^{†}

σ_{3 1^{'}}^{†}

ρ_{3}

ρ_{2}

σ_{3 1^{'}}

σ_{23}

a

, that is, the stochastic c-number variables corresponding to the operators

a

σ_{23}

σ_{3 1^{'}}

ρ_{3}

ρ_{2}

are replaced with their classical counterparts

A

p_{23}

p_{3 1^{'}}

n_{3}

n_{2}

. Then the c-number Langevin equations of Eqs. (3) can be written as

\frac{d A}{d t} = - \frac{κ}{2} A - i g^{*} p_{23} + F_{A} (t)

(6a)

\frac{d p_{23}}{d t} = - γ_{23} p_{23} + i g (n_{3} - n_{2}) A + F_{23} (t)

(6b)

\frac{d p_{3 1^{'}}}{d t} = - γ_{3 1^{'}} p_{3 1^{'}} - i \frac{Δ_{3 1^{'}}}{ℏ} p_{3 1^{'}} - \frac{i Δ_{0}}{2 ℏ} (n_{1^{'}} - n_{3}) + F_{3 1^{'}} (t)

(6c)

\frac{d n_{3}}{d t} = - γ_{s p} n_{3} - γ_{3} n_{3} + i (g^{*} A^{†} p_{23} - g p_{23}^{†} A) - \frac{i Δ_{0}}{2 ℏ} (p_{3 1^{'}}^{†} - p_{3 1^{'}}) + F_{3} (t)

(6d)

\frac{d n_{2}}{d t} = γ_{s p} n_{3} + (1 / τ_{32}) n_{3} - γ_{2} n_{2} + i (g p_{23}^{†} A - g^{*} A^{†} p_{23}) + F_{2} (t)

(6e)

Here, the Langevin forces

F_{μ}

have the following properties,

〈 F_{μ} (t) 〉 = 0

(7a)

〈 F_{u}^{†} (t) F_{v} (t^{'}) 〉 = 2 D_{u^{†} v} δ (t - t^{'})

(7b)

The diffusion coefficients

D_{u^{†} v}

of c-number Langevin forces may be different from the corresponding diffusion coefficients

d_{u^{†} v}

of the operator Langevin forces defined by Eqs. (A5). The diffusion coefficients

D_{u^{†} v}

are determined in such a way that the second moments calculated from the c-number equations agree with those calculated from the operator equations as explained in Appendix B.

3.2. Steady-state solution for above-threshold operation

The steady-state solutions for the mean values of the field and electron number variables above threshold are obtained by setting the time derivatives to be zero and dropping the noise terms in Eqs. (6). Neglecting the correlations between the electron populations and the photons, one then finds for the steady-state mean photon numbers, population inversion and population of the upper and lower laser levels

I_{0} = A_{0}^{2} = \frac{γ_{2} n_{2, 0} - γ_{s p} n_{3, 0} - (1 / τ_{32}) n_{3, 0}}{κ}

(8a)

n_{3, 0} - n_{2, 0} = \frac{γ_{23} κ}{2 g g^{*}}

(8b)

n_{3, 0} = [- \frac{γ_{2} γ_{23} κ}{2 g g^{*}} - {(\frac{Δ_{0}}{2 ℏ})}^{2} (N_{0} + \frac{γ_{23} κ}{2 g g^{*}}) \frac{2 γ_{3 1^{'}}}{{(γ_{3 1^{'}})}^{2} + {(Δ_{3 1^{'}} / ℏ)}^{2}}] / [1 / τ_{32} - γ_{3} - γ_{2} - {(\frac{Δ_{0}}{2 ℏ})}^{2} \frac{6 γ_{3 1^{'}}}{γ_{3 1^{'}}^{2} + {(Δ_{3 1^{'}} / ℏ)}^{2}}]

(8c)

n_{2, 0} = [- \frac{γ_{23} (1 / τ_{32} - γ_{3}) κ}{2 g g^{*}} + {(\frac{Δ_{0}}{2 ℏ})}^{2} (\frac{γ_{23} κ}{g g^{*}} - N_{0}) \frac{2 γ_{3 1^{'}}}{γ_{3 1^{'}}^{2} + {(Δ_{3 1^{'}} / ℏ)}^{2}}] / [1 / τ_{32} - γ_{3} - γ_{2} - {(\frac{Δ_{0}}{2 ℏ})}^{2} \frac{6 γ_{3 1^{'}}}{γ_{3 1^{'}}^{2} + {(Δ_{3 1^{'}} / ℏ)}^{2}}]

(8d)

where N₀ is the electron number of one module in the cavity.

On the other hand, the steady-state polarization can be expressed in terms of the mean value of the field as

p_{23, 0} = \frac{κ}{i 2 g} A_{0}

(9a)

p_{3 1^{'}, 0} = - \frac{i Δ_{0}}{2 ℏ} (N_{0} - n_{2, 0} - 2 n_{3, 0}) / (γ_{3 1^{'}} + i \frac{Δ_{1^{'} 3}}{ℏ})

(9b)

From Eq. (9a), we find that if

A_{0}

is real,

p_{23, 0}

is purely imaginary. Therefore, we can choose the mean value of the phase of the laser field to be zero for simplification.

The current at and above threshold can be derived as, respectively

I_{c, t h} = | e | (\frac{n_{3, t h}}{τ_{31}} + \frac{n_{2, t h}}{τ_{21}}) = | e | γ_{23} κ \frac{γ_{2} / 2 τ_{31} + (γ_{s p} + 1 / τ_{32}) / 2 τ_{21}}{(γ_{2} - γ_{s p} - 1 / τ_{32}) g g^{*}}

(10a)

I_{c} = | e | (\frac{n_{3, 0}}{τ_{31}} + \frac{n_{2, 0}}{τ_{21}}) = | e | \frac{{(\frac{Δ_{0}}{2 ℏ})}^{2} [- \frac{(τ_{31} + τ_{21}) N_{0}}{τ_{31} τ_{21}} + \frac{(2 τ_{31} - τ_{21}) γ_{23} κ}{2 τ_{31} τ_{21} g g^{*}}] \frac{2 γ_{3 1^{'}}}{{(γ_{3 1^{'}})}^{2} + {(Δ_{3 1^{'}} / ℏ)}^{2}}}{1 / τ_{32} - γ_{3} - γ_{2} - {(\frac{Δ_{0}}{2 ℏ})}^{2} \frac{6 γ_{3 1^{'}}}{γ_{3 1^{'}}^{2} + {(Δ_{3 1^{'}} / ℏ)}^{2}}}

(10b)

where

τ_{31}

and

τ_{21}

is the lifetime from the upper and lower laser level to injector level, respectively. It is noted that current contributions are made by electrons from the second phonon energy level to injector level in addition to the transition from the upper level to injector level for two-phonon resonance mid-IR QCLs.

3.3. Dynamics of fluctuations around steady state

To determine the frequency noise spectrum and linewidth we need first of all to investigate the small fluctuations of the field and electron number variables around the steady state. Neglecting terms of the second and higher order in the fluctuations, we set

A = A_{0} + δ A

(11a)

p_{23} = p_{23, 0} + δ p_{23}

(11b)

p_{3 1^{'}} = p_{3 1^{'}, 0} + δ p_{3 1^{'}}

(11c)

n_{3} = n_{3, 0} + δ n_{3}

(11d)

n_{2} = n_{2, 0} + δ n_{2}

(11e)

Thus, we get the following equations for fluctuations

\frac{d δ A}{d t} = - \frac{κ}{2} δ A - i g^{*} δ p_{23} + F_{A} (t)

(12a)

\frac{d δ p_{23}}{d t} = - γ_{23} δ p_{23} + i g (δ n_{3} - δ n_{2}) A_{0} + i g (n_{3, 0} - n_{2, 0}) δ A + F_{23} (t)

(12b)

\frac{d δ p_{3 1^{'}}}{d t} = - γ_{3 1^{'}} δ p_{3 1^{'}} - i \frac{Δ_{3 1^{'}}}{ℏ} δ p_{3 1^{'}} - \frac{i Δ_{0}}{2 ℏ} (δ n_{1^{'}} - δ n_{3}) + F_{3 1^{'}} (t)

(12c)

\begin{array}{l} \frac{d δ n_{3}}{d t} = - γ_{s p} δ n_{3} - γ_{3} δ n_{3} + i (g^{*} δ A^{†} p_{23, 0} - g p_{23, 0} δ A) + i (g^{*} A_{0} δ p_{23} - g δ p_{23}^{†} A_{0}) \\ - \frac{i Δ_{0}}{2 ℏ} (δ p_{3 1^{'}}^{†} - δ p_{3 1^{'}}) + F_{3} (t) \end{array}

(12e)

\frac{d δ n_{2}}{d t} = γ_{s p} δ n_{3} + (1 / τ_{32}) δ n_{3} - γ_{2} δ n_{2} + i (g p_{23, 0} δ A - g^{*} δ A^{†} p_{23, 0}) + i (g δ p_{23}^{†} A_{0} - g^{*} A_{0} δ p_{23}) + F_{2} (t)

(12f)

These equations can now be solved exactly by taking the Fourier transform of

δ n_{3} (ω)

δ n_{3} (ω) = \int_{- \infty}^{+ \infty} e^{- i ω t} δ n_{3} (t) d t

(13)

as well as all the other variables. Then, one gets

i ω δ A (ω) = - (κ / 2) δ A (ω) - i g^{*} δ p_{23} (ω) + F_{A} (ω)

(14a)

i ω δ p_{23} (ω) = - γ_{23} δ p_{23} (ω) + i g [δ n_{3} (ω) - δ n_{2} (ω)] A_{0} + i g (n_{3, 0} - n_{2, 0}) δ A (ω) + F_{23} (ω)

(14b)

i ω δ p_{3 1^{'}} (ω) = - γ_{3 1^{'}} δ p_{3 1^{'}} (ω) - i \frac{Δ_{1^{'} 3}}{ℏ} δ p_{3 1^{'}} (ω) - \frac{i Δ_{0}}{2 ℏ} [δ n_{1^{'}} (ω) - δ n_{3} (ω)] + F_{3 1^{'}} (ω)

(14c)

\begin{array}{l} i ω δ n_{3} (ω) = - γ_{s p} δ n_{3} (ω) - γ_{3} δ n_{3} (ω) + i [g^{*} δ A^{†} (- ω) p_{23, 0} - g p_{23, 0} δ A (ω)] \\ + i [g^{*} A_{0} δ p_{23} (ω) - g δ p_{23}^{†} (- ω) A_{0}] - \frac{i Δ_{0}}{2 ℏ} [δ p_{3 1^{'}}^{†} (- ω) - δ p_{3 1^{'}} (ω)] + F_{3} (ω) \end{array}

(14d)

\begin{array}{l} i ω δ n_{2} (ω) = γ_{s p} δ n_{3} (ω) + (1 / τ_{32}) δ n_{3} (ω) - γ_{2} δ n_{2} (ω) \\ + i [g p_{23, 0} δ A (ω) - g^{*} δ A^{†} (- ω) p_{23, 0}] + i [g δ p_{23}^{†} (- ω) A_{0} - g^{*} A_{0} δ p_{23} (ω)] + F_{2} (ω) \end{array}

(14e)

where the frequency-dependent fluctuation forces satisfy

〈 F_{u} (ω) F_{v} (ω^{'}) 〉 = 4 π D_{u v} δ (ω + ω^{'})

(15)

The solution of this linear system can then be easily obtained.

3.4. Noise spectra

To calculate noise spectra, we use the semiclassical expression of the laser field [13

L. Davidovich, “Sub-Poissonian processes in quantum optics,” Rev. Mod. Phys. 68(1), 127–173 (1996). [CrossRef]

]

A = \sqrt{I} \exp (i φ)

(16)

Then the dynamical equations of phase fluctuations are obtained as

δ φ (t) = \frac{i}{2 \sqrt{I_{0}}} [δ A^{†} (t) - δ A (t)]

(17)

By taking Fourier transform of Eq. (17) on both sides, one gets

δ φ (ω) = \frac{i}{2 \sqrt{I_{0}}} [δ A^{†} (- ω) - δ A (ω)]

(18)

According to Eq. (14a), the Fourier amplitudes of the field fluctuations can be expressed as

δ A (ω) = \frac{i g δ p_{23} (ω) + F_{A} (ω)}{i ω + κ / 2}

(19)

where

δ p_{23} (ω) = \frac{i g [δ n_{3} (ω) - δ n_{2} (ω)] A_{0} + i g (n_{3, 0} - n_{2, 0}) \frac{F_{A} (ω)}{i ω + κ / 2} + F_{23} (ω)}{i ω + γ_{23} - \frac{g g^{*} (n_{3, 0} - n_{2, 0})}{i ω + κ / 2}}

(20)

Therefore, with applying Eq. (19) and Eq. (20) into Eq. (18), the phase fluctuations can be explicitly written as

δ φ (ω) = \frac{i}{2 \sqrt{I_{0}}} \frac{(i ω + γ_{23}) [F_{A}^{†} (- ω) - F_{A} (ω)] + i g [F_{23}^{†} (- ω) - F_{23} (ω)]}{(i ω + γ_{23}) (i ω + κ / 2) - g g^{*} (n_{3, 0} - n_{2, 0})}

(21)

The autocorrelation function of the phase fluctuations can be expressed as

〈 δ φ (ω) δ φ^{*} (ω^{'}) 〉 = {(δ φ)}_{ω}^{2} δ (ω - ω^{'})

(22)

According to the Wiener-Khinchine theorem [17

K. Petermann, Laser diode modulation and noise (Kluwar Academic, 1988).

], the relations between autocorrelation function and power spectral density satisfy

〈 δ φ (ω) δ φ^{*} (ω^{'}) 〉 = 2 π S_{φ} δ (ω - ω^{'})

(23)

where the power spectral density is defined as the following Fourier transform

S_{φ} = \int_{- \infty}^{+ \infty} 〈 δ φ (t) δ φ (t - τ) 〉 \exp (- i ω τ) d τ

(24)

The frequency noise spectral density S_f, therefore, can be derived as

S_{f} = ω^{2} S_{φ} = \frac{κ (ω^{2} + γ_{23}^{2}) (2 n_{t h} + 1) + 2 g g^{*} γ_{23} (n_{3, 0} + n_{2, 0} - 2 n_{3, 0} n_{2, 0} / N_{0})}{4 I_{0} [ω^{2} + {(γ_{23} + κ / 2)}^{2}]}

(25)

Once we know the power spectrum density of the frequency noise, the linewidth can then be obtained by the following derivation. The autocorrelation function of the lasing field

A (t)

can be explicitly derived for single-mode QCLs when the phase fluctuations

δ φ (t)

are considered as a Gaussian distribution [17

K. Petermann, Laser diode modulation and noise (Kluwar Academic, 1988).

〈 A (t) A^{*} (t - τ) 〉 = I_{0} \exp (i ω_{λ} t) \exp (- 〈 δ φ^{2} 〉 / 2)

(26)

where we have neglected the amplitude fluctuations. The phase change induced by the amplitude change can be characterized by the α-parameter (linewidth enhancement factor) [18

C. Henry, “Theory of the linewidth of semiconductor lasers,” IEEE J. Quantum Electron. 18(2), 259–264 (1982). [CrossRef]

]. The overall linewidth should be multiplied by 1 + α² if this factor is considered. This parameter is not included in our model since the α-value is negligible in QCLs, which has been confirmed by experiments [19

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]

In general, if the spectral density of the frequency fluctuation is known, the mean squared value of the phase change can be obtained as [20

B. Daino, P. Spano, M. Tamburrini, and S. Piazzolla, “Phase noise and spectral line shape in semiconductor lasers,” IEEE J. Quantum Electron. 19(3), 266–270 (1983). [CrossRef]

]

〈 δ φ^{2} 〉 = \frac{2}{π} \int_{0}^{\infty} \frac{S_{f} (ω)}{ω^{2}} [1 - \cos (ω t)] d ω

(27)

namely

\begin{array}{l} 〈 δ φ^{2} 〉 = \frac{2 κ n_{t h} + κ + 2 (g g^{*} / γ_{23}) (n_{3, 0} + n_{2, 0} - 2 n_{3, 0} n_{2, 0} / N_{0})}{4 {(1 + κ / 2 γ_{23})}^{2} I_{0}} t \\ + \frac{- 2 κ^{2} n_{t h} / γ_{23} - κ^{2} / γ_{23} + 2 (g g^{*} / γ_{23}) (n_{3, 0} + n_{2, 0} - 2 n_{3, 0} n_{2, 0} / N_{0})}{4 {(1 + κ / 2 γ_{23})}^{2} I_{0}} [\frac{e^{- (γ_{23} + κ / 2) | t |} - 1}{γ_{23} + κ / 2}] \end{array}

(28)

The linewidth is then obtained from the following spectral density by taking the Fourier transform of Eq. (26)

S_{A} (ω) = I_{0} \int_{- \infty}^{+ \infty} \exp [- i (ω - ω_{λ}) t] \exp (- 〈 δ φ^{2} 〉 / 2) d τ

(29)

4. Results and discussions

As shown in Eq. (28), fluctuations of laser field can be attributed to three sources, i.e. thermal photons (2κn_th), vacuum fluctuations (κ), and spontaneous emission processes (2gg*/γ₂₃). These sources induce linewidth broadening in lasers, which cannot be overcome due to fundamental quantum limitations. It needs to be mentioned that the vacuum fluctuations as one of noise sources are not included in rate equation model. Although the fact that vacuum fluctuations cannot be detected directly, the interaction of this vacuum field and the laser field can result in a modulation of the photon flux, which causes noise in the cavity [21

C. H. Henry and R. F. Kazarinov, “Quantum noise in photonics,” Rev. Mod. Phys. 68(3), 801–853 (1996). [CrossRef]

]. As a result, our model based on quantum Langevin equations provides more information about the fundamental physical origins of the laser linewidth.

In contrast to the traditional average squared phase fluctuation expression (such as Eq. (B10) in Ref. 9

), there is an additional term (the second line in Eq. (28)). This term can finally determine the intrinsic linewidth, which depends on the time of measurement. If the time of measurement is shorter than the characteristic time τ = (γ₂₃ + κ/2)⁻¹, i.e. using ultrafast photodetectors, the phase fluctuation caused by spontaneous emission will disappear. This is attributed to suppression of spontaneous emission quantum noise caused by the electron memory effect associated with the transient behavior of the polarization. Therefore, for short-time measurements, the laser linewidth will reduce due to the disappearance of the spontaneous emission contribution. This property is not shown in the rate equation model. Since the rate equation model cannot include the dynamics of coherent interaction, the linewidth of the laser can only be based on the long-time measurement. Owing to the limit of bandwidth of the photodetector, we will not consider the cases of short-time measurements, and only investigate the laser linewidth of a long-time measurement in the following discussion.

For a long-time measurement, the exponential term becomes zero, and the phase fluctuation increases linearly with increasing time. The linewidth of the laser can then be expressed as

δ f = \frac{κ (2 n_{t h} + 1) + 2 (g g^{*} / γ_{23}) N_{0} n_{3, 0}}{8 π {(1 + κ / 2 γ_{23})}^{2} I_{0}}

(30)

According to above Eq. (30) and Eq. (B12) in Ref. 9

, we see that the linewidth deduced from our quantum mechanical Langevin model is different from the expression based on the rate equation model. The explicit difference between these two formulas is caused by the coefficient (1 + κ/2γ₂₃)². In the rate equation model,

γ_{23} ≫ κ

is assumed, and the polarization of the active medium is adiabatically eliminated. Hence, the information on dynamics of coherent interaction is missed. The inclusion of the coherent interaction can induce a smaller linewidth by a factor of (1 + κ/2γ₂₃)² if the coherent time is comparable to the cavity loss according to Eq. (30). The intrinsic difference between these two models is from the resonant-tunneling effects. Electrons in the ground injector level are injected into the upper laser level in QCLs through resonant tunneling. Since coherence plays an important role in the resonant-tunneling mechanism, which can significantly influence the electron transport [10

S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]

–12

H. Callebaut and Q. Hu, “Importance of coherence for electron transport in terahertz quantum cascade lasers,” J. Appl. Phys. 98(10), 104505 (2005). [CrossRef]

], it is necessary to include this factor in laser dynamics simulations of QCLs. Comparing these two formulas, the rate equation model overestimates the influence of thermal photon, which lead to a higher intrinsic linewidth calculation for THz QCLs during high temperature operation.

Figure 2 shows the comparisons of the laser linewidths derived from our quantum mechanical model and the classical rate equation model at different operation currents in both THz (Fig. 2(a)) and Mid-IR (Fig. 2(b)) QCLs. Comparing the two models, the apparent difference for THz QCLs and relatively smaller difference for Mid-IR QCLs are caused by including the resonant-tunneling transport effect in our model. These differences increase as the operation current increases. Because the effects of resonant tunneling increases when the injector level and upper laser level tend to be aligned as the increasing operation current. It is noted that the effects of dynamics of coherent interaction associated with laser transitions on the linewidth of QCLs are not significant according to our calculation with the present parameters.

Fig. 2 The comparisons of linewidths derived from our model and the classical rate equation model at different operation currents (a) in THz QCLs; The default parameters are chosen from the typical values for resonant-tunneling injection THz QCLs designs (Refs. 10

S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]

,23

S. Kumar, Q. Hu, and J. L. Reno, “186 K operation of terahertz quantum-cascade lasers based on a diagonal design,” Appl. Phys. Lett. 94(13), 131105 (2009). [CrossRef]

): γ₂₃ = 4.5 × 10¹² s⁻¹, γ_31′ = 2.5 × 10¹² s⁻¹, κ = 2 × 10¹¹ s⁻¹, 1/τ₃₂ = 3 × 10¹¹ s⁻¹, γ₃ = 5 × 10¹¹ s⁻¹, γ₂ = 2 × 10¹² s⁻¹, ε_31′ = 0, ∆₀ = 1.5 meV, μ = 3.7 nm, Γ ≈1, n = 3.6, T = 20 K, the doping sheet density is 3 × 10¹⁰ cm⁻². (b) in Mid-IR QCLs; The default parameters are chosen from two phonon resonance mid-IR QCLs designs (Ref. 24

D. Hofstetter, M. Beck, T. Aellen, and J. Faist, “High-temperature operation of distributed feedback quantum-cascade lasers at 5.3 μm,” Appl. Phys. Lett. 78(4), 396 (2001). [CrossRef]

): γ₂₃ = 20 × 10¹² s⁻¹, κ = 2 × 10¹¹ s⁻¹, 1/τ₃₂ = 3 × 10¹¹ s⁻¹, γ₃ = 8 × 10¹¹ s⁻¹, γ₂ = 4 × 10¹² s⁻¹, γ_31′ = 15 × 10¹² s⁻¹, ε_31′ = 0, ∆₀ = 8 meV, μ = 2.1 nm, Γ ≈0.5, n = 3.2, T = 300 K, the doping sheet density is 2 × 10¹¹ cm⁻².

The role of coherence in resonant-tunneling transport in QCLs has been investigated [10

S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]

–12

H. Callebaut and Q. Hu, “Importance of coherence for electron transport in terahertz quantum cascade lasers,” J. Appl. Phys. 98(10), 104505 (2005). [CrossRef]

]. It has been suggested that the inclusion of coherent transport and dephasing in calculations is essential especially when transport is dominated by transitions between weakly coupled states. We take THz QCLs structures [10

S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]

,12

H. Callebaut and Q. Hu, “Importance of coherence for electron transport in terahertz quantum cascade lasers,” J. Appl. Phys. 98(10), 104505 (2005). [CrossRef]

] based on resonant-tunneling injection to investigate how coupling strength and dephasing rate associated with resonant tunneling influence the intrinsic linewidth. Figure 3 shows the effects of coupling strength on the linewidth at resonance. As shown in Fig. 3(a), the linewidth decreases as the injection coupling strength increases. When the injection coupling strength is below 3 meV, the linewidth diminishes rapidly. But once the coupling strength exceeds 3 meV, there only exists slight changes. This can be attributed to the coherence of electrons transport across the injector barrier, and this coherence can be determined by the factor

{(Δ_{0} / ℏ)}^{2} {(γ_{3 1^{'}} γ_{3})}^{- 1}

(

≫ 1

coherent) and (

≪ 1

incoherent) [10

S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]

]. The small injection coupling strength corresponds to the regime of incoherent resonant-tunneling transport. But when the coupling becomes strong, the resonant-tunneling transport tends to be coherent. As a result, more excited electrons are injected into upper laser level, and then stimulated to the lower laser level by emitting photons, current density and photon number increases, as shown in Fig. 3(b). Hence, the noise associated with spontaneous emission is strongly suppressed, leading to a major reduction in the linewidth. Therefore, THz QCLs will have a much larger linewidth for structures with largely incoherent resonant-tunneling transport, but a smaller linewidth with coherent resonant-tunneling transport. Since the injector coupling strength of THz QCLs based on resonant-tunneling injection is typically 1~2 meV, the resonant tunneling can strongly influence the linewidth of THz QCLs.

Fig. 3 The effects of coupling strength on the linewidth of THz QCLs at resonance. (a) Linewidth as a function of coupling strength. The linewidth decreases as the coupling strength between the injector level and the upper laser level increases. The factor

{(Δ_{0} / ℏ)}^{2} {(γ_{3 1^{'}} γ_{3})}^{- 1}

is used to determine the transition of resonant tunneling from coherent (

≫ 1

) to incoherent (

≪ 1

) (Ref. 10

S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]

). Coupling strengths more than 3 meV lay in the coherent regime. (b) Current density and photon number as a function of coupling strength.

Since the resonant-tunneling transport can be strongly influenced by the scattering, the linewidth has a strong dependence on the damping of the coherent interaction. Figure 4(a) shows the effects of dephasing rate associated with resonant-tunneling transport on the linewidth. When the dephasing rate γ_31′ increases, only a small fraction of the electrons tunnels through the injector barrier into the upper laser level. As a result, the current density and photon number decreases (Fig. (4b)) and the noise associated with spontaneous emission becomes strong. Therefore, the linewidth increases as the dephasing rate increases.

Fig. 4 The effects of dephasing rate associated with resonant tunneling on the linewidth of THz QCLs at resonance. (a) Linewidth as a function of dephasing rate. The linewidth increases as the dephasing rate increases. (b) Current density and photon number as a function of dephasing rate.

As demonstrated in the above analysis, the characteristics of resonant-tunneling transport can strongly influence the linewidth of THz QCLs in the incoherent resonant-tunneling transport regime, but only have small effect on the linewidth in the coherent resonant-tunneling transport regime. For the mid-IR QCLs, the injector barrier is usually designed to be thinner, and the electrons are injected into the upper laser level by coherent resonant-tunneling processes. We can deduce that the resonant-tunneling transport can exhibit negligible effect on the linewidth of mid-IR QCLs. Figures 5 and 6 show the effects of coupling strength and dephasing rate associated with resonant tunneling on the linewidth. Owing to the much larger coupling strength of mid-IR QCLs, the electron can be injected into the upper laser level by coherent resonant tunneling. As a result, further increasing the coupling strength and dephasing rate γ_31′ can only induce few changes of photon numbers, and hence only introduce little influence on the linewidth of mid-IR QCLs.

Fig. 5 The effects of coupling strength on the linewidth of mid-IR QCLs at resonance. (a) Linewidth as a function of coupling strength. The linewidth slowly decreases as the coupling strength of the injector level and the upper laser level increases. (b) Current density and photon number as a function of coupling strength.

Fig. 6 The effects of dephasing rate associated with resonant tunneling on the linewidth of mid-IR QCLs at resonance. (a) Linewidth as a function of dephasing rate. The linewidth slowly increases as the dephasing rate increases. (b) Current density and photon number as a function of dephasing rate.

In the rate equation model, the operation current is taken as an input parameter. However, many other parameters e.g. coupling strength, dephasing rate and doping level can influence the current. Therefore quantum design with the optimization of the active region of a QCL to reduce the intrinsic linewidth is limited by using rate equation model, but it can be easily made according to our quantum mechanical Langevin model. To take the effect of doping on the linewidth as an example, we can optimize the doping density to achieve a reduced linewidth. Increasing doping density results in an increase in free-carrier absorption and waveguide loss. According to Ref. 22

A. M. Andrews, A. Benz, C. Deutsch, G. Fasching, K. Unterrainer, P. Klang, W. Schrenk, and G. Strasser, “Doping dependence of LO-phonon depletion scheme THz quantum-cascade lasers,” Mater. Sci. Eng. B 147(2-3), 152–155 (2008). [CrossRef]

, the loss can be assumed to be roughly proportional to the doping density. Doping not only affects the cavity loss (thus the intrinsic linewidth) but also the electron population distribution and photon number according to our model (we neglect the effects of doping density on relaxation rates in our chosen ranges of doping). For the rate equation model, doping affects the linewidth only by the cavity loss. Figure 7(a) shows the self-induced linewidth variation by doping, the linewidth variation due to the change of doping induced cavity loss, and the overall effects from these two factors. The linewidth is enhanced as the doping density increases. Doping density does not cause the change of photon number but induce the increasing of current density in this calculation, as shown in Fig. 7(b). This shows that more electrons populate the upper laser level and take part in the spontaneous emission, then induce the linewidth broadening. The simulation tells us that we can optimize the doping density to reduce the linewidth of QCLs. It is noted that the optimization of the injector barrier is complex, since the thickness of barrier not only influences the coupling strength, but also it affects the lifetimes of the upper and lower laser level. Figure 8 shows the effects of lifetimes of the upper and laser levels on the linewidth. The linewidth increases and decreases as increasing the relaxation rate of upper laser and lower laser level, respectively. This is due to the variation of ratio of current to threshold current and photon number caused by the change of their relaxation rates, as shown in Fig. (9) . Thus, a careful design of active region should be made to reduce the intrinsic linewidth of QCLs.

Fig. 7 Effect of doping density on the intrinsic linewidth of THz QCLs at resonance. (a) The self-induced linewidth variation by doping, the linewidth variation due to the change of doping induced cavity loss, and the overall effects from these two factors. (b) Current density and photon number as a function of doping density.

Fig. 8 (a) Linewidth as a function of relaxation rate of the upper laser level. (b) Linewidth as a function of the relaxation rate of the lower laser level

Fig. 9 (a) Ratio of current to threshold current and photon number as a function of relaxation rate of the upper laser level. (b) Ratio of current to threshold current and photon number as a function of relaxation rate of the lower laser level.

5. Conclusions

In conclusion, we have developed a new model for the calculation of intrinsic linewidth of QCLs based on quantum Langevin equations. It includes the effects of resonant-tunneling transport and dynamics of coherent interaction. We conclude that the linewidth of the laser field can be attributed to three sources, e.g. thermal photons, vacuum fluctuations and spontaneous emission processes. The results show that for the short time measurement, the effects of electron memory can lead to suppression of spontaneous emission quantum noise; and the intrinsic linewidth of QCLs are reduced with the consideration of dynamics of coherent interaction. We also demonstrate that the coupling strength and dephasing rate have significant effects on the linewidth of THz QCLs in the incoherent resonant-tunneling regime, but small effects on that of mid-IR QCLs due to their strong coherent resonant-tunneling. The linewidth decreases with the increase of the injection coupling strength and reduction of the dephasing rate associated with resonant tunneling. According to our model, a reduced intrinsic linewidth can be easily designed through optimization of the injector barrier and doping density of the active region of QCLs.

Appendices

Appendix A: operator diffusion coefficients

For the dynamical variation

o_{u}

, we have the following quantum Langevin equations

\frac{d o_{u}}{d t} = d_{u} + f_{u} (t)

(A1)

where

d_{u}

is the drift coefficients and

F_{u}

is assumed to have zero average and to be delta-correlated Langevin forces, i.e.

〈 f_{u} (t) 〉 = 0

(A2a)

〈 f_{u^{†}} (t) f_{v} (t^{'}) 〉 = 2 〈 d_{u^{†} v} 〉 δ (t - t^{'})

(A2b)

The associated diffusion coefficients

d_{u v}

can be calculated using the generalized Einstein relations [15

M. O. Scully and M. S. Zubairy, Quantum optics (Cambridge University Press, 1997).

,16

W. H. Louisell, Quantum statistical properties of radiation (Wiley, 1973).

]

2 〈 d_{u^{†} v} 〉 = \frac{d}{d t} 〈 o_{u}^{†} o_{v} 〉 - 〈 d_{u}^{†} o_{v} 〉 - 〈 o_{u}^{†} d_{v} 〉

(A3)

For example, the diffusion coefficient corresponding to the fluctuation force

f_{3} (t)

〈 f_{3} (t) f_{3} (t^{'}) 〉 = 2 d_{33} δ (t - t^{'})

(A4)

can be obtained from Eq. (A3)

\begin{array}{l} 〈 f_{3} (t) f_{3} (t^{'}) 〉 = 〈 \frac{d}{d t} {(ρ_{3})}^{2} - (\frac{d}{d t} ρ_{3}) ρ_{3} - ρ_{3} (\frac{d}{d t} ρ_{3}) 〉 δ (t - t^{'}) \\ = 〈 - γ_{s p} ρ_{3} - γ_{3} ρ_{3} + i (g^{*} a^{†} σ_{23} - g σ_{23}^{†} a) - \frac{i Δ_{0}}{2 ℏ} (σ_{3 1^{'}}^{†} - σ_{3 1^{'}}) 〉 \\ - 〈 - γ_{s p} ρ_{3} ρ_{3} - γ_{3} ρ_{3} ρ_{3} + i (g^{*} a^{†} σ_{23} - g σ_{23}^{†} a) ρ_{3} - \frac{i Δ_{0}}{2 ℏ} (σ_{3 1^{'}}^{†} - σ_{3 1^{'}}) ρ_{3} 〉 \\ - 〈 - γ_{s p} ρ_{3} ρ_{3} - γ_{3} ρ_{3} ρ_{3} + i ρ_{3} (g^{*} a^{†} σ_{23} - g σ_{23}^{†} a) - \frac{i Δ_{0}}{2 ℏ} n_{3} (σ_{3 1^{'}}^{†} - σ_{3 1^{'}}) 〉 \\ = 〈 (γ_{s p} + γ_{3}) ρ_{3} 〉 δ (t - t^{'}) \end{array}

(A5a)

Applying the same procedure to other Langevin forces, one can obtain

〈 f_{2} (t) f_{2} (t^{'}) 〉 = 〈 (γ_{s p} + 1 / τ_{32}) ρ_{3} + γ_{2} ρ_{2} 〉 δ (t - t^{'})

(A5b)

〈 f_{23}^{†} (t) f_{23} (t^{'}) 〉 = 〈 (2 γ_{23} - γ_{s p} - 1 / τ_{32}) ρ_{3} + γ_{2} ρ_{2} 〉 δ (t - t^{'})

(A5c)

〈 f_{23} (t) f_{23}^{†} (t^{'}) 〉 = 〈 (γ_{s p} + 1 / τ_{32}) ρ_{3} + (2 γ_{23} - γ_{2}) ρ_{2} 〉 δ (t - t^{'})

(A5d)

〈 f_{23} (t) f_{23} (t^{'}) 〉 = 0

(A5e)

〈 f_{3 1^{'}}^{†} (t) f_{3 1^{'}} (t^{'}) 〉 = 〈 γ_{2} n_{2} + (γ_{3} - 1 / τ_{32}) ρ_{3} + 2 γ_{3 1^{'}} ρ_{1^{'}} 〉 δ (t - t^{'})

(A5f)

〈 f_{3 1^{'}} (t) f_{3 1^{'}}^{†} (t^{'}) 〉 = 〈 (1 / τ_{32} - γ_{3} + 2 γ_{3 1^{'}}) ρ_{3} - γ_{2} ρ_{2} 〉 δ (t - t^{'})

(A5g)

Appendix B: c-number diffusion coefficients

For c-number

O_{u}

corresponding to the operator

o_{u}

, we have the following c-number quantum Langevin equations

\frac{d O_{u}}{d t} = D_{u} + F_{u} (t)

(B1)

The diffusion coefficients

D_{u v}

can be calculated using the generalized Einstein relations

2 〈 D_{u^{†} v} 〉 = \frac{d}{d t} 〈 O_{u}^{†} O_{v} 〉 - 〈 D_{u}^{†} O_{v} 〉 - 〈 O_{u}^{†} D_{v} 〉

(B2)

o_{u}^{†} o_{v}

is normally ordered, its expectation value is equal to the expectation value of the corresponding c-number product [15

M. O. Scully and M. S. Zubairy, Quantum optics (Cambridge University Press, 1997).

]. Therefore, we have

\frac{d}{d t} 〈 o_{u}^{†} o_{v} 〉 = \frac{d}{d t} 〈 O_{u}^{†} O_{v} 〉

(B3)

According to Eq. (B3) and (A3), we find that

〈 F_{u^{†}} F_{v} 〉 = 〈 f_{u^{†}} f_{v} 〉 + 〈 d_{u^{†}} o_{v} 〉 + 〈 o_{u^{†}} d_{v} 〉 - 〈 D_{u^{†}} O_{v} 〉 - 〈 O_{u^{†}} D_{v} 〉

(B4)

o_{u}^{†} o_{v}

is not normally ordered, we can use the commentator relationship of operator

o_{u}^{†} o_{v} = o_{v} o_{u}^{†} - [o_{v}, o_{u}^{†}]

to bring it into the chosen order.

Then the c-number second moments have

〈 F_{A} (t) F_{A} (t^{'}) 〉 = 〈 f_{a} (t) f_{a} (t^{'}) 〉

(B5a)

〈 F_{23}^{†} (t) F_{23} (t^{'}) 〉 = 〈 f_{23}^{†} (t) f_{23} (t^{'}) 〉

(B5b)

〈 F_{23} (t) F_{23} (t^{'}) 〉 = 〈 f_{23} (t) f_{23} (t^{'}) 〉 + 2 i g 〈 p_{23} A 〉 δ (t - t^{'})

(B5c)

Acknowledgments

We would like to acknowledge financial support from Nanyang Technological University (NTU) (M58040017), Defense 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.	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]
2.	T. Liu and Q. J. Wang, “Fundamental frequency noise and linewidth broadening caused by intrinsic temperature fluctuations in quantum cascade lasers,” Phys. Rev. B 84(12), 125322 (2011). [CrossRef]
3.	R. F. Curl, F. Capsso, C. Gmachl, A. A. Kosterev, B. McManus, R. Lewicki, M. Pusharsky, G. Wysocki, and F. K. Tittel, “Quantum cascade lasers in chemical physics,” Chem. Phys. Lett. 487(1-3), 1–18 (2010). [CrossRef]
4.	M. S. Vitiello, L. Consolino, S. Bartalini, A. Tredicucci, M. Inguscio, and P. De Natale, “The intrinsic linewidth of THz quantum cascade lasers,” CLEO: Science and Innovations, (Optical Society of America, 2012), paper CTu2B.2.
5.	C. Jirauschek, “Monte Carlo study of intrinsic linewidths in terahertz quantum cascade lasers,” Opt. Express 18(25), 25922–25927 (2010). [CrossRef] [PubMed]
6.	S. Bartalini, S. Borri, I. Galli, G. Giusfredi, D. Mazzotti, T. Edamura, N. Akikusa, M. Yamanishi, and P. De Natale, “Measuring frequency noise and intrinsic linewidth of a room-temperature DFB quantum cascade laser,” Opt. Express 19(19), 17996–18003 (2011). [CrossRef] [PubMed]
7.	L. Tombez, S. Schilt, J. Di Francesco, P. Thomann, and D. Hofstetter, “Temperature dependence of the frequency noise in a mid-IR DFB quantum cascade laser from cryogenic to room temperature,” Opt. Express 20(7), 6851–6859 (2012). [CrossRef] [PubMed]
8.	S. Bartalini, S. Borri, P. Cancio, A. Castrillo, I. Galli, G. Giusfredi, D. Mazzotti, L. Gianfrani, and P. De Natale, “Observing the intrinsic linewidth of a quantum-cascade laser: beyond the Schawlow-Townes limit,” Phys. Rev. Lett. 104(8), 083904 (2010). [CrossRef] [PubMed]
9.	M. Yamanishi, T. Edamura, K. Fujita, N. Akikusa, and H. Kan, “Theory of the intrinsic linewidth of quantum-cascade lasers: hidden reason for the narrow linewidth and line-broadening by thermal photons,” IEEE J. Quantum Electron. 44(1), 12–29 (2008). [CrossRef]
10.	S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]
11.	C. Sirtori, F. Capasso, J. Faist, A. L. Hutchinson, D. L. Sivco, and A. Y. Cho, “Resonant tunneling in quantum cascade lasers,” IEEE J. Quantum Electron. 34(9), 1722–1729 (1998). [CrossRef]
12.	H. Callebaut and Q. Hu, “Importance of coherence for electron transport in terahertz quantum cascade lasers,” J. Appl. Phys. 98(10), 104505 (2005). [CrossRef]
13.	L. Davidovich, “Sub-Poissonian processes in quantum optics,” Rev. Mod. Phys. 68(1), 127–173 (1996). [CrossRef]
14.	G. P. Agrawal and C. M. Bowden, “Concept of linewidth enhancement factor in semiconductor lasers: its usefulness and limitations,” IEEE Photon. Technol. Lett. 5(6), 640–642 (1993). [CrossRef]
15.	M. O. Scully and M. S. Zubairy, Quantum optics (Cambridge University Press, 1997).
16.	W. H. Louisell, Quantum statistical properties of radiation (Wiley, 1973).
17.	K. Petermann, Laser diode modulation and noise (Kluwar Academic, 1988).
18.	C. Henry, “Theory of the linewidth of semiconductor lasers,” IEEE J. Quantum Electron. 18(2), 259–264 (1982). [CrossRef]
19.	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]
20.	B. Daino, P. Spano, M. Tamburrini, and S. Piazzolla, “Phase noise and spectral line shape in semiconductor lasers,” IEEE J. Quantum Electron. 19(3), 266–270 (1983). [CrossRef]
21.	C. H. Henry and R. F. Kazarinov, “Quantum noise in photonics,” Rev. Mod. Phys. 68(3), 801–853 (1996). [CrossRef]
22.	A. M. Andrews, A. Benz, C. Deutsch, G. Fasching, K. Unterrainer, P. Klang, W. Schrenk, and G. Strasser, “Doping dependence of LO-phonon depletion scheme THz quantum-cascade lasers,” Mater. Sci. Eng. B 147(2-3), 152–155 (2008). [CrossRef]
23.	S. Kumar, Q. Hu, and J. L. Reno, “186 K operation of terahertz quantum-cascade lasers based on a diagonal design,” Appl. Phys. Lett. 94(13), 131105 (2009). [CrossRef]
24.	D. Hofstetter, M. Beck, T. Aellen, and J. Faist, “High-temperature operation of distributed feedback quantum-cascade lasers at 5.3 μm,” Appl. Phys. Lett. 78(4), 396 (2001). [CrossRef]

OCIS Codes
(140.3070) Lasers and laser optics : Infrared and far-infrared lasers
(270.2500) Quantum optics : Fluctuations, relaxations, and noise
(300.3700) Spectroscopy : Linewidth
(140.5965) Lasers and laser optics : Semiconductor lasers, quantum cascade

ToC Category:
Lasers and Laser Optics

History
Original Manuscript: May 23, 2012
Revised Manuscript: June 30, 2012
Manuscript Accepted: July 9, 2012
Published: July 12, 2012

Citation
Tao Liu, Kenneth E. Lee, and Qi Jie Wang, "Effects of resonant tunneling and dynamics of coherent interaction on intrinsic linewidth of quantum cascade lasers," Opt. Express 20, 17145-17159 (2012)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-20-15-17145

Sort: Author | Year | Journal | Reset

References

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]
T. Liu and Q. J. Wang, “Fundamental frequency noise and linewidth broadening caused by intrinsic temperature fluctuations in quantum cascade lasers,” Phys. Rev. B 84(12), 125322 (2011). [CrossRef]
R. F. Curl, F. Capsso, C. Gmachl, A. A. Kosterev, B. McManus, R. Lewicki, M. Pusharsky, G. Wysocki, and F. K. Tittel, “Quantum cascade lasers in chemical physics,” Chem. Phys. Lett. 487(1-3), 1–18 (2010). [CrossRef]
M. S. Vitiello, L. Consolino, S. Bartalini, A. Tredicucci, M. Inguscio, and P. De Natale, “The intrinsic linewidth of THz quantum cascade lasers,” CLEO: Science and Innovations, (Optical Society of America, 2012), paper CTu2B.2.
C. Jirauschek, “Monte Carlo study of intrinsic linewidths in terahertz quantum cascade lasers,” Opt. Express 18(25), 25922–25927 (2010). [CrossRef] [PubMed]
S. Bartalini, S. Borri, I. Galli, G. Giusfredi, D. Mazzotti, T. Edamura, N. Akikusa, M. Yamanishi, and P. De Natale, “Measuring frequency noise and intrinsic linewidth of a room-temperature DFB quantum cascade laser,” Opt. Express 19(19), 17996–18003 (2011). [CrossRef] [PubMed]
L. Tombez, S. Schilt, J. Di Francesco, P. Thomann, and D. Hofstetter, “Temperature dependence of the frequency noise in a mid-IR DFB quantum cascade laser from cryogenic to room temperature,” Opt. Express 20(7), 6851–6859 (2012). [CrossRef] [PubMed]
S. Bartalini, S. Borri, P. Cancio, A. Castrillo, I. Galli, G. Giusfredi, D. Mazzotti, L. Gianfrani, and P. De Natale, “Observing the intrinsic linewidth of a quantum-cascade laser: beyond the Schawlow-Townes limit,” Phys. Rev. Lett. 104(8), 083904 (2010). [CrossRef] [PubMed]
M. Yamanishi, T. Edamura, K. Fujita, N. Akikusa, and H. Kan, “Theory of the intrinsic linewidth of quantum-cascade lasers: hidden reason for the narrow linewidth and line-broadening by thermal photons,” IEEE J. Quantum Electron. 44(1), 12–29 (2008). [CrossRef]
S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]
C. Sirtori, F. Capasso, J. Faist, A. L. Hutchinson, D. L. Sivco, and A. Y. Cho, “Resonant tunneling in quantum cascade lasers,” IEEE J. Quantum Electron. 34(9), 1722–1729 (1998). [CrossRef]
H. Callebaut and Q. Hu, “Importance of coherence for electron transport in terahertz quantum cascade lasers,” J. Appl. Phys. 98(10), 104505 (2005). [CrossRef]
L. Davidovich, “Sub-Poissonian processes in quantum optics,” Rev. Mod. Phys. 68(1), 127–173 (1996). [CrossRef]
G. P. Agrawal and C. M. Bowden, “Concept of linewidth enhancement factor in semiconductor lasers: its usefulness and limitations,” IEEE Photon. Technol. Lett. 5(6), 640–642 (1993). [CrossRef]
M. O. Scully and M. S. Zubairy, Quantum optics (Cambridge University Press, 1997).
W. H. Louisell, Quantum statistical properties of radiation (Wiley, 1973).
K. Petermann, Laser diode modulation and noise (Kluwar Academic, 1988).
C. Henry, “Theory of the linewidth of semiconductor lasers,” IEEE J. Quantum Electron. 18(2), 259–264 (1982). [CrossRef]
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]
B. Daino, P. Spano, M. Tamburrini, and S. Piazzolla, “Phase noise and spectral line shape in semiconductor lasers,” IEEE J. Quantum Electron. 19(3), 266–270 (1983). [CrossRef]
C. H. Henry and R. F. Kazarinov, “Quantum noise in photonics,” Rev. Mod. Phys. 68(3), 801–853 (1996). [CrossRef]
A. M. Andrews, A. Benz, C. Deutsch, G. Fasching, K. Unterrainer, P. Klang, W. Schrenk, and G. Strasser, “Doping dependence of LO-phonon depletion scheme THz quantum-cascade lasers,” Mater. Sci. Eng. B 147(2-3), 152–155 (2008). [CrossRef]
S. Kumar, Q. Hu, and J. L. Reno, “186 K operation of terahertz quantum-cascade lasers based on a diagonal design,” Appl. Phys. Lett. 94(13), 131105 (2009). [CrossRef]
D. Hofstetter, M. Beck, T. Aellen, and J. Faist, “High-temperature operation of distributed feedback quantum-cascade lasers at 5.3 ?m,” Appl. Phys. Lett. 78(4), 396 (2001). [CrossRef]

Appl. Phys. Lett.

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]
S. Kumar, Q. Hu, and J. L. Reno, “186 K operation of terahertz quantum-cascade lasers based on a diagonal design,” Appl. Phys. Lett. 94(13), 131105 (2009). [CrossRef]
D. Hofstetter, M. Beck, T. Aellen, and J. Faist, “High-temperature operation of distributed feedback quantum-cascade lasers at 5.3 ?m,” Appl. Phys. Lett. 78(4), 396 (2001). [CrossRef]

Chem. Phys. Lett.

R. F. Curl, F. Capsso, C. Gmachl, A. A. Kosterev, B. McManus, R. Lewicki, M. Pusharsky, G. Wysocki, and F. K. Tittel, “Quantum cascade lasers in chemical physics,” Chem. Phys. Lett. 487(1-3), 1–18 (2010). [CrossRef]

IEEE J. Quantum Electron.

M. Yamanishi, T. Edamura, K. Fujita, N. Akikusa, and H. Kan, “Theory of the intrinsic linewidth of quantum-cascade lasers: hidden reason for the narrow linewidth and line-broadening by thermal photons,” IEEE J. Quantum Electron. 44(1), 12–29 (2008). [CrossRef]
C. Sirtori, F. Capasso, J. Faist, A. L. Hutchinson, D. L. Sivco, and A. Y. Cho, “Resonant tunneling in quantum cascade lasers,” IEEE J. Quantum Electron. 34(9), 1722–1729 (1998). [CrossRef]
B. Daino, P. Spano, M. Tamburrini, and S. Piazzolla, “Phase noise and spectral line shape in semiconductor lasers,” IEEE J. Quantum Electron. 19(3), 266–270 (1983). [CrossRef]
C. Henry, “Theory of the linewidth of semiconductor lasers,” IEEE J. Quantum Electron. 18(2), 259–264 (1982). [CrossRef]

IEEE Photon. Technol. Lett.

G. P. Agrawal and C. M. Bowden, “Concept of linewidth enhancement factor in semiconductor lasers: its usefulness and limitations,” IEEE Photon. Technol. Lett. 5(6), 640–642 (1993). [CrossRef]

J. Appl. Phys.

H. Callebaut and Q. Hu, “Importance of coherence for electron transport in terahertz quantum cascade lasers,” J. Appl. Phys. 98(10), 104505 (2005). [CrossRef]

Mater. Sci. Eng. B

A. M. Andrews, A. Benz, C. Deutsch, G. Fasching, K. Unterrainer, P. Klang, W. Schrenk, and G. Strasser, “Doping dependence of LO-phonon depletion scheme THz quantum-cascade lasers,” Mater. Sci. Eng. B 147(2-3), 152–155 (2008). [CrossRef]

Opt. Express

Phys. Rev. B

T. Liu and Q. J. Wang, “Fundamental frequency noise and linewidth broadening caused by intrinsic temperature fluctuations in quantum cascade lasers,” Phys. Rev. B 84(12), 125322 (2011). [CrossRef]
S. Kumar and Q. Hu, “Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers,” Phys. Rev. B 80(24), 245316 (2009). [CrossRef]

Phys. Rev. Lett.

S. Bartalini, S. Borri, P. Cancio, A. Castrillo, I. Galli, G. Giusfredi, D. Mazzotti, L. Gianfrani, and P. De Natale, “Observing the intrinsic linewidth of a quantum-cascade laser: beyond the Schawlow-Townes limit,” Phys. Rev. Lett. 104(8), 083904 (2010). [CrossRef] [PubMed]

Rev. Mod. Phys.

L. Davidovich, “Sub-Poissonian processes in quantum optics,” Rev. Mod. Phys. 68(1), 127–173 (1996). [CrossRef]
C. H. Henry and R. F. Kazarinov, “Quantum noise in photonics,” Rev. Mod. Phys. 68(3), 801–853 (1996). [CrossRef]

Science

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]

Other

M. S. Vitiello, L. Consolino, S. Bartalini, A. Tredicucci, M. Inguscio, and P. De Natale, “The intrinsic linewidth of THz quantum cascade lasers,” CLEO: Science and Innovations, (Optical Society of America, 2012), paper CTu2B.2.
M. O. Scully and M. S. Zubairy, Quantum optics (Cambridge University Press, 1997).
W. H. Louisell, Quantum statistical properties of radiation (Wiley, 1973).
K. Petermann, Laser diode modulation and noise (Kluwar Academic, 1988).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper