Vector mixed-gap surface solitons

We elucidate the properties of mixed-gap vector surface solitons supported by the interface between a uniform medium and an optical lattice imprinted in a Kerr-type nonlinear media. The components of such mixed-gap solitons emerge from different gaps of lattice spectrum and their mutual trapping results in the formation of stable vector states. The unstable soliton component is stabilized by the cross-coupling with the stable component. We show that vector mixed-gap surface solitons exhibit a new combination of properties of vectorial surface waves and gap solitons.

Cross-coupling between several light waves can considerably enrich the dynamics of their propagation. Such coupling provides a strong stabilizing mechanism, which results in the formation of light patterns that often do not even exist in scalar settings. Vector solitons forming due to the coupling of two waves have been encountered in a number of geometries, including uniform media with both coherent [1][2][3][4] and incoherent [5][6][7][8] interactions. On the other hand, several previous investigations have shown that a transverse periodic modulation of the refractive index strongly affects the properties of vector solitons. For example, twodimensional optical lattices [9] and arrays of weakly coupled waveguides support strongly localized vectorial modes that have no analogs in uniform media [10][11][12]. Similarly, mixedgap vector solitons whose components emerge from different gaps of the Floquet-Bloch spectrum of the periodic structure have been introduced recently [13,14]. Optical lattices support a variety of other stable vector solitons [15].
The presence of interfaces between different materials substantially modifies the properties of both scalar and vector solitons. It was shown recently that the interface between a waveguide array and a uniform material allows formation of new types of scalar surface solitons [16,17]. Prior to [17], optical surface waves were also observed at the nonlinear interface of a photorefractive crystal with diffusion nonlinearity [18], as well as at linear interfaces of periodic media [19]. Observation of highly nonlinear surface waves at interfaces of natural materials is a challenge due to the large powers required for their excitation (see [20][21][22][23] for reviews). Interfaces of periodic media support also gap scalar solitons that were predicted theoretically in [24] and observed experimentally in [25]. Gap solitons may also exist at surfaces that are periodically modulated in the direction of light propagation [26].
Besides scalar solitons, interfaces of periodic structures support vector surface waves. The simplest vector surface solitons were studied in [27], for the case of coherently interacting components that emerge from a semi-infinite gap. Nevertheless, the most interesting situation is encountered when the components of the vector surface soliton originate from different gaps. Several examples of such solitons were obtained recently at the interfaces of discrete binary waveguide arrays with focusing nonlinearity [28]. In this paper we study mixed-gap vector surface waves arising at the interface of a uniform medium and an optical lattice imprinted in a Kerr nonlinear medium, for focusing and defocusing nonlinearities. We reveal that cross-coupling results in the formation of stable vector states, where instability of one of the components is suppressed by the interaction with the stable component.
We consider the propagation of two mutually incoherent laser beams at the interface between a cubic uniform medium and a periodic optical lattice. For concreteness, here we assume that the strength of the cross-phase-modulation is equal to that of the self-phase-modulation. In this case the light propagation is described by two coupled nonlinear Schrödinger equations for the dimensionless complex amplitudes q : Here the transverse η and longitudinal ξ coordinates are scaled in terms of the beam width and diffraction length, respectively; the parameter σ stands for focusing/defocusing; is the lattice depth, and the function for and for describes the transverse refractive index profile. The quantity Ω is the lattice frequency. Eqs.
(1) conserve the total energy flow and . Focusing medium. Here and throughout the paper we set Ω .
Interfaces between uniform and periodic media can be fabricated technologically, e.g. by etching array of waveguides on a nonlinear substrate at η [11,12]. Optical induction offers a potentially powerful alternative. In this case, a periodic pattern is imprinted in a photorefractive crystal by using several interfering plane waves and subsequently erased at η with an intense green wave propagating in the direction orthogonal to and axes.
We search for vector soliton solutions of Eq. (1) in the form q w , where w are real functions and b are propagation constants. We are interested in solitons whose maxima seat on the lattice channel nearest to the interface. The periodicity of the lattice imposes restrictions on the available values of propagation constants of localized solitons, that should fall into gaps of the Floquet-Bloch lattice spectrum (white regions in Fig. 1(a)). In the gray regions (bands), lattices can support only delocalized Bloch waves. Surface solitons residing at the interface with the lattice penetrate into the uniform medium too, which imposes the conditions b on the propagation constants. When b , the components w have similar shapes and their properties are analogous to those for scalar gap solitons [24]. We assume that b and b fall into neighboring gaps of the Floquet-Bloch spectrum, so that w feature different structures. In a focusing medium ( the components of the simplest mixed-gap vector soliton emerge from the semi-infinite and the first finite gaps (see Fig. 1(a)). Representative profiles are shown in Figs. 2(a)-2(c). , . Panels (a)-(c) correspond to p and focusing medium. Panels (d)-(f) correspond to p and defocusing medium. In gray regions R , while in white regions R .
Physically, each component of the vector gap surface soliton forms by the nonlinear coupling between the incident and the reflected waves having equal wavevectors along the longitudinal coordinate (i.e. ) and opposite wavevectors in the transverse coordinate (i.e. η ), when both of them experience Bragg scattering from the periodic structure. The symmetry of the soliton components reflects the symmetry of the Bloch waves near the corresponding gap edges. Such symmetry is determined by the effective sign of the diffraction required for the formation of localized states in focusing (normal diffraction) and defocusing (anomalous  In the case σ the component w is always positive, while w exhibits multiple oscillations at η . Both components decay exponentially as η . For fixed p and , the energy flow U is a monotonically increasing function of b ( Fig. 1(b)). There exist a lower cutoff and an upper b cutoff on for the existence of composite vector solitons. Figure 1(c) shows the domain of existence on the ( , plane for . We found vector solitons for arbitrary large values of b . At large b the component w is narrow and its amplitude is much higher than that of w . The upper cutoff b coincides with the upper edge of the first gap ( b at ) and for b the second component deeply penetrates into the lattice region thus approaching a delocalized Bloch wave. In the upper cutoff w remains well localized (Fig. 2(c)). As b , the second component vanishes, while its transverse extent depends on position of b in the first gap. For , the upper cutoff b abruptly decreases when b also decreases. In this regime, vanishes in the upper cutoff. Figure 1(d) shows the energy sharing S U between the vector soliton components for this parameter range. With further decrease of b to , the lower cutoff b gradually approaches the lower edge of first gap, so that w at and b the soliton closely resembles a Bloch wave corresponding to this gap edge ( Fig. 2(a)). we did not find solitons for which w has a maximum in the first lattice channel. For b strong secondary peaks appear in w , while the maximum of shifts into the lattice so that the two components decouple. Similar domains of existence were encountered for various lattice depths. For small enough p , when the propagation constant corresponding to the lower edge of the first gap becomes negative, the domain of existence shown in Fig. 1(b) is bounded from below by the line b . also support mixed-gap solitons. In this case, vector soliton components emerge from the first and the second finite gaps of the Floquet-Bloch spectrum ( Fig. 3(a)). At fixed b and the energy flow is a monotonically decreasing function of b ( Fig. 3(b)). At p the domain of existence of composite states appears at b corresponding to the lower edge of the first gap (Fig. 3(c)). In the vicinity of this point, w turns out to be weakly localized and closely resembles a Bloch wave in the region . The lower cutoff is given by b , and as b , w deeply penetrates into the uniform medium ( Fig. 2(d) , since the maximum of w shifts into the lattice region and w decouple. The energy sharing between components of the vector soliton vs b is shown in Fig.  3(d). One concludes that the rich internal structure of the vector surface waves is determined by the location of the propagation constants b in the band-gap spectrum, locations which are affected by the presence of the interface. . We verified that in the absence of w , the 1,2 ( ) η 1 scalar surface solitons that emerge from the first gap in focusing media and from the second gap in defocusing media (i.e. those having internal structure similar to that of w in Fig. 2) are strongly unstable. In contrast, scalar solitons having internal structure of w can be stable in large parts of their existence domain. Intuitively, this suggests the possibility to stabilize composite vector solitons made of a relatively weak w and a strong w . This occurs in wide regions near the lower cutoff b in focusing media and near the upper cutoff b in defocussing media. Thus, at σ and vector solitons are stable for b at b , while at b the stability domain is given by b . We found that the stability domain substantially extends with b and the upper edge of the instability domain gradually approaches b . Stable propagation of the vector surface soliton depicted in Fig. 2(b) is shown in Figs. 4(a) and 4(b). At σ the stability domain at p is given by b for and b for b . The width of the stability domain is maximal for . Stable propagation of the vector soliton depicted in Fig. 2(e) is shown in Figs. 4(c) and 4(d). In all cases, the stable vector solitons retain their input structure over indefinitely long distances even in the presence of considerable input random perturbations. However, in both cases, increasing the weight of the w component results in the destabilization of the vector soliton. All instabilities encountered for vector surface solitons are oscillatory. They result in progressively increasing oscillations of amplitudes of the components of unstable soliton, emission of radiation into the lattice and decay into stable solitons. Here we present results of stability analysis for only one lattice depth, but we verified by a detailed analysis that different lattice depths yield qualitatively similar results. In summary, we reported a detailed analysis of the properties and stability of new vector solitons supported by the interface of a uniform medium and optical lattice. The components of such vector solitons emerge from different gaps of the Floquet-Bloch lattice spectrum, and exhibit different internal structures and stability properties. Nevertheless, we revealed that their cross-coupling leads to formation of new types of completely stable, mixed-gap surface vector solitons.