Spectral signature of nonlinear effects in semiconductor optical amplifiers

Optical spectra of signals at the output of semiconductor optical amplifiers (SOA) provide useful insight into amplifier nonlinearities. In this work, we determine the parameters of an analytical SOA model with a pump-probe experiment by evaluating the measured spectra of the pump and probe pulses at the SOA output. The analytical lumped SOA model considers carrier depletion, carrier recovery, spectral hole burning, two-photon absorption, and we include an additional effect termed ‘two-photon induced free-carrier absorption’, that is responsible for creating an identifiable blue-shifted component in the spectra. We are able to relate the underlying physical nonlinear effects to the spectral peculiarities of the output pump and probe spectra, and give guidelines for the exploitation of these nonlinear effects for optical signal processing. In addition, with a much-simplified SOA model and by replacing the pump pulse with modulated data we show that the output spectrum is altered in a manner consistent with phase patterning effects. ©2017 Optical Society of America OCIS codes: (250.5980) Semiconductor optical amplifiers; (190.4360) Nonlinear optics, devices. References and links 1. B. Filion, W. C. Ng, A. T. Nguyen, L. A. Rusch, and S. Larochelle, “Wideband wavelength conversion of 16 Gbaud 16-QAM and 5 Gbaud 64-QAM signals in a semiconductor optical amplifier,” Opt. Express 21(17), 19825–19833 (2013). 2. S. T. Naimi, S. P. Ó. Duill, and L. P. Barry, “All optical wavelength conversion of Nyquist-WDM superchannels using FWM in SOAs,” J. Lightwave Technol. 33(19), 3959–3967 (2015). 3. B. Filion, L. Jiachuan, A. T. Nguyen, X. Zhang, S. LaRochelle, and L. A. Rusch, “Semiconductor optical amplifier-based wavelength conversion of Nyquist-16QAM for flex-grid optical networks,” J. Lightwave Technol. 34(11), 2724–2729 (2016). 4. A. P. Anthur, R. Zhou, S. O’Duill, A. J. Walsh, E. Martin, D. Venkitesh, and L. P. Barry, “Polarization insensitive all-optical wavelength conversion of polarization multiplexed signals using co-polarized pumps,” Opt. Express 24(11), 11749–11761 (2016). 5. A. P. Anthur, R. T. Watts, K. Shi, J. O. Carroll, D. Venkitesh, and L. P. Barry, “Dual correlated pumping scheme for phase noise preservation in all-optical wavelength conversion,” Opt. Express 21(13), 15568–15579 (2013). 6. G. Contestabile, “Wavelength conversion of PAM signals by XGM in SOAs,” in Asia Communications and Photonics Conference, OSA Technical Digest (Optical Society of America, 2016), paper AF4H.4. 7. C. T. Hultgren and E. P. Ippen, “Ultrafast refractive index dynamics in AlGaAs diode laser amplifiers,” Appl. Phys. Lett. 59(6), 635–637 (1991). 8. T. Vallaitis, C. Koos, R. Bonk, W. Freude, M. Laemmlin, C. Meuer, D. Bimberg, and J. Leuthold, “Slow and fast dynamics of gain and phase in a quantum dot semiconductor optical amplifier,” Opt. Express 16(1), 170–178 (2008). 9. C. Dorrer and I. Kang, “Simultaneous temporal characterization of telecommunication optical pulses and modulators by use of spectrograms,” Opt. Lett. 27(15), 1315–1317 (2002). 10. M. G. Kane, I. Glesk, J. P. Sokoloff, and P. R. Prucnal, “Asymmetric optical loop mirror: analysis of an alloptical switch,” Appl. Opt. 33(29), 6833–6842 (1994). 11. J. Leuthold, P. A. Besse, J. Eckner, E. Gamper, M. Dulk, and H. Melchior, “All-optical space switches with gain and principally ideal extinction ratios,” IEEE J. Quantum Electron. 34(4), 622–633 (1998). 12. C. H. Henry, “Theory of the linewidth of semiconductor lasers,” IEEE J. Quantum Electron. 18(2), 259–264 (1982). 13. J. Wang, A. Maitra, C. G. Poulton, W. Freude, and J. Leuthold, “Temporal dynamics of the alpha factor in semiconductor optical amplifiers,” J. Lightwave Technol. 25(3), 891–900 (2007). Vol. 25, No. 24 | 27 Nov 2017 | OPTICS EXPRESS 29526 #300436 https://doi.org/10.1364/OE.25.029526 Journal © 2017 Received 27 Jun 2017; revised 23 Sep 2017; accepted 24 Sep 2017; published 13 Nov 2017 14. M. Willatzen, A. Uskov, J. Mørk, H. Olesen, B. Tromborg, and A.-P. Jauho, “Nonlinear gain suppression in semiconductor lasers due to carrier heating,” IEEE Photonics Technol. Lett. 3(7), 606–609 (1991). 15. J. Mørk and A. Mecozzi, “Response function for gain and refractive index dynamics in active semiconductor waveguides,” Appl. Phys. Lett. 65(14), 1736–1738 (1994). 16. S. Krishnamurthy, Z. G. Yu, L. P. Gonzalez, and S. Guha, “Temperature and wavelength-dependent two-photon and free-carrier absorption in GaAs, InP, GaInAs, and InAsP,” J. Appl. Phys. 109(3), 033102 (2011). 17. P. Y. Yu and M. Cardona, Fundamentals of Semiconductors (Springer, 2010), Ch. 5. 18. E. M. Conwell and M. O. Vassell, “High-field distribution function in GaAs,” IEEE Trans. Electron Dev. 13(1), 22–26 (1966). 19. G. Klimeck, R. C. Bowen, T. B. Boykin, and T. A. Cwik, “sp3s* tight-binding parameters for transport simulations in compound semiconductors,” Superlattices Microstruct. 27(5–6), 519–524 (2000). 20. P. Lunnemann, S. Ek, K. Yvind, R. Piron, and J. Mørk, “Nonlinear carrier dynamics in a quantum dash optical amplifier,” New J. Phys. 14(1), 013042 (2012). 21. J. Mørk and A. Mecozzi, “Theory of the ultrafast optical response of active semiconductor waveguides,” J. Opt. Soc. Am. B 13(8), 1803–1816 (1996). 22. G. 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). 23. A. Mecozzi and J. Mørk, “Theory of heterodyne pump–probe experiments with femtosecond pulses,” J. Opt. Soc. Am. B 13(11), 2437–2452 (1996). 24. 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). 25. M. Premaratne, D. Nesic, and G. P. Agrawal, “Pulse amplification and gain recovery in semiconductor optical amplifiers,” J. Lightwave Technol. 26(12), 1653–1660 (2008). 26. I. D. Rukhlenko, M. Premaratne, and G. P. Agrawal, “Nonlinear Silicon Photonics: Analytical Tools,” IEEE J. Sel. Top. Quantum Electron. 16(26), 200–215 (2010). 27. J. E. Carroll, J. Whiteaway, and D. Plumb, Distributed Feedback Semiconductor Lasers (Institute of Electrical Engineers, 1998). 28. J. M. Tang and K. A. Shore, “Strong picosecond optical pulse propagation in semiconductor optical amplifiers at transparency,” IEEE J. Quantum Electron. 34(7), 1263–1269 (1998). 29. A. Mecozzi and J. Mørk, “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). 30. M. Sheik-Bahae and E. W. Van Stryland, “Optical nonlinearities in the transparency region of bulk semiconductors” in Nonlinear Optics in Semiconductors I, E. Garmire and A. Kost, eds. (Academic, 1999). 31. D. Hutchings and E. Van Stryland, “Nondegenerate two-photon absorption in zinc blende semiconductors,” J. Opt. Soc. Am. B 9(11), 2065–2074 (1992). 32. B. Olson, M. Gehlsen, and T. Boggess, “Nondegenerate two-photon absorption in GaSb,” Opt. Commun. 304, 54–57 (2013). 33. C. H. Henry, R. A. Logan, and K. A. Bertness, “Spectral dependence of the change in refractive index due to carrier injection in GaAs lasers,” J. Appl. Phys. 52(7), 4457–4461 (1981). 34. R. Giller, R. J. Manning, and D. Cotter, “Gain and phase recovery of optically excited semiconductor optical amplifiers,” IEEE Photonics Technol. Lett. 18(9), 1061–1063 (2006). 35. I. Kang and C. Dorrer, “Measurements of gain and phase dynamics of a semiconductor optical amplifier using spectrograms,” in Optical Fiber Communication Conference, OSA Technical Digest Series (Optical Society of America, OFC 2004), paper MF43. 36. R. S. Grant and W. Sibbett, “Observations of ultrafast nonlinear refraction in an InGaAsP optical amplifier,” Appl. Phys. Lett. 58(11), 1119–1121 (1991). 37. Y. Ueno, S. Nakamura, and K. Tajima, “Nonlinear phase shifts induced by semiconductor optical amplifiers with control pulses at repetition frequencies in the 40-160-GHz range for use in ultrahigh-speed all-optical signal processing,” J. Opt. Soc. Am. B 19(11), 2573–2589 (2002). 38. O. Raz, J. Herrera, and H. J. Dorren, “Enhanced 40 and 80 Gb/s wavelength conversion using a rectangular shaped optical filter for both red and blue spectral slicing,” Opt. Express 17(3), 1184–1193 (2009). 39. J. Leuthold, R. Ryf, D. N. Maywar, S. Cabot, J. Jaques, and S. S. Patel, “Nonblocking all-optical cross connect based on regenerative all-optical wavelength converter in a transparent demonstration over 42 nodes and 16800 km,” J. Lightwave Technol. 21(11), 2863–2870 (2003). 40. M. L. Nielsen, B. Lavigne, and B. Dagens, “Polarity-preserving SOA-based wavelength conversion at 40 Gbit/s using bandpass filtering,” Electron. Lett. 39(18), 1334–1335 (2003). 41. Y. Liu, E. Tangdiongga, Z. Li, H. de Waardt, A. M. J. Koonen, G. D. Khoe, X. Shu, I. Bennion, and H. J. S. Dorren, “Error-free 320-Gb/s all-optical wavelength conversion using a single semiconductor optical amplifier,” J. Lightwave Technol. 25(1), 103–108 (2007). 42. J. Leuthold, D. M. Marom, S. Cabot, J. J. Jaques, R. Ryf, and C. R. Giles, “All-optical wavelength conversion using a pulse reformatting optical filter,” J. Lightwave Technol. 22(1), 186–192 (2004). 43. J. Wang, A. Marculescu, J. Li, P. Vorreau, S. Tzadok, S. Ben Ezra, S. Tsadka, W. Freude, and J. Leuthold, “Pattern effect removal technique for semiconductor-optical-amplifier-based wavelength conversion,” IEEE Photonics Technol. Lett. 19(24), 1955–1957 (2007). 44. X. Wei and J. Leuthold, “Relation between vestigial-sideband filtering and pi/2 progressive phase shift,” Opt. Lett. 29(14), 1599–1601 (2004). Vol. 25, No. 24 | 27 Nov 2017 | OPTICS EXPRESS 29527


Introduction
Semiconductor optical amplifiers (SOA) are key components for both, linear amplification in access networks, and nonlinear ultrafast all-optical signal processing such as wavelength conversion, or all-optical logic gating.Cross gain modulation (XGM) and cross phase modulation (XPM) are the techniques of choice for signals with intensity modulation; the recent deployment of advanced modulation formats such as quadrature phase shift keying (QPSK), 16-state quadrature amplitude modulation (16-QAM) and 64-QAM has motivated studies into using four-wave mixing (FWM) effects in SOAs to perform wavelength conversion [1][2][3][4].The power penalty for employing FWM in SOAs can be well below 1 dB [5].The information carrying capacity of intensity-modulated signals can be doubled by increasing the number of amplitude levels to four by using four-level pulse amplitude modulation (PAM-4), and XGM has been studied to perform wavelength conversion of PAM-4 signals in SOAs [6].Given the huge interest in all-optical signal processing using SOAs and the potential benefits arising thereof, a proper design of nonlinear SOAs requires models that allow to derive the slow and ultrafast gain and refractive index dynamics, where a multitude of nonlinear effects act on different timescales.
There are several methods to characterize nonlinear effects in an SOA.Unfortunately, all of them are experimentally quite demanding: Pump-probe sampling techniques together with heterodyne detection allow phase and amplitude resolutions that are only limited by the full width at half maximum (FWHM) of the pump pulse [7,8].Alternatively, indirect methods such as frequency-resolved electro-absorption gating (FREAG) [9] can be employed.The temporal resolution is limited by the rise time of the optical gate, and computationally expensive algorithms are required to resolve the gain and index dynamics.In addition, various other techniques, where devices are positioned in an interferometer and the amplitude and phase response are derived from the interferometer pattern, have been elaborated [10,11].
In this paper, we derive the temporal gain and phase dynamics of SOAs operated in the nonlinear regime from measured time-averaged optical spectra at the SOA output.We are including the following nonlinear effects in SOAs [12][13][14][15]: bandfilling that treats both carrier depletion and recovery, carrier heating, spectral hole burning, two-photon absorption, and we introduce an effect termed two-photon absorption induced free-carrier absorption (FCATPA) [16] that accounts for the large blue-shifted component especially in the probe spectrum.The measured output spectra reflect the temporal dependence of the output pulse and continuous wave (CW) probe that undergoes cross-gain (and cross-phase) modulation.By fitting a closed-form model to the measured spectra, we can separate the gain and refractive index dynamics of the SOA as well as extracting physically meaningful model parameters.We explain how the nonlinear processes have their own unique spectral signatures; we pay special attention to the blue-shifted components of the probe spectrum.We identify that a distinctive feature in the blue-shifted sideband of the probe spectrum indicates the occurrence of the FCA TPA process.This blue-shifted component is recognizable by mere inspection of the probe spectrum and is exploitable for ultrafast all-optical signal processing purposes.Only modest measurement apparatus is required comprising of a picosecond pump laser source, a CW probe laser, and an optical spectrum analyzer (OSA); this equipment is standard within optics and photonics research laboratories.Neither high-speed electronic modulation, nor detection, sampling, or all-optical probing with femtosecond pulses is required.Moreover, in a wavelength conversion experiment with a data signal and a probe light, the probe spectrum after the SOA provides information about phase patterning effects.
The paper is organized as follows: In Section 2, we derive the nonlinear model of the gain and refractive index dynamics.These dynamics are formulated in terms of differential equations with boundary or initial conditions that are known to have unique solutions.In the following Section 3, we compare the model to the measured spectra and extract the model parameters, which describe the relevant nonlinear dynamic processes within the SOA.In Section 4, we use this knowledge to design filters for all-optical wavelength conversion (AOWC).Finally, in Section 5, we show that patterning effects in AOWC can be predicted from a signature analysis of the converted signal spectrum.

Model of SOA gain and refractive index dynamics
The SOA model in this work relates to a typical experimental condition.We evaluate the SOA gain and refractive index dynamics by the mutual influence of a strong 2.8 ps wide input pump pulse and a weak continuous wave (CW) input probe light.Pump and probe wavelengths are closely neighbored within the gain bandwidth of the SOA.If the pump saturates the SOA gain, the local charge carrier density reduces, which in turn changes the temporal evolution of the pump signal, Fig. 1(a), and the probe light, Fig. 1(b), at the SOA's output.Changes of the modal gain g and of the refractive index n with respect to the carrier density N are interrelated by the Kramers-Kronig relations, which can be approximated by the linewidth enhancement factor (Henry factor, α-factor) [12] where 0 ω stands for the closely neighbored pump signal (subscript S) and probe (subscript ≈ , and c is the vacuum speed of light.The changes in refractive index and modal gain are denoted by Δn and Δg , respectively.

General description
Before explaining the details of the model, we review the most relevant gain dynamics and gain-saturating processes in SOAs which occur on vastly different time scales [13].Different α-factors can be attributed to these processes.For simplicity, carrier diffusion along the SOA is neglected.
(i) Band filling (BF).BF relates to the total number of carriers within the energy bands of the SOA.The dynamics comprises carrier depletion (CD) followed by subsequent carrier recovery (CR).The first effect is fast, because stimulated recombination is virtually instantaneous on the timescale of the pulse duration.Carrier depletion increases the refractive index.The recovery, which implies refilling of electrons and holes back to the concentration of the small-signal gain, takes place on a timescale of about 100 ps , which is much longer than the pump pulse duration.The α-factor associated with BF is denoted by BF α .In some publications, the symbol N α is used instead.
(ii) Spectral hole burning (SHB).The pump signal depletes the carriers virtually instantaneously at the pump photon energy so that a "hole" is burned into the gain spectrum.The empty states of the band are refilled through carrier-carrier scattering on a time scale of about 100 fs [14].The change in gain due to SHB has a symmetric spectral profile around the center frequency S ω of the pump and does not change the refractive index [15].Because the probe frequency ω × is close to S ω , the spectral symmetry applies also to the probe, and we therefore neglect the probe's SHB contribution to the refractive index.Consequently, the associated SHB α factor is assumed to be zero.
(iii) Carrier heating (CH).If carrier equilibrium inside a band is destroyed by rapid processes, the instantaneous carrier occupation probability cannot be described anymore with a Fermi function.However, after refilling the unoccupied states left by SHB, the carrier occupation probability follows again a Fermi distribution with an effective "hot" carrier temperature larger than the lattice temperature.The process is named carrier heating.By heating the carrier distribution, the gain decreases [14] and the refractive index increases.Subsequently, the hot carriers cool down (carrier cooling, CC), and the gain relaxes back on a timescale of up to a picosecond (carrier heating relaxation time) [14].The gain changes lead to refractive index changes, and the associated α -factor is denoted by CH α .
Two photons with the same or with different energies can be absorbed simultaneously.This process is rather unlikely, and therefore large powers (large photon numbers) are needed for this effect to be significant.Electrons are lifted virtually instantaneously from the valence band (VB) to the conduction band (CB).The TPA process can occur with simultaneous absorption of two pump photons (degenerate TPA), or with the simultaneous absorption of a pump and of a probe light photon (nondegenerate TPA).
For nondegenerate TPA, the polarizations of the optical fields play an important role, such that even for a weak probe light the corresponding TPA absorption coefficient can become twice as large as for degenerate TPA.Although a nondegenerate TPA-process may enhance the absorption coefficient for the probe light, the additional number of photoexcited electrons is too small for a significant change of the refractive index if the probe light is sufficiently weak.The associated increased optical loss with degenerate and nondegenerate TPA and the increased number of carriers lead to an instantaneous decrease of the refractive index [7].The α -factor associated with the TPA process is thus negative and will be denoted by TPA α .Carriers generated via TPA can become available for free-carrier absorption (FCA), we denote this process FCA TPA , and explain it as follows: For sufficiently strong pulses and semiconductors of suitably small bandgap, two-photon absorption generates an appreciable population of photoexcited electrons in the Γ -valley of the lowest conduction band (CB) [16].If this Γ-valley is filled or if the sum energy of the two photons is large enough, lower unoccupied states in side-valleys [17], e.g., in the X-valley, can be filled.In this case, however, an additional phonon has to deliver the required difference in crystal momentum.This is also true if electrons in Γ -valley states are to be scattered into the same-energy X-valley states [17,18].The transfer of the photoexcited electrons from the Γ -valley to the satellite Xvalley needs a finite amount of time 0 τ [18], usually termed intervalley scattering time.
Since the bottom of the X-valley lies close [19] to the next higher CB, electrons which reached this side-valley can absorb one more photon (direct absorption) and make a momentum-conserving upwards transition to the next higher conduction band without the assistance of a phonon.This direct, TPA-induced free-carrier absorption (FCA TPA ) has been demonstrated [16] for the binary compounds GaAs and InP.Similar to TPA, TPAinduced FCA increases the optical loss so that the refractive index decreases [20].However, this refractive index change occurs delayed by the intervalley scattering time 0 τ .An instantaneous FCA TPA -effect requires a simultaneous absorption of a photon to conserve both, energy and momentum, and can be neglected [16].A pictorial explanation of both instantaneous and delayed FCA TPA -processes is given in Appendix A. Intervalley scattering effect may well have other significance for SOAs.In previous work from our group [13] intervalley scattering would explain for the required higher-order correction term applied to the carrier density rate equation in [13], Eq. ( 13), whereby carrier leakage at high carrier energies was necessary to get a good fitting with experiments.
The aforementioned effects can be identified when looking at the evolution of the input pump signal (subscript S ) and the input CW probe light (subscript × ) inside the SOA.We we write the normalized root-mean-squared (RMS) electric fields as , , exp j , exp j .
The SOA length is denoted by L .We attach the superscript in (out) to the functions of time at the SOA input (output) at position 0 z = ( z L = ).Then the input and output normalized RMS fields are given by and respectively.We further assume that the probe light power is much smaller than the pump light power so that four-wave mixing between pump and probe remains insignificant.In this case coherent interaction is excluded so that the phase relation between pump and probe is unimportant, and the assignments ( ) in , 0 S t γ ϕ × = for both the pump and the probe are justified.
Pump wave and probe wave (influenced by the pump) have complex valued envelopes with a typical width of 2.8 ps , ( S S S z t P z t z t When pump and probe propagate simultaneously through the amplifier, the response of the nonlinear medium depends on the powers ( ) While the dominant coordinate-independent part is virtually the same for both signals, the coordinate-dependent contributions ( ) ( ) are in general different because the SOA nonlinearity couples the two optical fields through a nondegenerate TPA process.We also assume that both signals propagate with the same group velocity g v , have a negligible group-velocity dispersion [21], and the wavelengths of both pump and probe are close to the wavelength of the SOA gain peak so that SOA gain-dispersion can be ignored, i.e., , , 0 For simplicity, throughout this paper we neglect the internal loss coefficient int α .In the slowly-varying envelope approximation [22], the complex-valued envelope function It is convenient to formulate the propagation of the light signals in a retarded time frame Since the retarded time τ at the SOA input ( 0 z = ) coincides with the natural time t , the input pump and CW probe power dependencies as well as the input phase shifts are At the SOA output ( z L = ) we have the output pump and CW probe power dependencies as well as the output phase shifts We now perform the coordinate transformations ( ) =   , and obtain the propagation equations for the powers and the phase shifts, ( and respectively.Integrating Eq. ( 11) over the SOA length leads to The power gain coefficient ( ) in Eq. ( 11) is positive and takes care of the fact that the signal power is distributed over a cross-section so that the total power is larger by the field concentration factor [23].The highest value 0 g represents the unsaturated small-signal gain, see Eq. (6).The dimensionless integrated gain coefficient

( )
, S h τ × (logarithmic power gain) is also positive and determines the amplification of the respective optical signal The integrated unsaturated gain coefficient 0 0 h g L = corresponds to the largest amplification.
It is convenient to subtract 0 h from

( )
, S h τ × and introduce the integrated gain decrease From Eq. ( 12) we find the solution ( ) for the spatially integrated phase change of the output light with respect to the input light, ( ) The nonlinear output phase shift follows from Eqs. ( 9) and ( 16), In Fig. 1 we display a schematic of the dependencies listed in Eqs. ( 9), (10) and (16).The pump pulse ( ) in S P τ at the input has a symmetric shape, see red curve shown in Fig. 1(a).
While propagating through the SOA, the leading edge of the pump pulse experiences a larger gain than the trailing edge, due to gain saturation.The pump pulse emerging at the SOA output, which is depicted in blue color in Fig. 1(a), has therefore an asymmetric shape with its peak power shifted to smaller τ , this has been documented in [22].The gain dynamics which results from the previously mentioned effects is probed with a weak CW input light having constant power P × , which becomes modulated by the instantaneous gain.The power dependency of the modulated probe light ( ) out , ( ) 1(b).The gain decrease as induced by the pump pulse causes the falling slope of out ( ) P τ × .While CD occurs over the whole pump pulse duration, only the processes dominant in each temporal segment are explicitly marked.The probe power minimum is reached after the (input) pump pulse maximum at 0 τ = , because CD, TPA and SHB are still in effect.Eventually, the SOA gain recovers, and subsequently the probe power increases towards its initial value.
The time dependency of the change in the refractive index ( ) is shown in Fig. 1(c).The combined effects of CD (Section 2.1(i)) and TPA (Section 2.1(iv)) reduce the gain.However, the corresponding positive (CD) and negative changes (TPA) in the refractive index tend to balance, resulting in a time-independent region (CD + TPA) in Fig. 1(c).The onset of SHB in Fig. 1(b) reduces the gain further, but has no effect on the refractive index change, Fig. 1(c) (Section 2.1(ii)).While all nonlinear effects in Fig. 1(b) are of first order and FCA can be neglected, the index change in Fig. 1(c) due to FCA TPA decreases quadratically with the pump pulse power, so that TPA-induced free-carrier absorption (FCA TPA ) is responsible (Section 2.1(iv)) for the pronounced dip in the refractive index curve.During the ultrafast gain recovery (CH) the refilling of the spectral hole counteracts CD and increases the gain.Consequently, the SOA refractive index increases (Section 2.1(iii)).Once the pump signal leaves the SOA, the carrier distributions start resuming their intraband equilibrium states (CC).Finally, CR through the injection current begins and re-establishes the original carrier concentration.

Description of the analytical model
We now describe the gain and phase evolution by an analytical model in which the probe power can be considered negligible.The pump power ( ) , S P z τ influences the amplification and phase shift for both pump and probe light.According to the discussion of Section 2.1(iv), a weak probe light plays an important role only for TPA processes.
As a consequence, we will treat the other processes (BF, SHB and CH) as if they had occurred in the absence of the probe, and therefore omit the subscript S for all dynamical quantities except ( ) , S P z τ .We emphasize that the model described below applies only to bulk SOAs.
The local modal gain coefficient ( ) , g z τ comprises the small signal gain coefficient 0 g and the coefficients of the gain changes from the contributions of the processes labelled by " x ", i.e.,

L g z g g z h h h
In analogy to the general relation Eq. ( 1), the x α -factors associated with the above processes are denoted by The x α -factors do not depend on time or space, because the functional dependencies of We substitute Eq. ( 20) in the right-hand side of Eq. ( 17) in order to determine the output phase shift of the pump field.We denote the contribution of the various processes by Further, taking into account that ( ) ( ) We now start with the analytical description for band filling as developed by Agrawal and Olsson [22], and then we then include the ultrafast gain saturation as well as the TPA effect in a phenomenological fashion [21,24].In the following, we calculate the individual gain changes and the integrated coefficients ( ) x h τ ∆ of Eq. ( 18) as well as the quantities ∆ related by Eq. ( 21).Spatial integration techniques to account for nonlinear effects were previously described in [25] for gain recovery in SOAs, and in [26] for analyzing active and passive silicon waveguides.Such techniques are very different from the traditional discretization approach used in the context of distributed feedback semiconductor lasers [27].

Logarithmic gain
The description of the carrier dynamics in the active medium is based on rate equations [28,29].

Band filling
We assume a linear dependence of the band filling gain coefficient where Γ denotes the field confinement factor, a is the differential gain and st N represents the unsaturated value of the carrier density [24].In this paper we include Γ, a and st N in the unsaturated (small-signal) gain 0 g , the saturation power [24] sat P , and in the spontaneous carrier lifetime s τ .We now put down a differential equation that describes the temporal evolution of band filling during a pump-probe experiment with a negligible probe light power.From the differential equation for the charge carrier concentration we find the differential equation [28] for Band filling effects comprise carrier depletion and carrier recovery.As the unsaturated gain coefficient is always larger than the band filling gain coefficient, The first term on the right-hand side of Eq. ( 24) is negative and therefore stands for the power gain decrease which is due to induced carrier depletion (CD).The first term also includes TPA gain.According to [28], the TPA contribution to the carrier density increase has a weight factor 1 2 , therefore we have to subtract from the total gain ( ) , g z τ in Eq. ( 18).The second term is responsible for an increase of the gain and therefore for carrier recovery (CR) on a timescale of s τ .For picosecond pulses, CD and CR operate on two vastly different timescales (see Section 2.1(i)) and both effects can be treated separately.Over the duration of the pump pulse, CD is instantaneous and dominates over a timescale much shorter than s τ , therefore CR is negligible allowing for the CR term in Eq. ( 24) to be set to zero.
After the pump pulse has vanished, there is no more carrier depletion and now CR dominates, the first term in Eq. ( 24) can be neglected.We now solve Eq. ( 24) analytically for both cases separately.When CD dominates we use the notation

Carrier depletion
For the case of weak to moderate pump energies, i.e., when the effects of CH, SHB and TPA on pump-induced changes in carrier density are neglected and the conditions ( ) ( ) As a consequence, in determining we can approximate ≈ the right-hand side of Eq. ( 24) simplifies to , .
We integrate both sides of Eq. ( 27) over the whole SOA length, substitute With the integrated unsaturated small-signal gain as an initial condition, 28) can be integrated ( [22], see also Appendix B) resulting in For the case of Gaussian input pulses, the integral in Eq. ( 29) has a closed-form expression in terms of the error-function [22].After The energy content of the part up to R τ − in the leading edge represents thus the largest energy for which the saturation can be neglected in the CD process.We can assume that the amplification ceases when the energy content of the final part in the trailing edge attains the amount considered negligible.For symmetric pulses, i.e., for .Given a certain gain reduction q at the beginning of the saturation (at time R τ τ = − ), one can solve the resulting Eq. ( 30) for R τ .This procedure is particularly useful for a Gaussian power pulse of full-width half-maximum p τ , where the ratio R p τ τ can be expressed in closed-form through the complementary error-function involving q (and other independent parameters).A practicable criterion for determining R τ will be presented at the end of Section 3.

Carrier recovery
The differential equation for the gain recovery coefficient is found by setting the first term on the right-hand side of Eq. ( 24) to zero, .
The CR gain dynamics for ∆ is obtained by integrating both sides of Eq. ( 31) over the SOA length and by solving the resulting differential equation for The value of ( ) CD R h τ ∆ on the right-hand side has to be found from Eq. ( 29).

Overall band filling
The complete gain contribution to BF is then given by ( ) Our approximations lead to an analytical result, however, the time derivative of ( )

Spectral hole burning and carrier heating
The equations describing the intraband dynamics are stated in [29].If intraband relaxation times are much smaller than the pulse width, the time derivative in these equations can be adiabatically eliminated, i.e., set to zero, and the local gain changes ( ) can be immediately written as [24] ( ) ( ) ( ) ( ) Equation ( 35) describes the ultrafast gain reduction due only to the stimulated emission that removes the "cool" carriers from the distribution.We do not consider here the carrier heating via TPA and FCA.The ultrafast processes occur simultaneously with CD, however we calculate their contributions by only considering CD, see Eq. ( 25).For strong pulses the amount of CD would also be affected by the strength of the ultrafast gain compression, but for moderate levels of gain compression it is sufficient to replace in Eq. ( 35 The ultrafast compressions are obtained by integrating both sides of Eq. ( 36) over z from 0 z = to z L = , and by taking into account the boundary conditions 0, and , exp with the result The above results remain valid in the presence of a weak probe light.By contrast, as will be shown below, the probe light can sensibly influence the TPA-process, even when its power P × is small.

Two-photon absorption
We describe the local gain reduction due to a degenerate TPA-process (simultaneous absorption of two identical photons) by an absorption coefficient ( ) which is proportional to the power of the propagating light S P and to the inverse of an effective TPA cross-section TPA A [23,28], with a proportionality factor ( ) β ω called the TPA coefficient at the pump frequency S ω : ( ) ( ) ( ) .
The area TPA A is larger than the cross-section A of the active region [23], because it always includes some adjacent layers within the cladding region.
In a non-degenerate TPA-process a VB electron can reach the CB by absorbing one photon from the pump and another photon from the probe.The pump having power S P and angular frequency S ω is assumed to be dominantly linear polarized in the direction of the unit vector S e  , which lies in the cross-section area of the SOA waveguide.The probe light with power P × and angular frequency ω × is assumed to be linearly polarized in the same plane, but along a different direction e ×  , and its contribution has to be weighted by a factor In analogy, we define the local TPA absorption coefficient experienced by the probe light in the presence of the pump light as The functions β ω and simplify Eqs. ( 40) and ( 41) to The ratio TPA r , hereafter called the TPA non-degeneracy factor, is positive and depends on the relative polarization angle.For semiconductors with zinc blende crystal structure TPA 2 r = for parallel polarizations and takes a value 2-3 times smaller for orthogonal polarizations, depending on the probe frequency [31,32].Moreover, for a weak probe light only the first term on the right-hand side of Eqs. ( 43) and ( 44) contributes, and we are left with Under the assumptions stated above, the local absorption coefficient from the pump in the presence of the (weak) probe light coincides with the degenerate case, Eq. ( 39), while the local absorption coefficient from the (weak) probe light differs from the latter by the TPA non-degeneracy factor TPA r due to the pump-probe coupling.Proceeding in analogy to the fourth formula in Eq. ( 18) and substituting the subscript x by TPA,S and TPA,× , the integrated absorption coefficients Similarly, proceeding in analogy to Eq. ( 21), we find the TPA-induced phase change contributions .For the calculation of this contribution according to Eqs. ( 18) and ( 39), it is necessary to integrate ( ) , S P z τ over the SOA length.Since in Eq. ( 25) we assumed that ultrafast processes including TPA are negligible in comparison with CD, we replaced and by using the result, see Eq. ( 98), one finds the following TPA contribution ( ) We prove in Appendix B that all coefficients ( ),

FCATPA induced refractive index
In order to include the FCA TPA process (as described in [16]), we assume that many of the photoexcited electrons reach the X-valley of the lowest CB by absorbing optical phonons of suitable wave vectors.Owing to the low amount of energy carried by the phonon [17], the energy of the electron remains nearly constant, while its momentum increases accordingly.Also the transit time to the X-valley expressed by the intervalley scattering time 0 τ is of order of the electron-phonon scattering time.Once arrived in the new X-valley state, the electrons can be further excited into the next higher CB by absorption of a photon.The contribution of this free-carrier absorption to the index change can be estimated on the basis of the classical oscillator model of Drude and Lorentz.By assuming the electrons moving freely in the optical field according to Newton's second law, one finds the decrease of the refractive index as given by the Drude formula [33] (elementary charge e , vacuum permittivity 0 ε ) An instantaneous FCA TPA process in the Γ -valley where photoexcited electrons are produced needs phonon assistance and is unlikely.
In the following we consider photoexcited electrons with average lifetime ex τ arising in the active region of cross-section A by TPA from the pump (power S P , angular frequency S ω ) and probe light (power P × , angular frequency ω × ).We describe the dynamics of these carriers by the following rate equation (  is the reduced Planck constant): The first term on the right-hand side is positive and apart from the factor 1 2 whose origin has been already explained (see Eq. ( 24)), represents the local density of photons absorbed per second from both, pump and probe light.The second term is negative and accounts for the loss of photoexcited electrons.By setting 0 The dependency of FCA n ∆ on the power squared, qualifies the FCA TPA process as a fifth-order nonlinear optical effect.According to the assumption in Section 2.2 that TPA-induced freecarrier absorption has a negligible gain coefficient FCA g ∆ , the nonlinear susceptibility ( ) 5 χ [30] is real and is proportional to the product FCA TPA η ε .Although both factors in this product are measured in 1 W , FCA η cannot represent a compression coefficient and will be simply termed as nonlinear (optical) coefficient.The FCA TPA phase contribution, Eq. ( 21), can be found by integration over the SOA length.As before we use the approximation from Eq. ( 25) to replace The function ( ) C τ is related to the function ( ) C τ given in Eq. ( 49) by the expression ( ) As the FCA TPA effect is significant, most of the electrons produced by TPA in the conduction band will have an average lifetime close to 0 τ , i.e., one also expects ex 0 τ τ  .On the other hand, the validity of the steady-state solution, Eq. ( 55), requires ex τ to be shorter than the pulse width.The large value of FCA η (see Table 1) leads to a strong decrease of the refractive index.Notice that the corresponding phase shift contribution starts changing later, with a time delay 0 τ after the input pulse peak.

Summary
By definition the power amplification of a light signal is the ratio of the output to the input signal power.For a graphical representation of  33), ( 38) and ( 50), into the third formula in Eq. ( 18), one obtains In the pump-probe experiments the overall SOA power gain According to the discussion in Section 2.1(ii), we will assume that SHB does not contribute to the phase shifts ( ) out S ϕ τ , and set the corresponding α -factor SHB 0 α = .By substituting Eqs.
(33), ( 38), ( 50) and (57) in Eq. ( 22) we find The phase shift of the probe light is related to the change in the refractive index, ( ) , see Eq. ( 16), and is schematically depicted in Fig. 1(c).The quantity obtained in a similar way as ϕ τ by taking into account Eq. ( 47), As an example, we calculate the gain dynamics of an SOA, driven with a Gaussian input pump pulse having 3.8 mW peak power and a full width at half maximum (FWHM) of 2.8 ps , i.e., with a pump input energy of 11.3 fJ .The parameters given in Table 1 are typical for an SOA being operated at a wavelength of 1.55 µm corresponding to a frequency of about reported in [22].The α -factors are taken from [13].
We postpone the discussion of the other parameters to the end of Section 3 and first calculate the optical spectra, which we then compare to measured data.
TPA nondegeneracy factor The results of the calculated gain and phase shift dynamics are plotted in Fig. 2. They are in agreement with the measured gain and phase dynamics from pump-probe experiments [13,34,35].Starting with the basic BF contribution, we successively add the contributions of the other processes, and are able to see how each process contributes individually to the overall gain.An important feature of the plots is the delay between the gain minimum and the refractive index maximum.This was shown experimentally and theoretically by others [13,34,35] and is also confirmed by our model.
Since the FCA TPA process starts 0 1 ps τ = later than the TPA itself, the refractive index is stationary during this delay, see Fig. 2(c).Notice that due to the fast and ultrafast saturation effects the refractive index minimum is reached somewhat earlier than 1 ps.Neglecting FCA g ∆ can only be compatible with moderately high TPA effects in bulk SOAs.In such cases there are fewer free-carriers arising during the FCA TPA process than removed by carrier depletion.However, the absorption cross-sections calculated in [16] are in the order of the BF differential gain.This fact can satisfactorily explain the dramatic FCA TPA phase-shift shown in Fig. 2(c) by simultaneously keeping the negligible free-carrier absorption stipulated in Fig. 2(b).The topic is thoroughly discussed in Appendix A. The advantage of using this approach is the ability to analytically model the temporal gain dynamics with the few parameters of Table 1.Most of these parameters can be derived by fitting the calculated spectrum to the measured one.The model does not require details on the geometry and on the other SOA parameters which typically are inaccessible.

Pump and probe light spectra at SOA output: comparison of model and measurements with experiment
The closed form expressions for the gain and phase dynamics, in which only a few constant parameters are chosen independently, allows us to derive the spectral shapes and to compare them with the experiment.Moreover, one can determine the remaining constant parameters by fitting the spectral traces and hence get a better understanding of the dynamic behavior of an SOA at the respective operation point.The model parameters discussed at the end of this section are found to be within a 10% ± spread of their absolute extracted values.In the following we assume that the condition in 0 P P × > is met, and in 0 P denotes the peak power of the pump pulse ( ) in P τ , i.e., ( )  and P × is the CW probe light power.We also assume -as before -that the pump and CW frequencies are not too far apart.According to Eqs. ( 61) and (63), one needs an additional parameter, the TPA non-degeneracy factor TPA r , for describing the amplification and the phase shift of the probe light.Only if TPA 1 r ≈ holds, the propagation of a weak probe light does not influence the active medium.Since the SOA material is assumed to have a negligible dispersion, we can use the same values of the constant parameters other than TPA r , for both signals.For a comparison with the experiment, we calculate the one-sided power spectra , S × Θ for the real optical pump (subscript S) and probe (subscript × ) fields at the SOA output.Starting with the output normalized RMS electric fields as given by Eq. ( 4), we perform the Fourier transform of the real part . We then shift it from the carrier angular frequencies , S ω × to the baseband, observe a factor of 1 2 for the one-sided transform, take the modulus square and express the resulting power spectrum versus the frequency shift The powers and phase shifts of the output signals are obtained in terms of ( ) , see Eqs.
the retarded time dependencies from Eqs. ( 9) and ( 10) allow to express the optical spectra as ( The calculated and measured spectra of the pulse and probe light at the SOA output are shown in Fig. 3(a) and Fig. 3(b).Mode-locked pulses with a 2.8 ps FWHM are generated at the repetition rate of 10 GHz in order to allow the complete CR between two successive pulses.Pump and probe light are injected in co-propagating direction into the SOA, utilizing a 3 dB coupler.The spectra were recorded by an optical spectrum analyzer with a resolution bandwidth of 0.1 nm , and correspond to the envelope of the oscillatory structure related to the pulse repetition rate.
Pulse spectra help investigating the dynamics of the gain and refractive index.Due to the short pulse duration of a few picoseconds, only nonlinear effects that occur while the pump is present show up in the output spectrum of the pump.Whereas, the probe senses both the ultrafast effects that occur when the pump is present and also the slower effects like CR that do not influence the pump.Hence, the temporal limitation of the injected optical power acts as a gating mechanism which allows to get insight of the dynamics at the beginning, while suppressing later effects.The full dynamics of the SOA gain and refractive index can be concluded from the probe light spectra.Although possibly significant effects such as spectral artifacts [21] and amplified spontaneous emission (ASE) were not included in the analytical model, the calculated spectra (solid black lines) resemble the measured ones (light grey lines).For instance, one observes the pronounced side lobes in the blue-shifted spectral domain of calculated and measured pump and probe spectra.The origin of this spectral region can be traced back to the FCATPA nonlinear effect -as described by the parameters FCA η and 0 τ , see Fig. 3(c) and Fig. 3(d).
Both measured and simulated probe light spectra have broad and pronounced red-shifted spectral components extending beyond 200 GHz − with respect to the center frequency, see Fig. 3(b).They mainly originate from CD and CH effects.Towards higher frequencies with respect to the center frequency one can notice two blue-shifted spectral regions.First, there is a blue-shifted spectral range starting closely to the carrier frequency.This can be attributed to the slow carrier recovery process.We attribute the small blue-shifted shoulder with a peak around 100 GHz to the presence of both TPA and FCA TPA nonlinear effects, Fig. 3(d).Other researchers have found similar spectral characteristics, particularly when experiments have been performed with short pulses of 1 ps [36] and 2 ps [37].The discussion here has shown that the presence or absence of spectral components gives an indication for the occurrence and strength of a particular nonlinear effect.Specifically, one expects a distinct blue-shifted side-lobe far off from the center frequency in the presence of TPA and FCA TPA , and a red-shifted spectral component due to CD and CH as well as a blueshifted spectral component nearby the center frequency due to CR.
In Fig. 4 we present pump and probe spectra for different pulse input powers which allow to verify the differences in order of nonlinearity.Note the characteristic behavior of the blue sideband intensity associated with FCA TPA .For in 7 dBm P = , dip due to SHB can be seen at the input peak frequency.The method outlined here allows to easily recognize when TPA and FCA TPA play a role by just looking at the spectrum and by observing that the ultrafast decrease of the refractive index occurs during the trailing edge of the pump pulse.In addition, the method allows for an estimation of the nonlinear coefficients in a particular SOA at a certain operating point (see, e.g., Table 1).For instance, in our experiment we have employed an SOA with the structural parameters from [13] .Following [21], we take the TPA effective cross-section to be of the same size as the modal cross-section, TPA A A = Γ .By using Eq.(48), we estimate the nonlinear TPA coefficient 2 β from TPA ε which is given in Table 1, and find a value of 37 cm GW which is close to the value of 35 cm GW reported in [21] for InGaAsP-SOAs.We also express the ultrafast compression coefficients in 3  m facilitating the comparison with the values given in literature: We could not see any spectral change for 0.03 dB q < or, equivalently, for R 9.5% X < , i.e., for transition times from CD to CR larger then 2.9 ps .The value 0.017 dB q = reported in Table 1 corresponds to R 3.1 ps τ = and a ratio R 6% X = .In order to assign values to the remaining fit-parameters, a sensitivity analysis has been performed.For this purpose, we consider the relative error , and denote its averaged value over the whole measured spectrum with , S F × .We determine by calculation the changes of a single parameter leading, or corresponding to a minimum of , S F × .We then look for physical or model dependent reasons impeding the realization of this minimum.As a plausibility check, the parameter q is a useful indicator for the onset of carrier saturation.It also stands for a time τ R , which has to be longer than the pulse duration (FWHM of the pump pulse).Approximate values for 0 τ and TPA r result from their physical meanings and can be used as a first guess in the optimization algorithm.The value in Table 1 can be subsequently determined from the agreement of the calculated and measured spectra.The relative error averaged over the whole measured probe light spectrum F × has a minimum for 0.017 dB q = . By keeping this value for q constant, a sensitivity analysis has been performed for the parameters FCA η , 0 τ and TPA r .The following ranges have been established in this way: -1 14.5 15.5 W  for FCA η , 0.9 1.1 ps  for the X Γ → intervalley scattering time 0 τ and 0.9 1.1


for the TPA non-degeneracy factor TPA r .Recalling that TPA 1 r = describes the TPA-process for the pump alone, we conclude that in our measurements, the probe light has no measurable effect on the active medium.

Selecting spectral regions to perform optical filter based wavelength conversion
The previous section has shown that nonlinear effects in an SOA lead to unique spectral traces.With our model, we can estimate their qualitative and quantitative relevance.This information is useful to judge if a particular effect is present at all, and if it is sufficiently large to be exploited in a nonlinear experiment.
There are three distinct spectral ranges which feature strong nonlinear effects.We now discuss the SOA with respect to its capability to perform all-optical signal processing by exploiting the nonlinear cross-gain modulation (XGM) and cross-phase modulation (XPM) effect.
The basic features within the phase dynamics and the spectral components as seen by a probe light are schematically depicted in Fig. 5.Both figures are derived from the plots in Fig. 2(c) and Fig. 3(b).The colors indicate three temporal regions associated with three distinct nonlinear effects as discussed in Fig. 1(c).In the previous section we have shown that the dominant nonlinear effects in the three different temporal regions yield clear spectral traces in different spectral regions.So, for instance, the blue-shifted (BS) shoulder away from the central frequency (shown in dark blue) is due to the ultrafast sharp decrease of the refractive index due to TPA and FCATPA in the first time segment.Next, the strong redshifted (RS) spectral component (shown in red color) is mainly induced by the carrier depletion in the SOA.Finally, the spectral component close to the CW (in light blue color) results from the carrier recovery in the active medium.Note that a fast CR process causes a broad component.
Thus, SOA nonlinear effects create new spectral components, first in the far blue, then in the red and again in the near blue spectral domain [38].By combining an SOA with properly adjusted optical filters one can build simple and efficient all-optical wavelength converters.
A wavelength converter generating non-inverted signals is schematically depicted in Fig. 6.Using suitable red-shifted optical filters (RSOF) after an SOA, a wavelength converter can be realized.In [39], a 10 Gbit/s and a 40 Gbit/s all-optical regeneration experiment has been performed with a red-shifted optical filter offset by about 50 GHz from the center frequency.In another demonstration, the blue spectral range close to the CW carrier has been exploited for performing all-optical wavelength conversion, using blue-shifted (BSOF) optical filtering.In [40], a BSOF experiment where the filter is offset by approximately 60 GHz has been demonstrated at 40 Gbit/s .
Finally, the ultra-fast nonlinear TPA and FCA TPA effects can be used.While the power in the spectral component is lowest, this spectral component is created by the fastest nonlinear effect.Since this effect does not depend on the carrier density, there is no contribution from the carrier dynamics that suffer from patterning effects.The spectral contribution due to FCA TPA is very distinctive in the output probe spectrum and its presence in the spectrum is clearly recognizable.One would immediately know of the presence of the ultrafast FCA TPA effects in their SOA and be able to optimize the experimental conditions in order to maximize this FCA TPA contribution for their optical signal processing needs.Indeed, successful patterning free all-optical wavelength conversion at 320 Gbit/s has been demonstrated [41].In that experiment, the filter was offset by more than 150 GHz from the center frequency, where we previously identified the dominant nonlinear TPA and FCA TPA spectral component.By combining both filter types, RSOF and BSOF, a pulse reformatting optical filter (PROF) has been realized and tested with signals up to 40 Gbit / s [42,43].
The PROF works as follows: The signal after the SOA is split into two copies.One signal copy is routed through an RSOF, delayed and attenuated with respect to the second copy which passes a BSOF.The introduced delay compensates for the group delay difference between RS and BS components.Both modified signal copies, which contain only parts of the BS and RS components, are superimposed after being equalized in power and group delay such that a return-to-zero (RZ) signal at the novel wavelength will result.All-optical wavelength conversion supported by RSOF, BSOF or PROF leads to a signal consisting of Fourier transform limited pulses.By properly offsetting the filter, each of the three spectral regions in Fig. 5(b) can be used to induce a progressive phase shift between subsequent pulses.This way RZ, vestigialsideband (VSB) and carrier-suppressed return-to-zero (CSRZ) signals can be generated [44][45][46].
In this paper we have concentrated on the characterization and interpretation of the nonlinear responses of the active medium to Gaussian pulses.

Evidence of phase patterning effects in the spectrum
We now consider all-optical wavelength converted data signals rather than wavelength converted clock pulse-trains.The spectrum now changes, see Fig. 7.The spectrum of a 33% RZ signal at 40 Gbit/s, see Fig. 7(a), shows a central tone at the carrier frequency and two sidebands that contain the information data with tones located at multiples of the data rate.The lateral tones are accompanied by inverted spikes that are blue shifted with respect to each tone.These characteristics remain even after passing the signal through a filter for reshaping and pattern effect mitigation, e.g., a PROF as in Fig. 7  At high modulation frequencies, e.g., 40 GHz, the CR between consecutive signal pulses becomes incomplete leading thus to variations in the pulse amplitude and phase, i.e., to a patterning effect.For simplicity, we consider below only the incomplete phase restoration from bit-to-bit.To demonstrate that this phase patterning leads to inverted spikes and contiguous Sshaped shoulders we propose the following model: Every mark bit entering the SOA induces a decrease of the carrier density and thus an increase of the refractive index.As a consequence, the phase of each mark bit changes by a certain value 1 φ , while a space bit allows for a weak phase recovery 2 1 φ φ < .So each pulse entering the SOA experiences a phase accumulated from prior pulses.For simplicity, we assume constant phase shifts 1 φ and 2 φ .
They might indeed be realized for weak signals where patterning effect might add up linearly [47].The analysis can be extended to format and wavelength converted signals.In this case a progressive phase shift PPS is added to the total phase [46].To see the implications of phase patterning on the optical spectrum, we calculate the spectrum of an RZ and a CSRZ signal consisting of Gaussian shaped pulses with data encoded as pseudo-random bit sequence (PRBS) of 11  2 1 − .To each bit of the signal we then assign an accumulated phase.Table 2 shows the signal phase assignment for an exemplary bit sequence 110100.The time evolution of the amplitude and phase shift for the bit sequence is shown in Fig. 8 for (a) RZ and (b) a signal where the filter is offset such that a CSRZ signal is generated [46].Figure 9 shows the calculated spectra of the RZ and the CSRZ signals for The position of the tones does not change with 1 φ , 2 φ but the spectral region in-between tones shows the characteristic structure similar to the one observed in the experiment.The presence of phase patterning effect thus leaves traces in the signal spectrum.Higher order pseudo-random bit sequences (e.g. 31  2 1 − ) will show a slightly different output spectrum due to saturation and recovery effects.Nevertheless, since longer mark-andspace sequences contribute only to the spectral components at rather low frequencies, it will be difficult to observe the difference in the output spectrum.

Conclusion
We present an analytical model for the optical response of an SOA to short input pulses.By suitably adjusting the few parameters of the model we have successfully identified the nonlinear effects in the SOA from spectral traces.In particular, by comparing the calculated and measured output pump spectra, we clarified the dynamics of the refractive index during the pulse-SOA interaction.The gating mechanism of the pulse allows for the experiential to give information to conclude that the significant initial decrease of the refractive index is mainly due to an FCATPA process, which occurs as a result of an intervalley scattering.By fitting the calculated spectra for both pump pulse and CW, we determined the magnitude of the time delay and identified it with the X Γ → intervalley scattering time in the SOA.More generally, we know the full dynamics of the SOA gain and refractive index.On this basis, we relate the different spectral components of the modulated probe light to the respective physical processes in SOAs.We analyze the spectrum of a converted signal and estimate the importance of three dominant effects.
When converting signals with data rates of 40 Gbit/s, the carrier density cannot recover completely within a single bit slot.The incomplete carrier recovery between subsequent marks induces a progressive increase of the refractive index inside the SOA, and a corresponding decrease of the probe light phase, thereby leading to characteristic features in the optical spectra of a signal.This way, the spectral methods indicate the presence of patterning effects in SOA.The characteristic features observed in the spectra of the PRBS signals at high data powers are explained by assuming the constant progressive increase of the refractive index after each mark bit, and a smaller constant progressive decrease for each space bit.The regions between consecutive tones in the spectra of PRBS signals carry the information on the incomplete carrier recovery in the SOA allowing thus to see the patterning effect in the signal spectrum.

Appendix A.
The generation of TPA-induced free-carrier can occur either instantaneously (instantaneous FCATPA) by the absorption of a photon and a phonon, or via a two-stage process (delayed FCA TPA ), which we now describe.In the first stage, electrons are excited to the Γ -valley of the lowest conduction band (CB 1 ) at energies above the X-valley minima by the absorption of two photons.These photoexcited electrons are then scattered to the X-valley by short wavelength lattice phonons.The transfer to the X-valley (intervalley scattering X Γ → ) is favored by the large density of states in the X-valley as compared to the Γ -valley [16].
Both steps of this two-stage process do not occur simultaneously as there is a short delay of a few picoseconds between these two processes.Due to the specific band structure of the semiconductor material, once in the X-valley, those electrons can absorb a photon, without a phonon, to reach the next higher conduction band (CB 2 ).
Previous work from our group may have encountered an FCA TPA process delayed by the intervalley scattering X Γ → , because the carrier density rate equation in [13], Eq. ( 13), needed a correction term for carrier leakage at high carrier energies.Let X W be the energy separation between the X-valley of CB 1 and the top of the valence band for the quaternary compound InGaAsP (see Fig. assume single-mode propagation and represent the waveguide mode by a plane wave with an effective propagation constant 0 0 k n and an effective refractive index 0 n .The fields depend on the coordinates z (position) and t (time), and are represented by the square root of the respective powers.Denoting the powers by

Fig. 1 .
Fig. 1.Schematic input and output pump power, probe output power, and refractive index change as a function of retarded time g t z v τ = − be the same.For simplicity, we neglect the coefficient FCA g ∆ associated with the two-photon induced free-carrier absorption in our analysis and the corresponding local change index is derived by applying the Drude model, see Eq. (51).The total change in the refractive index τ by the solution ( ) CD , P z τ of the propagation equation

(
attained its minimum, it recovers back to the unsaturated value 0 h .As the model is based on a sharp separation of CD and CR, the temporal threshold R τ τ = needs to be introduced.A suitable time R τ − can be identified as the moment when the gain saturation begins.The gain reduction (in dB) up to this moment can be written as as defined in Eq.(26 pump in the presence of the probe light is related to each other, and become equal for S ω ω × =[31].In general, a TPA in pump-probe experiments consists of three kinds of absorption processes occurring at the same time: absorption processes of identical photon pairs from the pump and probe light, and of different photon pairs from the pump-probe coupling.For neighboring angular frequencies 0 S in Eqs.(46) and (47), associated to the TPA-process are all expressed in terms of the contribution weight and represent the three kinds of aforementioned TPA-processes, two degenerate and one non-degenerate TPA processes.For a sufficiently weak probe light, one can neglect the terms containing .(52), (55) and (39) into the right-hand side of the Drude formula (Eq.(51)), one gets τ , perform the integral with the methods developed in Appendix B, and obtain

(
. 1(a).Inserting Eqs. ( 60) is always given by the probe light amplification.It is represented in Fig. 1(b) for a CW probe light.The logarithmic integrated gain in agreement with the typical range of 5 10 pJ 

Fig. 2 .
Fig. 2. Calculated time dependence of SOA gain G = exp(h × ).The origin τ = 0 of the retarded time coincides with the maximum of the pump pulse P s (0,τ), see Fig. 1.With our model we calculate ln (G) = h × and the phase shift 0 n k L . (62) and (63).Writing

Fig. 3 .
Fig.3.Calculated spectra (log scale) using the parameters from Table1and comparison with measured ones.(a) Measured (light gray lines) and calculated (solid black lines) spectra after the SOA of the pump pulse and (b) measured (light gray lines) and calculated (solid black lines) spectra after the SOA of the probe light.The decomposition of the spectra is shown in (c) for the pump pulse and in (d) for the CW pulse light.(BF: Band filling; TPA: Two photon absorption; FCA TPA : Free-carrier absorption induced by TPA; SHB: Spectral hole burning; CH: Carrier heating).

FrequencyFig. 4 .
Fig. 4. Optical spectrum for P in = 0, 3, 7 dBm; P cw = 2 dBm; I = 400 mA of the (a) pump light and (b) probe light behind the SOA.The frequency origin corresponds to the input pump peak and to the probe light frequency, respectively.
now look for the parameters q , FCA η , 0 τ and TPA r , which do not change the calculated optical spectra.A reduction of the gain by 0.01 dB q = compared to the unsaturated gain marks the beginning of saturation and corresponds to a time R 3.3 ps τ = for a transition from CD to CR and to a ratio

Fig. 5 .
Fig. 5. Probe light phase shift and optical spectrum.(a) Phase shift dynamics and (b) corresponding induced spectral components in a measured probe light spectrum

Fig. 6 .
Fig.6.Scheme of an optical filter based all-optical wavelength conversion system.

Fig. 7 .
Fig. 7. Measured spectra of a 40 Gb/s wavelength converted signal.Optical spectrum (a) directly after the SOA and (b) after the PROF.

Table 2 .Fig. 8 .
Fig. 8. Schematic description of a signal evolution.Amplitude (green line) and phase (red line) evolution of (a) an RZ signal PPS 0 ( ) = and (b) a CSRZ signal PPS ( ) π = − after a filter supported wavelength conversion.A mark bit decreases the phase by 1 φ and a space bit

Fig. 9 .
Fig. 9. Spectral traces of a phase-patterning effect.Calculated spectra for standard (a) RZ and (b) CSRZ signals and for phase patterning affect (c) RZ and (d) CSRZ signals after wavelength conversion, respectively wavelength and format conversion.

10
), and X,cl W be the energy separation between the aforementioned X-valley and the top of the valence band for the cladding.Then X electrons can drift into the cladding, accounting thus for the carrier leakage. .
(18)e the first Eq.(18).By making use of this relation in Eq. (34), one obtains ).As a consequence, verify that ex τ is the photoexcited carrier lifetime.The double heterostructure confines all carriers in the active region, i.e., the carriers are distributed over the cross-section A .In the stationary regime