Modulational instability of discrete solitons in coupled waveguides with group velocity dispersion

We study temporal modulational instability of spatial disc rete solitons in waveguide arrays with group velocity dispersio n (GVD). For normal GVD we report existence of the strong ’neck’-type ins tability specific for the discrete solitons. For anomalous GVD the instab ili y leads to formation of the mixed discrete-continuous spatio-tempor al quasi-solitons. Feasibility of experimental observation of these effects i n the arrays of silicon-on-insulator waveguides is discussed. © 2006 Optical Society of America OCIS codes: (060.5530) Pulse propagation and solitons References and links 1. D.N. Christodoulides, F. Lederer, and Y. Silberberg, ”Di scretizing light behaviour in linear and nonlinear waveguide lattices,” Nature424, 817-823 (2003). 2. F. Wise and P. Di Trapani, ”The hunt for light bullets Spat io-temporal solitons,” Opt. Phot. News (February, 2002), pages 29-32. 3. D. Cheskis, S. Bar-Ad, R. Morandotti, J.S. Aitchison, H.S . Eisenberg, Y. Silberberg, and D. Ross, ”Strong Spatiotemporal Localization in a Silica Nonlinear Waveguide Ar ray,” Phys. Rev. Lett. 91, 223901 (2003). 4. P. St.J. Russell, ”Photonic crystal fibers,” Science 299, 358-362 (2003). 5. A. B. Aceves, C. De Angelis, A. M. Rubenchik, and S. K. Turit syn, ”Multidimensional solitons in fiber arrays,” Opt. Lett.19, 329-331 (1994) 6. A.B. Aceves, G. Fibich, and B. Ilan, ”Gap-soliton bullets in waveguide gratings,” Physica D 189, 277-286 (2004). 7. B.B. Baizakov, B.A. Malomed, and M. Salerno, ”Multidimensi onal solitons in a low-dimensional periodic potential,” Phys. Rev. A70 053613 (2004). 8. S. Droulias, K. Hizanidis, J. Meier, and D. Christodoulid es, ”X – waves in nonlinear normally dispersive waveguide arrays,” Opt. Express 13, 1827-1832 (2005). 9. A.A. Sukhorukov and Y.S. Kivshar, ”Slow-light optical bu llets in arrays of nonlinear Bragg-grating waveguides,” http://xxx.lanl.gov/abs/physics/0605194. 10. V.E. Zakharov and A.M. Rubenchik, ”Instability of waveg uides and solitons in nonlinear media,” Zh. Eksp. Teor. Fiz. 65, 997-1011 (1973) [Sov. Phys. JETP 38, 494 (1974)]. 11. Y.S. Kivshar, D.E. Pelinovsky, ”Self-focusing and tran sverse instabilities of solitary waves,” Phys. Rep. 331, 118-195 (2000). 12. A.V. Buryak, P. Di Trapani, D.V. Skryabin, S. Trillo, ”Op tical solitons due to quadratic nonlinearities: from basic physics to futuristic applications,” Phys. Rep. 370, 63-235 (2002). #77646 $15.00 USD Received 29 November 2006; accepted 30 November 2006 (C) 2006 OSA 11 December 2006 / Vol. 14, No. 25 / OPTICS EXPRESS 12347 13. S. Minardi, J. Yu, G. Blasi, A. Varanavicius, G. Valiulis , A. Berzanskis, A. Piskarskas, P. Di Trapani ”Red solitons: Evidence of spatiotemporal instability in χ(2) spatial soliton dynamics,” Phys. Rev. Lett. 91, 123901 (2003). 14. S.D. Jenkins, D. Salerno, S. Minardi, G. Tamosauskas, T.A .B. Kennedy, P. Di Trapani, ”Quantum-noise-initiated symmetry breaking of spatial solitons,” Phys. Rev. Lett. 95, 203902 (2005). 15. S.P. Gorza, P. Emplit, M. Haelterman, ”Observation of the s nake instability of a spatially extended temporal bright soliton,” Opt. Lett.31, 1280-1282 (2006). 16. P. Dumon, G. Priem, L.R. Nunes, W. Bogaerts, D. van Thourhou t, P. Bienstman, T.K. Liang, M. Tsuchiya, P. Jaenen, S. Beckx, J. Wouters and R. Baets, ”Linear and nonlin ear anophotonic devices based on silicon-oninsulator wire waveguides,” Jap. J. of Appl. Phys. 45, 6589-6602 (2006). 17. F.N. Xia, L. Sekaric, and Y.A. Vlasov, ”Mode conversion l osses in silicon-on-insulator photonic wire based racetrack resonators,” Opt. Express 14, 3872-3886 (2006). 18. M.A. Foster, A.C. Turner, J.E. Sharping, B.S. Schmidt, M. Lipson, and A.L. Gaeta, ”Broad-band optical parametric gain on a silicon photonic chip,” Nature 441, 960-963 (2006). 19. J. Meier, G.I. Stegeman, D.N. Christodoulides, Y. Silber berg, R. Morandotti, H. Yang, G. Salamo, M. Sorel, J.S. Aitchison, ”Experimental observation of discrete modulatio nal instability,” Phys. Rev. Lett. 92 163902 (2004). 20. D.N. Christodoulides and R.I. Joseph, ”Discrete self-f ocusing in nonlinear arrays of coupled waveguides,” Opt. Lett. 16, 446-448 (1991). 21. S. Darmanyan, I. Relke, and F. Lederer, ”Instability of co ntinuous waves and rotating solitons in waveguide arrays,” Phys. Rev. E55, 7662-7668 (1997). 22. M. Stepic, L. Hadzieski, M. Skoric, ”Stability of one-di mensional array solitons,” Phys. Rev. E 65, 026604 (2002). 23. Y.V. Kartashov and L. Torner, ”Surface soliton arrays,” Eur. Opt. Soc. Topical Meeting on Nonlinear Optics, TOM 6 (Paris, 2006). 24. E.W. Laedke, K.H. Spatschek, and S.K. Turitsyn, ”Stabil ity of Discrete Solitons and Quasicollapse to Intrinsicall y Localized Modes,” Phys. Rev. Lett. 73, 1055–1059 (1994). 25. N. Akhmediev and J.M. Soto-Crespo, ”Generation of a train of three-dimensional optical solitons in a selffocusing medium,” Phys. Rev. A 47, 1358-1364 (1993). 26. B. Deconinck, D.E. Pelinovsky, and J.D. Carter, ”Transv er e instabilities of deep-water solitary waves,” Proc. of the Royal Soc. A462, 2039-2061 (2006).


Introduction
Interplay between nonlinear effects and diffraction control in planar arrays of optical waveguides and in the light induced lattices in photorefractive materials have recently attracted substantial interest, see, e.g., [1] for a review.Spatiotemporal nonlinear effects relying on group velocity dispersion and, in particular, the topic of light bullets, is another active research area [2].Despite numerous interesting effects, which interplay of controllable diffraction and GVD can bring into nonlinear optics, the experimental efforts devoted to these problems in recent years have been limited [3].This is most likely because limitations imposed by the material and structure parameters.Photonic crystal fibers could be considered as a notable exception [4], which allows large degree of the GVD engineering by periodic variations of the refractive index, but on the other hand it most often enforces geometrical suppression of the diffraction, rather than allowing exploration of the functionality of the diffraction control as it has been done in planar structures [1].On the theoretical side the results on the spatio-temporal and soliton effects in nonlinear waveguide arrays and related band-gap systems are not that numerous either.They include the existence of the continuous in time and discrete in space bullets [5], multi-dimensional solitons in systems with a gap-like dispersion along one coordinate and the free space diffraction along the other dimensions [6,7], spatio-temporal X-waves in waveguide arrays [8] and solitons in coupled waveguides with Bragg gratings [9].
One of the theoretical methods of dealing with spatio-temporal problems, which has been well developed for the continuous systems, is to take a spatial soliton and consider its robustness with respect to the time dependent perturbations [10,11,12].It is known that presence of the GVD usually destabilizes spatial solitons in the continuous systems via growth of the frequency side-bands.This process is usually called modulational instability (MI).Spatio-temporal dynamics induced by solitonic MI in continuous systems has not only been studied theoretically [10,11,12], but more recently observed experimentally [13,14,15].Similar instabilities can be expected to occur for the discrete solitons in waveguide arrays.Our primary goal here is to consider temporal MI of the single peak spatial solitons in waveguides arrays with GVD.This problem becomes particularly relevant now because of the advances in fabrication of low loss planar silicon-on-insulator (SOI) structures with strong and controllable GVD and short coupling length, see, e.g., [16,17,18].Most importantly the strong spatial confinement on the nano-scale and large Kerr nonlinearity bring the power levels required for the MI and soliton related effects to show on the millimeter to centimeter lengths down to few watts.The estimates for parameters typical for SOI waveguides show that all the effects described below are within the experimental reach, while GVD induced modulational instability in a single SOI waveguide has already been reported in [18].
We should also mention here that the MIs of spatially extended, i.e. non-localized, supermodes of the waveguide arrays induced by the discrete diffraction have been recently observed experimentally [19] and previously studied theoretically [20].Papers [21,22] studied the same case of diffraction induced MI, but for the supermodes consisting from the temporal solitons.Spatial MI of spatial surface solitons in optical lattices has been recently reported [23].However, none of the known to us studies explored the problems of GVD induced instabilities of the spatial discrete solitons.

Model
We model an array of dielectric waveguides by a set of coupled NLS equations Here n enumerates the waveguides, τ and ζ are the time and coordinate along the waveguide, respectively.γ = 2πn 2 /(Sλ) is the nonlinearity parameter, where S is the effective mode area and n 2 is the Kerr coefficient.β 2 is the GVD coefficient.κ = π/(2l c ) is the coupling parameter and l c is the coupling length.In order to put Eqs. (1) into dimensionless form we divide them by some fixed length l.Then in terms of dimensionless propagation distance z = ζ/l and dimensionless time t = τ / |β 2 |l Eqs.(1) take the form

Stability analysis
Let us consider time-independent discrete soliton solutions of Eqs.(1) [20].These solutions having the form U n = a n e iqz can be found numerically, see Figs. 1(a,d,g).
It is known that the discrete solitons are dynamically stable for β 2 = 0 and the Vakhitov-Kolokolov stability criterion , is satisfied for them [24].As it is shown below the instabilities can arise as soon as GVD is taken into account.
To study stability with respect to time-dependent perturbations, we make the following ansatz Here ω is the perturbation frequency.The linearized equations for the amplitudes of small perturbations n,± can be transformed into the operator form where where q = q + 2C.MI of discrete solitons manifests itself through a growth of perturbations in a certain range of frequencies ω.This means that in this range there exists an eigenvalue of the problem (3) such that Im(λ) > 0. According to the classical results on MI of the bright solitons in the continuous NLS model [10] the eigenvectors of the corresponding eigenvalue problem can be either symmetric or antisymmetric with respect to the reflection about the soliton center.The instabilities associated with symmetric and antisymmetric eigenvectors are usually referred to as a 'neck'-instability and 'snake'-instability, respectively [10,11,12].We will adopt the same terminology in our stability analysis of the discrete solitons.If the soliton is centered at n = n 0 , then the eigenvectors of Eq. ( 3) can be either symmetric on the replacement of n 0 + m with n 0 − m (as the soliton itself) or antisymmetric.
It follows from Eq. ( 3) that the stability of a discrete soliton with respect to the perturbations with ω = 0 is determined by the eigenvalues of the matrix L0 .Since the soliton is stable in the absence of GVD, all the eigenvalues of L0 are located on the real axis of the complex λ-plane.The zero eigenvalues of L0 , which are related to the continuous symmetries of Eqs.(1) can be used to get some analytical results.The symmetry transformation U n → U n e iθ (where θ is an arbitrary phase) together with the Hamiltonian structure of these equations imply the existence of the two zero eigenvalues of the matrix L0 with the same eigenvector x 0 .This means that the following two identities: L0 x 0 = 0 and L0 x 1 = x 0 , which can be rewritten Figure 1: MI of discrete solitons.The first row corresponds to C = 7, the second to C = 15 and the third one to C = 30, respectively.q = 10 and N = 51 for all the panels.The right column shows transverse profiles of the discrete solitons.The middle column presents the frequency dependence of the MI growth rate (Imλ > 0) in the anomalous GVD regime (s < 0).The right column shows all the unstable eigenvalues in the case of the normal GVD (s > 0).Letters 'N' and 'S' mark the 'neck' and 'snake' instabilities, respectively.in the form L2 0 x 1 = 0, are satisfied.Here x 0 = (a 1 , −a 1 , . . ., a N , −a N ) T , x 1 = −∂ q (a 1 , a 1, . . ., a N , a N ) T .When the frequency ω deviates from zero, see Eq. ( 3), these two eigenvalues move away from the origin in the complex λ-plane.They move in the opposite directions either along the real or along the imaginary axis.The latter case corresponds to the 'neck'-instability associated with the vector x 0 that has the same symmetry as the soliton itself.Note, that in the discrete case the continuous symmetry with respect to the lateral shifts of the soliton position is absent, and therefore the degeneracy of the zero eigenvalue of L0 in the discrete model is half less than in the continuous one [10,11,12].
Assuming that ω 1 we write the following asymptotic expansions for the eigenvalue and corresponding eigenvector: λ = ωλ 1 + ω 2 λ 2 + . . .and x = x 0 + iλ 1 ωx 1 + ω 2 x 2 + . ... Substituting these expansions into Eqs.(3), we get in the third order of the perturbation theory the solvability condition λ 2 1 = sQ/(2∂ q Q).According to this condition the long wavelength instability of the 'neck'-type takes place for s < 0 (anomalous GVD).This result is not surprising because essentially the same instability persists even for κ = C = 0.The middle column in Fig. 1 shows the 'neck' instability growth rates for different values of C. The low frequency part of the instability growth rates is approximately described by the above analytical expression.Direct numerical modeling of Eqs.(1) with the initial condition corresponding to a discrete soliton shows in this case formation of the regular trains of spatio-temporal quasi-solitons, which are discrete in space and continuous in time, see the left column in Fig. 2. Thus development of MI for the anomalous GVD case is qualitatively similar to the MI of the spatial solitons in the continuous NLS equation with saturable nonlinearity [25].
In the case of normal GVD (s > 0) we have found complex instability spectra consisting from multiple sidebands, see the rightmost column in Fig. 1.These type of spectra appear to be specific to the discrete solitons.For the relatively small coupling strength, i.e. sufficiently far from the continuous limit, the dominant instability is of the 'neck'-type.This instability leads to the break-up of the initial soliton to the localized lumps of light, which disperse with further propagation, see the middle column in Fig. 2. Contrary, for the continuous 2D NLS with the normal GVD the anti-symmetric 'snake'-like instability dominates dynamics of 1D bright solitons [10,26].More, recent studies [26] have demonstrated that the symmetric 'neck'-type eigenvectors also can be unstable in the continuous hyperbolic 2D NLS, with their growth rate been below the one for the 'snake' MI.As we have found in our model the 'neck' instability dominates the dynamics of the discrete solitons for small to moderate values of the coupling coefficient, see Fig. 1.Only for strong coupling, when the system becomes quasi-continuous, the 'snake' instability starts to be dominant over the 'neck' one.This leads to the break up of the discrete soliton in the snake-like fashion, see the right most column in Fig. 2. Note, that the dominant MI band of the 'neck'-type found for small frequencies ω close to zero is associated with a pair of complex eigenvalues λ.At the same time the dominant band of the 'snake' MI and the strongest peak of the 'neck'-MI at the relatively large frequencies ω have purely imaginary eigenvalues λ.

Physical estimates
As a guideline for physical estimates we consider parameters typical for SOI waveguides.In particular for a channel waveguide with width 480nm and thickness 220nm [16] GVD at 1.5μm is anomalous and its value is 580ps/nm/km.For width below 400nm or above 640nm GVD becomes normal.The coupling length for spacing around 400nm can be estimated at 200μm [17].The kerr coefficient n 2 for silicon is ∼ 6 × 10 −14 cm 2 /W.For our power estimates we take the effective area S 0.3μm 2 .Then one can show that C = 7, 15, 30 used in numerical modeling give the following values for the dimensionless unit of z: 0.88mm, 2mm, 3.8mm.One unit of the time t corresponds to 22fs, 34fs, 47fs and unit of the peak power -to 1.7W, 0.8W, 0.4W, respectively.Considering that in the best SOI waveguides the loss is few dB/cm and remembering about two-photon and free carrier absorption, our power estimates should be scaled up.In particular, MI in a single SOI waveguide reported in [18] has been observed for 10W of pump power.Note, that though more detailed account of the above absorption mechanisms is desirable in future research, it can be forecasted, that they will simply proportionally suppress the instability growth rate, without qualitative changes in the effect itself.The above time estimates show that MI can be observed already with pico-second pump pulses, when the role of the free carrier absorbtion and dispersion is negligible.

Summary
We have analyzed modulational instability of bright discrete solitons in the waveguide arrays with group velocity dispersion.In the case of normal GVD we have found multiple instability bands.For weak to moderate strength of coupling the discrete solitons exhibit the 'neck' instability leading to breakup of the solitons into a train of dispersive pulses.Only for strong coupling, i.e., in the quasi-continuous limit, this instability is getting gradually suppressed by the 'snake' instability known for the 2D continuous NLS model [10].In the case of anomalous GVD the expected neck type instability leads to formation of composite discrete-continuous spatio-temporal quasi-solitons.

Figure 2 :
Figure 2: The left column shows patterns of the 'neck' instability for anomalous GVD (s < 0) for 3 consequential values of the propagation distance z: C = 7.The middle column shows patterns of the 'neck' instability for normal GVD (s > 0): C = 7.The right column shows patterns of the 'snake' instability for normal GVD (s > 0): C = 30.q = 10 and N = 51 for all the panels.