Interlaced linear-nonlinear optical waveguide arrays

The system of coupled discrete equations describing a two-component superlattice with interlaced linear and nonlinear constituents is revisited as a basis for investigating binary waveguide arrays, such as ribbed AlGaAs structures, among others. Compared to the single nonlinear lattice, the interlaced system exhibits an extra band-gap controlled by the, suitably chosen by design, relative detuning. In more general physics settings, this system represents a discretization scheme for the single-equation-based continuous models in media with transversely modulated linear and nonlinear properties. Continuous wave solutions and the associated modulational instability are fully analytically investigated and numerically tested for focusing and defocusing nonlinearity. The propagation dynamics and the stability of periodic modes are also analytically investigated for the case of zero Bloch momentum. In the band-gaps a variety of stable discrete solitary modes, dipole or otherwise, in-phase or of staggered type are found and discussed.


Introduction
Periodic photonic structures in nonlinear dielectric media are a subject of intense theoretical and experimental research in our days. The formation of self-trapped localized solitary modes, among others, is of major importance. 1,2,3 Such solitons have eigenvalues inside the gaps of the band structure and their existence is due to the interplay between the effective lattice diffraction and nonlinearity. A major attribute of theirs, that facilitates their experimental observation 3 , is their robustness under propagation. This attribute makes the periodic photonic structures ideal for applications in integrated photonic devices and waveguide arrays, such as multiport beam coupling, steering, and switching. [5][6][7][8][9][10] Furthermore, the related field of Bose-Einstein condensates loaded in optical lattices [11][12][13][14][15] is of increasing scientific interest and, thus, both fields progress in parallel as far as the theoretical investigation is concerned. The latter and, more specifically, the study of the formation and propagation of spatially localized modes in photonic structures, has been mostly based of either the tight-binding approximation or the coupled-mode theory, thus, rendering simplified discrete models. [16][17][18][19][20] However, these approximations provide accurate modeling only under the corresponding assumptions. A more general model is the nonlinear Schrödinger equation (NLSE) with spatially periodic coefficients, where the longest transverse dimension, y, has been integrated out by appropriate averaging, namely, ( ) ( ) where E is the y-averaged envelope electric field intensity, k 0 =2π/λ (λ is the wavelength of the carrier), β 0 = k 0 <n(x)> (<n(x)> is the mean value of the linear refractive index), n 2 (x) models the periodic nonlinear refractive index in (m/V) 2 and Δn(x) models the periodic modulation of the linear refractive index. Such photonic structures where both linear and nonlinear refractive indices are transversely modulated are of increasing interest lately due to special properties they possess. 21 The particular cases of periodic functions in the form of periodic sequences of Dirac functions 2, 22-24 or piecewise-constant coefficients 25-27 (nonlinear Kronig-Penney model) have previously been considered. Apart of particular cases as the latter, Eq. (1) is not amenable to a straightforward and conclusive analytical treatment.
In the present work the discrete re-casting of Eq. (1) is revisited. The investigation is aiming towards a model of a binary optical waveguide array 28 amenable to a thorough analytical treatment and, thus, predictable and controllable as far as the functions of the conceptual device it models are concerned. The main incentive is the experimental observation of discrete gap solitons in such AlGaAs binary arrays 28 . Periodic modulation of both refractive indices can be conceptually achieved in ribbed AlGaAs waveguides. In Fig. 1 such a configuration is illustrated: The mole fraction of the aluminum in the substrate (S) is higher than the respective one in the guiding ribbed region (W); that leads to a relatively higher refractive index in the latter. The modulated rib structure, on the other hand, renders the effective linear refractive index in the ribs piece-wise constant with higher value (n L ) in the transversely (i.e., in y) and longitudinally (i.e., in x) wider segment and lower (n N ) in the transversely narrower one. However, there exist a geometric refractive index contrast between each rib and the region that separates them from the neighboring uneven ribs. Thus, light can be guided in the region under both the families of ribs. On the other hand, the effective nonlinear refractive index is also piece-wise constant attaining a considerably higher value in the narrower segment. Eventually, focusing will enable light to be efficiently guided in the small ribs while, at the same time, the much higher linear refractive index contrast of the large ribs will ensure the capability of guiding light by them too. However, relatively speaking, the wider segment can be considered almost linear as we will show shortly. Similar binary arrays can, conceptually at least, be thought in silica (of considerably lower nonlinearity than AlGaAs): A planar arrangement of two interlaced silica waveguides of relatively large and small diameters embedded in a silica cladding of slightly lower refractive index.
For a system such as the one illustrated in Fig. 1 and for optical propagation inside the first band of the band structure, Eq. (1) can readily be replaced by a pair of discrete nonlinear Schrödinger equations (DNLSEs) via a tight-binding approximation performed on the binary structure involved. By choosing even (odd) node-indices for the set of the relatively small (large) segments one readily obtains, where C is the coupling coefficient (in m -1 ), the magnitude of A n (related to the power of the beam) is in W 1/2 , G even and G odd are the respective nonlinear coefficients in W -1 m -1 and V even where, the nonlinear term in the second DNLSE can be ignored since G odd /G even ≈1.8% for the aforementioned example. The latter ratio is usually less than 1% to 2%, depending upon the geometrical details of the optical arrays in hand such as the ones based on ribbed AlGaAs configurations. Equations of the type of Eq. (3) are the subject matter of this work. This system of coupled discrete equations corresponds to a two-component superlattice that has also been previously proposed. 24 : It provides a generalization of the DNLS theory with a wider applicability than the tight-binding approximation and its importance has been explicitly stressed. It represents, therefore, a suitable discretization scheme for the continuous NLS model in media with transversely modulated linear and nonlinear. However, its apparent complexity (compared with previous single DNLS models) was considered in the past as a drawback for further pursuing it.
In the present work it is shown that, due to the effectively linear character of the second DNLSE (practically, a discrete linear Schrödinger equation), the system is actually amenable to analytical investigation of its behavior to a considerable extent. Thus, useful conclusions could be obtained as far as properties, functionality and controllability of devices of the type discussed are concerned. It must be emphasized that the investigation presented is not exhaustive since the problem in hand has a multitude of aspects which worth further investigation. However, the findings and the conclusions of the present study provide a firm analytical basis for this endeavor.
The paper is organized as follows: The model and the corresponding nonlinear diffraction relation along with the comparison of the transmission spectrum of propagating continuous wave modes with the single DNLS are presented in Section II. In Section III, the equationcondition for the modulational instability (MI) of these modes is derived and the latter are investigated at the base, the edge and inside the Brillouin zone as far as their stability and their characteristic features are concerned. In Section IV, the propagation dynamics of the continuous wave solutions of zero transverse momentum are analytically investigated and discussed. In Section V, numerically found solitary solutions are presented. Finally, in Section VI the main conclusions are summarized.

The model
In the following we consider Eq.(3) within a more general perspective, namely we allow for defocusing nonlinearity as well: where σ=±1. The phase shift parameters ε even and ε odd , from the point of view of optical configurations as discussed in the Introduction are respectively the normalized detuning parameters V even /2C and V odd /2C. In a broader perspective, however, they may vary in a wide range of values. The nonlinear spatial dispersion (diffraction) condition for the existence of propagating modes can be found by assigning to ψ n the form ψ n =v n exp(ik z ζ) This system renders the following expression for the even sites: for continuous (plain) wave solutions (CW) for the even-site lattice of the form v 2m =Aexp (imq), where the transverse (Bloch) momentum q is q=2k x D with D being the common spatial period of the lattice. Then, the odd-site lattice also supports CW solution of the form v 2m±1 = v 2m (1+1/e ±iq )/[2(κ+Δε)], where Δε=ε even -ε odd is the normalized relative detuning (for optical configurations as such as in Fig. 1, usually negative) and κ = k z +1-ε even is related to the propagation constant. In order for such periodic solutions to exist |A| 2 must be positive. Therefore, the following condition must hold: Note that for a single nonlinear lattice of spatial period 2D the respective condition is fundamentally different, In Fig. 2 various examples for both cases and for both signs of σ and Δε are presented. At the Diffraction condition for the existence of propagating modes DD(shaded areas) for self-focusing (first row) and self-de focusing (second row) nonlinearity, for interlaced lattices (a, b, d, e) and single nonlinear (c, f). base (q=0) and the edge (q=π) of the Brillouin zone the group velocity, dκ/dq, vanishes. The respective Δε-dependent pairs of κ's (characterizing diatomic lattices) are Δε/2± (Δε 2 /4+1) 1/2 and Δε/2± | Δε/2|. In contrast, the single nonlinear lattice of spatial period 2D the respective single κ's are 1 and -1. More importantly, there is a fundamental difference between the diatomic lattice and the single nonlinear lattice: There exist two transmission bands, namely: max(-Δε,0)0≤κ≤(Δε 2 /4+1) 1/2 -Δε/2 and -(Δε 2 /4+1) 1/2 -Δε/2≤κ≤min(-Δε,0). For Δε= 0 these transmission bands are reduced to the familiar one of the DNLS, -1≤κ≤1.
In the shades areas in Fig. 2, either inside the transmission bands or in the band-gaps, the aforementioned nonlinear diatomic CW solutions can be supported. Furthermore, solitary solutions outside the transmission bands are also expected as we are going to show in Section V. Equation (6) can also be exploited for the investigation of various families of stationary solutions of transverse even-wise periodicity 1, 2, 3 and so on. This is an extensively investigated area in the scientific literature since Eq. (6) is the stationary version of a DNLS for the even sites. However, the propagation dynamics of the interlaced system, Eq.(4), is qualitatively and quantitatively different than a system posed by a single z-depended DNLS.

Modulational instability of the CW solution
Perturbing the amplitude and the phase of a CW solution v 2m =Aexp (-imq where u and w are functionally independent small amplitude perturbations, Γ the propagation constant and Q is their respective spatial frequency. One can easily obtain the following expression for the respective linear response by solving the second linear equation in Eq.(4), Utilizing Eqs. (9) and (10) in Eq.(4) and linearizing, one readily obtains, The MI-condition equation can be easily obtained from Eq.(11), On the basis of Eq. (12), a general remark one can make is that the stability character is unique for fixed values of σκ and σΔε: That is, |Im[Γ(κ, Δε)]| for focusing media (σ=1) coincides with |Im[Γ(-κ, -Δε)]| for the defocusing ones (σ=-1). Therefore, without loss of generality, one may fix the signs of σ (say, σ=1) and Δε and investigate the associated growth rates of the MI over the entire κ-region (-∞<κ<∞). In the following, the investigation concerns the upper and the lower edges of the Brillouin zone as well as an indicative case within the that zone for the sake of stressing morphological differences among the edges and the interior. In Fig. 3, the growth rate of MI, Im(Γ), is shown for the upper edge of the Brillouin zone (q=π) in gray scale diagram as function of the effective propagation constant of the original modulationally unstable CW, κ, and the spatial frequency of the modulation Q. It is worth mentioning that, at this edge, the instability region in κ is bounded as shown in Fig. 3. The solid gray curve traces this region which consists of a main lobe extending over the entire range of Q and of two secondary lobes in the vicinity of Q=0, π. In Fig. 4, the growth rates for two, manifestly unstable choices of κ, that is, κ=0.5, 2 (respectively within the secondary and the main lobes) are shown [ Fig. 4(a,c)] along with their respective propagation constants, Re(Γ) [ Fig. 4(b,d)]. In the latter, the unstable modes are shown in gray. The unstable modes within the secondary lobes [ Fig. 4(a,b)] are co-moving with the CW, that is, they are characterized Re(Γ)=0 for MI. In Fig. 5a, the propagation of a perturbed mode corresponding to Fig. 4(c,d) (with κ=2, that is, within the instability region) is shown. In contrast, in Fig. 5b, the propagation of a perturbed mode with κ=4 (that is, within the stability region) is also presented. The initial (ζ=0) amplitudes in both cases satisfy the condition Eq. (7). The perturbations on the amplitude and phase ζ=0 are random within ±5%. The maximum propagation distance (ζ=50) corresponds to about z max = 2 cm for an AlGaAs type of configuration as in Fig. 1  At the lower edge of the Brillouin zone (q=0), MI exhibits quite different qualitative features: The instability occurs within a semi-bounded region in κ. For self-focusing media CW modes exist for -Δε/2-[1+(Δε/2) 2 ] 1/2 ≤κ <1 and -Δε/2+[1+(Δε/2) 2 ] 1/2 ≤κ<∞. The latter semi-bounded region exhibits MI, while the former one is MI-stable. For self-defocusing media, on the other hand, the bounded and the semi-bounded regions become respectively, 1<κ≤-Δε/2+[1+(Δε/2) 2 ] 1/2 which is MI-stable, and -∞<κ≤-Δε/2-[1+(Δε/2) 2 ] 1/ 2 where MI is expected to develop. In Fig. 6 the region of instability is shown for self-focusing media and Δε=-1. In Fig. 7, the growth rates for two, manifestly unstable choices of κ, that is, κ=1.7, 3.5 (respectively near the edge and in the bulk of the instability region) are shown [ Fig. 7(a,c)] along with their respective propagation constants, Re(Γ) [ Fig. 7(b,d)] where the unstable modes are shown in gray.
It is worth mentioning, at this point, that in the case of a single DNLS lattice the steady state solutions at the lower edge of the Brillouin zone (the so-called unstaggered modes) are always MI-unstable (MI-stable) for self-focusing (self-defocusing) media. However, in the  interlaced model in hand, at the lower edge of the Brillouin zone, there always co-exist MIstable (MI-unstable) CW modes for self-focusing (self-defocusing) media within an isolated region in κ-space. In Figs. 8a-8b, the propagation of two perturbed modes corresponding to Figs. 7a-7d.are shown The perturbations on the amplitude and phase at the launching point (ζ=0) are random within ±5% as previously. In Fig. 8c the propagation of a perturbed mode with κ=0.4, that is within the stability region, is also shown. The short scale amplitude variations are intrinsic to the particular choice of the amplitude at ζ=0 and they are going to be elucidated in the next Section. The sustained longer scale amplitude variations are due to the random amplitude and phase perturbations initially imposed.
Within the Brillouin zone MI attains different characteristics. To qualitatively illustrate this point two examples are provided: For self-focusing media, Δε=-1 and q=2π/3 (that is near the upper edge), the instability remains bounded (as in the case q=π) and it is shown in Fig. 9. The solid gray curve traces this region which consists of a main lobe extending over the entire range of Q while the two secondary lobes in the vicinity of Q=0, π, which previously characterized the upper edge of the Brillouin zone (in Fig. 3), now coalesce in a single disjoint region of instability. In Fig. 10, the growth rates for three manifestly unstable choices of κ, that is, κ=0, 1.5, 2.0 (one within the secondary region and two within the main one) are shown (first row) along with their respective propagation constants, Re(Γ) (second row). In the latter, the unstable modes are shown in gray. In Figs. 11a-11b, the propagation of two perturbed modes corresponding to the first two columns of Fig. 10 (κ=0, 1.5) are shown. In contrast, in Fig. 11c the propagation of a perturbed mode with κ=5 (that is, well inside the Fig. 10. Growth rates for, κ=0, 1.5, 2 (one value within the secondary region shown in Fig. 9 and two within the main one) are shown (first row) along with their respective propagation constants (second row). In the latter, the unstable modes are shown in gray.
stability region) is shown. Again, in all these cases, the initial (ζ=0) amplitudes satisfy the condition Eq. (7) Fig. 10; (c) stable propagation of perturbed mode with κ=5 (that is, beyond the instability region shown in Fig. 9).
Δε=-1 is shown (Fig. 12). This case is equivalent (upon κ -> -κ) to a case of a self-focusing medium, as in the previous example, with, Δε=1. In Fig. 13, the growth rates an unstable choice of κ, (κ=1.12) is shown along with the respective propagation constants, Re(Γ). Additionally, in Fig.14, the propagation of MI-stable mode (for a negative value of κ in this particular case) is contrasted with an unstable one (right sub-figure). Since the instability Fig. 12. Gray scale diagram of the MI growth rate within the Brillouin zone (q=2π/3) as function of the effective propagation constant of the original modulationally unstable CW, κ, and the spatial frequency of the modulation Q for self-defocusing nonlinearity (σ=-1) and Δε=-1.
is weak, the propagation distance is set to a very large value in order to allow the instability to develop. The random perturbations on the amplitude and phase at ζ=0 are also set within the quite wide margin of ±25%. The short scale amplitude variations of the stable mode are both due to this quite strong initial perturbations and to the particular choice of the amplitude at ζ=0 as well.

Propagation dynamics of the CW solutions at the lower edge of the Brillouin zone
Starting from the original system, Eq.
It is clear that the combined intensity is a constant of the motion (S 0 ). By introducing a new combined dynamical variable D= ρ 2 even -ρ 2 odd , the previous system becomes, By eliminating ρ odd and ρ even one finally gets, ( ) Notice that, for the case of a single lattice, this system reduces to the following one: The For fixed values of Δε and S 0 (constant of propagation) one or three fixed points exist (i.e., one or three distinct values of D) as shown in Fig. 15. The transition from one to three occurs only for σcos[Δφ(ζ=0)]>0 [that is, for focusing (defocusing) nonlinearity and phase difference 0 (π), as in Fig. 15a where beyond some value of S 0 , D(S 0 ) becomes multi-valued function]. In Fig. 15 contours of constant Δε are shown. The red line corresponds to Δε=0. One may easily show that the value of D at the transition point can then be expressed in terms of S 0 as follows: The dynamical system exhibits Hamiltonian structure: ( ) Therefore, the system can be integrated in terms of elliptic functions, where the sign in front of the integral is determined by Δφ(ζ=0). Only periodic or bounded Jacobi elliptic functions are acceptable. Therefore, proper conditions on the choice of Δφ(ζ=0), ρ even (ζ=0), ρ odd (ζ=0) must be met. However there is plenty of flexibility available because one has three parameters to adjust. The stability of the fixed points involved for a particular choice of the parameters Δε and S 0 is controlled via the sign of the parameter Δ (Δ>0: stable, Δ<0: unstable) at the fixed points, In the Figs. 16a-16f, the parameter regions the nature and the number of fixed points are depicted. For zero phase difference between the odd and the even sites (Δφ=0, first column) and for focusing nonlinearity there exist both stable (elliptic) and unstable (hyperbolic) fixed points while for defocusing nonlinearity there exist only one elliptic fixed point. On the other hand, for phase difference Δφ=π (second column) for focusing nonlinearity there exist only one elliptic point, while for defocusing nonlinearity there exist both stable (elliptic) and unstable (hyperbolic) fixed points. In the vicinity of the elliptic points stable and periodic in ζ CW solutions exist possessing not necessarily the same amplitude in the odd and even sites. This periodic in ζ behavior is visible, for example, in Fig. 8c: Due to the random perturbation imposed at ζ=0, the CW modes exhibit short scale stable variations in ζ since the system has been slightly displaced from its stable fixed point configuration.
In Fig. 17  separates the latter family for the innermost ζ-periodic modes. It should be emphasized that the separatrix corresponds to a family of unstable modes with infinite period in ζ: If they were stable, they could preserve indefinitely the initially chosen Δφ and D or, equivalently, the same amplitudes in even and odd lattice sides and the same difference in their respective phases. One final important remark: There are no amplitude-only-stationary CW solutions with q=0, that is, solutions of the form as in Eq. (13) with constant amplitude and periodic phase variation. That is, the single lattice respective modes, as in Eq. (17) cannot be extended to the interlaced lattice model. For CW solutions with nonzero transverse momentum (q≠0) a multitude of possibilities opens up among which there exist modes of transverse period 1, 2, etc in each sub-lattice. However, as far as the particular underlying dynamical system is concerned several degrees of freedom may now be involved that render the system non-integrable in the general case. This, of course, goes beyond the scope of the present work.

Solitons
We have shown in Section II that the band structure of the interlaced lattice system is different than the one for the single nonlinear lattice case: There is a third bounded band-gap between the two semi-infinite band-gaps, that is, located between the acoustic and the optical branches. Because of this extra gap, the comparison with the DNLS of the single lattice becomes rather tedious and, therefore, in this Section we will limit ourselves in a few basic remarks. Characteristic cases of solitary modes with propagation constants inside these three band gaps will be presented. All solitary solutions have been found via the Newton's iteration method In Fig. 18 two cases with propagation constant located in the upper semi-infinite band gap (that is, above the acoustic branch) are shown: The case in Fig. 18(a,c) corresponds to an inphase fundamental solution, while the one in Fig. 18(b,d) represents a bound state with dipole-like initial profile. Both cases refer to a self-focusing nonlinearity with κ=0. 8 and Δε=1.15. The transmission bands for this choice of Δε, are -1.73<κ<-1.15and 0<κ<0.58). In both cases, the initial (at ζ=0) profile is perturbed randomly in amplitude and phase within 10%.
In the lower semi-infinite band gap (that is, below the optical branch) one expects localized solutions only for self-defocusing nonlinearity (σ=-1). Indeed, in Fig. 19 two such cases with propagation constant, κ=-3, located in the lower semi-infinite band gap (for the same choice of Δε as in Fig. 18) are shown. The case in Fig. 19(a,c) corresponds to an inphase fundamental solution, while the one in Fig. 19(b,d) represents a bound state with dipole-like staggered initial profile. In both cases, the initial (at ζ=0) profile is perturbed randomly in amplitude and phase within 5%.  58<κ<0 and 1.15<κ<1.73). The propagation constant is κ=0.46. The case in Fig.  20(a,c) corresponds to a staggered type of solution (apart of the central sites, the sign flips for even and odd ones separately), while the one in Fig. 20(b,d) represents a bound state with an initial profile which is dipole-like and staggered too. In both cases, the initial (at ζ=0) profile is perturbed randomly in amplitude and phase within 5%. Comparatively speaking (Fig. 20 vs.  Fig. 19, though for different sign in the nonlinearity) the numerical investigation have shown that the extra gap is populated by more structurally rich localized modes, that is, modes with several sites appreciably active

Summary and conclusions
In this Section the main observations and new features that the interlaced lattice model introduces are summarized. The model is basically a diatomic one and it is characterized by two coupled DNLSEs. However, due to the effectively linear character of the second DNLSE, it was possible to analytically investigate the behavior of the system to a considerable extent. However, the present investigation was not exhaustive since the problem in hand has a multitude of aspects which worth further investigation. However, the findings and the conclusions of this study provide a firm analytical basis for further investigation.
An important design parameter for the lattice which can be suitably chosen is the normalized relative detuning Δε which dictates the existence of two distinct CW modes at each constituent lattice. As far as the modulational instability of the CW modes is concerned, several observations can be made: The stability character is unique for fixed values of the σκ and σΔε (σ is the sign of the cubic nonlinearity). Therefore, without loss of generality, one may fix the signs of the cubic nonlinearity and the sign of the design parameter Δε and investigate the associated growth rates of the MI over the entire κ-region (-∞<κ<∞). At the upper edge of the Brillouin zone, namely for Bloch momentum q=π, the instability region in κ is bounded. This region consists of a main lobe extending over the entire range of the values of modulation frequency, Q, and of two secondary lobes in the vicinity of Q=0, π. The unstable modes within the secondary lobes are co-moving with the CW, that is, their respective real part of the exponent Γ for MI is zero. On the other hand, at the lower edge of the Brillouin zone (q=0), MI exhibits quite different qualitative features: For self-focusing media, the instability occurs only within the semi-bounded region in κ. while the bounded one is MI-stable. The same applies for selfdefocusing media where MI is expected to develop within the respective semi-bounded region. It is worth mentioning that, in the case of a single DNLS lattice, the steady state solutions at the lower edge of the Brillouin zone (the so-called unstaggered modes) are always MI-unstable (MI-stable) for self-focusing (self-defocusing) media. In contrast, in the interlaced model in hand, at the lower edge of the Brillouin zone, there always co-exist MIstable (MI-unstable) CW modes for self-focusing (self-defocusing) media within an isolated region in κ-space. Within the Brillouin zone MI is modified accordingly: For instance, for self-focusing media and near the upper edge, the instability remains bounded in the main lobe, while the two aforementioned secondary lobes coalesce in a single disjoint region of instability. However, near the same edge, for self defocusing media, MI spreads now inside the bounded region as well.
The propagation dynamics of CW modes at the lower edge of the Brillouin zone has also been investigated since, at this edge, the behavior is amenable to a rather straightforward analytical treatment. Additionally, from the practical point of view, these modes are fundamental: That is, they represent CW modes launched transversely onto the face of the lattice sides. It is worth mentioning that there are no amplitude-only-stationary CW solutions in this case, that is, solutions with constant amplitude and periodic phase variation. The corresponding dynamical system has two evolving in z dynamical variables, namely, the common phase difference among the sites of the constituent lattices and the difference of their respective squared amplitudes (i.e., their intensities). The system exhibits Hamiltonian structure, it possesses two constants of the motion (propagation invariants), namely the Hamiltonian function itself [Eq.(19)] and the sum of the squares of the respective amplitudes. The equations of motion are readily integrated and, therefore the dynamical variables are expressed in terms of elliptic functions. The stability of these CW modes is tightly associated with the existence and the character of the fixed point of this underlying dynamical system. The character and the number of these points is controlled by the design parameter Δε, which dictates the MI as we have already shown, and the common sum of the intensities of adjacent sites (as illustrated in Fig. 16). It is worth mentioning that, for focusing (defocusing) nonlinearity and phase difference 0 (π) [that is, for in-phase (staggered-type)] among adjacent sites, there exist three fixed points beyond a particular value of this common sum (which is conserved upon propagation). In such a case, two of the fixed points are elliptic and one is hyperbolic. In the vicinity of the elliptic points there exist stable and periodic in z CW modes possessing not necessarily the same amplitude in the odd and even sites. On the other hand, through the hyperbolic point passes a separatrix. The latter separates three families of periodic in z modes: Periodic CWs which encircle all three fixed points and CWs in the vicinity of each elliptic point. The former one represents modes with respective amplitudes and phases in adjacent sites varying in a considerably wider range.
The existence and the form of stable localized solutions were also investigated. However, the presence of an extra gap (controlled by the design parameter Δε), in a region which is basically a transmission region for the single DNLS lattice, renders the comparison with the latter quite tedious. Therefore, the investigation was limited to a few characteristic cases in all three band-gaps. In the upper semi-infinite band-gap, that is, beyond the acoustic branch, inphase discrete solitons were found for self-focusing materials. In this band-gap, there exist a variety of both fundamental (center-peaked) localized solutions as well as dipole ones. Inside the extra gap, there also exist both previous types. They are, however of a staggered sort of type. None in phase stable localized modes were found in this gap. Finally, in the semiinfinite band-gap below the optical branch, localized modes were found for self-defocusing media, as expected. They are, too, of staggered sort of type, fundamental or dipole ones. In conclusion, by choosing the propagation constant κ, and accordingly adjusting the design parameter Δε, a multitude of localized modes can, thus, be generated. There is evidence that the extra gap is populated by more rich structurally (several sites are active) localized modes.
The appearance and the controllability, by design, of the extra gap in the interlaced lattice model presented in this work, lead to a richness in both the form and the dynamics of the stable modes which are supported. These attributes, the interlaced lattice system possesses, may decisively contribute to an enhanced functionality of devices of the type discussed. .