## Abstract

We investigate the nonlinear propagation of an ultra-short, 150 fs, optical pulse along the waveguide of a quantum dot (QD) laser operating above threshold. We demonstrate that among the various nonlinear processes experienced by the propagating pulse, four-wave mixing (FWM) between the pulse and the two oscillating counter-propagating cw fields of the laser is the dominant one. FWM has two important consequences. One is the creation of a spectral hole located in the vicinity of the cw oscillating frequency. The width of the spectral hole is determined by an effective carrier and gain relaxation time. The second is a modification of the shape of the trailing edge of the pulse. The wave mixing involves first and second order processes which result in a complicated interaction among several fields inside the cavity, some of which are cw while the others are time varying, all propagating in both directions. The nonlinear pulse propagation is analyzed using two complementary theoretical approaches. One is a semi-analytical model which considers only the wave mixing interaction between six field components, three of which propagate in each direction (two cw fields and four time-varying signals). This model predicts the deformation of the tail of the output signal by a secondary idler wave, produced in a cascaded FWM process, which co-propagates with the original injected pulse. The second approach is a finite-difference time-domain simulation, which considers also additional nonlinear effects, such as gain saturation and self–phase modulation. The theoretical results are confirmed by a series of experiments in which the time dependent amplitude and phase of the pulse after propagation are measured using the cross-frequency-resolved optical gating technique.

© 2013 OSA

## 1. Introduction

Ultra-fast nonlinear dynamics in quantum dot (QD) edge emitting lasers depend on a number of phenomena which are often hard to distinguish and identify experimentally. The dynamics are usually studied using time [1–3] or frequency [4,5] domain techniques to probe semiconductor optical amplifiers (SOAs), rather than lasers operating above threshold. However, since the rate of steady state stimulated emission is much higher in lasers compared to SOAs, there are important differences in the carrier and gain dynamics of the two structures. One way to examine directly the ultra-fast dynamics of a laser is by mapping the nonlinear propagation of an ultra-short pulse along the waveguide, while operating above threshold [6,7]. The nonlinear interaction between the pulse and the two counter-propagating cw fields of the laser imprints a distinct signature on the amplitude and spectral content of the pulse. The nature of this signature is determined by a variety of nonlinear processes, the most important ones being four-wave mixing (FWM) [8–11] and self-phase modulation (SPM) [12].

This paper describes results of modeling and experimental characterization of these nonlinear interactions for InAs/InP QD lasers. We observe the creation of a hole in the spectrum of the pulse resulting from FWM among all the fields propagating within the cavity. The spectral hole is located in the vicinity of the cw frequency and its width is determined by an effective relaxation time of the carriers. The time-domain manifestation of this signature is an amplitude tail. In addition, second order FWM processes are predicted to take place within the cavity generating two counter-propagating idlers. One of these secondary idlers co-propagates with the original pulse and contributes to the modification of the trailing edge of the output signal. Similar results were obtained for quantum dash lasers (such as the ones described in [11]) but here we present only QD laser results.

The nonlinear interactions are analyzed using two approaches. First, a semi-analytical model is presented which considers only FWM interaction among six waves which propagate in the cavity. This model isolates the wave mixing process neglecting saturation and SPM. Second, the complete interaction is analyzed using a numerical model based on a finite-difference time-domain (FDTD) simulation of the Maxwell-Schrödinger equations. The QD gain medium is described in the simulation as a sequence of two-level systems, which are fed, under an electrical bias, from a high energy carrier reservoir.

Experimental verification of the model requires the pulse to be characterized using ultra-fast coherent techniques. We have used cross-frequency-resolved optical gating (X-FROG) measurements [13,14] which reveal all the predicted modifications induced on the complex field of the pulse, such as: a distinct hole in the pulse spectrum (in the vicinity of the cw oscillation wavelength) and deformation of the pulse trailing edge. Detection of such complex ultra-fast effects requires separating the pulse response from any steady state background signal using FROG-like sampling techniques [15].

The paper is organized as follows. Section 2 describes two theoretical models. The experimental results are presented in Section 3 together with an analysis and a comparison to the predictions, while Section 4 is devoted to conclusions.

## 2. Theory

The interaction between a propagating pulse and the two counter-propagating cw oscillating fields is analyzed in this section. A semi-analytical FWM model is described in sub-section 2.1 and a comprehensive numerical simulation in 2.2.

#### 2.1. Semi-analytical FWM model

The propagation of a short pulse along a laser diode driven above threshold initiates a complicated FWM process between the pulse and the two counter-propagating cw fields. The formalism we present for describing this interaction is an extension of the models by Agrawal for counter-propagating cw signals [8] and for pulse amplification in a saturated medium [12] as well as of Shtaif for FWM between two co-propagating pulses [9].

Assuming propagation along the waveguide of the laser, parallel to the z-axis, the wave equation for the electric field is:

^{9}

In order to properly account for all the possible FWM products, one needs to consider six field components (as illustrated in Fig. 1): Two counter-propagating cw fields, representing the oscillating mode of the cavity, and serving as the pumps of the FWM interaction; The injected pulse, whose carrier wavelength is defined as the probe wavelength; Two counter-propagating idler signals, carried by a conjugate wavelength; Another signal having the probe wavelength, propagating in the opposite direction relative to the input pulse.

The fields are assumed to propagate adiabatically, with negligible modifications to their (single) lateral mode profile and to their TE polarization. For simplicity, we model the laser fields as two external cw signals, which are injected from both sides of the waveguide with equal frequencies, amplitudes and phases, therefore:

Substitution of the total electric field Eq. (4) into Eq. (3), and applying the small-signal approximation, yields the evolution of several time and space dependent components of the susceptibility:

Using Eq. (5) we can separate Eq. (3) into individual evolution equations:

Eliminating the lateral profile of the fields in the wave Eq. (4), using

The susceptibility is related to the gain and refractive index by:

Equations (6)-(12), together with the matching equations for the down-converting susceptibility components, and Eqs. (14)-(19), are solved numerically in time and length of propagation. As a first stage, the steady-state susceptibility is calculated together with the cw fields. A transform limited Gaussian pulse is introduced after laser oscillations have stabilized. At each time step the pulse and the other wave mixing products are advanced in space and the susceptibility components are evaluated at each point along the waveguide.

The parameters used in the calculations are listed in Table 1. They are set so that the net average modal gain under steady-state conditions (which represents the threshold modal gain of the laser) is about 11 dB, while matching the basic assumptions of the theory.

The evolution of the signals along the waveguide is presented in Fig. 2(a), 2(b) and 2(c). Shown in logarithmic scale are the amplitudes of the injected pulse and the wave mixing products (relative to the pulse peak), calculated when the pulse is located near the input facet, in the middle, and close to the output facet, respectively. The injected pulse (shown in blue) creates, via FWM with the two cw fields, both counter and co-propagating idlers (shown in red and green, respectively). Upon propagation, the idlers interact with the cw fields in a cascaded FWM process which produces additional idlers. One of these (shown in black) propagates towards the left (input) facet, and is marked in the figure as the left propagating pulse. The other secondary idler co-propagates with the injected pulse and contributes to the modification of the trailing edge of the output signal in addition to the pulse deformation originating from the first order FWM process.

Figure 3 presents the total output signal and the various wave mixing products, as obtained at both facets of the cavity, as function of propagation time. The output signal emerges about 12 ps after the original pulse was injected to the cavity. The low level tail is observed once more. The peak of the left propagating idler (shown in red) is calculated at the input (left) facet after it has been amplified along the entire length of the cavity. After it reaches its peak, this idler is not immediately suppressed because of the secondary FWM processes, which mediate continuous transfer of energy between the idler and the pulse signals (in both propagation directions). This is also inferred by the inset of Fig. 3(a), which shows the trailing edge of the total output signal at the right facet, as well as the pulse and idler signals of which it is comprised. In this image, the signals show a cyclic transfer of energy between pulse and idler after the main part of the pulse has passed. Finally, the existence of the secondary FWM process is also implied by the presence of the pulse signal that propagates towards the left facet, and plotted in black in Fig. 3(b). This signal could only be created by FWM between the idler signals and the cw fields, which is a cascaded process.

Figure 4 shows the time-dependent phase of the total output signal. The figure actually shows an inverted phase profile (with respect to the semi-analytical formalism) in order to match the notations used in the FDTD simulation in sub-section 2.2, and in the measurements described in section 3. The phase profile shows a distinct pattern on the trailing edge, followed by an almost linear decrease (folded to fit the $[-\pi ,\pi ]$ domain), indicating a frequency shift of the FWM products relative to the pulse carrier frequency.

The spectral modifications of the output signal are shown in Fig. 5. The output spectra, normalized to an arbitrary background level, are shown for effective gain recovery times of 1 ps, 2 ps and 5 ps. The spectral hole due to FWM mediated wave mixing is consistent with previous results [5,8,17] and is clearly visible at about 1500 nm. It is asymmetric relative to the cw frequency due to the well-known Bogatov effect [18]. The hole width is approximately proportional to the inverse of the gain recovery time, as expected [17,19]. This signature is consistent with the time-domain phase profile, since it represents a relatively narrow-banded signal, which is also shifted in wavelength from the pulse carrier wavelength. The inset shows a spectrum for the case of 5 ps gain recovery time indicating that the spectral hole is created on the short wavelength side of the pulse peak. At the wavelength of the cw signal, the pulse amplitude is reduced from its value at its peak wavelength by about 5dB. The assumption of no spectral overlap between the pulse spectrum and the cw field is not fully satisfied, but since the spectral width of the wave mixing signature is much smaller than 10 nm, the approximation is reasonably accurate.

A calculated X-FROG spectrogram of the output signal (gated by the input pulse) is shown in Fig. 6 (enhanced by power of 1/5 for improved visualization). The long pattern on the right is a trace of the temporal tail, and the distortion at a delay of 400 fs represents the spectral hole due to wave mixing.

#### 2.2. Numerical Model

In order to widen the scope of the modeling and to include additional effects, we developed a comprehensive numerical model, which solves the Maxwell-Schrödinger equations. It uses the density matrix formalism under the dipole approximation together with the rate equations, all of which are solved self-consistently using the FDTD technique [20,21].

Our model calculates the light-matter interaction along the propagation (z) axis and therefore we evaluate the x-component of the electrical field, E_{x}, and the y-component of the magnetic field, H_{y}. Each z-interval contains one QD with a single excited state for the electrons, which is coupled to a carrier reservoir and is biased electrically. We neglect the fact that the field addresses actually more than one QD simultaneously. The carrier reservoir supplies both electrons and holes. The valance band contains several closely placed states which are highly occupied, and from which holes relax to the ground state at a high rate. For describing the dynamical behavior, it is sufficient therefore to consider for holes only the ground state. Spontaneous electron-hole recombination processes are assumed to exist for electrons in the reservoir, in the excited state and in the ground state of the QD. A schematic diagram of the energetic arrangement of the QD is shown in Fig. 7.

The interaction is expressed by calculating the co-evolution of the coupled electromagnetic field and the two-level system as given by the Maxwell-Schrödinger equations; the two-level system interacts with the field through the dipole-moment driven polarization, while the field affects the two-level system through the potential term in the Schrödinger’s equations.

The evaluation of the ground state electron and hole level occupation probability-amplitudes, ${\rho}_{11}\left(t,z\right)$ and${\rho}_{22}\left(t,z\right)$, respectively, and the off-diagonal elements, ${\rho}_{12}\left(t,z\right)$ and${\rho}_{21}\left(t,z\right)$, is calculated using Eqs. (21)–(23), which are written here for simplicity without the time and space notations.

P_{res} and N_{es} are the hole reservoir and the electron excited state carrier densities, respectively. N_{D} is the spatial dot density and D_{res} is the electron reservoir density of states. The total QD density of electrons in the ground state, which is two-fold degenerate, equals$2{N}_{D}{\rho}_{11}$.

${\mu}_{x}$ is the x-component of the QD dipole moment and ω is the resonance angular frequency. ${\gamma}_{c}$, ${\gamma}_{v}$ are the upper and lower level decay rates, while ${\gamma}_{in}$ is the off-diagonal element decay rate which determines the linewidth of the resonant transition.

The additional terms in Eq. (21) are responsible for the interaction with the excited-state and for the coupling to the host reservoir. ${\tau}_{es\to 11}$ and${\tau}_{11\to es}$, are the relaxation and emission time constants between the excited state and the ground state of the electrons. Similarly, ${\tau}_{esc}^{h}$ and ${\tau}_{cap}^{h}$ are the escape and capture time constants to and from the hole reservoir.

The dynamical behavior of the QD and the reservoir are evaluated using a set of rate equations:

_{res}being the reservoir carrier density.

*I*,

*η*and

_{i}*V*are the bias current, injection efficiency and the active layer volume, respectively. ${\tau}_{res}$, is the reservoir lifetime. ${\tau}_{cap}$ and ${\tau}_{esc}$ are the reservoir to dot capture and escape lifetimes which are related, for electrons and holes separately, according to the principle of detailed balance [19].

Applying Maxwell's equations to calculate E_{x}(t,z) and H_{y}(t,z), completes the description of the light-matter interaction.

We distinguish between the polarization which originates from the background material refractive index and the dot polarization as expressed in Eqs. (30)–(32). The background electrical susceptibility is χ_{e}. The confinement factor Г accounts for the difference between the volumes of the optical mode and the gain medium. The internal losses are included through the conductance, σ, which is assumed to be uniform throughout the structure.

The linewidth enhancement factor α is accounted for by formulating a parametric dependence of the refractive index on carrier population:

_{res}and C

_{QD}are chosen so that the index depends only on the reservoir population [22].

The numerical calculation space is arranged in a way that the QD laser cavity is placed in the center of the z-axis computational space and is surrounded by free-space where the pulse excitation occurs. At the end points of the calculation space (which is longer than the laser cavity), absorbing boundary conditions are artificially applied in order to prevent the outgoing *E _{x}* and

*H*fields from being reflected back into the calculation space. In order to achieve both high precision and rapid calculations, the Schrödinger and Maxwell equations were calculated simultaneously using a central-difference method while the rate equations, which are sampled at a significantly higher rate than actually needed, were evaluated using a forward-difference scheme.

_{y}The model does not account for the inhomogeneous broadening of the self-assembled QDs although additional energy levels may be added easily at the expense of increasing the computational effort. In order to accommodate the entire spectrum of the short, 150 fs pulses, we increase the linewidth of the resonant transition to 16.5 meV (30 nm).

Dielectric coatings were added to the laser facets yielding reflectivities of ~1%. The increased mirror losses yield a large threshold gain as in the experiment, which is presented in Section 3.

The simulation proceeded in two stages. Steady-state laser oscillations were first built-up from a random electrical field noise. Following the stabilization of the various transient effects, a short Gaussian pulse was applied which experienced the various nonlinear effects upon propagation; its complex properties were analyzed at the output. The parameters used for the calculation are listed in Table 2.

In the second stage we propagated a 150 fs wide unchirped Gaussian pulse with an energy of 14 pJ through the oscillating medium. Figure 8 shows the calculated multimode laser oscillations (in red, and enlarged in the inset), and the pulse spectrum before (green) and after propagation (blue). The figure shows a shift of the pulse peak wavelength by approximately 2 nm; such a shift was indeed observed in previous experiments [14]. As in the analytical model of sub-section 2.1, the output pulse spectrum exhibits a spectral hole near the wavelength of the oscillating field. The width of the hole is estimated to be about 10 nm which corresponds to an effective response time of about 1 ps, consistent with the QD capture time used in the simulation. Nevertheless, since other nonlinear effects such as SPM take place simultaneously, the width of the hole is not determined solely by the capture rate.

The complex temporal field envelopes (phase and amplitude) of the pulse before and after propagation are shown in Fig. 9. Both the phase and amplitude are presented without the fast oscillating carrier frequency and have been extracted according to the relation:

Another feature of the interaction is that the spectral hole vanishes for high peak pulse peak powers, due to an increased spectral modification, caused by SPM which masks the effects of FWM, as seen in Fig. 10. This behavior is also observed in the experiments described in Section 3.

Finally, in order to be able to compare the numerical simulation to the experiments and the semi-analytical calculation, we plot in Fig. 11 the calculated spectrogram (X-FROG trace) of the emitted pulse, where the trace is enhanced by power of 1/5 to improve the visualization of the low power features. The trace shows a somewhat symmetric shift at delays between −0.2 ps to 0.2 ps which is a well-known signature of SPM [15]. In addition, the trace shows the trailing excitation at delays longer than 0.4 ps which originate from the wave mixing process.

## 3. Experimental results

Time dependent amplitudes and phases, as well as spectra of a 150 fs FWHM pulse, generated by a tunable optical parametric oscillator propagating through a QD laser, were measured by a X-FROG system under various operating conditions. The X-FROG measurement generates a spectrogram (X-FROG trace) of the sum-frequency product (in a thin LiNbO_{3} crystal) of the sampled pulse with a pre-characterized pulse, at different temporal delays between them. The bandwidth of the nonlinear sampling and detection is sufficiently large to easily support the 150 fs pulse namely, it is about 80 nm at the doubled frequency. A phase retrieval technique enables to extract the precise temporal and spectral characteristics of the sampled pulse, which otherwise could not be separated from the background amplified spontaneous emission (ASE).

The laser we tested was a 1.5-mm-long, buried-channel, InAs/InP QD laser with reduced facet reflectivity yielding a high threshold gain. Its bias dependent ASE spectra are depicted in Fig. 12. The wide gain spectrum below threshold results from the inhomogeneously broadened gain. The oscillation threshold is reached at around 75 mA. Short (150 fs) pulses were launched into the laser and propagated along its waveguide in the presence of the cw oscillating fields.

Figure 13(a) and 13(b) show a typical sum-frequency X-FROG trace, and the corresponding X-FROG trace after retrieval, respectively, obtained for a 0.25 pJ input pulse at 1510 nm (which is longer than the cw oscillation wavelength) and a bias of 100 mA. Both traces are enhanced by power of 1/3 for improved visualization of details. The retrieved trace, which is clearer than the measured one, helps in matching the experimental results to the theory. The traces are similar to those predicted by both the semi-analytical calculation and the FDTD simulation (Fig. 6 and 11). The tail and dip in its base, which originate from the FWM process, are circled in red. The weak, trailing signal and the dip are observable at wavelengths longer than the marked tail, and are consistent with the spectrogram obtained by the FDTD simulation. The differences between the measured (retrieved) and the predicted spectrograms originate from the fact that the simulations assume a perfect input pulse while the experiment used an imperfect input pulse which was not purely Gaussian and also contained some chirp. Using the real pulse as an input to the simulation can improve the accuracy of the model. Similar results were also obtained for an input pulse centered at 1470 nm.

The temporal profiles of the retrieved pulses, at various bias levels are presented in Fig. 14(a) and 14(b), for pulses centered at 1470 nm and 1510 nm, respectively. A pulse asymmetry rises with bias level, as the medium gets saturated more deeply. The saturation evolution, with respect to the pulse wavelength, follows the spectral shape of the ASE. Furthermore, for above-threshold bias levels, the predicted tail is visible (magnified in the right inset of Fig. 14(b)), and grows with bias in proportion to the amount of ASE in the spectral domain of the pulse. The time-dependent phase of the pulses is plotted in the insets, on the left of Fig. 14(a) and 14(b). Once more, the predictions of the simulations are confirmed, as the phase rises or falls linearly according to the shift of the oscillating mode frequency relative to the pulse carrier frequency. Thus, we conclude that FWM originated wave mixing is indeed dominant, together with some SPM, in the nonlinear interactions of these experiments.

The retrieved output pulse spectra for a pulse at 1470 nm are shown in Fig. 14(c). For bias levels above threshold, the spectra exhibit the predicted hole in the vicinity of the cw oscillation wavelength. As the pulse is tuned to 1510 nm, the spectral location of the hole remains approximately unchanged (Fig. 14(d)) proving, unequivocally, that it originates from the FWM process with the oscillating cw fields. The width of the spectral hole is 3-4 nm, (~500 GHz), consistent with the reciprocal time constant for carrier relaxation in QD gain media [19]. The nature of the X-FROG system is such that the cw spectral line is not observed; its wavelength is determined separately by an optical spectrum analyzer.

The spectral hole is most significant for a pulse energy of 0.25 pJ, while at larger energies other effects, such as SPM, mask it. This is demonstrated in Fig. 15, which shows the retrieved spectra of 1510 nm pulses, with 50 pJ energies, at various bias levels.

The wave mixing spectral dip is seen; however, it is either shifted or masked, and in general it is significantly less clear compared to the case of the 0.25 pJ pulses. This observation is consistent with the numerical model, Fig. 10.

## 4. Conclusion

To conclude, we analyzed theoretically the properties of an ultra-short pulse, which propagates through a QD laser operating above threshold, and interacts with the oscillating modes. A semi-analytic calculation of the FWM process between the pulse and the oscillating field predicts a dip in the spectrum of the output signal close the cw wavelength. The width of the spectral hole is determined by an effective gain recovery time. It also predicts a tail in the time dependent amplitude which is spectrally narrow and is centered near the oscillating mode frequency; These patterns are predicted as well, and in more detail, by a FDTD simulation, which considers charge carrier dynamics as well as two-level interactions with the electromagnetic fields and includes additional nonlinear effects other than FWM such as saturation induced SPM. However, the FWM interactions are dominant and are easily observed in both the time and frequency domains, as well as in the calculated spectrograms.

Finally, a X-FROG characterization was used to examine experimentally the properties of a 150 fs pulse which propagated along an InAs/InP QD laser. The experimental results are consistent with all the numerical predictions. An effective gain recovery time of about 2 ps was measured, in agreement with known values for such lasers. Many fine details, including small trailing pulses as well as spectral modifications predicted by the models, were confirmed by the X-FROG experiments, leading to the conclusion that FWM, accompanied by SPM, is the dominant nonlinear interaction in these experiments.

## Acknowledgments

This work was partially supported by the Israeli Science Foundation and by the European Commission through the project GOSPEL. The work of the authors from HHI and TUB was funded by DFG in the framework of SFB 787. G.E. would like to thank A. von Humboldt Foundation and DFG/SFB 787 for generous support. O.K acknowledges support of RBNI and A.C thanks the Clore Foundation.

## References and links

**1. **K. L. Hall, J. Mark, E. P. Ippen, and G. Eisenstein, “Femtosecond gain dynamics in InGaAsP optical amplifiers,” Appl. Phys. Lett. **56**(18), 1740–1742 (1990). [CrossRef]

**2. **P. Borri, W. Langbein, J. M. Hvam, F. Heinrichsdorff, M.-H. Mao, and D. Bimberg, “Ultrafast gain dynamics in InAs-InGaAs quantum-dot amplifiers,” IEEE Photon. Technol. Lett. **12**(6), 594–596 (2000). [CrossRef]

**3. **A. Capua, G. Eisenstein, and J. P. Reithmaier, “A nearly instantaneous gain response in quantum dash based optical amplifiers,” Appl. Phys. Lett. **97**(13), 131108 (2010). [CrossRef]

**4. **A. Uskov, J. Mork, and J. Mark, “Wave mixing in semiconductor laser amplifiers due to carrier heating and spectral-hole burning,” IEEE J. Quantum Electron. **30**(8), 1769–1781 (1994). [CrossRef]

**5. **A. Capua, V. Mikhelashvili, G. Eisenstein, J. P. Reithmaier, A. Somers, A. Forchel, M. Calligaro, O. Parillaud, and M. Krakowski, “Direct observation of the coherent spectral hole in the noise spectrum of a saturated InAs/InP quantum dash amplifier operating near 1550 nm,” Opt. Express **16**(3), 2141–2146 (2008). [CrossRef]

**6. **J. Zimmermann, S. T. Cundiff, G. von Plessen, J. Feldmann, M. Arzberger, G. Böhm, M.-C. Amann, and G. Abstreiter, “Dark pulse formation in a quantum-dot laser,” Appl. Phys. Lett. **79**(1), 18–20 (2001). [CrossRef]

**7. **C. Sun, B. Golubovic, H. Choi, C. A. Wang, and J. G. Fujimoto, “Femtosecond investigations of spectral hole burning in semiconductor lasers,” Appl. Phys. Lett. **66**(13), 1650–1652 (1995). [CrossRef]

**8. **G. P. Agrawal, “Population pulsations and nondegenerate four-wave mixing in semiconductor lasers and amplifiers,” J. Opt. Soc. Am. B **5**(1), 147–159 (1988). [CrossRef]

**9. **M. Shtaif and G. Eisenstein, “Analytical solution of wave mixing between short optical pulses in a semiconductor optical amplifier,” Appl. Phys. Lett. **66**(12), 1458–1460 (1995). [CrossRef]

**10. **A. Mecozzi and J. Mork, “Saturation effects in nondegenerate four-wave mixing between short optical pulses in semiconductor laser amplifiers,” IEEE J. Sel. Top. Quantum Electron. **3**(5), 1190–1207 (1997). [CrossRef]

**11. **A. Bilenca, R. Alizon, V. Mikhelashvili, D. Dahan, G. Eisenstein, R. Schwertberger, D. Gold, J. P. Reithmaier, and A. Forchel, “Broad-band wavelength conversion based on cross-gain modulation and four-wave mixing in InAs-InP quantum-dash semiconductor optical amplifiers operating at 1550 nm,” IEEE Photon. Technol. Lett. **15**(4), 563–565 (2003).

**12. **P. Agrawal and N. A. Olsson, “Self-phase modulation and spectral broadening of optical pulses in semiconductor laser amplifiers,” IEEE J. Quantum Electron. **25**(11), 2297–2306 (1989). [CrossRef]

**13. **S. Linden, H. Giessen, and J. Kuhl, “XFROG — A new method for amplitude and phase characterization of weak ultrashort pulses,” Phys. Status Solidi, B Basic Res. **206**(1), 119–124 (1998). [CrossRef]

**14. **A. Capua, A. Saal, O. Karni, G. Eisenstein, J. P. Reithmaier, and K. Yvind, “Complex characterization of short-pulse propagation through InAs/InP quantum-dash optical amplifiers: from the quasi-linear to the two-photon-dominated regime,” Opt. Express **20**(1), 347–353 (2012). [CrossRef]

**15. **R. Trebino, *Frequency-Resolved Optical Gating: The Measurement of Ultrashort Laser Pulses* (Kluwer Academic Publishers, Norwell, 2002).

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

**17. **M. Shtaif and G. Eisenstein, “Noise properties of nonlinear semiconductor optical amplifiers,” Opt. Lett. **21**(22), 1851–1853 (1996). [CrossRef]

**18. **A. P. Bogatov, P. G. Eliseev, and B. N. Sverdlov, “Anomalous interaction of spectral modes in a semiconductor laser,” IEEE J. Quantum Electron. **11**(7), 510–515 (1975). [CrossRef]

**19. **A. Bilenca and G. Eisenstein, “On the noise properties of linear and nonlinear quantum-dot semiconductor optical amplifiers: the impact of inhomogeneously broadened gain and fast carrier dynamics,” IEEE J. Quantum Electron. **40**(6), 690–702 (2004). [CrossRef]

**20. **G. M. Slavcheva, J. M. Arnold, and R. W. Ziolkowski, “FDTD simulation of the nonlinear gain dynamics in active optical waveguides and semiconductor microcavities,” IEEE J. Sel. Top. Quantum Electron. **10**(5), 1052–1062 (2004). [CrossRef]

**21. **J. E. Kim, E. Malic, M. Richter, A. Wilms, and A. Knorr, “Maxwell–bloch equation approach for describing the microscopic dynamics of quantum-dot surface-emitting structures,” IEEE J. Quantum Electron. **46**(7), 1115–1126 (2010). [CrossRef]

**22. **J. Kim, C. Meuer, D. Bimberg, and G. Eisenstein, “Role of carrier reservoirs on the slow phase recovery of quantum dot semiconductor optical amplifiers,” Appl. Phys. Lett. **94**(4), 041112 (2009). [CrossRef]