We introduce soliton percolation phenomena in the nonlinear transport of light packets in suitable optical lattices with random properties. Specifically, we address lattices with a gradient of the refractive index in the transverse plane, featuring stochastic phase or amplitude fluctuations, and we discover the existence of a disorder-induced transition between soliton-insulator and soliton-conductor regimes. The soliton current is found to reach its maximal value at intermediate disorder levels and to drastically decrease in both, almost regular and strongly disordered lattices.
©2007 Optical Society of America
Analogies between the electron dynamics in perfect crystals and light propagation in periodic optical media guide the elucidation of a variety of new physical phenomena and related applications [1–4]. Bloch oscillations and Zener tunneling [5–7] are just salient examples of effects that arise in the linear regime of light propagation in periodic optical media, while discrete and lattice solitons [8–13] as well as complex soliton trains  are examples of the possibilities afforded by the addition of nonlinearity. All these phenomena occur in regular, periodic or weakly modulated lattices.
Nevertheless, disordered solid-state periodic materials are known to exhibit a wealth of unique electron dynamics phenomena. In optics, the interplay between disorder and nonlinearity was studied, in particular, in systems that may be modeled by the nonlinear Schrödinger equation with random-point impurities [15,16] and in discrete systems [17–19]. Soliton propagation in random potentials and the combined effect of periodic and random potentials on the transmission of moving and formation of stable stationary excitations has been addressed in a number of previous studies (see, for example, Refs. [20–23] and reviews [24–26]). Different regimes of light localization in linear disordered photonic crystals were addressed in Ref. . Band theory of light localization in one-dimensional linear disordered systems was developed in Ref. , where it was illustrated that Bragg reflection is responsible for Anderson localization of light. The phenomenon of Anderson-type localization of walking spatial solitons in the nonlinear optical lattices with random frequency modulation was studied in Ref. . Anderson localization in random optical lattices imprinted in photorefractive crystal has been recently observed in a landmark experiment by Segev and co-workers . The phenomena predicted in disordered optical lattices also occur in Bose-Einstein condensates [31–35] and vice versa.
A universal feature of wave packet and particle dynamics in disordered media in different areas of physics is percolation [36,37]. Percolation occurs in all types of physical settings, including high-mobility electron systems , Josephson-junction arrays , two-dimensional GaAs structures near the metal-insulator transition , or charge transfer between superconductor and hoping insulator , to cite a few.
In this paper we introduce, for the first time to our knowledge, the nonlinear optical analog of biased percolation that is related to disorder-induced soliton transport in randomly modulated optical lattices with Kerr-type focusing nonlinearity in the presence of linear variation of the refractive index in the transverse plane, thus generating a constant deflecting force for light beams entering the medium. When such a force is too small, solitons in perfectly periodic lattice are trapped in the vicinity of the launching point due to Peierls-Nabarro potential barriers, provided that the launching angle is smaller than a critical value . Under such conditions soliton transport is suppressed, and thus the lattice acts as a soliton insulator. However, random modulations of the lattice parameters turns soliton transport possible again, with the key parameter determining the soliton current being the standard deviation of phase/amplitude fluctuations. Hereby we discover that the soliton current in lattices with amplitude and phase fluctuations reaches its maximal value at intermediate disorder levels and that it drastically reduces in both, almost regular and strongly disordered lattices. This suggests the possibility of a disorder-induced transition between soliton insulator and soliton conductor lattice states.
Our analysis is based on the nonlinear Schrödinger equation describing propagation of a laser beam in a medium with focusing Kerr-type nonlinearity and spatial modulation of the refractive index in the transverse direction, namely
Here q is the dimensionless complex amplitude of light field; the transverse Laplacian writes Δ⊥=∂2/∂η2 or Δ⊥=∂2/∂η2+∂2/∂ζ2 in the case of one- or two-dimensional geometries, respectively; the transverse η, ζ and longitudinal ξ coordinates are scaled in terms of beam radius and diffraction length, respectively; the function R describes the lattice profile, while the parameter α characterizes the rate of linear growth of the refractive index in the direction of η axis. Notice that such a refractive index gradient produces a constant deflecting “force” for the light beams entering the medium. Here we concentrate in rather shallow, high-frequency lattices featuring a harmonic regular component with small additive random amplitude or phase fluctuations, with the general functional form R(η)=pcos(Ωη)+σaρ(η), R(η)=pcos[Ωη+σpρ(η)], respectively. Here p is the depth of the refractive index modulation, Ω is the modulation frequency, ρ(η) is a random function with zero mean value 〈ρ(η)〉=0, and unit variance 〈ρ 2(η)〉=1(the angular brackets stand for statistical averaging), while the parameters σ a and σ p define the depth of the amplitude or phase stochastic modulations. We assume the correlation function of the random field ρ(η) to be Gaussian 〈ρ(η1)ρ(η2)〉exp[-(η2-η1)2/L 2 cor] with a correlation length L cor larger than the regular lattice period L cor>2π/Ω. In the two-dimensional case, we considered lattices with amplitude and phase fluctuations having the functional shapes a R=pcos(Ωη)cos(Ωζ)+σaρ(η,ζ) and R=pcos[Ωη+σpρ(η,ζ)]cos[Ωζ+σpρ(η,ζ)], respectively, where random field ρ(η,ζ) has the same statistical properties as its one-dimensional counterpart. Such refractive index landscapes can be induced in photorefractive crystals (see Ref. ) by combining regular periodic lattices with random nondiffracting patterns. Random component can be generated, e.g., by illuminating a narrow annular slit with a random transmission function placed at the focal plane of a lens. Modification of the transmission function, e.g., by rotating a suitable light diffuser placed after the slit, results in the formation of different random lattice realizations. In the particular case of optical lattices imprinted in SBN crystals biased with a dc electric field ~5kV/cm and laser beam with widths 5µm at λ=532 nm, a length ξ~1 corresponds to some 0.7 mm, Ω=6 sets a lattice period ~5µm, the parameter p=1 corresponds to a refractive index variation ~10-4, and q=1 corresponds to a peak intensity of the order of 100 mW/cm2.
To gain physical insight into the full nonlinear wave dynamics we start our analysis with an effective particle approach , which uses the integral coordinates of the beam center
is the total energy flow that remains constant upon propagation. In the deterministic case σ a,p→0, the effective particle approach yields coupled equations for the beam center coordinates, which in a two-dimensional geometry write d 2 ηc/dξ 2+WΩcos(Ωζc)sin(Ωη c)=α and d 2ζc/dξ 2+WΩcos(Ωηc) sin(Ωζc)=0. Here we take a trial function |q(η,ζ,ξ)|=q 0×exp[-χ 2(η-η c)2]exp[-χ 2(ζ-ζ c)2], with the amplitude q 0 and form-factor χ. The parameter W=pexp(-Ω2/2χ 2) characterizes the height of the Peierls-Nabarro potential barrier, which grows linearly with the lattice modulation depth p and diminishes rapidly with increasing carrying frequency Ω. Other trial functions provide qualitatively similar results. The evolution equation for the beam center in one-dimensional geometries is obtained by setting ζ c≡0. The critical value of rate α of linear increase of refractive index is directly related to height of the barrier, and is given by α cr=WΩ. At α>α cr the beam starts to drift across the lattice in the positive direction of η axis, while at α<α cr it is trapped in the vicinity of the launching point. However, in the presence of fluctuations the height of potential barrier becomes a random function of the transverse coordinates, thereby generating a non-vanishing probability of soliton transport, or percolation, even at cr α<α cr. This phenomenon may be viewed as a disorder-induced transition between soliton-insulator and soliton-conductor regimes. Series of comprehensive numerical simulations for different regions of the parameter space revealed that this transition appears most clearly in shallow, high-frequency lattices, where the soliton mobility is relatively high and the radiative losses are small in unperturbed lattices.
To elucidate the exact implications of this result, we conducted a comprehensive Monte- Carlo numerical investigation, by integrating Eq. (1) numerically with a split-step Fourier method. We propagated solitons up to a given large distance, termed ξ end, for different series of random realizations of the refractive index profiles Rk(η,ζ)1≤k≤Nr, with Nr=103. We first focus on a detailed quantitative analysis of one-dimensional geometries. We set a lattice frequency at Ω=6 and launching point at η 0=-10π/Ω(coinciding with one of local lattice maxima). We used a wide integration domain η∈[-100,100] to eliminate any boundary effects. An input beam with the shape q|ξ=0=sech(η-η 0) and an input energy flow U=2 was selected. Such input would propagate undistorted in a uniform cubic medium. The critical refractive index slope for such beam in the optical lattice of depth p=0.4 amounts to α cr=0.0027, and such beam would be trapped in regular lattice with α=0.0025. Notice that the shapes of truly stationary solitons supported by the lattice with a refractive index gradient may be slightly asymmetric and that they are fairly similar to those existing in lattices with a linear amplitude modulation [see, e.g., Fig. 1(d) of Ref. ]. Such stationary states exist above a power threshold that depends on the refractive index gradient. For the seek of generality and experimental relevance, in what follows we consider symmetric inputs with a bell-shape in the form q|ξ=0=sech(η-η 0).
Figure 1 shows an illustrative set of soliton intensity distributions in a lattice with random phase and amplitude fluctuations for different standard deviation values obtained by direct numerical integration of Eq. (1). When the local restoring force acting on the soliton in the random lattice compensates the random-induced deflecting force, the soliton remains trapped in the input lattice channel. Otherwise, it is accelerated and moves along a complex path. The kinetic energy of the associated effective particle gradually increases and the probability of its trapping in one of subsequent lattice channels or backward reflection in the vicinity of regions with locally decreased refractive index decreases. Yet, it does not vanish completely [see, e.g., Fig. 1(b)]. The radiation, that appears unavoidably when solitons cross the lattice channels, is weak until the local propagation angle does not approach the Bragg angle θ B=Ω/2. Note that the soliton propagation trajectories approach parabolic ones in lattices with weak fluctuations σ a,p and that they can be quite complex already when σ a,p~0.2. In the case of quasi-parabolic trajectory the overall displacement of soliton can reach values ~αξ2/2, which for a propagation distance ξ~150 can exceed the input soliton width by more than one order of magnitude, as readily visible in Fig. 1(a).
It is worth stressing that the dynamics depicted in Fig. 1 cannot be considered as an initial stage of Bloch oscillations that take place in linear systems [5–7]. In our model the strongly nonlinear excitations experience only small modulations due to the presence of a shallow lattice and their diffraction spreading is balanced mostly by nonlinearity. Importantly, the corresponding localized excitations move across the regular lattice over the significant distances without changing their shape, provided that their velocities are sufficient to overcome the so-called Peierls- Nabarro potential barrier. If such nonlinear wave packet is accelerated by a constant force and the deflection angle approaches the Bragg one, it spreads rapidly due to strong radiative losses, as shown in Refs. [44,45] instead of being reflected as a single object, as it occurs upon Bloch oscillations in arrays of evanescently-coupled waveguides described by discrete model. Note that when the propagation angle is smaller than the Bragg angle, solitons always preserve their internal structure even upon reflection on the lattice defects [an effect visible in Fig. 1(b)], in contrast to the phenomena generated by Bloch oscillations.
A key percolation feature to be monitored is the dependence of the soliton current on the level of randomness present in the system. We calculated the total number of solitons N c launched with zero velocity at η 0=-10π/Ω, that reach the point |η 0|=10π/Ω at any random distance ξ c for a set of N r=103 realizations of optical lattices with amplitude or phase fluctuations. This enables to introduce soliton current, defined as J=N c/N r. Figure 2(a) shows the dependence of such quantity on the standard deviation of phase fluctuations σ p for different values of the lattice depth p. Since the value of deflecting force was below the critical value (α<α cr), a rapid growth of the soliton current with increasing σ p occurs. This is entirely due to the formation of continuous cluster of opened lattice sites connecting the input and output points. This phenomenon is intuitively analogous to the disorder-induced phase transition that arises between insulating and conducting states. However, further increasing σ p results in a decrease of soliton current due to growing probability of soliton trapping or reflection in the vicinity of lattice sites that are sufficiently strongly distorted.
Therefore, one of the central results of this work is that soliton current in the lattices with fluctuations reaches its maximal value at intermediate disorder level and substantially decreases in almost regular (σ p→0) and strongly disordered (σ p→1) lattices. Growth of the lattice depth p results in a drop off of the maximal value of the soliton current due to the reduction of soliton mobility in deeper lattices. The soliton current maximum is achieved at slightly larger disorder levels in deeper optical lattices. Figure 2(b) illustrates the dependence of the soliton current on the rate α of the linear refractive index variation (i.e., the equivalent to the applied voltage in electro-dynamics). Instead of the sharp jump from zero to unity at α=α cr that takes place in regular lattices, the current is growing smoothly in disordered lattices. Remarkably, the current may be nonzero even for very small slopes α that are far below the critical value and can be less than unity for α>α cr. Growth of the disorder level results also in an increase of the transition region of the J(α) dependence that can be defined as a region where the current increases from 0.1 to 0.9. We have found qualitatively similar dependencies of the current on the disorder level and slope α in lattices with amplitude fluctuations [Figs. 2(c) and 2(d)]. Notice, however, a remarkably higher sensitivity of the soliton current on the standard deviation in the case of amplitude fluctuations. Note that the results depicted in Fig. 2 remain almost unchanged for much longer integration distances, exceeding ξ=150. This is so because solitons that become trapped at a certain lattice realization typically remain in the vicinity of the launching point for distances that far exceed 150 diffraction lengths, while the probability of trapping for solitons that start walking and being accelerated typically decreases with the propagation distance, so that such solitons will contribute to soliton current on larger distances too.
As mentioned above, the distance ξ c at which soliton reaches the point η=|η 0|, is a random value because soliton follows a random path to reach it. Figure 3 shows histograms of this random value for lattices with either phase or amplitude fluctuations. The histograms show ratio between the number of lattice realizations where ξ c falls into fixed intervals, and the total number of events when soliton eventually reaches the point |η 0|=10π/Ω. The maximum of the histogram corresponds to the most probable distance at which soliton reaches the point η=|η 0|(one has 〈ξ c〉≈78 in the case of phase fluctuations and 〈ξ c〉≈72 in the case of amplitude fluctuations for selected disorder levels), while the histogram width characterizes the standard deviation of ξ c. We found that in the regime corresponding to a growth of the soliton current with σ a,p [see Figs. 2(a) and 2(c)] an increase of the fluctuation level causes increase of the histogram width (standard deviation of ξ c) and an increase of the distance 〈ξ c〉 mainly due the expansion of the histogram tail toward the region of higher ξ c values. Note the analogy of this phenomenon with the scattering of ultra-short pulses in turbid media, where an increase of the turbidity results in larger group delays and pulse broadening. Note also that the sensitivity of ξ c on the standard deviation is remarkably higher in the case of amplitude fluctuations.
We also investigated the possibility of soliton percolations in two-dimensional geometries. In these geometries, soliton mobility in a medium with cubic nonlinearity is strongly restricted because of the possibility of development of collapse or decay of the input beams. Thus, soliton moving across two-dimensional lattices may drastically broaden in between sites of the lattice and can hardly jump between more than two sites. To avoid such type of behavior, we added phenomenological nonlinearity saturation into Eq. (1) by rewriting the nonlinear term in the form -q|q|2/(1+S|q|2). Notice that nonlinearity saturation is typical for photorefractive crystals, where disordered lattices might be imprinted optically. It is worth mentioning that in the previous cases we addressed the canonical model of Kerr nonlinearity in a (1+1)- dimensional geometry because of its generality. The model is directly applicable to Bose-Einstein condensates and to light propagation in photorefractive crystals (with rather intensive background illumination), as well as in suitable semiconductor lattices. In (2+1)-dimensions Kerr models exhibit collapse, which may be arrested by saturation, e.g., in photorefractive media with a suitable background illumination.
Figure 4(a) shows a typical profile of a two-dimensional lattice with random phase modulation. The fundamental soliton corresponding to the saturation parameter S=0.2 and energy flow U=8 was launched in the point η 0=-10π/Ω, ζ=0 into such lattice in the presence of linear refractive index variation with slope α=0.0015. Figure 4(b) shows an illustrative set of random trajectories corresponding to different lattice realizations. One of them (a small almost closed loop) illustrates the possibility of soliton trapping in the vicinity of the launching point, while other trajectories confirm that soliton percolation (namely, drift in the positive direction of η axis) does occur in two-dimensional settings. One of the characteristic features of two-dimensional soliton percolation, in comparison with one-dimensional case, is a faster increase of soliton current with disorder level due to much richer set of possible propagation paths in the two-dimensional geometry.
We thus conclude by stressing that we have introduced the phenomenon of percolation of spatial solitons in optical lattices with random phase or amplitude fluctuations. The central discovery reported is the possibility of disorder-induced transition between soliton insulator and soliton conductor lattice states. Such phenomenon opens up further analogies with solid-state and statistical physics, whose experimental exploration has become feasible after the recent observation of Anderson localization in random optically induced lattices .
This work has been supported in part by the Government of Spain through the Ramon-y-Cajal program and the grant TEC2005-07815, and by CONACYT through the grant 46552.
References and links
3. E. Yablonovitch, “Photonic crystals: Semiconductors of light,” Scientific American 285, 47 (2001). [CrossRef]
4. J. D. Joannopoulos, “Photonic-bandgap microcavities in optical waveguides,” Nature 386, 143 (1997). [CrossRef]
5. T. Pertsch, P. Dannberg, W. Elflein, A. Bräuer, and F. Lederer, “Optical Bloch oscillations in temperature tuned waveguide arrays,” Phys. Rev. Lett. 83, 4752 (1999). [CrossRef]
6. H. Trompeter, W. Krolikowski, D. N. Neshev, A. S. Desyatnikov, A. A. Sukhorukov, Y. S. Kivshar, T. Pertsch, U. Peschel, and F. Lederer, “Bloch oscillations and Zener tunneling in two-dimensional photonic lattices,” Phys. Rev. Lett. 96, 053903 (2006). [CrossRef] [PubMed]
9. V. V. Konotop, O. A. Chubykalo, and L. Vazquez, “Dynamics and interaction of solitons on an integrable inhomogeneous lattice,” Phys. Rev. E 48, 563 (1993). [CrossRef]
10. J. W. Fleischer, T. Carmon, M. Segev, N. K. Efremidis, and D. N. Christodoulides, “Observation of discrete solitons in optically induced real time waveguide arrays,” Phys. Rev. Lett. 90, 023902 (2003). [CrossRef] [PubMed]
11. J. W. Fleischer, M. Segev, N. K. Efremidis, and D. N. Christodoulides, “Observation of two-dimensional discrete solitons in optically induced nonlinear photonic lattices,” Nature 422, 147 (2003). [CrossRef] [PubMed]
13. H. Martin, E. D. Eugenieva, Z. Chen, and D. N. Christodoulides, “Discrete solitons and soliton-induced dislocations in partially coherent photonic lattices” Phys. Rev. Lett. 92, 123902 (2004). [CrossRef] [PubMed]
16. V. A. Hopkins, J. Keat, G. D. Meegan, T. M. Zhang, and J. D. Maynard, “Observation of the predicted behavior of nonlinear pulse propagation in disordered media,” Phys. Rev. Lett. 76, 1102 (1996). [CrossRef] [PubMed]
17. V. V. Konotop, D. Cai, M. Salerno, A. R. Bishop, and N. Grønbech-Jensen, “Interaction of a soliton with point impurities in an inhomogeneous discrete nonlinear Schrödinger system,” Phys. Rev. E 53, 6476 (1996). [CrossRef]
18. P G. Kevrekidis, Y. S. Kivshar, and A. S. Kovalev, “Instabilities and bifurcations of nonlinear impurity modes,” Phys. Rev. E 67, 046604 (2003). [CrossRef]
19. T. Pertsch, U. Peschel, J. Kobelke, K. Schuster, H. Bartelt, S. Nolte, A. Tunnermann, and F. Lederer, “Nonlinearity and disorder in fiber arrays,” Phys. Rev. Lett. 93, 053901 (2004). [CrossRef] [PubMed]
20. R. Knapp, “Transmission of solitons through random-media,” Physica D 85, 496 (1995). [CrossRef]
21. P. L. Christiansen, Y. B. Gaididei, M. Johansson, K. O. Rasmussen, D. Usero, and L. Vazquez, “Stabilization of nonlinear excitations by disorder,” Phys. Rev. B 56, 14407 (1997). [CrossRef]
22. Y. B. Gaididei, D. Hendriksen, P. L. Christiansen, and K. O. Rasmussen, “Stationary states of the two-dimensional nonlinear Schrodinger model with disorder,” Phys. Rev. B 58, 3075 (1998). [CrossRef]
23. J. Garnier, “Propagation of solitons in a randomly perturbed Ablowitz-Ladik chain,” Phys. Rev. E 62, 026608 (2001). [CrossRef]
24. Fluctuation Phenomena: Disorder and Nonlinearity, Ed. by A. R. Bishop, S. Jimenez, and L. Vazquez (World Scientific, Singapore, 1995).
25. Nonlinearity with Disorder, Ed. by F. Kh. Abdullaev, A. R. Bishop, S. Pnevmatikos, and E. N. Economou (Springer, Berlin, 1992). [CrossRef]
26. F. Abdullaev, Theory of Solitons in Inhomogeneous Media (Wiley, New York, 1994).
27. Y. A. Vlasov, M. A. Kaliteevski, and V. V. Nikolaev, “Different regimes of light localization in a disordered photonic crystal,” Phys. Rev. B 60, 1555 (1999). [CrossRef]
28. A. P. Vinogradov and A. M. Merzlikin, “Band theory of light localization in one-dimensional disordered systems,” Phys. Rev. E 70, 026610 (2004). [CrossRef]
29. Y. V. Kartashov and V. A. Vysloukh, “Anderson localization of solitons in optical lattices with random frequency modulation,” Phys. Rev. E 72, 026606 (2005). [CrossRef]
30. T. Schwartz, G. Bartal, S. Fishman, and M. Segev, “Transport and Anderson localization in disordered two-dimensional photonic lattices,” Nature (London) 446, 52 (2007). [CrossRef]
32. R. C. Kuhn, C. Miniatura, D. Delande, O. Sigwarth, and C. A. Muller, “Localization of matter waves in twodimensional disordered optical potentials,” Phys. Rev. Lett. 95, 250403 (2005). [CrossRef] [PubMed]
33. T. Schulte, S. Drenkelforth, J. Kruse, W. Ertmer, J. Arlt, K. Sacha, J. Zakrzewski, and M. Lewenstein, “Routes towards Anderson-like localization of Bose-Einstein condensates in disordered optical lattices,” Phys. Rev. Lett. 95, 170411 (2005). [CrossRef] [PubMed]
35. D. Clement, A. F. Varon, J. A. Retter, L. Sanchez-Palencia, A. Aspect, and P. Bouyer, “Experimental study of the transport of coherent interacting matter-waves in a 1D random potential induced by laser speckle,” New J. Phys. 8, 165 (2006). [CrossRef]
36. G. R. Grimmett, Percolation (Springer, Berlin, 1999).
37. B. I. Shklovskii and A. L. Efros, Electronic Properties of Doped Semiconductors, Springer Series in Solid-State Sciences (Springer, Berlin, 1984).
38. S. Das Sarma, M. P. Lilly, E. H. Hwang, L. N. Pfeiffer, K. W. West, and J. L. Reno, “Two-dimensional metal-insulator transition as a percolation transition in a high-mobility electron system,” Phys. Rev. Lett. 94, 136401 (2005). [CrossRef] [PubMed]
40. G. Allison, E. A. Galaktionov, A. K. Savchenko, S. S. Safonov, M. M. Fogler, M. Y. Simmons, and D. A. Ritchie, “Thermodynamic density of states of two-dimensional GaAs systems near the apparent metal-insulator transition,” Phys. Rev. Lett. 96, 216407 (2006). [CrossRef] [PubMed]
42. Y. S. Kivshar and D. K. Campbell, “Peierls-Nabarro potential barrier for highly localized nonlinear modes,” Phys. Rev. E 48, 3077 (1993). [CrossRef]
43. Y. V. Kartashov, V. A. Vysloukh, and L. Torner, “Oscillations of two-dimensional solitons in harmonic and Bessel optical lattices,” Phys. Rev. E 71, 036621 (2005). [CrossRef]
44. A. V. Yulin, D. V. Skryabin, and P. St. J. Russell, “Transition radiation by matter-wave solitons in optical lattices,” Phys. Rev. Lett. 91, 260402 (2003). [CrossRef]
46. Y. V. Kartashov, V. A. Vysloukh, and L. Torner, “Soliton control in chirped photonic lattices,” J. Opt. Soc. Am. B 22, 1356 (2005). [CrossRef]