The modulation bandwidth of quantum well nanoLED and nanolaser devices is calculated from the laser rate equations using a detailed model for the Purcell enhanced spontaneous emission. It is found that the Purcell enhancement saturates when the cavity quality-factor is increased, which limits the maximum achievable spontaneous recombination rate. The modulation bandwidth is thereby limited to a few tens of GHz for realistic devices.
©2010 Optical Society of America
In the emerging field of photonic signal processing there is an expressed need for small and efficient light emitters that can be integrated in photonic circuits and which have large modulation bandwidths. Today’s fastest diode lasers exhibit modulation speeds on the order of a few tens of GHz and nanolasers and nanoLEDs have been suggested as interesting alternatives. In the past two decades process technology has matured sufficiently to enable fabrication of active cavities with high quality-factors and low mode volumes , where the Purcell enhancement of spontaneous emission  can become large enough to affect the dynamics of the device [3–6]. Early predictions of the properties of nanolasers by Altug et al. , based on an analysis of the laser rate equations with a phenomenological Purcell enhancement of the spontaneous emission rate, suggested modulation speeds exceeding 100 GHz. Recently, the same system was investigated numerically in a paper by Lau et al.  where the importance of gain suppression was pointed out and it was indicated that the ultrahigh modulation speeds reported by Altug et al. were the result of the measurement method and are in reality only achievable in non-lasing devices with ultralow mode volume. Because of the potentially very high Purcell factor, the dynamics of the device is very dependent on the details of the spontaneous emission and a rigorous treatment of spontaneous emission is needed. In this paper we therefore calculate the spontaneous emission from fundamental principles and study the effect on the dynamic properties of nanoLEDs and nanolasers. We find that the high-speed properties predicted in ref . for nanoLEDs are deteriorated due to a reduction and saturation of the Purcell enhancement factor, which is intrinsic to quantum well and bulk active materials.
First, we will state the basic laser rate equations and a general expression for the 3dB-bandwidth. Then we will outline the model for the Purcell enhanced spontaneous emission rate and discuss the general behavior of the calculated rates. Next, the steady-state carrier and photon densities, the spontaneous and stimulated emission rates and the 3dB-bandwidth are calculated for two characteristic devices: one in the LED regime and one in the lasing regime and the results are discussed in the context of previous models. Finally, the 3dB-bandwidth is calculated for a wide range of nanoLEDs and nanolasers.
General rate equations and modulation response8] (see Table 1 for parameter values used in this paper). The total carrier recombination rate has been separated into contributions from stimulated emission (Rst), spontaneous emission into the cavity (Rc), spontaneous emission into all other modes (Rb) and non-radiative losses (Rnr). For the latter we approximate Rnr = N/τnr with τnr being the non-radiative life time. For the stimulated emission we use7]. Here vg = c/n is the group velocity, G0 is the material gain, Ntr is the transparency carrier density and Ns is a gain parameter . The gain, G, in Eq. (3) is inversely proportional to 1 + εS  to include gain suppression due to spectral hole burning and carrier heating at high photon densities.
An important measure for the performance of laser diodes is the modulation bandwidth which can be estimated from a small-signal analysis of Eq. (1) and (2). Following a standard derivation (ref . p. 195-201) the general form of the small-signal response function can be written asEq. (1) and (2), ωR and γR are given byEq. (5) - (7) to calculate the modulation response.
Purcell enhanced spontaneous emission
For nanocavity devices with high Q and low mode volume (V) the Purcell effect should be taken into account when calculating the spontaneous emission [2,5,6,11]. The Purcell factor reflects the increased spectral density of a single mode centered in a cavity and can be written asEq. (8) to describe its magnitude. We consider two models for spontaneous emission, which include the Purcell enhancement.
In conventional semiconductor diodes operated near threshold the spontaneous emission can be approximated as being linear in the carrier density, Rsp = N/τsp, where τsp is the spontaneous lifetime in bulk. In the linear model the spontaneous emission is often split so that Rc = βRsp and Rb = (1-β)Rsp, where β = Rc/Rsp is the spontaneous emission factor.
The Purcell effect is taken into account phenomenologically by multiplying F onto Rc [5,6,12,13] to include the relatively larger emission into the cavity. We use the spontaneous emission terms from ref . as an example of the linear model, i.e.Eq. (9) and (10) into Eq. (6) and (7) above, we thus recover Eq. (9) and (10) of ref . For systems with β ~1 this model assumes that the Purcell effect acts on all the carriers in the system as it does not take into account the details in the optical density-of-states (DOS), which are important when discussing devices with large Purcell factors. Furthermore, the linear model is only valid close to its Taylor expansion point.
In the full model we follow ref . and write the differential spontaneous recombination rate for a semiconductor device in terms of the differential carrier population, dN2, as dRsp = A21dN2, where the Einstein coefficient (A21) is proportional to the optical DOS (ρop) and the B21 coefficient, i.e. A21 = hνρopB21. dN2 can be written in terms of the electron (f2) and hole (f1) Fermi functions and the reduced electronic DOS (ρel) as dN2 = ρelf2(1-f1)dE21 and homogeneous broadening is taken into account by convoluting with the lineshape function L(E - hν), which is usually taken to be a Lorentzian. Integrating over energy we arrive at an expression for spontaneous emission rate (ref . p. 459-472)
The optical DOS is modeled as the bulk DOS plus a Lorentzian shaped cavity spectrum  placed in the center of a photonic bandgap, i.e.
Here H is the Heaviside function and the cavity is centered at hν0 = (hνU - hνL)/2, with hνL (hνU) being the lower (upper) edge of the photonic band gap and the width of the Lorentzian (c.f. Figure 1 ). The Einstein B21 coefficient is found by setting A21(hν0) = 1/τ21 = hν0B21ρop(hν0), where the differential recombination time (τ21) is chosen so that the bulk spontaneous emission time (τsp) is recovered in the absence of the cavity, i.e. so that Rsp(Ntr) = Ntr/τsp in bulk. We then have17] and follows a sum rule , is not taken into account.
The reduced electronic DOS is given by 
In this model (Eq. (11) - (14)), which we will refer to as the “full model”, Rb is the integral over the bulk part of Eq. (13) while the Lorentzian part gives Rc. Although the full model specifically describes quantum well devices it is also suitable for bulk devices if an appropriate electronic DOS is chosen. Quantum dot devices, however, are governed by a different set of laser rate equations than Eq. (1) and (2) and should account for the dynamics in both the wetting layer and quantum dot levels. Furthermore, the electronic DOS for quantum dots is markedly different from that of bulk and quantum wells and including quantum dots in the treatment would make the discussion less clear. Although the use of quantum dot structures may be very promising, we therefore only treat quantum well (and bulk) devices in this paper.
Saturation of Purcell enhancement
In Fig. 2a the spontaneous emission calculated in the linear and the full model are compared. It is seen that the spontaneous emission in the full model is significantly lower than in the linear model and does not have a linear dependence on N. To see clearly the contribution from the electronic part of Eq. (11) we also plot (for a Q of 100)Fig. 2a, Rel approaches the slope of the linear model for high carrier density, i.e. for the quasi-Fermi level separation high above the electronic band edge, where the integral in Eq. (15) is almost equal to the carrier density. For low carrier density Rel follows Rc, although Rc is somewhat smaller due to the limiting effect of the cavity DOS. The saturation of Rc for N > Ntr is not seen in Rel and originates from the finite quantum well DOS and associated band filling effects at the cavity resonance E0, i.e. the quasi-Fermi level separation becomes much larger than E0 so that f2 (1 – f1) ~1 near E0.
Perhaps the most interesting feature in Fig. 2a is that Rc does not increase with the Q-factor in the full model. This can be explained by evaluating the integral over frequency in Eq. (11). This is usually done by assuming the homogeneous broadening term, L (E - hν), to be sharply peaked at hν = E compared to the other terms in Eq. (11) so that it can be replaced by a Dirac delta-function. However, for high-Q cavities the cavity linewidth can easily become smaller than the homogeneous broadening and in this case the integral over frequency in Eq. (11) should be evaluated exactly. The results in the present model, however, are not significantly changed when evaluating the frequency integral exactly using a Lorentzian for the homogeneous broadening as the electronic DOS is much wider than the homogeneous broadening. Therefore, we will use L(E – hν) = δ(E - hν) from here on.
The cavity part of the integral over frequency in Eq. (11) becomesEq. (11) becomesEq. (17), which is the case for high Q, the remaining integral is effectively over the Lorentzian and we end up withEq. (17), i.e. for Q larger than a few hundreds. The reduction thus originates from the mismatch between cavity bandwidth and the effective inhomogeneous broadening expressed by the electronic DOS. This agrees with the qualitative discussion given in ref .
Notice that the expression in Eq. (18) is still inversely proportional to the mode volume, so that the rate can be enhanced by lowering Vn. Note also that we here neglect features in the electronic DOS such as the exciton peak near the band edge, which may challenge the requirement that the electronic DOS be slowly varying. Deviations of this type from the assumed smooth electronic DOS do not affect the conclusions, but rather increase the Q-factor at which the Purcell enhancement saturates. However, a rigorous treatment of the electronic DOS is beyond the scope of this work.
Finally, we note that the emission into the background (Rb) shown in Fig. 2a is negligible for N < 10 Ntr. Reducing the size of the photonic bandgap increases Rb, but the change is insignificant compared to Rc.
Effective Purcell factor
In order to make the subsequent discussion more transparent we now introduce an effective Purcell factor, which we define as
Using the definition in Eq. (19) and the expression for Rc in Eq. (9) the effective Purcell factor reduces to Feff = βF in the linear model. In the full model the effective Purcell factor is a function of the optical and electronic DOS and thus depends on Q, Vn and N and is therefore markedly different from the Feff in the linear model.
The effective Purcell factors in the two models are plotted in Fig. 2b for the devices in Fig. 2a. For N < Ntr we have Feff in the full model almost constant, but lower than Feff in the linear model due to the reduction originating from the limited cavity bandwidth. The decrease for N > Ntr in the full model is a combination of the bandfilling effect and Rb growing large.
We now proceed to calculate the steady-state carrier and photon densities, the spontaneous and stimulated emission and the 3dB-bandwidth for two specific devices in the linear and full model. To ease comparison, the β-factor calculated from the full model is used in the linear model. Device A has a Q-factor of 104 and a mode volume of 10 Vn and device B has a Q-factor of 102 and a mode volume of 0.1 Vn. The results for device A are plotted in Fig. 3a-c , and for device B in Fig. 3d-f.
In Fig. 3a the carrier and photon densities are shown as a function of the pump (in units of J0 = Ntr/τsp). For J < J0 the carrier density in the linear model is much lower than in the full model. This is an effect of the reduction of the effective Purcell factor for high Q-factors in the full model, which limits Feff to ~1 for low carrier density, whereas Feff ~610 in the linear model. The higher spontaneous emission rate in the linear model explains the lower carrier density.
The carrier density in the full model clamps for J ~J0, when lasing sets in, and remains constant until J ~100 J0, where it again begins to increase. This increase is due to the gain suppression that becomes a significant process for S0 ~1/ε ~ 0.05 Ntr and must be compensated by an increase of the linear gain. The β-factor is close to unity for the entire pumping range, which explains why the photon densities for the two models are equal for J > 10 J0 even though the linear model is dominated by spontaneous emission, while the full model is dominated by stimulated emission.
Figure 3b compares the spontaneous and stimulated emission for the two models. In the full model the spontaneous emission increases strongly until the carrier density clamps, where after the stimulated emission becomes dominant. The spontaneous emission remains constant until J ~100 J0, where the carrier density begins to increase again due to the gain suppression effect discussed above. In the linear model the carrier density is much lower than in the full model, giving a lower stimulated emission. This pushes the threshold pump up to around J ~1000 J0.
Several features observed in Fig. 3a and b are also found in the 3dB-bandwidth in Fig. 3c. In the full model the 3dB-bandwidth is dominated by spontaneous emission until J ~J0, after which stimulated emission dominates until gain suppression becomes important around J ~100 J0. Over the entire pumping range the 3dB-bandwidth does not exceed 20 GHz. In the linear model the high 3dB-bandwidth for J < 100 J0 is due to the high spontaneous emission and the drop-off for J >100 J0 is due to the damping rate (γR) increasing more rapidly than the resonance frequency (ωR) as explained in ref .
The results for device B are shown in Fig. 3d-f, where both the Q-factor and mode volume, Vn, are 100 times lower than in device A, so that Feff remains ~610 in the linear model, but changes to ~70 in the full model (for low carrier density). For this device the photon loss is 100 times larger than for device A and therefore the photon densities for the two models in Fig. 3d are 100 times lower. In the linear model, the relatively lower photon density is reflected in the carrier density, which is also lowered to balance the photon loss.
In the full model the carrier density is only slightly lower than for device A for J < J0. This is because the spontaneous emission rate in the full model follows another dependence on N (approximately N2) than in the linear model and therefore a smaller adjustment of N is necessary to compensate the lower photon density.
In Fig. 3e the spontaneous emission rate in the linear model is lower by a factor corresponding to the lowering of carrier density compared with device A. In the full model the spontaneous emission almost follows the increase in Feff, but is also modified by the lower carrier density. In neither models the gain becomes large enough to initiate lasing and this is reflected in the carrier densities, which do not clamp in this device.
Figure 3f shows that the 3dB-bandwidth calculated in the two models differs significantly. The 3dB-bandwidth is almost an order of magnitude lower in the full model compared to the linear model. This can be explained by studying Eq. (5) in the LED regime, i.e. for dominating spontaneous emission. For the linear model Eq. (5) reduces toFig. 3f is due to the effective Purcell factor growing large. In the full model the 3dB-bandwidth is also given by Eq. (21), but in this model we haveEq. (21), Eq. (23) and Fig. 2b the behavior of the 3dB-bandwidth in the full model now becomes clear.
For J < 50 J0 the effective Purcell factor is almost constant and Rb,N ≈0, so that the first term in Eq. (23) determines τeff. For 50 J0 < J < 800 J0 the term Rb,N is still low, while Feff starts to decrease so that the first term of Eq. (23) becomes smaller and the second term becomes negative, leading to a decrease in the 3dB-bandwidth. For J > 800 J0 the background emission increases sharply, making τeff decrease and leading to the final increase in the 3dB-bandwidth. Thus the 3dB-bandwidth in the full model is roughly an order of magnitude lower than in the linear model and this clearly underlines the necessity for a detailed description of the spontaneous emission.
The same analysis can be made for Fig. 3c for J < 2J0, as spontaneous emission is also dominant in this pumping range. Here the effective Purcell factor in the full model is 100 times lower than for device B, making 1/τnr the dominating term in Eq. (23) and giving the low 3dB-bandwidth compared to device B. Another way of expressing the behavior is that the Purcell enhancement only affects the carriers associated with transitions within the bandwidth of the cavity resonance. When the cavity resonance is much narrower than the electronic bandwidth the influence from the enhanced spontaneous emission only has a small effect on the total carrier density life time and thus also the speed.
In order to examine the modulation bandwidth for a large range of devices, the 3dB-bandwidth has been calculated in the full model versus Q-factor and mode volume for a range of devices and the result is plotted in Fig. 4 for J = J0 and J = 100 J0. The white line in the figure shows, for the given pump current, the device parameters for which stimulated and spontaneous emission are equal in magnitude. Devices to the right of this boundary are dominated by stimulated emission. The black line shows the limiting value of this boundary for large pump current, i.e. devices to the left of this line will always be dominated by spontaneous emission, independently of the strength of the pump, and are thus always in the LED regime.
All the devices in Fig. 4a are dominated by spontaneous emission for the given pump and are therefore in the LED regime. It is seen that the 3dB-bandwidth is below 3 GHz for most devices, except at extremely low mode volumes, which are probably unattainable in practice, where the bandwidth starts to increase. This is due to the effective Purcell factor, which is large at low mode volume, so that the first term in Eq. (23) dominates. The slight dependence on Q is due to the photon life time that becomes smaller at low Q and thereby increases the 3dB-bandwidth in Eq. (21).
The same dependence on Q is seen in Fig. 4b, where the 3dB-bandwidth exceeds 200 GHz in the lower left corner, i.e. in the LED regime. In the top right area, which corresponds to conventional laser structures, the effective Purcell factor saturates at a Q of a few hundreds and the stimulated emission therefore becomes the dominant recombination process so that the term aS0/τp in Eq. (6) becomes large, giving the large modulation speed in this area. In the lower right corner the photon loss is too large to meet the lasing condition and the effective Purcell factor is low, giving the lower 3dB-bandwidth. In general the ultrahigh modulation speeds previously reported  are not seen, neither at low or high pump, because the effective Purcell factor saturates and becomes independent of the Q-factor as discussed above.
We note that for two level structures, such as quantum dot devices the upper left corner of Fig. 4a and b may correspond to devices operating in the strong coupling regime , where the calculated 3dB-bandwidth would be invalid. However, for quantum well structures and operation at room temperature, as considered here, this is not considered to be an issue.
Discussion and conclusions
In this paper we have discussed how to treat Purcell enhancement of spontaneous emission in quantum well nanoLEDs and nanolasers in the laser rate equation model. In the present model for the spontaneous emission rate, which incorporates details about the optical density-of-states, it was shown that the spontaneous emission rate saturates for high carrier density and that it does not follow the usual linear enhancement with the cavity Q-factor. The independence of the Q-factor was found to originate from the cavity bandwidth being small compared to the electronic density-of-states and this limits the magnitude of the enhancement of the spontaneous emission in nanocavity systems. We introduced an effective Purcell enhancement that incorporates the important features due to the optical and electronic density-of-states and can be used to describe the spontaneous emission in a simple way. The limited effective Purcell enhancement entails that for devices with mode volumes attainable with today’s technology the highest modulation bandwidth is found above threshold and limited by well-known damping effects due to gain nonlinearities [6,18]. The ultrafast modulation speeds previously reported in ref . are only found in devices with extremely low mode volume and even for these devices the 3dB-bandwidth is below a few hundred GHz and the spectral bandwidth is very broad making wavelength division multiplexing systems difficult.
The limited modulation bandwidth is due to the wide electronic density-of-states as compared to the cavity spectrum. However, it can also be seen from Eq. (17) that if the electronic density-of-states function is narrower than the cavity bandwidth the saturation of the effective Purcell factor is lifted. This indicates that ultra-high bandwidths may be obtainable in quantum dot systems. However, it will be necessary to include inhomogeneous broadening in the analysis of quantum dot systems as this effect will widen the electronic density-of-states significantly. The same is true for the homogeneous broadening, which will be present even if the inhomogeneous broadening is made insignificant. A rigorous treatment of quantum dot systems is beyond the scope of this paper and whether the ultra-high modulation speeds are attainable in quantum dot systems remains an open question.
The authors would like to thank Michael Lorke for helpful discussions and Villum Fonden for financial support via the NATEC (Nanophotonics for TErabit Communications) center.
References and links
1. Y. Akahane, T. Asano, B.-S. Song, and S. Noda, “High-Q photonic nanocavity in a two-dimensional photonic crystal,” Letters to Nature 425(6961), 944–947 (2003). [CrossRef]
2. E. M. Purcell, “Spontaneous emission probabilities at radio frequencies,” Phys. Rev. 69, 681 (1946).
4. O. Painter, R. K. Lee, A. Scherer, A. Yariv, J. D. O’Brien, P. D. Dapkus, and I. Kim I, “Two-dimensional photonic band-Gap defect mode laser,” Science 284(5421), 1819–1821 (1999). [CrossRef] [PubMed]
5. H. Altug, D. Englund, and J. Vučković, “Ultrafast photonic crystal nanocavity laser,” Nat. Phys. 2(7), 484–488 (2006). [CrossRef]
7. L. A. Coldren, and S. W. Corzine, Diode Lasers and Photonic Integrated Circuits (John Wiley & Sons, inc., New York, 1995).
8. J. D. Jackson, Classical Electrodynamics (John Wiley & Sons, inc., New York, 1998)
9. A. Mecozzi and J. Mørk, “Saturation induced by picosecond pulses in semiconductor optical amplifiers,” J. Opt. Soc. Am. B 14(4), 761–770 (1997). [CrossRef]
10. J. D. Joannopoulos, S. G. Johnson, and J. N. Winn, Photonic Crystals: Molding the Flow of Light (Princeton University Press, New Jersey, 2008).
11. J.-M. Gerard, “Solid-state cavity-quantum electrodynamics with self-assembled quantum dots”, in Single Quantum Dots, Fundamentals, Applications and New Concepts, P. Michler (Springer, Berlin, 2003), pp. 269–314.
12. H. Yokoyama and S. D. Brorson, “Rate equation analysis of microcavity lasers,” J. Appl. Phys. 66(10), 4801 (1989). [CrossRef]
13. T. Baba, “Photonic crystals and microdisk cavities based on GaInAsP-InP system,” IEEE J. Sel. Top. Quantum Electron. 3(3), 808 (1997). [CrossRef]
14. D. Englund, D. Fattal, E. Waks, G. Solomon, B. Zhang, T. Nakaoka, Y. Arakawa, Y. Yamamoto, and J. Vucković, “Controlling the spontaneous emission rate of single quantum dots in a two-dimensional photonic crystal,” Phys. Rev. Lett. 95(1), 013904 (2005). [CrossRef] [PubMed]
16. T. Baba, D. Sano, K. Nozaki, K. Inoshita, Y. Kuroki, and F. Koyama, “Observation of fast spontaneous emission decay in GaInAsP photonic crystal point defect nanocavity at room temperature,” Appl. Phys. Lett. 85(18), 3989–3991 (2004). [CrossRef]
17. M. P. Marder, Condensed Matter Physics (John Wiley & Sons, inc., New York, 2000). [PubMed]
18. R. S. Tucker, “High-speed modulation of semiconductor lasers,” J. Lightwave Technol. 3(6), 1180–1192 (1985). [CrossRef]