Compacton matter waves in binary Bose gases under strong nonlinear management
Abstract
The existence of compacton matter waves in binary mixtures of quasi onedimensional BoseEinstein condensates in deep optical lattices and in the presence of nonlinearity management, is first demonstrated. For this, we derive an averaged vector discrete nonlinear Schrödinger equation (DNLSE) and show that compacton solutions of different types can exist as stable excitations. Stability properties are studied by linear analysis and by direct numerical integrations of the DNLSE system and their dependence on the inter and intraspecies scattering lengths, investigated. We show that under proper management conditions, compactons can be very robust excitations that can emerge spontaneously from generic initial conditions. A possible experimental setting for compacton observation is also discussed.
I Introduction
It has been recently demonstrated that BoseEinstein condensates (BEC) in deep optical lattices (OL) when exposed to strong and rapid periodic time modulations of the scattering lengths can support matter wave compactons, i.e. localized excitations with a compact support (1). Similarly to discrete breathers, compactons are intrinsically localized stable excitations. In contrast with them, however, compactons have no exponential tails, the lacking being due to the effective nonlinear dispersion induced by the modulation, which permits the vanishing of the tunneling just at the compacton edges. In particular, it was shown (see (1) for details) that after averaging out the fast time scale dynamics, the tunneling rate of the resulting averaged system depends not only on modulation parameters but also on local field density differences (atom numbers in the BEC case) between neighboring sites. Thus, the periodic modulation in time of the scattering length (nonlinear management) can be used to change the original dispersion (e.g. the one in absence of management) into an effective nonlinear dispersion essential for the compacton existence.
Field dependent tunneling suppression in the presence of strongly and rapidly modulated interactions was demonstrated in (2) for the case of a BEC trapped in a double well potential (compacton formation being severely restricted in this case by the size of the system) and in (1) for a onedimensional BEC array modeled by the discrete nonlinear Schrödinger equation (DNLSE). The phenomenon of compacton formation, however, is of general validity and, as we demonstrate in this paper, can occur also in more complicated BEC systems.
From an experimental point of view, periodic time modulations of the scattering length can be achieved by the Feshbach resonance (FR) technique (3), e.g. by varying the external magnetic field near a resonant value. Besides compacton excitations, periodic modulations of the scattering lengths can be used to create density dependent gauge fields (4), this being presently a field of rapidly growing interest, connected with interesting physical phenomena, including pair superfluidity, exactlydefect free Mottinsulator states etc. (5). Properties of the superfluid Mott transition in a 2D square and a 3D cubic optical lattice with periodic modulation of the atomic scattering length have been investigated in (6). Modulations of the interactions were used also to design new correlatedhopping models for fermions in optical lattices (7) and, in combination with OL shaking, for engineering unconventional BoseHubbard models (8).
Spatial and temporal periodic changes of the scattering lengths have been shown to be effective tools to change stability properties of nonlinear excitations, leading to the existence of two dimensional bright solitons in one and two component attractive condensates (9); (10); (11); (12); (13). Moreover, they were used to induce long lived Bloch oscillations (14); (15); (16), dynamical localization (17), Rabi oscillations (18) of BEC gap solitons in optical lattices, Faraday waves (19); (20) etc. Nonlinear management techniques were considered also in nonlinear optics to stabilize 2D and 3D solitons and to reduce collapse in optically layered media with selffocusing interaction (9), to improve communication capacities via soliton dispersion management in optical fibers (21), to create linear superpositions of gap solitons in periodic Kerr media (22), etc.
All these studies refer to the single component case, e.g. BEC made by a a single atomic species. An interesting question to ask, however, is wether compacton excitations could exist also in multicomponent systems of interest both for BEC and nonlinear optics. The aim of the present paper is just to provide an answer to this question.
In this respect we introduce an averaged vector DNLSE with effective nonlinear discrete dispersion terms which describes the dynamics of a BEC array in the presence of strong interspecies and intraspecies scattering length modulations. The existence of different compacton states in the form of brightbright (BB), brightdark (BD), and darkdark (DD) pairs, is first demonstrated and stability properties of these states investigated both by means of a linear stability analysis and by direct numerical integrations of the model equations. As a result we find that while single site and two sites out of phase BB compactons are stable in the whole parameter range, for the other modes there exists thresholds in the tunneling constant rate below which they cannot exist as stable excitations. The dependence of the compacton stability on interspecies interaction is also investigated. For DD compactons we find that when the interspecies scattering length is detuned to zero (e.g. in the uncoupled single component limit) the stability becomes narrower compared to that of the pure two components case. The predictions of our analysis are shown to be in good agreement with the results obtained direct numerical simulations of the model. The emergence of compact excitations from generic initial excitations and possible experimental settings for their observation are also discussed.
The paper is organized as follows. In Sec. II we introduce the model equations and discuss the averaged equations. In Section III we derive the conditions for existence of two component of BB, BD and DD compactons and discuss stability properties. In section IV we consider the case of DD compactons in the uncoupled limit corresponding to the single component dark compacton case. In section V the emergence of compact excitations from generic initial excitations is investigated and possible experimental parameter design provided. Finally, in Sec. VI the results of the paper are briefly summarized.
Ii Model equations and averaging
Twocomponent atomic BEC in a deep OL can be described in the tight binding approximation by the following vector DNLSE (23)
(1)  
where the overdot stands for time derivative, the coefficients are related to the tunneling rates of atoms between neighboring wells of the optical lattice and , are nonlinear coefficients related to the interspecies () and intraspecies () scattering lengths, respectively. Eq. (1) arises also in nonlinear optics where they model the propagation of the electric field in array of optical waveguides with variable Kerr nonlinearity. In this context the role of time is played by the longitudinal propagation distance along the optical fiber and the nonlinear coefficients correspond to self and crossphase modulations of the electric field components, respectively (25); (24). Notice that the above two component DNLSE has the Hamiltonian form with and the Hamiltonian given by
(2) 
Here the stands for the complex conjugation and denotes the complex conjugate of the expression in the parenthesis. Also notice that the number of atoms are conserved for each component.
Taking into account that Eq. (1) arises in the tight binding approximation (26); (27), the nonlinear coefficients can be expressed in terms of the overlap integrals of Wannier functions of the corresponding continuous periodic GrossPitaevskii mean field model (28) as:
with denoting the atomic masses and the reduced mass (functions are normalized on the whole line: ). The modulations of the nonlinear coefficients is assumed to be of the form
(3) 
with a constant, a rapidly varying periodic function of time and a small parameter that controls the strength and the frequency of the modulation (strong nonlinearity management corresponds to ). For we consider periodic modulations of the type
(4) 
of period in the fast time variable . Although in general one could consider independent modulations of the nonlinear coefficients, in practical contexts one usually deals with the simpler settings

i) , i=1,2

ii) , i=1,2.
In the BEC context, case i corresponds to a modulation in time of the intraspecies scattering lengths, keeping interspecies scattering length constant, while in case ii it is done just the opposite. Although these settings are both feasible for BEC mixtures, case ii) is not easy to implement in nonlinear optics because in this context crossphase modulations are usually difficult to control. Moreover, the averaged equations obtained in the case ii) are mathematically more involved due to the presence of a complicated nonlinear dispersion (see below). In the following we concentrate mainly on case i) and discuss case ii) only briefly at the end of the section.
Let us consider then a fixed interspecies scattering length () and assume the intraspecies scattering lengths modulated as follows:
(5) 
To find effective nonlinear evolution equations, we use averaging method to eliminate the fast time, , dependence. In this respect it is convenient to perform the following transformation
(6) 
where are the antiderivatives of ,
where is the period. By substituting Eq.(6) into Eqs.(1) we obtain:
(7) 
(8) 
where , and
(9) 
The average over the rapid modulation can be done with the help of the relations
(10) 
where , are Bessel functions (29) of first kind of zeroth and first order, respectively, and with given by
(11) 
Here the angular bracket denotes average with respect to the fast time variable, e.g. . The system of averaged equations is then obtained as:
(12) 
(13) 
Note that Eqs. (12, 13) have Hamiltonian form with averaged Hamiltonian, , given by:
(14) 
By comparing Eq. (14) with the corresponding unperturbed Hamiltonian Eq. (2), one can see that the effect of the scattering lengths modulation simply reflects in the following nonlinear rescaling of the tunneling constants:
(15) 
From this equation and from Eq. (11) it is clear that the interspecies interaction play no role in determining the lattice sites where the zero tunneling condition occur (compacton boundaries) when the nonlinear management is made with respect to the intraspecies interactions (see below for explicit examples). Thus, in analogy with the single component case considered in (1) the tunneling constants depend on the atom difference between neighboring sites (2); (1). This introduces an effective nonlinear dispersion in the system which leads to the existence of twocomponent compactons.
In closing this section it is interesting to discuss the changes in the above derivation implied by a nonlinear management of the interspecies scattering length
(16) 
with constant (not necessarily equal) intraspecies parameters , . In this case, to remove the explicit time dependence from Eq. (1), the transformation in (6) must be replaced by
(17) 
with the antiderivative of . Following the same approach as before, it is not difficult to show that one arrives at the same averaged Hamiltonian as in Eq. (14) but with the interchanges , , operated in the tunneling terms and the obvious replacements , . This implies a different rescaling of the tunneling constants
(18) 
and leads to complicated nonlinear dispersion terms in the averaged equations. In view of this complexity, in the following we restrict only to compactons induced by intraspecies management. A detailed study of the interspecies management requires more investigations and will be discussed elsewhere (30).
Iii Existence and stability of compactons in binary BEC mixtures
Exact compacton solutions of the averaged system can be searched as stationary states of the form:
(19) 
with chemical potentials of the two atomic species. Substituting these expressions into Eq. (12,13) one gets the following stationary equations:
(20)  
(21)  
to be solved for the chemical potentials and amplitudes of the compacton modes. The compact nature of the solution ( outside a finite (small) range of sites), allows to truncate the above infinite system into a finite number of relations between the above variables, which can be solved exactly. In the following subsections this is shown explicitly for the different compacton types.
iii.1 BrightBright compactons
To search compacton solutions of the BB type we need to look for the last sites of vanishing amplitude, say , where the vanishing of the tunneling rate is realized. For a single site BB compacton we assume for all . Substituting this ansatz in Eqs.(20,21) we obtain the corresponding condition for the compacton existence as
(22) 
where is a zero of the Bessel function (in all numerical calculations below we take the first zero of : ). This condition together with
(23) 
gives us the single site BB compacton pair.
Typical examples of single site BB compactons are depicted in top panels of Figs. 1,2. Stability properties of the solution have been investigated by standard linear analysis. Denoting by the eigenfrequencies of the linearized averaged equations, associated to growing perturbations of the form , we have that linear stability is granted if all have zero imaginary parts. For the considered single site BB compacton, this condition is well satisfied, as one can see from the left panel in the second row of Fig. 1. Actually we find that single site BB compactons are generically stable for wide range of parameters. This can be seen from the second top right panel of Fig. 1 where the dependence of the lowest squared eigenfrequencies, , is reported as a function of (to reduce number of parameters we fix, here and in the following , ) . Extensive numerical linear stability studies show that is always non negative, meaning that the solution is linearly stable.
Nonlinear stability properties have been investigated by direct numerical integrations of Eqs. (1213) taking as initial conditions exact compactons perturbed with a random uniformly distributed noise field of amplitude . Results are found in full agreement with the linear stability results discussed above. Excellent agreement is also obtained from direct numerical integrations of the original vector DNLSE (1), as one can see from third and fourth row panels of Fig. 1. In particular notice from the bottom row panels that the deviation of the original dynamics from the exact averaged dynamics is very small (four orders of magnitude small for ) and can be made even smaller by further decreasing .
In Fig. 2 stability properties of single site BB compactons have been investigated as a function of the interspecies interaction, , for a particular choice of the remaining parameters. As one can see, the obtained behavior is very similar to the one obtained for the dependence in Fig. 1, this further confirming the robustness of the one site BB solution with respect to wide parameters variations.
Twosite BB compactons can also be found and, similarly to discrete breathers, can be of two types: inphase and of outofphase. Twosite inphase compactons follow from the ansatz: . By substituting into Eqs.(20,21) one gets Eq. (22) as before with chemical potentials given by:
First two row panels of Fig. 3 show a typical twosite inphase BB stationary compacton. Notice from the second row panels that the linear eigenfrequencies are all positive in the interval this implying a wide stability rate also in this case. Direct numerical integrations of the averaged equation is shown in last two rows of Fig. 3.
Outphase BB compactons can be obtained from the ansatz . In this case the chemical potentials are given by
with fixed as in Eq. (22).
Typical twosite outofphase compacton and related linear stability analysis and time evolutions, are reported in Fig. 4.
Quite remarkably it is also possible to find exact threesite BB compactons. In this case assuming , we obtain chemical potentials as
where and .
The constraint equations for the amplitudes are Eq.(22) and
Results are depicted in Fig. 5 for typical example. Notice from second row right panel, that stability range in this case is strongly reduced and stability is possible only for . Also notice from the bottom panels that the deviation of the original dynamics relatively higher than the previous cases of BB. The amplitude profile differs mainly at the middle point but this discrepancy, however, stays bounded in time and is quite small if compared to the maximum amplitudes of the exact solution for each component (for the considered case it never exceed the of the exact amplitudes).
We also find that threesite BB solutions with small values display higher discrepancy compared to the respective solution with large (by increasing , however, the amplitudes of the compacton also become smaller). These discrepancies may be ascribed to the management functions being not sufficiently strong. It is worth to remark here that in the limit (infinitely strong management) exact compactons of the averaged system should be exact also for Eq. (1) (limiting cases of very small , however, are quite difficult to investigate numerically due to the very high accuracy required).
iii.2 BrightDark compactons
It is interesting to look for novel type of solutions of the averaged DNLSE, which are more characteristic of the twocomponent systems, and in particular to BD (or darkbright) onsite compactons. Existence of such solutions follows from the general stationary Eqs. (20,21) by letting , and for . In this case one finds that exact solutions exist if chemical potentials and parameters of the BD compacton are related by the following equations
(24) 
(25) 
with a zero of . The last equations can be solved numerically to determine the amplitude and the background of the dark component.
Typical one site BD compactons are depicted in Figs. 6,7 (see top panels) together with their stability property (see left middle panels of these Fig.s). Stability properties versus and parameters are displayed in the second row right panels of Figs. 6,7, respectively, from which we see that in spite of the presence of some instability region, the range of stability is still quite large. Numerical time evolution obtained from Eq. (1) and deviation from the exact stationary dynamics, are investigated for a BD compacton taken in the left and right bottom panels of Fig. 6. The figure refers only to the dark component, being the bright component in very good agreement with the stationary solution. Notice from the right bottom panel the presence of fine periodical ripples on the background and small oscillations of the deep amplitude in the middle. These deviations, however, are relatively small compared to dark amplitude and background of the exact solution, remaining bounded on a long time scale.
iii.3 DarkDark Compactons
Single site darkdark (DD) compactons can be searched of the form and for . The equations for chemical potential and amplitudes in this case are:
(26) 
with .
The first row panels of Fig. 8 show the amplitude profile for and respective numerical linear stability as function of . In the second row panels we report stability properties for two increasing values of the interspecies interaction: (left and right panel, respectively). As one can see, the region of stability becomes wider as is increased, this indicating a tendency of the coupling to stabilize darkdark pairs. Time evolutions with respect to the averaged and original systems and their deviations from the exact solution are shown in the bottom two rows left and right panels, respectively. In general we find that DD system is more sensitive to disturbance compared to all previous cases and for higher value of and (such as and for example) much smaller and numerical time step are needed to numerically simulate stable dynamics on a long time scale (exact initial solutions can be perturbed only by very small noise of the order ). A possible reason for this weak stability is discussed in the next section.
Iv Single component dark compacton
From the above analysis it appears evident that stability properties become more critical for compactons that involve dark components. This is true for BD and even more for DD pairs. This fact may result from a possible instability of the single component dark compacton. We remark that for single component BEC, compactons were investigated only for the bright case (1), so the question of whether stable single component dark compactons can exist, is open. In this section we provide the answer to the question by studying the uncoupling limit of the DD compacton discussed before.
In this respect the averaged equations in (12,13) for with reduce to the same equations considered in (1) for the single field variable . Considering then the equation becomes
(27) 
with . Dark modes are possible only for repulsive interactions so that in the above equation . A dark single site compacton located at the site can be determined by setting if and for to yield
(28) 
where . It follows that to satisfy the equation at the nonvanishing sites, one must have
(29) 
then Eq.(28) can be numerically solved to obtain .
Numerical linear stability analysis for parameters values and (see second row, right panel of Fig. 9) shows that the solutions are unstable until . A comparisons with the first two rows of Fig. 8 indicates that the stability regime can be improved by interspecies coupling in the two component DNLS system. Nevertheless, when , stable dark compactons exist as shown in the third row of Fig. 9 for case .
From this it appears evident that, in contrast with the bright compactons case, the strong nonlinear management does not provide (at least for the simple real amplitudes ansatz assumed) stable single component dark compacton solutions for wide range of parameters. It is interesting that this situation is slightly improved when the interspecies interaction is switched on, as discussed before.
V Generation of binary BEC compactons and experimental setting
It is also interesting to discuss physical conditions for which compactons could be experimentally observed in ultracold BEC mixtures. From a first sight one could think that, in order to keep the zero tunneling condition satisfied at the edges one needs a very precise control of the number of atoms inside the compacton, a fact that could be difficult to arrange in a real experiments. On the other hand, if the excitation is very stable, it should appear also in the presence of generic fluctuations that are unavoidable in any real experiment. In this respect is of interest to address the problem of compacton generation from generic initial excitations, this providing evidence of their robust emergence even under unfavorable conditions. To this end, let us concentrate on the most stable compact excitations, which, as we have seen, are the ones of BB type. Without loss of generality, we take initial excitations of Gaussiantype for both components , with fixed as in the following numerical simulations (similar results can be obtained with other types of initial conditions).
In Fig. 10 we show the generation of a single site (see top two row panels) and a three sites (see bottom two row panels) BB compacton from a Gaussian initial condition (same for both cases). The first and third row panels refer the density profiles (at time and , respectively) while second and forth row panels refer to time evolutions of the two components. From top two row panels we see that after expelling some excess matter, a single site compacton on site (plus some background noise), emerge. The background noise is unavoidable because the number of atoms (squared amplitudes) in the compacton for the chosen parameter is fixed by the first zero of and is smaller than the number of atoms in the initial Gaussian. In spite of this, the densities of the two components in the neighboring sites, , keep very small (although not exactly zero) during their time evolution, irrespectively from the excess radiation. For the considered case, was and therefore the management was not very strong. One can expect that by further decreasing one can better and better approximate a true compacton solution with vanishing densities at (this becoming exactly true in the limit ). The effect of a stronger management can be see from the bottom two row panels of Fig. 10. In this case epsilon is one order of magnitude smaller (), and the mismatch between numbers of atoms in the initial and final states is quite small. As a result, we see that a three site compacton, practically indistinguishable from an exact solution based on the first zero of , is formed. This clearly demonstrates that, with proper parameter design and proper management conditions, compactons can be very robust excitations that can emerge spontaneously from generic initial conditions.
Further results shown in Fig. 11 address the problem of the influence of the interspecies interaction on the compacton formation . Here same parameters and same Gaussian initial excitations as in bottom two row panels in Fig. 10 are used, but with lower value of the interspecies interaction, . In this case we see that instead of a three site compacton, an intrinsic localized mode with exponential tail is formed, in full agreement with the fact that by lowering the interspecies interaction the compacton stability may be lost.
Let us now discuss possible experimental setting for the observation of compacton modes. In this respect we remark that the existence bound pairs of BD solitons (localized excitation with tails) has been experimentally demonstrated in (23) for the continuous case. A possible experimental setting to observe the corresponding compacton modes could be a BEC mixture of K and Rb atoms loaded in a deep optical lattice considered in the experiment (31) and subjected to interspecies scattering length modulations. An equivalent setting could be implemented also in nonlinear optics with an array of optical waveguides with varying Kerr nonlinearity along the propagation (32); (24). For the case of binary BEC mixtures, the time modulations of the scattering lengths can be easily implemented by means of the FR technique by varying the external magnetic field, , near a resonant value:
(30) 
Here denotes the background value of the respective scattering lengths and the resonant value of the magnetic field. In the case of Rb atoms, for example, there are several narrow FR which could be used, the broadest one lying at G. The tightbinding limit, appropriate for a vector DNLSE description, could be reached by considering optical lattices of amplitude , where is the recoil energy. Thus, by changing periodically and rapidly in time the magnetic field around a FR of the intraspecies scattering lengths it should be possible to observe twocomponent matter wave compactons to emerge from generic initial conditions in real experiments.
Vi Conclusions
In conclusion, we have investigated the existence and stability of binary mixtures matter waves in arrays of BEC subjected to time dependent periodic variations of the scattering length, by means of an averaged two component DNLSE. In addition of BB compactons we showed that BD and DD compactons are also possible. The stability of these modes has been investigated both by linear spectral analysis and by direct numerical integrations. We found that the single site and the twosite (outofphase) BB compactons are always very stable in the whole parameter range, while for the other modes there exist thresholds in the tunneling constant rate below which they cannot exist as stable excitations. The stability resulted in general to be more critical for pairs involving one (or both) component of dark type. In all cases, however, the predictions of the averaged system were found in good agreement (in some cases; excellent) with the results of the numerical simulations.
Acknowledgements
M. S. acknowledges partial support from the Ministero dell’Istruzione, dell’Universitá e della Ricerca (MIUR) through a PRIN (Programmi di Ricerca Scientifica di Rilevante Interesse Nazionale) 20102011 initiative.
References
 F. Kh. Abdullaev, P. G. Kevrekidis, and M. Salerno, Phys.Rev.Lett. 105, 113901 (2010).
 J. Gong, L. MoralesMolina, and P. Hänggi, Phys. Rev. Lett. 103, (133002) (2009).
 S. Inouye, M. R. Andrews, J. Stenger, H.J. Miesner, D. M. StamperKurn and W. Ketterle, Nature 392, 151 (1998).
 S. Greschner, G. Sun, D. Poletti, and L. Santos, Phys. Rev. Lett. 113, 215303 (2014).
 A. Rapp, X, Deng, and L. Santos, Phys.Rev.Lett. 109, 203005 (2012).
 T. Wang, X.F. Zhang, F. E. A. dos Santos, S. Eggert, and A. Pelster, Phys. Rev. A 90 013633 (2014).
 M. Di Liberto, C. E. Creffield, G. I. Japaridze, and C. M. Smith, Phys. Rev. A 89, 013624 (2014).
 S. Greschner, L. Santos, and D. Poletti, Phys. Rev. Lett. 113, 183002 (2014).
 H. Saito, M. Ueda, Phys. Rev. Lett. 90, 040403 (2003).
 F. Kh. Abdullaev, J.G. Caputo, R.A. Kraenkel, and B.A. Malomed, Phys. Rev. A 67, 013605 (2003).
 G. D. Montesinos, V. M. PerezGarcía and P. Torres, Physica D 191, 193 (2004).
 P. Kevrekidis, D. Pelinovsky, A. Stefanov, J. Phys. A 39, 479 (2006).
 F. Kh. Abdullaev, J. Garnier, Phys. Rev. E 72, 035603(R) (2005).
 M. Salerno, V.V. Konotop and Y.V. Bludov, Phys. Rev. Lett. 101, 030405 (2008).
 C. Gaul, R.P.A. Lima, E. Díaz, C.A. Müller, and F. DominguezAdame, Phys.Rev.Lett. 102, 255303 (2009).
 E. Díaz, C. Gaul, R. P. A. Lima, F. DominguezAdame, and C. A. Müller, Phys.Rev. A 81, 051607 (2010).
 Yu. Bludov, V.V. Konotop and M. Salerno, Europhys. Lett. 87 (2009) 20004.
 Yu. Bludov, V.V. Konotop and M. Salerno, Phys.Rev.A. 80 023623, (2009)
 K. Staliunas, S. Longhi, and G. J. de Valcarcel, Phys. Rev. Lett. 89, 210406 (2003).
 F.Kh. Abdullaev, M. Ogren, and M.P. Soerensen, Phys. Rev. A 87, 023616 (2013).
 I. Gabitov and S. K. Turytsin, Opt. Lett. 21, 327 (1996); N. J. Smith et al., Electron. Lett. 32, 54 (1996).
 Yu. Bludov, V.V. Konotop and M. Salerno, Europhys.Lett. 93, 30003 (2011).
 U. Shrestha, J. Javanainen, and J. Ruostekoski, Phys.Rev.Lett. 103, 190401 (2009).
 Gaetano Assanto, Luis A. Cisneros, Antonmaria A. Minzoni, Benjamin D. Skuse, Noel F. Smyth, and Annette L. Worthy, Phys.Rev.Lett. 104, 053903 (2010).
 A. Kobyakov, S. Darmanyan, F. Lederer, E. Schmidt, Optical and Quantum Electronics 30, 795 (1998).
 A. Trombettoni and A. Smerzi, Phys. Rev. Lett. 86, 2353 (2001).
 F.Kh. Abdullaev, B. B. Baizakov, S. A. Darmanyan, V. V. Konotop, and M. Salerno, Phys. Rev. A 64, 043606 (2001).
 G.L. Alfimov, P.G. Kevrekidis, V.V. Konotop, and M. Salerno, Phys. Rev E 66 046608 (2002).
 M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, (National Bureau of Standards, Washington, 1964).
 F.Kh. Abdullaev, M.S.A. Hadi, M. Salerno, and B. Umarov, unpublished.
 J. Catani, L. De Sarlo, G. Barontini, F. Minardi, and M. Inguscio, Phys. Rev. A 77, 011603(R) (2008).
 F. Kh. Abdullaev, E. N. Tsoy, B. A. Malomed, and R. A. Kraenkel, Phys. Rev. A 68, 053606 (2003).