we present a numerical study of exact breather dynamics in two - dimensional nonlinear schrdinger lattices .
this issue fits well with the enduring scientific interests of serge aubry in the many faceted subject of localization and transport in nonlinear macroscopic discrete systems , where serge s outstanding contributions are widely recognized .
many ideas and lines of study followed in the investigation reported here have found a source of inspiration in early works of aubry and collaborators , among which the doctoral works of th .
cretegny and j.l .
marn deserve special mention .
the existence and properties of localized solutions in extended discrete systems have atracted interest in a broad range of physical fields @xcite .
discrete breathers , sometimes referred to as intrinsic localized modes , are spatially localized and time periodic solutions .
these solutions arise in the context of nonlinear discrete systems and are of fundamental interest for varied physical applications such as pulse propagation in nonlinear optics , energy storing and transport in biomolecules , plasma physics , etc ... the existence of discrete breathers in these systems has been proved rigorously @xcite for a number of equations with physical relevance and , contrary to continuous nonlinear equations , their existence can be regarded as a generic feature of these systems .
one of the most important class of equations are the so - called _ discrete nonlinear schrdinger lattices _
the existence of discrete breathers has been proven for a wide range of systems belonging to this class of nonlinear difference - differential equations .
in particular , the most important example of wide aplicability is the standard nonlinear schrdinger equation .
for instance , this equation was employed in @xcite for describing the propagation of localized beams in an array of nonlinear ( kerr type ) waveguides , having experimental validation subsequently reported in @xcite .
the study of two - dimensional nonlinear schrdinger lattices has atracted much atention @xcite in recent years due to the new phenomena emerging when the dimensionality of the lattice is increased .
some examples of these new features are the existence of vortex - breathers @xcite which supports energy flux , the appearance of an energy threshold for the creation of discrete breathers @xcite and the ubiquity of an instability ( the quasi - collapse ) of some discrete breather solutions leading to a highly localized pulson state @xcite .
these theoretical efforts have their counterpart in recent advances in the field of nonlinear optics .
the studies of two - dimensional arrays of coupled nonlinear waveguides allow the experimental observation of those effects studied theoretically .
specially relevant is the recent experimental breakthrough ( theoretically designed in @xcite ) by fleisher _ et al _ @xcite , where a two - dimensional array of nonlinear waveguides is induced in a photosensitive material .
this technique provides a clear experimental verification of the two - dimensional discrete breather existence in this system .
in particular , besides the observation of standard discrete breathers , these works reported the first observations of staggered discrete breathers .
our study here will focus on the computation of numerically exact dicrete breathers in two - dimensional anisotropic nonlinear schrdinger lattices , _
i.e. _ where the couplings in the two spatial directions are different .
the use of shooting methods allow us to find these solutions and analyze their structural and stability properties . both pinned and mobile discrete breathers are studied . in the latter case
we will study only the ones whose motion is along one axis of the lattice .
the analysis of the numerically exact solutions shed light on some features of the properties and stability of localized solutions reported in previous works .
the plan of the paper is as follows : in section [ sec : numerical ] we provide the technical background and details needed for self - contained purposes .
first we summarize in [ subsec:1d ] the main conclusions on the dynamics of 1d schrdinger discrete breathers reported in @xcite
. a detailed account of the numerical methods ( svd - regularized newton continuation of operator fixed points ) that we have used can be found in that reference . also in [ subsec:1d ]
, we discuss briefly the most relevant formal differences with respect to alternative approaches to ( one - dimensional ) exact mobility of discrete breathers , _
e.g. _ those in refs @xcite and @xcite . in subsection
[ sec : numerics ] we introduce the two - dimensional anisotropic salerno lattice and provide explanations on the implementation of the numerical procedures used to study the dynamics of 2d discrete breathers .
the analysis of the results of our numerics on pinned discrete breathers for anisotropic nonlinear schrdinger lattices is reported in section [ sec : pinned ] .
we present the numerical computations of the fixed point norm , as a function of three parameters : breather frequency , transversal coupling , and nonlinearity ( see below ) .
they show , as anticipated , the so - called _ quasi - collapse _
transition , associated with the ( well - known ) existence of thresholds for the breather norm in two - dimensional lattices .
we present numerically computed sectors of the bifurcation surface .
we conclude this section with a brief look at the nonlinear dynamics on the unstable manifold , whose typical trajectories have been called _ pulson _ states
. early numerical work on the 2d quasi - collapse phenomena in isotropic lattices was reported in @xcite and @xcite .
a two - year - old account of the `` state of knowledge '' on 2d schrdinger lattices can be found in section six of @xcite .    in section
[ secc : mobile ] we show results on a type of mobile breather , namely those moving along the direction of stronger lattice coupling constant .
the structure of each of these mobile exact discrete breathers is that of a localized moving _ core _ superimpossed on a specific extended state of resonant small amplitude radiation , the _ background_. we present here the results of an extensive floquet stability analysis of this type of solutions in two sectors of the three - dimensional parameter space , which clearly show the existence of two different transitions .
the tangent space eigenvectors associated to each of the transitions are presented , and the relation of the unstable manifold trajectories to pulson states is analyzed afterwards .
we conclude with section [ secc : conclusions ] , where we briefly review the results obtained and illustrate their possible implications for mobility of 2d discrete breathers .
the use of numerical tools for the continuation of discrete breather solutions has been widely employed since their existence proof was reported ( see _ e.g _ @xcite ) . in particular , the design of numerical techniques for finding exact mobile breathers based on those employed for the pinned ones has been explored in the recent years @xcite and paves the way to resolving the still open question about discrete breather mobility . here , after reviewing the most important features of 1d mobile breathers in nonlinear schrdinger lattices , we briefly explain how the continuation method is implemented in our system .      as was introduced in section [ sec : introduction ] ,
exact mobility of dicrete breathers in 1d schrdinger nonlinear lattices were numerically studied by the authors in previous works @xcite .
in particular , a numerical continuation from the integrable _ ablowitz - ladik lattice _ , a - l , ( where exact mobile breathers can be calculated analytically @xcite ) , @xmath1 to the standard _ discrete nonlinear schrdinger equation _ , dnls , @xmath2 was performed within the so - called _ salerno model _
@xcite , @xmath3 where @xmath4 , @xmath5 and @xmath6 are the parameters accounting for the strength of the nonlinear terms .
the above equation includes the former two relevant equations ( [ eq : al ] ) and ( [ eq : dnls ] ) , for ( @xmath7 , @xmath8 ) and ( @xmath9 , @xmath10 ) respectively , providing the desired interpolation needed to develop the continuation scheme .
the integrability of equation ( [ eq : al ] ) provides for @xmath11 a two - parametric family of discrete breathers @xmath12\;{\mbox{exp}}[{\mbox{i}}(\alpha(n - x(t))-\omega(t))]\;,\ ] ] where the two parameters @xmath13 $ ] and @xmath14 describe the velocity and internal frequency of the solution @xmath15 the salerno lattice ( [ eq : salerno ] ) , possess two dynamical invariants , namely the hamiltonian , @xmath16 where @xmath17 denotes the complex conjugate of @xmath18 , and the norm @xmath19    in order to find exact mobile discrete breathers in the salerno model we define a ( @xmath20)-resonant solution @xmath21 referred to some time scale @xmath22 such that @xmath23 within the above definition , a mobile breather that translates @xmath24 sites after @xmath25 periods of the internal oscillation will satisfy equation ( [ eq : p - q_resonant ] ) when @xmath26 . in this sense ,
the continuation focuses on the families of ( @xmath20)-resonant discrete breathers , that is breather solutions with the two characteristic time scales ( corresponding to the breather velocity , @xmath27 , and the internal frequency , @xmath28 ) being commensurate . in all the computations we have used finite lattices with periodic boundary conditions ( pbc ) so that @xmath29 and @xmath30 ( with @xmath31 being the lattice size ) .    in the previous works @xcite the authors start from those a - l solutions which are ( @xmath20)-resonant , we discretize ( fine grid ) the path @xmath32 ( with @xmath5 and @xmath6 being positive ) along the salerno model , and go through it computing the corresponding ( @xmath20)-resonant solutions for the pairs @xmath33 .
we can choose the path without loss of generality because of the scaling property of equation ( [ eq : salerno ] ) .
then , each solution is numerically computed as a fixed point of the map @xmath34 where @xmath35 is the lattice translation operator @xmath36 , and @xmath37 is the @xmath38-evolution map ( @xmath39 ) following the dynamics dictated by equation ( [ eq : salerno ] ) for the corresponding value of @xmath6 ( @xmath40 ) ; _ i.e. _ , @xmath41=\{\phi_{n}(t_{0}+t_b)\}$ ] .
the continuation was then performed for a fine grid of frequencies belonging to the family of ( @xmath42)-resonant discrete breathers .
the most important conclusion about discrete breathers in these nonlinear schrdinger lattices is that mobility in the non - integrable regime ( @xmath43 ) demands the presence of an extended _ background _ to which the fixed point solution is spatially asymptotic ( @xmath44 ) , _ i.e. _ the solution is exactly written as @xmath45 this expression defines the purely localized component @xmath46 of the solution .
the background is a finite linear combination of nonlinear plane waves , @xmath47 $ ] .
these plane waves are exact solutions of the salerno model ( [ eq : salerno ] ) being the nonlinear dispersion relation @xmath48\cos k_{j } - 2\nu |a_{j}|^{2}$ ] .
the results concerning the characterization of the background are discussed in detail in @xcite , here we briefly summarize the most important features :    * the set of `` @xmath49 '' plane waves which take part in the background of a @xmath50resonant discrete breather with internal frequency @xmath51 is derived by the simple selection rule for the wave - numbers @xmath52 @xmath53 _ i.e. _ only the plane waves which are ( @xmath20)-resonant with the internal period of the breather can be components of @xmath54 . the number of solutions of ( [ resonant ] ) fixes `` @xmath49 '' . *
the amplitudes @xmath55 of the nonlinear plane waves differ by orders of magnitude .
* there exist a strong positive correlation between the amplitude of the background and the strength of the peierls - nabarro barrier arising from the periodic lattice .
this correlation is particularly clear when symmetry breaking transitions occur for the also studied case of @xmath56 and @xmath57 ( see reference @xcite ) , and reflects the link between non - integrability and the existence of the background dressing of the mobile core . *
finally , the interpretation of the correlation described in _ ( iii ) _ is reinforced from a study of the energy evolution of the mobile core : _ there is an energy balance brought by the background when the core moves along the lattice_. in particular , it can be observed how the core energy oscillates periodically so that it takes the maximum energy value when the core visits the inter - site configuration .
this extra energy periodically obtained by the core is provided by the background , with the energy maximum clearly related to the background amplitude .    before concluding this subsection
, it is worth commenting on some of the differences between the newton continuation of fixed points that we use in this paper , and other important recent approaches to breather numerics .
the work by ablowitz _
et al _ @xcite uses discrete fourier analysis to obtain a nonlinear nonlocal integral equation , from where the `` ... soliton is thus viewed as a fixed point of a nonlinear functional '' ( sic ) in the fourier transformed space of functions .
following these authors , their results seem to differ from those of early pioneering work @xcite ( nowadays textbook material @xcite ) `` in which a _ continuous _ travelling solitary waves were reported using fourier series expansions with finite period @xmath58 while assuming convergence as @xmath59 '' ( sic ) .
et al _ term _ continuous _ a solution that can be defined off the lattice points , which they see as `` necessary when discussing travelling waves in lattices '' ( sic ) , and disagree with some conclusions reported in the earlier works .
the ( `` orthodoxy matters '' ) discussion above helps us to clarify how differently our numerical approaches `` sees '' the discrete schrdinger breather problem : the very concept of a variable defined off the lattice points is intrinsically alien to our discrete approach , which neither needs of it nor excludes its eventual consideration . in contrast to those views ( but not at all in logical opposition ) , we consistently view the thermodynamical limit ( @xmath60 ) in lattice space , much in the sense used _
e.g. _ by serge aubry in his celebrated work on the frenkel - kontorova ground state problem @xcite : the infinite size limit is built up from a subsequence of pbc ( finite ) lattices for which the limit is well defined .
this will make the fourier - transformed @xmath61-space continuum .
closer to our approach in some respects , though technically different in many others , is the formal approach purposed recently by james and collaborators @xcite .
this approach , which uses recent central manifold theorems , was brought to our attention very recently and unfortunately , we have to defer to a future publication a comparison of our 1d computations with theirs , this is beyond the scope of this paper , namely , the anisotropic 2d schrodinger lattices , which we introduce in the next subsection .
motivated by the method described above , we can extend the continuation scheme for calculating exact discrete breathers in higher dimensional systems .
in particular we focus on the two - dimensional nonlinear schrdinger lattice @xmath62    this lattice can be viewed as the _ two - dimensional salerno model_. the two coupling parameters @xmath63 and @xmath64 provide a technical advantage for numerics ( see below ) , but they are also introduced for theoretical and experimental interest .
the possibility of controlling the ratio between the two linear couplings of the two transversal directions has been studied in various works as a way of analysing how the intrinsic 2d phenomena ( such as the quasi - collapse ) emerge .
in fact , for @xmath65 , @xmath9 and @xmath10 equation ( [ eq:2d - dnls ] ) describes a set of weakly coupled nonlinear wavequide arrays and can be considered as a case of `` _ intermediate dimensionality _ '' .
this extreme has been studied experimentally in @xcite and using perturbative methods in @xcite .
on the other hand , this equation incorporates , as two particular limits , the physically relevant standard two - dimensional dnls equation ( @xmath9 , @xmath10 ) and the two - dimensional a - l lattice ( @xmath7 , @xmath8 ) .
the continuation between these two limits provides a useful tool for studying the interplay between the on - site and inter - site nonlinearity in the 2d anisotropic lattice .
moreover , the anisotropy ( or freedom in the values of the coupling parameters @xmath66 and @xmath67 ) allows to include an integrable model among the members of the family of nonlinear lattices described by ( [ eq:2d - dnls ] ) . that is , for @xmath8 , @xmath68 and
@xmath69 one obtains a set of integrable a - l chains . in this way , every 2d model included in ( [ eq:2d - dnls ] ) is connected with an integrable model where analytic discrete breathers are available . in what follows
we set @xmath70 and @xmath71 with @xmath72 and @xmath57 .
the dynamics ( [ eq:2d - dnls ] ) can be derived from the salerno poisson structure @xmath73 with the hamiltonian @xmath74 as in the 1d salerno model there is also a second conserved quantity , the norm , due to the phase invariance of the equations of motion ( [ eq:2d - dnls ] ) @xmath75    in the same manner as in the 1d case we will focus on a special set of 2d discrete breathers . for this , we have to generalize the definition ( [ eq : p - q_resonant ] ) of a resonant solution in the two - dimensional case . in this context discrete breathers solutions
are characterized by three time scales .
namely , one associated with the internal oscillation @xmath28 and the other two derived from the translation of the localization center , _
i.e. _ its velocity @xmath76 .
the subset of @xmath77-tuples ( @xmath28,@xmath78 ) that fulfil the ( @xmath79)-resonance condition @xmath80 ( where @xmath81 , @xmath82 and @xmath25 are integers ) denote the breather solution that can be obtained with our continuation method .
these solutions are those that after @xmath25 periods of the internal frequency , @xmath83 , translates @xmath81 and @xmath82 lattice sites in the @xmath84 and @xmath85 direction of the square lattice , respectively , _
i.e. _ @xmath86 where again pbc are applied @xmath87 , @xmath88 , @xmath89 and @xmath90 ( with @xmath91 and @xmath92 being the lattice size in the @xmath84 and @xmath85 direction respectively ) . consequently , a ( @xmath79)-resonant state will be a solution of the following set of equations @xmath93 & = & { l}_{y}^{p_{y}}{l}_{x}^{p_{x } } { t}_{(\omega_{b } , \nu , c_{1 } , c_{2})}^q[\{\hat{\phi}_{n , m}(t_{0})\}]= \nonumber \\ & = & \{\hat{\phi}_{n , m}(t)\}\ ; , \label{eq:2d - map}\end{aligned}\ ] ] where the operators @xmath94 are the lattice translation in the @xmath95-direction , @xmath96&=&\{\phi_{n+1,m}(t_{0})\}\ ; , \\
l_{y}[\{\phi_{n , m}(t_{0})\}]&=&\{\phi_{n , m+1}(t_{0})\}\;.\end{aligned}\ ] ] besides , @xmath97 is the time evolution operator given by equation ( [ eq:2d - dnls ] ) over one period @xmath98 , @xmath99= \{\phi_{n , m}(t_{0}+t_{b})\}\;.\ ] ]    in order to illustrate the 2d time scales resonance we consider the plane wave solutions of equation ( [ eq:2d - dnls ] ) : @xmath100 $ ] .
these solutions possess the following nonlinear dispersion relation @xmath101 hence , we can obtain the subset of plane waves which are ( @xmath79)-resonant with some time scale @xmath22 ( _ i.e. _ after a time @xmath102 they have translated @xmath81 and @xmath82 sites in the @xmath84 and @xmath85 direction , respectively ) .
each member of these subsets will be labelled by the pair @xmath103(@xmath104 , @xmath105 ) and from the condition ( [ eq:2d - map ] ) it follows that the corresponding set of values of @xmath106 for each family will satisfy the relation @xmath107 where @xmath108 is an integer and @xmath109(@xmath81 , @xmath82 ) . in figure [ fig : resonances ] the corresponding values of @xmath106 are represented for two resonances of the type @xmath110 and @xmath111 .
cc    the method used for solving equation ( [ eq:2d - map ] ) for each resonant @xmath77-tuple ( @xmath28,@xmath78 ) is the same as in the 1d case , which is extensively described in @xcite .
the implicit function theorem assures that a fixed point solution of a map ( [ eq:2d - map ] ) given by @xmath112(@xmath81 , @xmath82 , @xmath25 , @xmath28 , @xmath6 , @xmath63 , @xmath64 ) can be obtained provided that _ ( i ) _ the jacobian of the operator @xmath113-{i}$ ] is invertible and _ ( ii ) _ we know a fixed point of a map corresponding to an infinitesimally close set of parameters , @xmath114 ( @xmath81 , @xmath82 , @xmath25 , @xmath115 , @xmath116 , @xmath117 , @xmath118 ) .
the first demand can be satisfied using a _ singular value decomposition _ ( svd ) @xcite of the jacobian in order to obtain the pseudo - inverse operator . on the other hand , when the second condition is fulfilled convergence of the newton - raphson iterative scheme is guaranteed . for this , we start with a sufficiently good trial solution , @xmath119 and solve the equation @xmath120^{-1}\cdot f_{\vec{\xi}}[\{\phi^{0}_{n , m}(t_{0})\}]\ ; , \label{eq : iteration}\ ] ] in order to obtain @xmath121 .
we iterate this calculations to the desired convergence , and then the solution @xmath122 , is obtained . in our numerics
this is the case when @xmath123<n\cdot10^{-16}\;,\ ] ] ( where @xmath31 is the total number of sites in the square lattice ) is fulfilled .
once the solution is found we use it as the following trial solution , @xmath119 , for solving the map ( [ eq:2d - map ] ) corresponding to the next set of parameters @xmath124 .    as an additional benefit ,
this method provides the linear stability of the computed solution . for this
, we only have to compute the eigenvalues of the jacobian of the operator around the solution @xmath125 $ ] .
in fact , this jacobian is the extended floquet matrix of the solution .
then , a solution is linearly stable if all the eigenvalues of the corresponding floquet matrix are inside the unit circle .
moreover , the symplectic character of the dynamics ( [ eq:2d - dnls ] ) implies that all the eigenvalues appear in quadruplets ( @xmath126 , @xmath127 , @xmath128 , @xmath129 ) and thus for a linearly stable solutions all the eigenvalues of its extended floquet matrix lie on the unit circle .
there are two possible paths for developing the continuation method depending on the choice of the starting point of the continuation .
one possibility is to start from the full anti - continuum limit , @xmath130 , where a pinned breather solution of frequency @xmath28 is written as @xmath131 starting from the above solution , we can perform the continuation increasing the parameters @xmath63 and @xmath64 as usual , and so obtain the whole family of ( @xmath132)-resonant discrete breathers .
an alternative path starts from the one - dimensional limit , @xmath133 .
the choice of this second limit ( which implies taking as the very initial trial solution of the continuation the whole set of 1d solutions described in section [ subsec:1d ] ) is justified when seeking mobile solutions .
as stated above , this limit offers the possibility of studying strongly anisotropic lattices as a controlled interpolating situation between one and two dimensions . on the other hand , employing this strategy we can only obtain those solutions which are ( @xmath134)-resonant , _ i.e. _ the two - dimensional continuation of those one - dimensional ( @xmath135)-resonant discrete breathers .
hence , the solution from which we start is @xmath136 where @xmath137 is the corresponding ( @xmath135)-resonant one - dimensional solution .
in what follows we will employ both continuation paths when we study the case of pinned breathers ( section [ sec : pinned ] ) , and we will show that the results obtained are the same when approaching the same limit ( the standard two - dimensional dnls ) .
as we have discussed , we can choose two different starting points for the continuation of @xmath138-resonant fixed points ( pinned breathers ) of equation ( [ eq:2d - map ] ) : _ ( i ) _ the full anti - continuum ( ac ) limit ( @xmath139 ) , or _ ( ii ) _ the ( one - dimensional , 1d ) limit of uncoupled chains ( @xmath140 ) , where they can be obtained from continuation along increasing values of the parameter @xmath6 from the 1d a - l lattice ( [ eq : al ] ) . as a test for our codes
, we have checked that both paths arrive to the same solution .
in fact , unique continuations can proceed along any path on the plane of parameters ( @xmath141 ) that we have explored .
early works @xcite on the isotropic two - dimensional standard dnls equation analyzed the so - called _ quasi - collapse _ instability of pinned discrete breathers , _
i.e. _ the condensation of all the energy into a few modes in discrete nonlinear systems , which corresponds to the onset of a singularity ( wave collapse ) @xcite in multidimensional continuum models .
subsequent numerical works @xcite extended these studies to the isotropic 2d salerno lattice and addressed the question of how the instability is affected by the presence of impurity lattice sites .
as expected , our results further corroborate the existence of quasi - collapse instabilities in the anisotropic case : the phase diagram in parameter space ( @xmath51 , @xmath67 , @xmath6 ) consists of two regions ( stable and unstable ) separated by the surface of transition . as we perform the continuation of breather solutions across the parameter space we scan the floquet stability of the computed solution . in figure
[ fig : ctresholds ] we present the two stability transition curves in the plane ( @xmath28 , @xmath64 , @xmath142 ) , _ i.e. _ the function @xmath143 , corresponding to the two different continuation starts .
the continuation from the ac limit is made through the path @xmath144 and the one from the 1d limit is made at @xmath145 . the convergence of the two paths at @xmath146 is clearly seen .
cc    the criterion for stability of the pinned discrete breather solution derived in @xcite , @xmath147 is of a very general character and our numerics illustrate it clearly . on the other hand ,
the floquet stability analysis detects the dimensionality ( and a basis in tangent space ) of the unstable linear manifold associated with the quasi - collapse instability that these exact discrete breathers experience for some parameter values .
we have computed numerically , for a fine grid of @xmath51 values and a coarser grid of @xmath67 and @xmath6 , the function @xmath148 , from which we show some sectors in figures [ fig : norma - c2 ] and [ fig : norma - nu ] .    in figure
[ fig : norma - c2 ] we show the numerically computed norm ( [ eq:2dnorm ] ) as a function of the breather frequency @xmath149 , for three different values of the transversal coupling @xmath67 , and a fixed value of @xmath150 ( anisotropic dnls limit ) .
we observe the existence of a minimum value , @xmath151 , which is thus seen as an _ excitation threshold _ for the creation of these solutions .
the position of the minimum @xmath152 , which _ naturally _ increases with @xmath67 , separates the stable and unstable branches of pinned breathers .
breathers corresponding to values of @xmath51 where @xmath149 has a negative slope are unstable : this is shown in the insets , where the floquet spectra of two representative examples of pinned discrete schrdinger breathers are plotted in the complex plane .
note that the high accuracy of the numerical solution allows an unprecedented detailed floquet analysis of the instability , paving the way to rigorous analytical characterizations of the quasi - collapse unstable manifold .
this is a one - dimensional manifold , as our numerical results unambiguously confirm .
then , in the regime of small time scales , the unstable manifold is fully characterized by a single floquet eigenvector .
cc    figure [ fig : norma - nu ] shows the ( surface ) function @xmath153 for the volume sector of constant @xmath67(@xmath154 ) .
most noticeably , the critical ( threshold ) line of bifurcation points ( @xmath155 ) , as seen in the inset , does not define a monotone function @xmath156 .
in fact , in the whole interval of @xmath157 values , the range of values of @xmath158 is quite short , indicating the insensitibity of the gross features of the quasi - collapse transition to the value of @xmath6 .
however , considering finer details , one sees that the threshold curve @xmath156 smoothly reaches its slightly larger values around midway between the dnls and the a - l limits . in other words , intermediate values of the interpolating ( salerno )
parameter @xmath6 somewhat favour the enhancement of the quasi - collapse unstable region .
these conclusions are in contrast with the stated conclusion ( for isotropic lattices ) in @xcite that the ablowitz - ladik term increases the stability regime .
cc    when instability is allowed to develop beyond the fixed point tangent space into the nonlinear realm of perturbations , the trajectory obtained by direct integration of the equations of motion invariably ends after a transient ( of time scale given by the real floquet exponent larger than @xmath159 ) in a localized solution with complex dynamics , the _
pulson states_. these states were characterized in @xcite in the following terms `` ... where the peak intensity @xmath160 oscillates between the central site and its four nearest neighbours ...
it is not known whether these pulson states represent true quasiperiodic solutions to the dnls equation '' .
what makes these trajectories on the unstable nonlinear quasi - collapse manifold of much practical relevance and interest is their ubiquity : they appear as persistent localized states in the hamiltonian dynamical evolution from a wide variety of initial conditions .
their description requires at least two frequencies , namely the internal ( genuine breather - like frequency ) and the frequency of the oscillations of the _ breather width _ around a mean width value , which turns out to be less than the width of the unstable exact discrete breather .
second and outer shells of neighbours ( in both lattice axes ) also participate in the width oscillations .
cc    though a more detailed characterization of the pulson states will be presented elsewhere , it is illustrative to consider ( figure [ fig : qc - fft ] ) the power spectrum of the field at the central site of a typical trajectory on the unstable nonlinear manifold of a quasi - collapsing pinned discrete breather .
this shows peaks at the combinations @xmath161 ( @xmath162 ) , where @xmath163 is the frequency of the width oscillations characterizing the pulson state , while @xmath164 is a frequency higher than the ( initial condition ) fixed point frequency @xmath51 .
the new frequency @xmath165 turns out to be very close to the breather frequency of the same ( initial ) norm on the stable branch . in other words
, the instability drives a shift of breathing frequency towards the stable branch , while the excess energy is transferred to the oscillatory motion of the observable width .
this behavior seems to be the essence of the physical characterization of the nonlinear quasi - collapse manifold dynamics .    the numerical observation of a two - frequency power spectrum for a typical pulson state points towards an eventual positive answer to the question ( on true quasiperiodicity ) arised in @xcite , though , for more details we have to refer to the phd dissertation @xcite .
this point serves to illustrate how the high accuracy of the fixed point numerical solution provides detailed clues on many still unsolved ( from a mathematical and physical point of view ) questions on two - dimensional schrdinger localization , which are of prospective experimental interest in nonlinear ( photonic , josephson , ... ) physics technologies .
early and current attempts to explore straightaway discrete breather mobility in isotropic 2d schrdinger lattices seem to agree @xcite that `` kicking '' procedures meet huge difficulties in delivering good mobile solutions , contrary to the numerical experiences in 1d lattices .
we note here that the formal basis for those methods @xcite takes advantage of the floquet spectra analysis of exact pinned breathers , where the so - called _ depinning _ ( symmetry - breaking ) mode is identified .
this allows , _ provided peierls - nabarro barriers are small enough _ , to obtain nice numerical 1d mobile discrete breathers , by computing trajectories from perturbations of the exact pinned breather along the tangent space direction specified by the depinning eigenvector .
the presence of symmetry - breaking instabilities leading to exchange of stability between one - site and two - site centered pinned breathers @xcite and the associated lowering of the peierls - nabarro barriers to breather displacements , hugely facilitating the success of these procedures when applied to ( both hamiltonian and dissipative ) one - dimensional lattices .
@xcite    in contrast , our `` anisotropic lattice '' continuation approach takes advantage of the availability of exact 1d mobile solutions by monitoring the parameter @xmath67 of transversal coupling , and then does not rely on how easily one promotes clean mobility from pinned localization .
in this way we obtain accurate numerical ( @xmath166 ) fixed points , that is schrdinger discrete breathers moving along the strong coupling direction . in a forthcoming paper
we will address the ( much more difficult ) question for arbitrary direction of motion .      in figure
[ fig : mob - structure ] we visualize the instantaneous real and imaginary components of the 2d discrete field profile of a typical ( 1,0,1 ) schrdinger breather .
its structure can be seen as the natural extension to two - dimensional lattices of the structure of mobile schrdinger breathers analyzed in @xcite .
the numerical solution is spatially asymptotic to a _ finely tuned _ small - amplitude extended ( delocalized ) radiation state @xmath167 when @xmath168 .
the fixed point solution can be thus decomposed as    @xmath169    which defines @xmath170 , the spatially _ localized _ component of the solution .
it turns out that the spatially _ delocalized _ component is a highly localized state in the ( continuum , in the thermodynamic limit ) @xmath61-space of wavevectors .
more precisely , @xmath171 is a finite linear combination of ( 1,0,1 ) -resonant nonlinear ( _ i.e. _ amplitude - dependent frequency @xmath172 ) 2d planewaves .
it can be said that , as might be expected , 1d schrdinger breather mobility smoothly persists when ( strong @xmath66-coupling ) 1d chains are coupled transversally .
importantly , the numerical continuation for increasing values of the transversal coupling @xmath67 proceeds far from the weak coupling regime into where the genuine two - dimensional effects start to be manifest , as we will see below .
cc    most noticeable , the svd - regularized newton procedure invariably selects the values @xmath173 for all values of @xmath67 and @xmath6 , and thus the values of @xmath174 for the 2d resonant planewave are independent of @xmath67 ( so it remains equal to the @xmath61 values of the 1d ( @xmath175 ) fixed point for the uncoupled chain ) .
the appearance of an extended background modulation in the transversal direction of @xmath176 appears naturally as the _ best _ choice to take advantage of approximately 1d breather propagation along strong coupling direction , for it keeps the value of @xmath174 favoured by the strong coupling @xmath66 value : any other value of @xmath177 would entail a different @xmath174 value .
note however that this provides only a _
plausibility argument _ for the interpretation of the numerical observation ( @xmath176 ) .
the high accuracy of the computed solutions allows a detailed analysis of many issues concerning 2d schrdinger breather exact mobility along the strong coupling direction .
we leave aside in this paper many of them , and focus here on how the existence of quasi - collapse instabilities of pinned schrdinger breathers , for increasing @xmath67-coupling values , influences the stability properties of moving ( @xmath178 ) breathers . in other words , we search here for genuine 2d effects on these `` strong - coupling -direction '' ( quasi-1d ) moving breathers .    we have performed an exhaustive exploration of two sectors of the parameter space ( @xmath179 ) , corresponding to the breather frequency values @xmath180 , and @xmath181 , by computing the continued ( @xmath178 ) fixed point .
these values of @xmath51 were chosen low enough to allow the analysis of pinned breather quasi - collapse effects on mobility , which occurs at relatively low values of @xmath182 for these values of @xmath172 .
cc     +    the floquet analysis of the computed solutions provides the stability diagrams represented in figures [ fig : mob - stability ] .
both show no qualitative differences : there are two regions in the ( @xmath183 ) plane where the ( @xmath178 ) mobile breather is linearly unstable .
the figures are not `` schematic '' : every point of the plane in a fine grid of values of @xmath67 and @xmath6 has been analyzed , _
i.e. _ the floquet spectrum of the computed ( @xmath178 ) fixed point is scrutinized , as shown in figures [ fig : mob - stability].b and [ fig : mob - stability].d , where the modulus of the floquet eigenvalues is shown as a function of either @xmath6 ( figure [ fig : mob - stability].b ) or @xmath67 ( figure [ fig : mob - stability].d ) .
the first unstable region appears at low values of @xmath67 and intermediate to high values of the salerno parameter @xmath6 , _
i.e. _ it does not occur close to the a - l limit .
this unstable region is also bounded above in the direction of @xmath67 : the variation of the modulus of the unstable floquet eigenvalue versus the transversal coupling parameter @xmath67 shows that the mobile breather becomes stable again at larger values of @xmath67 , before the second instability at even higher coupling takes place .
an important observation is that the pinned discrete breather of the same frequency is _ linearly stable _ at the points in this unstable region for ( 1,0,1 ) mobile breathers .
thus this instability can not be ascribed to pinned quasi - collapse effects .
the second transition occurs for values of @xmath67 close to , but slightly higher than , the values @xmath184 of the quasi - collapse of the pinned breather of the same frequency .
we had already seen in the previous section that the quasi - collapse transition @xmath182 is only very weakly dependent on the value of @xmath6 , and note that the same is true for this mobile breather bifurcation .
these results suggest that this second transition is related to quasi - collapsing phenomena .
significantly , the stability of the ( @xmath178 ) mobile breather persists for a small interval of coupling values above the pinned breather quasi - collapse .
this should be regarded as natural , for the mobile breather is a different solution .
note in figure [ fig : mob - stability].b that the modulus of the unstable floquet eigenvalue , in the interior of the unstable region , reaches much higher values than those typical for the first type of instability , and decreases for larger values of @xmath67 , before the breather solution ceases to exist and only plane wave solutions are obtained by our numerical method .
note that this behaviour of the unstable floquet eigenvalue also fits well to the main features of the pinned quasi - collapse instability strength , as described by the slope @xmath185 . from now on we will refer to this instability of mobile breathers as the quasi - collapse instability .    in the next subsection
we characterize both generic types of instability , by looking at the details of the unstable manifold associated with each type .
as we will see , pulson states turn out to play a role in the description of typical trajectories on the unstable nonlinear manifolds .
first , we analyze the quasi - collapse instability of ( @xmath178 ) mobile breathers .
the unstable linear subspace in the tangent space of the fixed point is one - dimensional .
the typical instantaneous profile of the ( modulus ) unstable floquet eigenvector driving the instability is shown in figure [ fig : unstabledyn].d .
it is an exponentially localized 2d profile which decays asymptotically to zero as @xmath186 , _
i.e. _ it does not excite radiation .
these characteristics are shared by the quasi - collapse unstable eigenvector of the pinned breathers , which further reinforce the previous considerations leading us to consider this instability as the mobile counterpart of the pinned quasi - collapse transition .    in figure
[ fig : unstabledyn].c we have visualized the time evolution of the field modulus contour plot for a typical trajectory on the unstable manifold . this is obtained by direct numerical integration of the equations of motion , from an initial condition in which a small perturbation along the quasi - collapse eigenvector has been added to the unstable fixed point solution .
one sees that the breather translational motion slows down , and the energy is transferred to width oscillations .
these oscillations turn out to be more irregular , see figure [ fig : qc - evolution ] , than those observed in section [ sec : pinned ] when we inspected typical trajectories on the unstable nonlinear manifold of pinned breathers
.    1    cc     +    the difference in the character of the width oscillations in both ( pinned and mobile ) cases may be ascribed to the presence of an extended background component in the mobile breather solution , which naturally enters into the energy transfer taking place during temporal evolution .
the slowing down of the translational motion continues and eventually the breather pins into a convulsive pulson state surrounded by the remaining radiation .
now we pay attention to the `` low @xmath67 '' instability of ( @xmath178 ) mobile breathers .
the modulus profile of the unstable floquet eigenvector that drives this instability is m - shaped ( bimodal ) , as shown in figure [ fig : unstabledyn].b , and is asymptotic to an extended planewave - like profile as @xmath168 , _ i.e. _ it is not a purely localized perturbation .
it is indeed rather different from the quasi - collapse unstable eigenvector analyzed above , which is consistent with the fact that the pinned breather of the same frequency is linearly stable in this region of parameter space .
as argued above , this instability is not related to quasi - collapse phenomena , and it does not appear in the region of small values of the salerno parameter @xmath6 , close to the a - l limit .
cc    a typical trajectory on the unstable manifold associated with this instability is shown in figure [ fig : unstabledyn].a , where we have plotted the time evolution of the field modulus contour plot .
we can see there that the mobile breather pins quickly while a small pulse moving backwards is ejected , which spreads and finally mixes with the remaining delocalized background .
however some energy is transferred to width oscillations of the pinned breather so that also in this case we observe the formation of pulson states surrounded by the remaining radiation .
as the main difference of this behaviour , with respect to the evolution observed on the quasi - collapse unstable manifold , is the ejection of the small moving pulse , we refer to this instability as _
fission_.    by increasing the strength of the initial perturbation along the direction of the unstable eigenvector , one observes that the size of the ejected pulse increases .
this observation is consistent with the results reported in @xcite , where the evolution of initial moving gaussian pulses in isotropic 2d schrdinger lattices was studied .
these numerical experiences lead the authors to conclude that `` the characteristic feature of the discrete quasi - collapse of a moving pulse is the splitting of the initially moving broad pulse into a track of the standing narrow structures ... '' ( sic ) .
however , we see from our study of the stability of exact moving discrete breathers that the fission and the quasi - collapse instabilities have different origins and they appear in different regions of parameter space . on the other hand , the ubiquitous phenomenon of width oscillations of pinned localized structures ( pulson states ) can not be ascribed to quasi - collapse .
they also appear as the preferred way to allocate excess of ( localization ) energy in regions of parameter space far from the quasi - collapse unstable region .
we have studied here the dynamics of exact numerical discrete breathers , both pinned and mobile , in two - dimensional anisotropic nonlinear schrdinger lattices .
these solutions are computed using a svd - regularized newton method by continuation from a set of uncoupled 1d chains into increasing non - zero values of the coupling in the transversal direction .
we have performed an extensive exploration in the parameter space ( @xmath28 , @xmath64 , @xmath6 ) of breather frequency , transversal coupling and salerno parameter , by computing the floquet spectra of the numerical solutions .
we have also computed the breather norm @xmath187 and further corroborate the general validity of the criterion found in @xcite , namely that the partial derivative @xmath188 is positive for stable pinned breathers .
furthermore , we have analyzed the dynamics on the quasi - collapse unstable manifold , where the unstable breather experiences a shift in frequency towards the ( higher ) value of the stable breather with the same norm .
the excess of energy is coherently transferred to oscillations of the breather width , so that the resulting pulson state is characterized by two frequencies .
we have studied discrete breathers moving along the strong coupling direction .
these solutions are composed of an exponentially localized core on top of an extended background which is itself the finite sum of a finite set of nonlinear 2d plane waves .
the time scales asociated with these plane waves are resonant with the core internal frequency as happens in the 1d case . in particular
, the background chooses a finite set of plane waves from a continuous family of resonant solutions .
the floquet analysis of these mobile discrete breathers reveals the existence of two distinct types of instability .
one is the counterpart , for mobile breathers , of the quasi - collapse experienced by pinned breathers .
the other instability occurs in a region of parameter space where pinned breathers are linearly stable .
the analysis of the dynamics on the unstable manifold show that the excess of energy is partly transferred to a small moving pulse , ejected from the center of localization , which justifies the designation of a fission instability .
however , part of the energy excess is also transferred to width oscillations .
the appearance of pulson states far from the quasi - collapse regime indicates that the tendency to allocate energy in the form of width oscillations is a general 2d feature , not exclusively associated to quasi - collapse instabilities .    in a future work
we will focus on mobility of 2d discrete breathers in an arbitrary lattice direction .
the results obtained here shed light about how this mobility can be obtained . in fact , our experiences show that mobility of pinned breathers can be induced based on the existence of the extended background in the numerically exact mobile solution . on the other hand , the results obtained here and the aforementioned future work may help to design and better understand recent numerical experiments reported in @xcite , concerning the interaction between high amplitude pinned breathers and mobile ones .
these experiments provides a possible way for routing and blocking mobile discrete breathers via the interaction with the high amplitude pinned ones , resulting in a plausible implementation of logical functions .
the authors acknowledge f. falo , yu .
kivshar , r.s . mackay and m. peyrard for sharing thoughts , and pointing out some important references to us .
jg - g and lmf are grateful to m. johansson and b. malomed for discussions on some issues regarding the `` travelling wave '' ( orthodox ) perspective on discrete breathers .
financial support came from mcyt ( project no .
bfm2002 00113 ) , dga and bifi .
jg - g akcnowledges financial support from the mecyd through a fpu grant .
work at los alamos performed under the auspices of the us doe .
j. c. eilbeck and m. johansson , in proceedings of _ localization and energy transfer in nonlinear systems _ , edited by l. vzquez , r. s. mackay , and m. p. zorzano ( world scientific , singapore , 2003),p.44 .
t. peschel and f. lederer , proceedings of the confererence on nonlinear coherent structures in physics and biology , d.b .
duncan and j.c .
eilbeck eds . , published on www : http://www.ma.hw.ac.uk/solitons/procs/