Nonlinear diffusion and beam self-trapping in diffraction-managed waveguide arrays

We study nonlinear propagation of light in diffractionmanaged photonic lattices created by periodically-curved arrays of optical waveguides. We identify different regimes of the nonlinear propagation of light in such structures depending on the input power. We start from the regime of self-collimation at low powers and demonstrate that, as the beam power increases, nonlinearity destroys the beam self-imaging and leads to nonlinear diffusion. At higher powers, we observe a sharp transition to the self-trapping and the formation of discrete diffraction-managed solitons. © 2007 Optical Society of America OCIS codes: (190.4420) Nonlinear optics, transverse effects in; (190.5940) Self-action effects. References and links 1. D. N. Christodoulides, F. Lederer, and Y. Silberberg, “Discretizing light behaviour in linear and nonlinear waveguide lattices,” Nature 424, 817–823 (2003). 2. H. S. Eisenberg, Y. Silberberg, R. Morandotti, and J. S. Aitchison, “Diffraction management,” Phys. Rev. Lett. 85, 1863–1866 (2000). 3. T. Pertsch, T. Zentgraf, U. Peschel, A. Brauer, and F. Lederer, “Anomalous refraction and diffraction in discrete optical systems,” Phys. Rev. Lett. 88, 093901–4 (2002). 4. P. T. Rakich, M. S. Dahlem, S. Tandon, M. Ibanescu, M. Soljacic, G. S. Petrich, J. D. Joannopoulos, L. A. Kolodziejski, and E. P. Ippen, “Achieving centimetre-scale supercollimation in a large-area two-dimensional photonic crystal,” Nature Materials 5, 93–96 (2006). 5. S. Longhi, M. Marangoni, M. Lobino, R. Ramponi, and P. Laporta, “Observation of Dynamic Localization in Periodically Curved Waveguide Arrays,” Phys. Rev. Lett. 96, 243901–4 (2006). 6. S. Longhi, “Self-imaging and modulational instability in an array of periodically curved waveguides,” Opt. Lett. 30, 2137–2139 (2005). 7. R. Iyer, J. S. Aitchison, J. Wan, M. M. Dignam, C. M. de Sterke, “Exact dynamic localization in curved AlGaAs optical waveguide arrays,” Opt. Express 15, 3212–3223 (2007), http://www.opticsinfobase.org/abstract.cfm?URI=oe-15-6-3212. 8. I. L. Garanovich, A. A. Sukhorukov, and Yu. S. Kivshar, “Broadband diffraction management and selfcollimation of white light in photonic lattices,” Phys. Rev. E 74, 066609–4 (2006). 9. D. N. Christodoulides and R. I. Joseph, “Discrete self-focusing in nonlinear arrays of coupled wave-guides,” Opt. Lett. 13, 794–796 (1988). 10. H. S. Eisenberg, Y. Silberberg, R. Morandotti, A. R. Boyd, and J. S. Aitchison, “Discrete spatial optical solitons in waveguide arrays,” Phys. Rev. Lett. 81, 3383–3386 (1998). 11. Yu. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic Press, San Diego, 2003). 12. M. J. Ablowitz and Z. H. Musslimani, “Discrete diffraction managed spatial solitons,” Phys. Rev. Lett. 87, 254102–4 (2001). 13. K. Staliunas, R. Herrero, and G. J. Valcarcel, de, “Arresting soliton collapse in two-dimensional nonlinear Schrodinger systems via spatiotemporal modulation of the external potential,” Phys. Rev. A 75, 011604–4 (2007). #82632 $15.00 USD Received 1 May 2007; revised 12 Jul 2007; accepted 15 Jul 2007; published 17 Jul 2007 (C) 2007 OSA 23 July 2007 / Vol. 15, No. 15 / OPTICS EXPRESS 9547 14. S. Somekh, E. Garmire, A. Yariv, H. L. Garvin, and R. G. Hunsperger, “Channel optical waveguide directional couplers,” Appl. Phys. Lett. 22, 46–47 (1973). 15. Y. V. Kartashov, L. Torner, and D. N. Christodoulides, “Soliton dragging by dynamic optical lattices,” Opt. Lett. 30, 1378–1380 (2005). 16. I. L. Garanovich, A. A. Sukhorukov, and Yu. S. Kivshar, “Soliton control in modulated optically-induced photonic lattices,” Opt. Express 13, 5704–5710 (2005), http://www.opticsinfobase.org/abstract.cfm?URI=oe-13-15-5704. 17. C. R. Rosberg, I. L. Garanovich, A. A. Sukhorukov, D. N. Neshev, W. Krolikowski, and Yu. S. Kivshar, “Demonstration of all-optical beam steering in modulated photonic lattices,” Opt. Lett. 31, 1498–1500 (2006). 18. Y. Kominis and K. Hizanidis, “Continuous-wave-controlled steering of spatial solitons,” J.Opt. Soc. Am. B 21, 562–567 (2004). 19. I. Tsopelas, Y. Kominis, and K. Hizanidis, “Soliton dynamics and interactions in dynamically photoinduced lattices,” Phys. Rev. E 74, 0366139–8 (2006). 20. A. C. Scott and L. Macneil, “Binding energy versus nonlinearity for a “small” stationary soliton,” Phys. Lett. A 98, 87–88 (1983).


Introduction
Propagation of light in dielectric media with a periodically-varying refractive index is known to demonstrate many novel features in both linear and nonlinear regimes [1].Over the recent years, different types of periodic photonic structures, including arrays of evanescently coupled optical waveguides, optically-induced lattices in photorefractive materials, and photonic crystals, have been employed to engineer and control many fundamental properties of light propagation.In particular, the idea to control the light spreading through diffraction management [2] has attracted a special attention.It was shown that both magnitude and sign of the beam diffraction can be controlled in periodic photonic structures.For example, diffraction can be made negative allowing for focusing of diverging beams [3], or it can even be suppressed leading to the self-collimation effect when the beam width does not change over hundreds of the free-space diffraction lengths [4].Self-collimation of light beams was also realized experimentally in periodically curved waveguide arrays [5][6][7].Recently, the concept of the broadband diffraction management of polychromatic light has been introduced for the periodically-curved waveguide arrays with special bending profiles [8].
The combination of tailored diffraction characteristics and light self-action opens new possibilities for the power-controlled beam shaping and switching in nonlinear photonic structures.Various schemes for active beam control based on the special properties of narrow selflocalized beams in straight waveguide arrays, called discrete spatial solitons, have been suggested and demonstrated [9][10][11].It was shown that solitons can exist in diffraction-managed lattices [12,13] in the regime when beam exhibits effectively averaged diffraction.However, the analysis of narrow beam self-action in periodically-curved waveguide arrays, beyond the applicability of averaging procedures, remained largely unexplored.An especially intriguing problem is the nonlinear beam self-action under the condition of linear self-collimation, where diffraction is suppressed in all orders.In the latter case, the modulational instability is suppressed [6], suggesting that discrete soliton formation may demonstrate unusual features.
In this work, we investigate nonlinear propagation of light beams in the diffraction-managed periodically curved photonic lattices [such as sketched in Fig. 1(a)], and identify different regimes of the light propagation depending on the input power.In particular, when linear discrete diffraction is fully suppressed, we observe the transition from the regime of selfcollimation, for low powers, to that of the discrete self-trapping and the formation of lattice solitons, for high powers.This occurs through the intermediate regime of the nonlinear diffusion, where limited nonlinear beam broadening takes place.We also observe the similar regime of nonlinear diffusion in other types of periodically-curved photonic lattices; this regime has no analogies with the nonlinear beam self-focusing or self-defocusing in a bulk medium [11], or discrete self-trapping of light in arrays of straight waveguides [1].We also show, that the crit- ical power required for the formation of lattice solitons in curved waveguide arrays is higher than that in the straight waveguide arrays, and it can be effectively controlled by changing the amplitude of the waveguide bending.

Propagation of light beams in arrays of curved waveguides
We study the propagation of light beams in a one-dimensional array of coupled nonlinear optical waveguides with the transverse period d in the x direction, where the waveguide axes are periodically curved in the propagation direction z with the period L d [see Fig. 1(a)].When the tilt of beams and waveguides at the input facet is less than the Bragg angle, the beam propagation is primarily characterized by coupling between the fundamental modes of the individual waveguides, and it can be described by the tight-binding equations taking into account the periodic waveguide bending [5,6], idΨ n /dz+C 0 (Ψ n+1 + Ψ n−1 ) = ω ẍ0 (z)nΨ n − γ|Ψ n | 2 Ψ n , where Ψ n (z) is the amplitude of the n-th waveguide, ω = 2πn 0 d/λ is the dimensionless frequency, λ is the vacuum wavelength, n 0 is the average refractive index of the medium, γ is an effective nonlinear coefficient which accounts for the Kerr-type nonlinear response of the waveguide material, and the dots stand for the derivatives.Transverse shift x 0 (z) ≡ x 0 (z + L) defines the periodic longitudinal bending profile of the waveguide axis.Coefficient C 0 defines the coupling strength between the neighboring waveguides, and it characterizes diffraction in a straight waveguide array with x 0 ≡ 0 [14].Then, the total electric field envelope E(x, z) is represented as a superposition of the modes E 0 (x) of the individual waveguides, The distinctive features of discrete beam dynamics become most evident when only one waveguide is excited at the input.Then the light evolution for both positive and negative nonlinearities is fully equivalent in the framework of the tight-binding model [5,6], and we take γ ≥ 0. In our simulations presented below, we use the following values which are typical for the experiments with optical waveguide arrays: d = 9 μm, n 0 = 2.35, λ = 532 nm, C 0 = 0.13 mm −1 , γ = 1.9.Normalization chosen is such that x is measured in μm and z is measured in mm.We use the discrete model for the calculations presented in this paper, however we have confirmed the validity of our results by simulating the full parabolic equations for the continuous electric field envelopes.We note that in case of the strong coupling between the waveguides, the long range coupling between non-nearest neighbours becomes important.In this case, exact dynamical localization can also be ralized in arrays of curved waveguides [7], and effects similar to the presented in this paper may be expected to take place in the nonlinear regime.

Discrete diffraction and self-collimation
In order to specifically distinguish the effects due to diffraction management, we consider the light propagation in the waveguide arrays with symmetric bending profiles, i.e. x 0 (z − z a ) = x 0 (z a − z) for a given coordinate shift z a , since asymmetry may introduce other effects due to the modification of refraction [15][16][17][18][19].Then, after the full bending period (z → z + L) the beam diffraction remains the same as in a straight waveguide array with the effective coupling coefficient [5,6] As a specific example, we consider a curved waveguide array with the sinusoidal bending of the form x 0 (z) = A{cos[2πz/L]− 1}, for which the effective coupling coefficient is C eff = C 0 J 0 (2πωA/L), where J 0 is the Bessel function of the first kind of zero order.The effective diffraction can be made either normal, zero, or anomalous depending on the value of the bending amplitude [see Fig. 1(b)].
Cancellation of the effective coupling and periodic beam self-collimation takes place at low powers when A = A sc = ξ L[2πω] −1 , where ξ 2.405 is the first root of the function J 0 .For example, for the bending period L = 15 mm self-collimation occurs when A sc = 23.0 μm [see Fig. 2(a), where we use the modal representation described above in order to reconstruct the optical field distribution].The beam width is determined as the width of the transverse crosssection function centered at the current center of mass of the beam where 75% of the beam power is concentrated.When the bending amplitude differs from the self-collimation value, the beam experiences discrete diffraction at low powers, similar to the effect observed in straight waveguide arrays [1] [see Fig. 2(b)].

Nonlinear beam propagation and control
When the power of the input beam increases, nonlinearity destroys the self-collimation condition of light by changing the refractive index of the waveguide material.Initially, we observe that the beam shape experiences irregular distortion, such that the periodicity of the selfcollimation is lost [see Fig. 3(a)].However, the beam does not broaden significantly, and it still experiences approximate self-restoration at some points.
When the input power increases further, the beam no longer experiences self-restoration,  and the nonlinear diffusion takes place, where the beam experiences significant broadening and self-defocusing, as shown in Fig. 3(b).This self-defocusing is intrinsically limited due to the diffraction cancellation in the waveguide array.After propagation over some distance the beam broadens and its intensity is reduced accordingly.Therefore, the further beam spreading stops when the average beam width achieves a certain value.Such a peculiar nonlinear beam dynamics has no analogies in bulk media [11] or discrete systems [1] analyzed before.At even higher input powers, transitional self-trapping of the beam is observed.The beam initially becomes self-trapped upon the launch into the array, but after propagation for some distance (which somewhat depends on the input power), it broadens rapidly and experiences again nonlinear diffusion [see Fig. 3(c)].Finally, at some critical power we observe a sharp transition from the nonlinear diffusion to the discrete self-trapping over the whole length of the array, and the discrete lattice soliton is formed [see Fig. 3(d)].
From the practical point of view, the output width of the beam which exits the waveguide array is of the main interest.In Figs.4(a) and 4(b) the beam width w out at the output facet of the array is shown as a function of the input power P normalized to the power P 0 required for the formation of one-site discrete lattice soliton in the straight array, where P 0 5C 0 /γ [20].Shading shows the range between the minimum and the maximum beam widths during the propagation over the last two periods of the array [note that the maximum width in Fig. 4(b) is the same as the output width].All the output, the minimum, and the maximum widths are normalized to the input beam width.In Fig. 4(a), one can clearly identify all the discussed above different nonlinear propagation regimes which happen in diffraction free waveguide arrays [compare this figure with the animation in Fig. 3].In Fig. 4(b), where we take the bending amplitude to be half of that required for the self-collimation, the regime of the intermediate nonlinear beam diffusion can also be identified, however the critical power required for the formation of the lattice soliton is different.We find that this critical power is considerably higher in arrays of periodically curved waveguides than in straight waveguide arrays.It grows substantially with the magnitude of the bending amplitude, as shown in Fig. 4(c).

Fig. 1 .
Fig. 1.(a) Schematic of an array of periodically curved waveguides.(b) Effective coupling C eff in the curved waveguide array normalized to the coupling C 0 in the straight array, as a function of the bending amplitude A normalized to the amplitude A sc which corresponds to the self-collimation condition.Point 'a' corresponds to the self-collimation, when the bending amplitude is A = A sc , while for the point 'b' we take A = 1/2A sc .

Fig. 2 .
Fig. 2. Linear propagation of light in a curved waveguide array.Shown are the absolute values of the amplitudes of the modes of individual waveguides (left) and corresponding optical field distribution (center).Right: beam width w [normalized to the input width] as a function of the propagation distance.(a) Bending amplitude corresponds to the selfcollimation condition [point 'a' in Fig. 1(b)].(b) Bending amplitude is half of that in (a) [point 'b' in Fig. 1(b)].Bending period is L = 15 mm.Waveguide array length is 150 mm.

Fig. 3 .Fig. 4 .
Fig. 3. (1.5MB) Nonlinear propagation of light in an array of periodically curved waveguides.Shown are the absolute values of the amplitudes of the modes of individual waveguides (left) and corresponding optical field patterns (center).Right: beam width w [normalized to the input width] as function of the propagation distance.Waveguide array is the same as in Fig. 2(a).The input power is (a) P/P 0 = 0.70, (b) P/P 0 = 1.7,(c) P/P 0 = 2.7, and (d) P/P 0 = 3.4,where P 0 is the power required for the formation of one-site discrete soliton in the straight array.The animation shows the beam propagation dynamics and the output beam width as the input power increases from P/P 0 = 0.00029 to P/P 0 = 3.4.