Institut f¨ur Theoretische Physik, Technische Universit¨at Berlin, Germany 2 Sandia National Laboratories, Albuquerque, NM, USA

Abstract: We investigate the dependence of the amplitude-phase coupling in quantum-dot (QD) lasers on the charge-carrier scattering timescales. The carrier scattering processes influence the relaxation oscillation parameters, as well as the frequency chirp, which are both important parameters when determining the modulation performance of the laser device and its reaction to optical perturbations. We find that the FM/AM response exhibits a strong dependence on the modulation frequency, which leads to a modified optical response of QD lasers when compared to conventional laser devices. Furthermore, the frequency response curve changes with the scattering time scales, which can allow for an optimization of the laser stability towards optical perturbations. © 2014 Optical Society of America OCIS codes: (140.5960) Semiconductor lasers; (230.5590) Quantum-well, -wire and -dot devices; (270.3100) Instabilities and chaos.

References and links 1. W. W. Chow and F. Jahnke, “On the physics of semiconductor quantum dots for applications in lasers and quantum optics,” Prog. Quantum Electron. 37, 109–184 (2013). 2. G. H. M. van Tartwijk and D. Lenstra, “Semiconductor laser with optical injection and feedback,” Quantum Semiclass. Opt. 7, 87–143 (1995). 3. T. B. Simpson, J. M. Liu, K. F. Huang, and K. Tai, “Nonlinear dynamics induced by external optical injection in semiconductor lasers,” Quantum Semiclass. Opt. 9, 765–784 (1997). 4. S. Wieczorek, B. Krauskopf, T. Simpson, and D. Lenstra, “The dynamical complexity of optically injected semiconductor lasers,” Phys. Rep. 416, 1–128 (2005). 5. N. A. Naderi, M. Pochet, F. Grillot, N. B. Terry, V. Kovanis, and L. F. Lester, “Modeling the injection-locked behavior of a quantum dash semiconductor laser,” IEEE J. Sel. Top. Quantum Electron. 15, 563–571 (2009). 6. B. Kelleher, D. Goulding, S. P. Hegarty, G. Huyet, D. Y. Cong, A. Martinez, A. Lemaitre, A. Ramdane, M. Fischer, F. Gersch¨utz, and J. Koeth, “Excitable phase slips in an injection-locked single-mode quantum-dot laser,” Opt. Lett. 34, 440–442 (2009). 7. J. Pausch, C. Otto, E. Tylaite, N. Majer, E. Sch¨oll, and K. L¨udge, “Optically injected quantum dot lasers - impact of nonlinear carrier lifetimes on frequency locking dynamics,” New J. Phys. 14, 053018 (2012). 8. D. Ziemann, R. Aust, B. Lingnau, E. Sch¨oll, and K. L¨udge, “Optical injection enables coherence resonance in quantum-dot lasers,” Europhys. Lett. 103, 14002 (2013). 9. B. Lingnau, K. L¨udge, W. W. Chow, and E. Sch¨oll, “Failure of the α-factor in describing dynamical instabilities and chaos in quantum-dot lasers,” Phys. Rev. E 86, 065201(R) (2012). 10. B. Lingnau, W. W. Chow, E. Sch¨oll, and K. L¨udge, “Feedback and injection locking instabilities in quantum-dot lasers: a microscopically based bifurcation analysis,” New J. Phys. 15, 093031 (2013). 11. R. Lang and K. Kobayashi, “External optical feedback effects on semiconductor injection laser properties,” IEEE J. Quantum Electron. 16, 347–355 (1980). 12. C. Masoller and M. S. Torre, “Influence of optical feedback on the polarization switching of vertical-cavity surface-emitting lasers,” IEEE J. Quantum Electron. 41, 483–489 (2005).

#204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4867

13. S. Azouigui, B. Dagens, F. Lelarge, J. G. Provost, A. Accard, F. Grillot, A. Martinez, Q. Zou, and A. Ramdane, “Tolerance to optical feedback of 10-gb/s quantum-dash-based lasers emitting at 1.51µ m,” IEEE Photon. Technol. Lett. 19, 1181–1183 (2007). 14. H. Erzgr¨aber, B. Krauskopf, and D. Lenstra, “Bifurcation analysis of a semiconductor laser with filtered optical feedback,” SIAM J. Appl. Dyn. Syst. 6, 1–28 (2007). 15. M. Gioannini, G. A. P. The, and I. Montrosset, “Multi-population rate equation simulation of quantum dot semiconductor lasers with feedback,” in “Numerical Simulation of Optoelectronic Devices, 2008. NUSOD ’08. International Conference on,” (2008), pp. 101–102. 16. C. Otto, K. L¨udge, and E. Sch¨oll, “Modeling quantum dot lasers with optical feedback: sensitivity of bifurcation scenarios,” phys. stat. sol. (b) 247, 829–845 (2010). 17. F. Albert, C. Hopfmann, S. Reitzenstein, C. Schneider, S. H¨ofling, L. Worschech, M. Kamp, W. Kinzel, A. Forchel, and I. Kanter, “Observing chaos for quantum-dot microlasers with external feedback,” Nature Communications 2, 366 (2011). 18. S. Osborne, P. Heinricht, N. Brandonisio, A. Amann, and S. O’Brien, “Wavelength switching dynamics of twocolour semiconductor lasers with optical injection and feedback,” Semicond. Sci. Technol. 27, 094001 (2012). 19. C. Otto, B. Globisch, K. L¨udge, E. Sch¨oll, and T. Erneux, “Complex dynamics of semiconductor quantum dot lasers subject to delayed optical feedback,” Int. J. Bif. Chaos 22, 1250246 (2012). 20. B. Globisch, C. Otto, E. Sch¨oll, and K. L¨udge, “Influence of carrier lifetimes on the dynamical behavior of quantum-dot lasers subject to optical feedback,” Phys. Rev. E 86, 046201 (2012). 21. B. Dagens, A. Markus, J. Chen, J. G. Provost, D. Make, O. Le Gouezigou, J. Landreau, A. Fiore, and B. Thedrez, “Giant linewidth enhancement factor and purely frequency modulated emission from quantum dot laser,” Electron. Lett. 41, 323–324 (2005). 22. H. Su and L. F. Lester, “Dynamic properties of quantum dot distributed feedback lasers: high speed, linewidth and chirp,” J. Phys. D: Appl. Phys. 38, 2112–2118 (2005). 23. K. L¨udge, R. Aust, G. Fiol, M. Stubenrauch, D. Arsenijevi´c, D. Bimberg, and E. Sch¨oll, “Large signal response of semiconductor quantum-dot lasers,” IEEE J. Quantum Electron. 46, 1755 –1762 (2010). 24. L. V. Asryan, Y. Wu, and R. A. Suris, “Carrier capture delay and modulation bandwidth in an edge-emitting quantum dot laser,” Appl. Phys. Lett. 98, 131108 (2011). 25. B. Lingnau, K. L¨udge, W. W. Chow, and E. Sch¨oll, “Influencing modulation properties of quantum-dot semiconductor lasers by electron lifetime engineering,” Appl. Phys. Lett. 101, 131107 (2012). 26. C. Wang, F. Grillot, and J. Even, “Impacts of wetting layer and excited state on the modulation response of quantum-dot lasers,” IEEE J. Quantum Electron. 48, 1144–1150 (2012). 27. D. Gready, G. Eisenstein, M. Gioannini, I. Montrosset, D. Arsenijevi´c, H. Schmeckebier, M. Stubenrauch, and D. Bimberg, “On the relationship between small and large signal modulation capabilities in highly nonlinear quantum dot lasers,” Appl. Phys. Lett. 102, 101107 (2013). 28. T. C. Newell, D. J. Bossert, A. Stintz, B. Fuchs, K. J. Malloy, and L. F. Lester, “Gain and linewidth enhancement factor in InAs quantum-dot laser diodes,” IEEE Photonics Technol. Lett. 11, 1527–1529 (1999). 29. F. Grillot, B. Dagens, J. G. Provost, H. Su, and L. F. Lester, “Gain compression and above-threshold linewidth enhancement factor in1.3µm InAs/GaAs quantum-dot lasers,” IEEE J. Quantum Electron. 44, 946–951 (2008). 30. S. Melnik, G. Huyet, and A. V. Uskov, “The linewidth enhancement factor α of quantum dot semiconductor lasers,” Opt. Express 14, 2950–2955 (2006). 31. M. Gioannini and I. Montrosset, “Numerical analysis of the frequency chirp in quantum-dot semiconductor lasers,” IEEE J. Quantum Electron. 43, 941–949 (2007). 32. E. J. Doedel, H. B. Keller, and J. P. Kervenez, “Numerical analysis and control of bifurcation problems. (I) Bifurcation in finite dimensions,” Int. J. Bif. Chaos 1, 493–520 (1991). 33. J. Kim, C. Meuer, D. Bimberg, and G. Eisenstein, “Numerical simulation of temporal and spectral variation of gain and phase recovery in quantum-dot semiconductor optical amplifiers,” IEEE J. Quantum Electron. 46, 405–413 (2010). 34. N. Majer, S. Dommers-V¨olkel, J. Gomis-Bresco, U. Woggon, K. L¨udge, and E. Sch¨oll, “Impact of carrier-carrier scattering and carrier heating on pulse train dynamics of quantum dot semiconductor optical amplifiers,” Appl. Phys. Lett. 99, 131102 (2011). 35. K. L¨udge, “Modeling quantum dot based laser devices,” in “Nonlinear Laser Dynamics - From Quantum Dots to Cryptography,” , K. L¨udge, ed. (WILEY-VCH Weinheim, Weinheim, 2012), chap. 1, pp. 3–34. 36. K. L¨udge, E. Sch¨oll, E. A. Viktorov, and T. Erneux, “Analytic approach to modulation properties of quantum dot lasers,” J. Appl. Phys. 109, 103112 (2011). 37. K. L¨udge, B. Lingnau, C. Otto, and E. Sch¨oll, “Understanding electrical and optical modulation properties of semiconductor quantum-dot lasers in terms of their turn-on dynamics,” Nonlinear Phenom. Complex Syst. 15, 350–359 (2012). 38. B. Kelleher, C. Bonatto, G. Huyet, and S. P. Hegarty, “Excitability in optically injected semiconductor lasers: Contrasting quantum-well- and quantum-dot-based devices,” Phys. Rev. E 83, 026207 (2011). 39. E. J. Doedel and B. E. Oldeman, Auto-07P: Continuation and bifurcation software for ordinary differential equations, Concordia University, Montreal, Canada (2009).

#204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4868

40. K. L¨udge and E. Sch¨oll, “Temperature dependent two-state lasing in quantum dot lasers,” in “Laser Dynamics and Nonlinear Photonics, Fifth Rio De La Plata Workshop 6-9 Dec. 2011,” (IEEE Publishing Services, New York, 2012), IEEE Conf. Proc., pp. 1–6. 41. B. Lingnau, K. L¨udge, W. W. Chow, and E. Sch¨oll, “Many-body effects and self-contained phase dynamics in an optically injected quantum-dot laser,” in “Semiconductor Lasers and Laser Dynamics V, Brussels,” , vol. 8432 of Proceedings of SPIE, K. Panajotov, M. Sciamanna, A. A. Valle, and R. Michalzik, eds. (SPIE, 2012), pp. 84321J 42. C. Harder, K. Vahala, and A. Yariv, “Measurement of the linewidth enhancement factor alpha of semiconductor lasers,” Appl. Phys. Lett. 42, 328–330 (1983). 43. M. Osinski and J. Buus, “Linewidth broadening factor in semiconductor lasers – an overview,” IEEE J. Quantum Electron. 23, 9–29 (1987). 44. T. Fordell and A. M. Lindberg, “Experiments on the linewidth-enhancement factor of a vertical-cavity surfaceemitting laser,” IEEE J. Quantum Electron. 43, 6–15 (2007). 45. F. Grillot, C. Wang, N. Naderi, and J. Even, “Modulation properties of self-injected quantum dot semiconductor diode lasers,” IEEE J. Quantum Electron. 19, 1900812 (2013).

1.

Introduction

Self-organized quantum-dot (QD) lasers have been the subject of extensive experimental and theoretical studies [1]. Their low threshold currents, high temperature stability, and low sensitivity to unwanted optical perturbations in comparison to conventional lasers make them potential candidates for many types of applications. Of particular interest in recent studies are the dynamical response of lasers to the deliberate injection of an optical signal into the cavity [2–10], their sensitivity to time-delayed optical self-feedback [11–20], as well as their modulation properties [21–27]. The unique dynamical behavior of QD lasers stems from the existence of localized electronic QD states. The charge carriers in these states are coupled to each other and to the surrounding charge-carrier reservoir via nonlinear scattering processes. This scattering can occur on different timescales, depending on, e.g., the material system and temperature of the device. If they are in the order of the relevant dynamical processes, they introduce additional dynamic degrees of freedom. Furthermore, the amplitude-phase coupling in QD lasers, i.e. the dynamic change of the active medium refractive index in relation to changes of the gain, commonly characterized by the linewidth enhancement factor α, has been shown to lie within a wide range of possible values [21, 28] and to exhibit strongly nonlinear dependences on the operating conditions [29]. A description of the amplitude-phase coupling that goes beyond the α-factor has been previously considered [9, 10, 30, 31], showing that in QD lasers the use of an α-factor is not always justified. Nevertheless, the characterization of QD laser devices in terms of an α-factor is still widely employed [29]. In this manuscript, we explore the influence of the charge-carrier dynamics on the dynamic response of QD lasers. Specifically, we employ a rate equation model to determine the dependence of the amplitude-phase coupling on the charge-carrier scattering processes. The model is based on the microscopically based balance equation (MBBE) model used in [10] which was further simplified to allow for numerical path continuation of the occurring bifurcations [32]. We evaluate the laser dynamics under injection of an external optical signal and compare our results with a model that treats the amplitude-phase coupling as a fixed parameter. We highlight the modifications to the laser behavior arising due to the complex carrier dynamics in QD lasers and clarify under which circumstances the approximation of the amplitude-phase coupling by an α-factor is valid. We then proceed by analyzing the amplitude and frequency response of the laser device under small-signal modulation conditions, and show that it exhibits a strong dependence on the frequency of the modulation as well as on the origin of the modulation.

#204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4869

2.

Model

We consider a 1.2mm long ridge waveguide edge-emitting dot-in-a-well (DWELL) singlemode laser device, consisting of a number of aL stacked InGaAs quantum-well (QW), each embedding a density of N QD InAs QDs. The QDs are assumed to have localized bound electron and hole states, with a ground state (GS) centered at a transition frequency corresponding to an emission linewidth of λ = 1.3µm, and the twofold degenerate first excited state (ES). Due to the inhomogeneous broadening of the QD states not all QDs contribute equally to the gain. An intuitive approach would be to calculate the dynamics of individual QD groups, each with a common transition energy, separately [10,33,34]. However, such an approach greatly increases the dimensionality of the dynamical system. A much simpler approach is to distinguish only between the resonant (active) and off-resonant (inactive) part of the QD distribution, where only the resonant part contributes to stimulated emission, as we have previously used in [23]. We extend the previous approach by explicitly tracking the occupation of inactive QDs. We assume a fraction of f act of QDs to be resonant ( f inact ≡ 1 − f act ). This approach was found to exhibit excellent quantitative agreement with the predictions of the multi-subgroup model. The QD laser is described by dynamic equations for the electric field inside the cavity E, the (in)act occupation probabilities of active and inactive QDs in the GS ρGS,b and the ES ρES,b , as well as the 2D charge-carrier density in each QW layer wb (b ∈ {e, h}): ∂ E =[(g − iδ ω) − κ]E + KκE 0 exp(−iωinjt) ∂t ∂ act 2g cap,act rel,act |E|2 − Ract ρ = − act sp,GS + SGS,b + Sb ∂t GS,b f aL 2NQD ∂ inact cap,inact ρ = − Rinact + Sbrel,inact sp,GS + SGS,b ∂t GS,b 1 act rel,act ∂ cap f Sb ρES,b = − Rsp,ES + SES,b − + f inact Sbrel,inact ∂t 2 ∂ cap,act cap,inact cap w wb =2N QD f act SGS,b + f inact SGS,b − 4N QD SES,b − rloss +J. ∂t

(1) (2) (3) (4) (5)

Here, κ is the cavity loss rate, K is the optical injection strength, E 0 is the steady-state electricfield amplitude, ωinj is the injection frequency. The number of active layers is aL , N QD is the 2D (in)act

(in)act

(in)act

QD density per QW layer, Rsp,m = Wm ρm,e · ρm,h describes spontaneous recombination w = Rw w · w is an effective QW loss rate, taking losses in the QD states (m ∈ {GS, ES}), rloss h loss e into account spontaneous emission losses as well as defect and Auger recombination as one phenomenological rate, J is the pump current density per QW layer. In the description of the stimulated emission, only the active QD GS states are assumed to contribute to the gain, while all off-resonant states only lead to an index change [31]. The magnitude of the gain and index change coefficients were calculated with a model treating the QD subgroups and the QW kstates separately, each contributing to the total gain and index change [9]. Using the complete microscopically based laser equations, the complex gain was found to be very well described by a linear dependence on the respective carrier occupations and densities over the operating parameter range typically encountered in the laser device. We can thus write the amplitude gain g and change of instantaneous frequency δ ω due to refractive index changes as follows: act act g =gGS (ρGS,e + ρGS,h − 1) ,

δω

e h wh . =δ ωES (ρES,e + ρES,h ) + δ ωQW we + δ ωQW

(6) (7)

The linear frequency change coefficients are extracted from our microscopic model [9]. #204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4870

Table 1. Fit parameters for charge-carrier scattering processes, extracted from microscopic calculations for a GS confinement energy of 64(35)meV and a GS-ES separation of 50(20)meV for electrons (holes), and T = 300K.

b: m: A (10−11 cm2 ns−1 ) B (10−11 cm2 ) C (ns−1 ) D (10−11 cm2 )

electrons GS ES 18.5 48.3 1.9 0.48 1014 1.4

holes GS ES 10.5 21.4 5.3 1.8 2272 2.3

The scattering processes between the different charge-carrier states are given by cap,(in)act (in)act cap,in cap,out (in)act Sm,b =Sm,b 1 − ρm,b − Sm,b ρm,b (in)act rel,(in)act (in)act Sb =Sbrel,in ρES,b 1 − ρGS,b − Sbrel,out 1 − ρES,b ρGS,b .

(8) (9)

We approximate the nonlinear in-scattering rates by the following expressions, derived from fits to the full microscopically calculated rates [34]: Am,b wb 2 Bm,b + wb Cb wb , Sbrel,in (wb ) = Db + wb

cap,in Sm,b (wb ) =

(10) (11)

with the fit parameters given in Table 2. Such an approach greatly reduces the numerical complexity involved to describe the charge-carrier scattering processes, while still maintaining the most important features of the carrier-density dependence of the scattering rates. The corresponding out-scattering rates are given by the detailed balance relationships [35]: eq EF − εm,b cap,out cap,in Sm,b =Sm,b exp − b (12) kB T εES,b − εGS,b , (13) Sbrel,out =Sbrel,in exp − kB T eq

where EF b is the quasi-Fermi level of the QW, εm,b is the mean energy of the corresponding localized QD state, kB is the Boltzmann constant, and T is the charge-carrier temperature. The numerical model treats the carrier dynamics within the QD laser on a microscopic level and thus allows for quantitative predictions of the QD laser dynamics if the bandstructure of the device is known. Thus, with the correct energy structure as an input parameter, it can be used for device optimization. In our simulations, we use scattering rates that were microscopically calculated for a given QD structure with confinement energies as given in Tab. 2. Since the scattering rates depend crucially on, e.g., the material system and device temperature, real devices can exhibit scattering rates that may differ from our calculated rates. Thus, in the following investigation, we will vary the rates and present results in dependence on their magnitude, in order to make predictions for a variety of possible QD laser devices. 3. 3.1.

Results Solitary laser dynamics

We simulate the laser device with parameters as given in Table 1. One of the most important properties of the QD laser defining the behavior under optical perturbations as well as its mod#204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4871

Table 2. Parameters used in the simulations unless noted otherwise.

Symbol f act N QD aL δ ωES e δ ωQW h δ ωQW

35

ωRO, ΓRO (ns−1 )

ΓRO

Symbol κ Rw loss WGS WES gGS

Value 50ns−1 0.54cm2 ns−1 0.44ns−1 0.55ns−1 230ns−1

rates ×1

25

ωRO

rates ×0.02

30

Value 0.5 1011 cm−2 15 125ns−1 11.3 × 10−11 cm2 ns−1 5.5 × 10−11 cm2 ns−1

20 15 10 5 0 109

1010 1011 1012 1013 −1 −1 inverse electron scattering lifetime τe (s )

Fig. 1. Dependence of the QD laser dynamics on the effective electron scattering lifetimes, as defined in Eq. (14). Shown are the relaxation oscillation angular frequency ωRO (red) and damping ΓRO (blue). The gray dashed vertical lines show the inverse electron lifetimes for the full rates as in Tab. 2 and 50-fold reduced rates, as will be used in Figs. 2 and 3, respectively. The gray shaded area shows the overdamped regime. The pump current was adjusted at each data point to twice the respective threshold current.

ulation capabilities are its relaxation oscillation angular frequency ωRO and damping ΓRO . We therefore determine the dependence of these quantities on the charge-carrier scattering lifetime. We thus define an effective GS scattering lifetime, quantifying the average time needed for the QD GS to be filled with a carrier: −1 cap,in cap,out rel,in rel,out τb = SGS,b + SGS,b + SGS ρES,b + SGS 1 − ρES,b . (14) Figure 1 shows the relaxation oscillation parameters, determined from fits to the intensity timeseries I(t) ≡ |E(t)|2 : I(t) = I 0 + a · exp − (ΓRO + iωRO )t + b · exp(−γt) , (15) where I 0 , a, b, γ are fit parameters, with I 0 the steady state intensity and γ ΓRO describing an underlying slow change of the intensity. The resulting relaxation oscillation parameters are plotted as a function of the inverse effective electron scattering lifetime τe −1 (the scattering lifetime is evaluated at steady-state). The QD laser exhibits pronounced relaxation oscillations both for very slow and very fast scattering rates, which become increasingly damped when the #204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4872

Fig. 2. Bifurcation diagram of the QD laser under optical injection, depending on injection strength K and detuning ∆νinj . The color code denotes the different regimes of dynamics from a phase-locked state (white) over periodic orbits of different periodicity (gray) to irregular or chaotic behavior (black) obtained by direct integration. Thick solid and dashed lines are obtained by numerical path continuation and denote saddle-node on an invariant cycle (SNIC) and Hopf-bifurcations, where the black parts denote bifurcations of the stable phase-locked solution, and light gray denotes bifurcations of an unstable solution. Blue (gray) solid thin lines denote period doubling (PD) bifurcations. (a) Results using the full model. (b) Results using an α-factor. J = 2Jth . The solitary laser effective QW recombiw ≡ [Rw (w + w )]−1 = 0.21ns, and τ = 2ps nation and QD scattering lifetimes are τloss e h loss e (τh = 1ps), respectively.

rates approach the overdamped region around τe −1 ≈ 1011 s−1 (gray shaded area) [25, 36, 37]. The relaxation oscillation parameters for the microscopic rates given in Tab. 2 are denoted by the vertical line labeled “rates×1”. Apart from the relaxation oscillation parameters, another important property determining the dynamics is the amplitude-phase coupling, i.e., the change of the instantaneous electric-field frequency under changes of the charge-carrier number due to refractive index variations. This effect becomes important as soon as electric-field interference effects arise, i.e., under optical perturbations, such as optical feedback [16] and optical injection. Since in our model the amplitude-phase coupling arises directly from the carrier dynamics, we will in the following investigate its dependence on the charge-carrier scattering timescales in an optical injection setup. 3.2.

Optical injection

In an optical injection setup, the laser dynamics depend on the strength of the optical sig1 nal K, as well as the frequency detuning ∆νinj ≡ 2π (ωinj − ω 0 ), where ω 0 denotes the optical frequency of the free-running laser. The dynamics of semiconductor lasers under optical injection is already well understood [4, 5, 7, 10, 37, 38], but now our microscopically based MBBE model [10] that is reduced to a 10-variable system allows for a thorough investigation of an optically injected QD laser by path continuation methods and opens possibilities for a deeper analytical understanding of the system. By direct integration as well as numerical continuation of occurring bifurcations using the

#204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4873

program AUTO07p [32, 39], we create bifurcation diagrams of the optically injected QD laser, which is plotted in Fig. 2(a). Around zero detuning there exists a parameter region for which the laser is phase-locked to the injected signal. This phase-locked solution is delimited by either saddle-node on an invariant cycle (SNIC) or Hopf-bifurcations, both of which lead to periodic solutions outside of that region. These periodic solutions can undergo a series of bifurcations [4], which may lead to solutions of higher periodicity and to irregular or chaotic behavior, denoted by the darker color shading in Fig. 2. The specific solutions within the different regions have been discussed elsewhere [4, 7, 10]. Our focus in this paper is to determine how these regions change if different quantum-dots are used and to understand the conditions that allow the description of the QD laser with an α-factor. The dynamics of the QD laser under optical injection depends crucially on the frequency and damping of its relaxation oscillations [7, 40] and the amplitude-phase coupling [4, 7]. A common way to describe this change of instantaneous laser frequency is the α-factor, which relates the changes in gain and refractive index by a proportionality factor. Such an approach may, however, lead to inaccurate predictions by neglecting the independent dynamics of gain and index change [9, 10]. Therefore, we compare our predictions using Eq. (7) with those obtained when using an α-factor, which we calculate from the changes of gain and frequency under a small injection strength. For the QD structure used here, we obtain αE = 1.42. The index E is used to denote the α-factor stemming from the variation of the electric field due to the injection [41]. The dynamic frequency change is then given by δ ω = αE g, and the calculated bifurcation diagram is shown in Fig. 2(b). A generally good agreement between the two approaches is visible. The period doubling curves span a slightly larger parameter region and are shifted towards higher injection strengths when using an α-factor, while the Hopf and SNIC bifurcations closely resemble those obtained using the full model equations. Compared to our previous results published in [10], the agreement between both descriptions is improved, which is related to the higher relaxation scattering rate in the used QDs, as discussed in the following. Given the carrier-dependent refractive index change, it is now interesting to see how the charge-carrier dynamics influence the bifurcation structure of the QD laser. We therefore simulate the laser with the scattering rates reduced by a factor of 50. Lower scattering rates can occur, e.g., in devices with stronger electronic confinement or at reduced temperature. The chosen slower carrier scattering leads to a regime where the QD states are only weakly coupled to the carrier reservoir, such that while the output power is decreased compared to the faster rates, the relaxation oscillation frequency and damping remain comparable [25,40] (compare the two cases indicated by the vertical lines in Fig. 1). The resulting bifurcation diagram is shown in Fig. 3(a) and 3(b). The dynamics surprisingly exhibits a drastically different behavior, with a more symmetric shape of the bifurcations with respect to the frequency detuning. These modifications can be explained by the reduced variation of off-resonant carrier occupation under optical injection due to the weaker coupling to the resonant QD GS carriers, which then leads to a lower refractive index change. This can be expressed by a lower alpha factor, which evaluated for the reduced scattering rates yields αE = 0.29. Comparing Fig. 3(a) and 3(b) reveals again a shift of the bifurcation structure of the periodic orbit outside of the phase-locked area when assuming an α-factor, which introduces an asymmetry in the bifurcation diagram that is not observed when using the full model. The bifurcations of the locked steady-state solutions, however, again show a very good agreement. The changes in the bifurcation structure introduced by altering the scattering processes reveals a strong dependence of the amplitude-phase coupling on the internal charge-carrier scattering timescales, as evident from the varying asymmetry of the locking regions. Even though two laser devices can have near-identical steady-state gain spectra, their response to optical perturbations can be very different (compare Figs. 2 and 3).

#204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4874

Fig. 3. Bifurcation diagram of the QD laser under optical injection, with capture and relaxation scattering rates reduced by a factor 50 compared to Tab. 2. Cf. Fig. 2. (a), (b) show the results for full losses, (c), (d) for reduced losses. J = 2Jth . Solitary laser full w = 0.15ns; τ = 76ps; τ = 44ps. Solitary laser reduced losses in the reservoir: losses: τloss e h w τloss = 1.5ns; τe = 80ps; τh = 45ps.

Apart from the carrier scattering processes, also the non-radiative losses in the surrounding w in Eq. (5) play an important role for the carrier dynamics. Lower reservoir indicated with rloss losses are associated with reduced threshold currents and improved laser efficiency. In addition to the reduced scattering rates, we now calculate the bifurcation diagram for reduced nonradiative losses in the reservoir, which is shown in Fig. 3(c) 3(d). Here, the modifications to the bifurcation diagram arising from reduced scattering rates become even more pronounced, with a strongly asymmetric locking region and symmetric bifurcations of periodic solutions. Clearly, in this case, the discrepancy between the full modeling and the α-factor approach is very strong. A simple fit of the shape of the locking area in order to retrieve a value for α can therefore not be expected to accurately reproduce the dynamics of periodic solutions and the bifurcations of the phase-locked solution simultaneously. The above discussion shows a strong dependence of the amplitude-phase coupling on the #204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4875

w ×0.1 Rloss

rates ×1

static αE

full losses rates ×0.02

1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 9 10

1010 1011 1012 1013 −1 −1 inverse electron scattering lifetime τe (s )

Fig. 4. Dependence of the static α-factor under optical injection on the scattering lifetimes. Solid and dashed lines show the dependence of αE for full and reduced losses, respectively. The gray shaded area shows the overdamped regime, cf. Fig. 1.

charge-carrier scattering dynamics. We thus evaluate the dependence of the static αE -factor under optical injection, i.e., the alpha-factor that is able to describe the phase locking boundary, on the charge-carrier scattering lifetime. We evaluate αE from the gain and index variation under injection of a small resonant optical signal. The resulting αE is plotted in Fig. 4. A strong dependence of the evaluated αE on the scattering timescales can be seen. Its value approaches ≈ 1.7 for very fast scattering and steadily decreases with lower scattering rates. As seen from the optical injection bifurcation diagrams before, a slower scattering of the charge carriers can help to increase the stability of the QD laser by reducing the amplitude-phase coupling. With decreased charge-carrier losses in the reservoir, on the other hand, αE increases over the whole scattering rate interval. Nevertheless, as previously seen, the increased value of αE for reduced losses does not lead to a strong amplitude-phase coupling for periodic solutions in the optical injection setup, as evident from the symmetric shape of the period doubling regions in Fig. 3. Therefore, in order to fully understand the phase response of the QD laser, we perform a frequency resolved analysis of the laser modulation dynamics in the following. 3.3.

Modulation dynamics

The previous results imply a complex dependence of the amplitude-phase coupling on the charge-carrier dynamic timescales. Furthermore, different types of dynamics vary in their corresponding refractive index dynamics, i.e., the amplitude-phase coupling for steady-state and periodic solutions under optical injection behaves differently. We thus calculate the laser response to different types of external modulation (either optical or electrical), in order to reveal the underlying effects of the changes in the amplitude-phase coupling. The small-signal modulation with frequency f of either the pump current or the electric field leads to a variation of the charge-carrier occupations in the system. These will in turn induce a variation of the electric-field gain and refractive index. The instantaneous electricfield frequency is then modulated as δ ω(t) = δ ω 0 + ∆ω( f ) cos 2π f t + φω0 , (16) with the steady-state frequency shift δ ω 0 , the frequency chirp ∆ω( f ), and a constant phase #204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4876

αX (f) ≡ ||∆∆ωg||

100

(a)microscopic scattering rates current modulation field modulation

10 1

(b) scattering rates ×0.02

αJ (f) α E (f)

0.1

0.010.1

full losses losses ×0.1 1

f (GHz)

10

0.1

1

f (GHz)

10

Fig. 5. Frequency dependence of the amplitude-phase coupling. Shown is the ratio of the frequency chirp to gain modulation depth |∆ω( f )|/|∆g( f )|, under either a direct modulation of the electric field (red), or pump current (blue). The results are plotted for full carrier reservoir losses (solid lines) and for reduced losses (dashed lines), with (a) full scattering, and (b) reduced scattering. J = 2Jth .

shift φω0 . Similarly, the amplitude gain also undergoes a modulation, with ∆g( f ) describing its modulation depth. When assuming the amplitude-phase coupling to be described by an αfactor, the frequency chirp of the laser would be given by ∆ω = α∆g .

(17)

Therefore, we define the frequency response as αX ( f ) ≡ |∆ω( f )|/|∆g( f )|, where f is the modulation frequency. The index X distinguishes between a modulation of either the pump current (X ≡ J), J(t) = J 0 + ∆J cos(2π f t)

(18)

with small ∆J (chosen equal to 0.01 × J 0 ), or a modulation of the the electric field (X ≡ E), modeled by ∂ E(t) = [(g(t) − iδ ω(t)) − κ]E(t) + kE 0 cos(2π f t) , ∂t

(19)

with a small injection rate k (chosen equal to 10−4 ns−1 ). The modulation amplitudes were chosen such that nonlinearities in the response of the laser are suppressed, i.e., only the smallsignal response is evaluated. Figure 5(a) shows αX ( f ) of the QD laser device for either a direct modulation of the cavity field (solid red line) or a modulation of the pump current (solid blue line). For the pump current modulation, this value corresponds to the evaluation of 2β /m in experimental FM/AM , m ≡ 2 |∆E| are the frequency and amplitude modulation measurements [42], where β ≡ ∆ω ω0 E0 coefficients, respectively. It can be seen that the frequency response curve exhibits a strong frequency dependence. Under pump current modulation αJ is very high at low modulation frequencies and decreases with increasing f , until it reaches a minimum around 10GHz before it increases again at high f & 50GHz [30]. The decrease at low frequencies is related to nonlinear gain suppression [43], while the increase at fast modulation is a direct result of the finite charge #204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4877

carrier scattering rate, which limits the maximum frequency up to which the gain modulation can follow the pump current modulation. In terms of α( f ) this means that the amplitude-phase coupling itself is frequency dependent, and modulations at low or very high frequencies will induce a chirp of higher magnitude than at intermediate frequencies. A direct modulation of the electric field, as occurs when considering optical perturbations either by feedback or injection of an external optical signal, leads to a different response of the QD laser (red solid curves, αE , in Fig. 5). The frequency response under such modulation displays a plateau for low modulation frequencies and a steady decrease for faster modulation. The plateau value of αE was used previously in Figs. 2 and 3 to approximate the locking regions in the optical injection setup. It is important to note that up until the cutoff frequency at ≈ 4GHz the gain and frequency modulations are nearly in phase, i.e., gain and frequency variations are proportional. Thus it is possible to accurately characterize the response of the laser under optical perturbations with frequencies below this cutoff value with an α-factor. This explains the good agreement of the bifurcations of steady-state solutions under optical injection between the two descriptions of the amplitude-phase coupling seen before. Modulations at higher frequency, however, will introduce a more complex phase response. The decrease of the phase response at high modulation frequencies is due to the finite scattering lifetime. Only the case of instantaneous charge-carrier scattering would lead to a frequency independent amplitude-phase coupling. Furthermore, αE ( f ) would then coincide with the value of αJ ( fmin ) at the frequency response minimum under pump current modulation. It was suggested to use this minimal value as the device’s α-factor [30]. Generally, however, an evaluation of the frequency response under pump current modulation will only give an upper limit for αE ( f ), and the amplitude-phase coupling under optical perturbations will be weaker than αJ ( fmin ). Thus trying to characterize the laser response with a single α will fail [9, 10, 30, 44]. As before, we want to clarify the impact of the scattering timescales and therefore the frequency dependent analysis is repeated for the reduced scattering rates that we already used in Fig. 3. Figure 5(b) shows the modulation behavior when considering slower scattering rates. The modulation capabilities of the QD laser become worse [25, 45], as the frequency chirp under pump current modulation is increased. This is due to a less efficient modulation of the resonant charge carriers, which increases the effect of the off-resonant carriers on the refractive index. Also, now a very pronounced minimum around f ≈ 4GHz is visible, with a sharp increase for higher frequencies. Under a modulation of the electric field, however, a lower frequency chirp is observed, since here the resonant carriers are directly modulated, and a slower scattering thus leads to a decreased modulation of off-resonant carrier occupation. Slow scattering processes can therefore greatly reduce the frequency chirp of QD lasers under optical perturbations and increase their stability. The modulation dynamics of lasers is also crucially affected by charge-carrier losses. Under modulations of the pump current, a slight reduction of the amplitude-phase coupling can be observed, which is more pronounced for lower modulation frequencies. Under optical modulations, however, the frequency variation is increased for low frequencies, and approaches the value of αE in the full loss case for higher f . While reduced non-radiative losses improve the efficiency of the laser device, they can therefore lead to a stronger amplitude-phase coupling and thus a potentially higher sensitivity to optical perturbations. The observed change of the bifurcation diagram under optical injection with decreased carrier losses can be explained from the changes to the frequency response curve. The low frequency αE is higher at lower losses, which explains the increased asymmetry of the locking region, while it remains the same at higher frequencies, thus leading to a similar bifurcation structure of the periodic solutions.

#204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4878

4.

Summary

We have presented a rate-equation model based on a microscopic approach, that implements the dynamic refractive index changes explicitly from off-resonant charge carriers in a QD laser. This model allows the evaluation of the amplitude-phase coupling without the need to assume a value for the α-factor. The model is thus suited for predicting the chirp of QD lasers under different kinds of modulation and for comparisons with experimentally determined α-factors. We have investigated the dynamic response of QD lasers to optical perturbations, specifically external optical injection, and to small-signal modulation of the pump current and cavity optical field, in dependence of the charge-carrier scattering and loss timescales. We calculated the bifurcation structure of the QD laser under optical injection, and compared the results with predictions when using a fixed α-factor, and found that in the case of fast carrier scattering the QD laser behavior can be reasonably well characterized by an α-factor. For slower scattering timescales, as, e.g., present at low temperatures, the bifurcation structure exhibits differences between the two approaches that become even more pronounced when carrier losses are low. In these cases, the bifurcations of periodic solutions and of the phase-locked solution can then no longer be simultaneously described by a single α. Instead, the bifurcation structure outside the locking region is found to be more symmetric, which is due to a reduced amplitude-phasecoupling at modulation frequencies of the periodic oscillations. The stability of the QD laser towards optical perturbations is thus increased for slower intrinsic time scales of the carrier dynamics, i.e., the minimum strength of optical perturbations required to induce instabilities is increased and the extent of chaotic regions in the parameter space is reduced. Simulations of the small-signal response of the QD laser reveal a strong dependence of the FM/AM response on the modulation frequency. Under pump current modulation, the QD laser reproduces the well-known behavior of QW lasers with a sharp decrease of the FM/AM response at first for increasing modulation frequencies until it reaches a constant plateau. For even higher frequencies, however, an increase is predicted that is especially pronounced for slow scattering rates. For direct modulations of the electric field, as occurring under optical perturbation, the FM/AM response is always smaller than under a modulation of the pump current. It shows a plateau for low modulation frequencies, and a steady decrease for modulation frequencies above a cutoff-frequency that depends on the charge-carrier scattering rates. Thus, for slower scattering between the QD and reservoir states, the overall chirp of the laser is reduced and cuts off at lower frequency, such that the amplitude-phase coupling for high-frequency modulations is significantly reduced. A reduction of charge-carrier losses, on the other hand, leads to an increase of the amplitude-phase coupling under optical perturbations and thus to a higher chirp, while it leads to a slightly reduced frequency variation under pump current modulation. The laser stability under optical perturbations can thus be improved by slow carrier scattering and sufficiently large carrier losses. Acknowledgments This work was supported by Deutsche Forschungsgemeinschaft within SFB 787, and by Sandia’s Solid-State Lighting Science Center, an Energy Frontier Research Center (EFRC) funded by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences. The authors thank E. Sch¨oll for fruitful discussions and careful reading of the manuscript.

#204443 - $15.00 USD (C) 2014 OSA

Received 8 Jan 2014; revised 12 Feb 2014; accepted 13 Feb 2014; published 24 Feb 2014 10 March 2014 | Vol. 22, No. 5 | DOI:10.1364/OE.22.004867 | OPTICS EXPRESS 4879