the concept of nonlinear self - localization is of importance for many physical phenomena , and has appeared in a number of different contexts since the pioneering work by landau @xcite on the polaron problem in the 1930s . in recent years , much attention has been devoted to studies of spatially localized and time - periodic vibrational modes in anharmonic lattices ( see e.g. @xcite , @xcite for recent reviews ) . the general existence of such modes , which have been termed _ discrete breathers _ , or _
intrinsic localized modes _ , as robust solutions to nonlinear ( and in general non - integrable ) lattice - equations was suggested in 1988 by takeno _
et al _ @xcite .
later , their existence was rigorously proven under rather general conditions by mackay and aubry @xcite by considering the limit of uncoupled oscillators ( the so called _ anticontinuous _ or _ anti - integrable _ limit ) . by means of the implicit function theorem
, they showed that the trivial solution of a single - site localized vibration at the uncoupled limit could be continued into a localized breather solution for non - zero coupling between the oscillators , provided that the individual oscillators are anharmonic , and that no multiples of the breather frequency resonate with the bands of linear excitations ( phonons ) . as was demonstrated first in @xcite
, the ideas of the rigorous proof can be turned into an efficient numerical scheme to calculate breather solutions to any desired accuracy .
since the discrete breathers appear under very general conditions in anharmonic lattices and provide efficient means of energy localization , they have been proposed as candidates to explain experimentally observed localization of energy in many different physical areas , e.g. dna dynamics @xcite .
although , from a fundamental and mathematical viewpoint , the existence theorems for discrete breathers provide an important cornerstone for understanding the dynamics of anharmonic lattices , it is probably of even larger physical importance to understand the behaviour of a system close to an exact breather solution . by linearizing the lattice - equations around the exact solution , one can obtain an approximate description of the dynamics of weakly perturbed breathers , and in particular the linear stability properties determining whether small perturbations will grow exponentially or not .
it was shown in @xcite,@xcite that the simplest , single - site , breathers are generally linearly stable close to the uncoupled limit , and numerical investigations using standard floquet analysis ( see e.g. @xcite ) have shown that linearly stable breathers typically exist also for rather large values of the inter - site coupling . however , when considering time - scales large compared to the breather period , the mere linear stability of a breather does no longer guarantee the eternal existence of the breather in the presence of small perturbations , and there are still many questions remaining concerning the different mechanisms by which breathers may grow or decay , or possibly finally be destroyed . if the breathers have a finite life - time , the determination of this life - time is of large importance for understanding the role of breathers in real systems .
it is the purpose of this paper to investigate in more detail some mechanisms for breather growth and decay in a simple model system , the discrete nonlinear schrdinger ( dnls ) equation .
the dnls equation is generic in the sense that it describes slowly ( in time ) varying modulational waves in discrete systems in a rotating - wave approximation ( see e.g. @xcite ) ; however , due to its extra symmetry properties ( see sec .  [ sec : perturbation ] ) it exhibits some nongeneric features among discrete systems , e.g. exact quasiperiodic breathers @xcite .
the single - site breathers of the dnls equation are stationary states which are linearly stable for all inter - site coupling , and which reduce to the nls soliton in the continuum limit ( see e.g. @xcite ) .
an important application appears in nonlinear optics , where the single - site dnls - breather describes a discrete spatial soliton in an array of weakly coupled waveguides @xcite ; recent experimental observations @xcite confirm the successful use of the dnls - model in this context .
some recent numerical investigations @xcite have shown that dnls - breathers can be spontaneously created from noisy backgrounds , in a similar manner as was previously observed for klein - gordon @xcite and fpu @xcite lattices .
typically , this spontaneous energy localization was observed to occur in two steps . in the first step ,
a large number of small breathers are created as a result of the modulational instability @xcite of travelling plane waves occurring for certain wave number regimes .
the second step proceeds by inelastic collisions between the breathers , in which systematically the big breathers grow at expense of the smaller ones .
thus , the outcome will be a small number of large breathers , together with some remaining background of small - amplitude ( phonon ) oscillations .
however , it can generally not be concluded from the numerical simulations that this is the true final state of the system , and actually long - time simulations for fpu - chains @xcite revealed also a third step , in which the interaction with the phonon oscillations leads to the final destruction of the breather and equipartition of energy .
thus , to elucidate the nature of the final states for typical initial conditions in anharmonic chains , it is necessary to obtain a better understanding of the mechanisms for interactions between breathers and small - amplitude perturbations .    in this paper , we take the following approach . as initial state , we consider an exact single - site breather solution , and add a small perturbation corresponding to a time - periodic eigensolution to the equations of motion linearized around the breather .
these solutions , which can be either localized or extended in space , constitute a complete set in which an arbitrary initial perturbation can be expanded .
the localized solutions correspond to internal modes of the breather @xcite , while the excitation of an extended solution corresponds to a standing - wave ( i.e. , non - propagating ) phonon interacting with the breather . in sec .
[ sec : perturbation ] we describe the model and outline the perturbational approach which forms the analytical backbone for the interpretation of the numerical results presented in secs .
[ sec : internal ] and [ sec : phonons ] .
[ sec : internal ] discusses the long - time consequences of the interaction between the breather and its internal modes , while sec .
[ sec : phonons ] concerns the interaction between the breather and small - amplitude standing - wave phonons of different wave vectors .
we will find that in both cases , a simple argument based on the conservation laws can be used to obtain a sufficient condition for breather growth . in sec .
[ sec : instabilities ] we discuss another type of mechanism for breather growth and destruction which becomes appreciable when the amplitude of the standing wave is non - negligible ( and consequently the perturbational approach can be expected to fail ) , and which has its origin in the recently discovered oscillatory instabilities of the nonlinear standing - wave phonon themselves @xcite . finally , we make some concluding remarks in sec .  [
sec : concluding ] .    concerning the numerical simulations of the dynamics presented in this paper , unless otherwise stated they always apply for a system of infinite size ( finite size systems are considered only in sec .
[ sec : instabilities ] ) .
the simulations have been performed either by using very large system sizes or by appending damping regions of various sizes to the boundaries ; in all cases we have carefully checked that the boundary conditions have no essential influence on our results .
we consider the following form of the dnls hamiltonian with canonical conjugated variables @xmath0 : @xmath1 this yields the dnls equation @xmath2 which , in addition to the hamiltonian ( [ dnlsham ] ) , also conserves the total _ excitation norm _ ( or _ power _ in nonlinear optics applications ) , @xmath3 the conservation laws for the norm and hamiltonian are , through noether s theorem , related to the invariance of the dnls equation ( or , more precisely , of its corresponding action integral ) under infinitesimal transformations in phase ( @xmath4 ) and time ( @xmath5 ) , respectively . defining the norm density @xmath6 and hamiltonian density @xmath7 as in eqs .
( [ norm ] ) resp .
( [ dnlsham ] ) , the conservation laws can be expressed in terms of continuity equations as @xmath8 @xmath9 with the _ ( norm ) current density _
@xmath10\ ] ] and the _ hamiltonian flux density _ @xmath11 , % \ ] ] respectively .
these conservation laws are discrete analogs to those existing for the continuous nls equations with general nonlinearities ( see e.g. @xcite ) , however there is no discrete counterpart to the momentum conservation law since the discrete equation has no continuous translational symmetry in space .
furthermore , we note that the transformation @xmath12 in ( [ dnls ] ) is equivalent to @xmath13 , and thus we will for the rest of this paper only consider @xmath14 without loss of generality .
the single - site dnls - breather is a stationary - state solution to eq .
( [ dnls ] ) of the form @xmath15 where the time - independent shape @xmath16 depends on the frequency @xmath17 and is spatially localized with a single maximum at a lattice site .
the breather exists for all @xmath18 ; the limit @xmath19 ( or @xmath20 ) corresponding to the anticontinuous limit , where @xmath16 is localized at a single lattice - site , while the limit @xmath21 corresponds to the continuous limit , where @xmath16 approaches the nls soliton . the single - site breather is a ground state solution to eq .
( [ dnls ] ) in the sense that it minimizes the hamiltonian ( [ dnlsham ] ) for a fixed value of the norm ( [ norm ] ) , i.e. , @xmath22 , where the frequency @xmath17 appears as the lagrange multiplier ( see e.g. @xcite ) . the norm for the single - site breather , @xmath23 ,
is known to be a monotonously increasing function of @xmath17 , while the hamiltonian , @xmath24 , is negative and monotonously decreasing ( see e.g. @xcite ) . from the minimization condition
, these functions will be related as @xmath25    to describe the dynamics close to the breather ( [ stationary ] ) , we introduce the following perturbation expansion : @xmath26 where @xmath27 is the initial perturbation and @xmath28 . thus , as for the usual stability analysis of stationary states ( see e.g. @xcite ) , the perturbation is applied in a frame rotating with the breather frequency @xmath17 . substituting into eq .
( [ dnls ] ) and identifying coefficients for consecutive powers of the small parameter @xmath29 yields an infinite set of equations , which from 0th to 4th order read : @xmath30 where the operator @xmath31 ( which is linear over the field of real numbers ) is defined from the first equality in eq .
( [ 1st ] ) . the 0th order equation ( [ 0th ] )
gives the breather shape @xmath16 ( which for the single - site breather can be assumed real and positive without loss of generality ) , while the 1st order equation ( [ 1st ] ) is the linearization of the dnls equation around the breather .      to obtain the solutions to the linearized equations ( [ 1st ] ) ,
we proceed in a similar way as is usually done for continuous generalized nls models ( see e.g.@xcite ) and introduce a substitution of the form @xmath32 so that @xmath33 substituting ( [ subst ] ) into eq .  ( [ 1st ] ) and assuming @xmath34 real yields @xmath35 where the operators @xmath36 and @xmath37 are hermitian .
thus , we can obtain the eigenfrequencies @xmath38 and the corresponding eigenvectors @xmath39 from matrix diagonalization , @xmath40 to make connection to earlier work @xcite , we remark that the vector @xmath41 is an eigenvector of the floquet matrix with eigenvalue @xmath42 , where the time - period @xmath43 here is arbitrary since the operators @xmath36 and @xmath37 are time - independent .
( the symplectic floquet matrix is @xmath44 , where @xmath45 is obtained from @xmath46 by changing @xmath37 into @xmath47 . )
thus , ( [ subst ] ) is the linear combination of two complex conjugated floquet eigensolutions which makes @xmath48 and @xmath49 real . for the single - site breather ,
all eigenvalues @xmath38 of @xmath46 are always real , implying the linear stability of the breather for all parameter values @xmath50 @xcite .
accordingly , we can also choose the eigenvectors @xmath39 of @xmath46 to be real and normalized , in which case the phase of the amplitude @xmath51 describes the symmetry properties of the solution ( [ subst ] ) under time reversal : choosing @xmath51 real yields a time - symmetric solution , @xmath52 , while choosing @xmath51 purely imaginary yields a time - antisymmetric solution , @xmath53 .
for an infinite system , the spectrum of the ( non - hermitian ) matrix @xmath54 can generally be divided into a continuous ( phonon ) part , corresponding to extended eigenvectors , and a point spectrum corresponding to localized eigenvectors . the phonon spectrum for any localized solution @xmath16
is easily obtained from the limit @xmath55 , since the condition @xmath56 implies that the operators @xmath57 and @xmath58 become identical and eqs .
( [ w_n ] ) - ( [ u_n ] ) reduce into two uncoupled equations for the linear combinations @xmath59 resp .
assuming @xmath61 and @xmath62 , respectively , yields the dispersion relations @xmath63 from eqs .
( [ w_n ] ) - ( [ u_n ] ) .
thus , the continuous spectrum of the matrix @xmath54 consists of two branches , symmetrically located around @xmath64 , and since @xmath65 for the single - site breather these two branches never overlap .
note also that two eigenvectors with eigenvalues @xmath66 correspond to the same solution to eq .
( [ 1st ] ) ( changing the sign of @xmath38 in eq .
( [ subst ] ) is equivalent to changing @xmath67 , @xmath68 , @xmath69 ) , and therefore it is enough to consider e.g. @xmath70 , in which case @xmath71 always vanishes exponentially as @xmath72 .
when @xmath73 is not too large , the linear spectrum around the single - site breather contains also two pairs of nonzero isolated eigenvalues @xmath38 , which correspond to the two internal modes of the breather @xcite .
one of these modes is a spatially symmetric , breathing , mode , while the other is a spatially antisymmetric , translational or pinning mode .
numerically , it has been found that the breathing mode exists for @xmath74 , while the pinning mode exists for @xmath75 .
the numerically calculated internal mode frequencies @xmath38 as a function of breather frequency @xmath17 are shown in fig .
note that as @xmath76 , the breathing mode frequency approaches the lower edge of the phonon band ( but always stays outside the band @xcite ) , while the pinning mode frequency approaches zero ( but always stays nonzero ) .
this is consistent with the fact that the soliton solution of the continuous nls equation has no breathing mode ( due to its exact integrability ) , and has a translational mode with zero frequency due to the translational symmetry of the nls equation .    to obtain a complete set of solutions to eq .
( [ 1st ] ) in which an arbitrary initial perturbation @xmath77 can be expanded , we must include also the zero - frequency solutions , which generally can be written as a superposition of two fundamental modes .
one of these modes ( phase mode @xcite ) is the solution @xmath78 to the homogeneous equation ( [ w_n ] ) , @xmath79 . the corresponding perturbation @xmath80 describes a rotation of the overall phase of the breather .
the second mode ( growth mode @xcite ) is obtained by solving the inhomogeneous equation @xmath81 , which has the solution @xmath82 . the corresponding solution to eq .
( [ 1st ] ) is @xmath83 , and corresponds to a time - linear growth of the perturbation representing a small change in the breather frequency .
although the set of eigensolutions ( [ subst ] ) together with the two zero - frequency modes forms a basis for the space of solutions to eq .
( [ 1st ] ) ( there are no bifurcations , which could result in additional marginal modes @xcite with time - linear growth at degenerate eigenvalues ) , this basis is in general not orthogonal using the ordinary scalar product , and typically there is a considerable overlap between the solution corresponding to the internal breathing mode and the zero - frequency modes .
however , in analogy with e.g. @xcite , we can define a pseudo - scalar product between any two vectors @xmath84 and @xmath85 by @xmath86 this product is formally not a true scalar product , since for the general case the product of a vector with itself as defined by ( [ scalar ] ) is not necessarily positive .
however , when @xmath39 is a real eigenvector of @xmath54 we have @xmath87 from ( [ w_n ] ) , and the operator @xmath57 is positive definite for all @xmath88 @xcite . with this product , it follows from ( [ w_n ] ) - ( [ u_n ] ) that all eigensolutions with different ( real ) eigenfrequencies @xmath38 are orthogonal in the sense that @xmath89 and the only nonzero product involving the zero - frequency modes is the cross - product between the phase mode and the growth mode @xcite : @xmath90 where @xmath23 is the norm ( [ norm ] ) of the breather with frequency @xmath17 .
we also remark that , since the product ( [ scalar ] ) multiplied by a factor @xmath91 is just the symplectic product between the two vectors @xmath92 and @xmath93 , the sign of the product of an eigenvector with itself can be interpreted as the negative of the krein signature of the corresponding pair of floquet eigenvalues @xcite .
we now consider as initial state a single - site breather perturbed in the direction of a single eigenmode ( [ subst ] ) ( localized or extended ) of the linearized equations ( [ 1st ] ) , and wish to describe qualitatively the long - time effects of this perturbation using the expansion ( [ perturbation ] ) . in general , taking into account terms up to order @xmath94 in this expansion yields a solution to the dnls equation which is correct to @xmath95 , i.e. , for long but finite time - scales for small initial perturbations ( note that the expansion parameter @xmath29 plays the same role as the mode amplitude @xmath51 . ) . as is wellknown however , this kind of expansion in general diverges due to resonances between solutions to the homogeneous equation ( [ 1st ] ) and the inhomogeneous terms appearing in the right - hand sides of eqs .
( [ 2nd])-([4th ] ) and the corresponding higher - order equations . a resonance with a solution belonging to the continuous spectrum results in a bounded but nonlocalized solution corresponding to outgoing radiation , while a resonance with an eigenfunction belonging to the discrete spectrum gives a spatially localized response which diverges linearly in time .
however , up to any finite order these divergences can be systematically removed by allowing a slow time dependence of the independent variables , which in our case are taken to be the mode amplitude @xmath51 and the breather frequency @xmath17 .
this procedure adds additional terms to the equations , which can be tuned so that the divergent parts of the response disappear . in other words ,
these two quantities are used as collective variables which , together with the outgoing radiation fields , are expected to describe the main features of the asymptotic dynamics if the initial perturbation is sufficiently small .
the second order correction is given by the inhomogeneous eq .
( [ 2nd ] ) , which with the substitution ( [ subst ] ) becomes ( choosing @xmath34 , @xmath96 , and @xmath97 real without loss of generality ) : @xmath98 .\end{aligned}\ ] ] thus , the right - hand side contains one static part and one part involving the frequencies @xmath99 .
it acts as a periodic force with frequencies 0 and @xmath100 , and since all terms contain the factor @xmath34 this force is localized at the breather region .
the response to this force will remain bounded and localized unless the corresponding homogeneous equation ( [ 1st ] ) has a solution with frequency 0 resp .
@xmath101 which is non - orthogonal to the corresponding part of the right - hand side in ( [ 2nd2 ] ) . as will be shown in sec .
[ subsec : higher - order ] , a non - zero overlap between the static part of ( [ 2nd2 ] ) and the zero - frequency solutions of ( [ 1st ] ) yields a ( time - independent ) shift of the breather frequency .
moreover , if @xmath102 , so that @xmath100 is inside the phonon band of the homogeneous equation , a resonance will generally occur , resulting in radiation with frequency @xmath103 emitted from the breather region . the strength of the radiation field is determined by the ( generally nonzero ) overlap between the @xmath101-part of ( [ 2nd2 ] ) and the corresponding homogeneous solution ( see sec .
[ subsec : higher - order ] ) . in a similar way
, we obtain that the right - hand - side of the third order equation ( [ 3rd ] ) contains the frequencies @xmath38 and @xmath104 , the fourth order equation ( [ 4th ] ) contains the frequencies 0 , @xmath105 , and @xmath106 , and in general the _
p_th order equation contains as its highest harmonic the frequency @xmath107 .
accordingly , we conclude that if @xmath108 so that @xmath107 belongs to the phonon band , the perturbed breather will radiate to _
p_th order .
the consequences of this radiation for the breather itself will be discussed in secs .
[ sec : internal ] and [ sec : phonons ] for the cases of localized and extended perturbations @xmath109 , respectively .
with the initial perturbation @xmath77 of the single - site breather corresponding to a spatially localized eigenmode of the linearized equations ( [ 1st ] ) , it is clear from the discussion in sec .
[ subsec : strategy ] that higher order radiation always will be created , since for any internal mode frequency @xmath38 there is always an integer @xmath94 such that @xmath107 belongs to the phonon band and ( [ pradiation ] ) is fulfilled .
moreover , from the numerical results presented in fig .
[ fig1 ] we find that the spatially symmetric breathing mode always radiates to second order , since ( [ pradiation ] ) always is fulfilled for @xmath110 , while the antisymmetric pinning mode radiates to second order only when @xmath111 .
thus , due to this radiation from the breather the total norm contained in any finite region around the breather(i.e . , the total norm of breather + internal mode ) will always decrease with time .
however , the main concern here is the long - time effect of the internal - mode excitation on the breather itself , and thus we must investigate whether there will also be some transfer of energy between the breather and its internal mode . we will first ( sec . [
subsec : numerics ] ) show results from direct numerical integration of eq .
( [ dnls ] ) ; then we will give two alternative approaches to the analytical interpretation of these results based on the higher - order equations ( [ 2nd])-([4th ] ) ( sec .
[ subsec : higher - order ] ) resp . the conservation laws ( [ normcons])-([hamcons ] ) ( sec .
[ subsec : conservation ] ) .      in fig .
[ fig2 ] , we show a typical example on the long - time evolution of a breather when the initial perturbation is taken in the direction of its internal breathing mode . as is seen from fig .  [ fig2 ] ( a ) , the amplitude of the breathing mode decays slowly with time as a consequence of the losses due to generation of second order radiation , and a careful study of its envelope @xmath112 indicates that it decays as @xmath113 where @xmath114 is a constant .
this is consistent with a similar result obtained for the continuum nls equation with generalized ( non - cubic ) nonlinearity @xcite ; the analytical motivation for this result ( which is analogous to that of the continuum model given in ref .
@xcite ) is given in the following subsections .
however , the main result of this section is illustrated in fig .
[ fig2 ] ( b ) and ( c ) . fig .
[ fig2 ] ( b ) is obtained by calculating the time - average of the central - site intensity as @xmath115 where @xmath116 is a set of closely spaced time - instants .
it is clear that the interaction between the breather and its internal mode asymptotically leads to an _ increase _ of the average peak intensity , i.e. , to _ breather growth_. the same phenomenon is illustrated also in fig .  [ fig2 ] ( c ) , where we have plotted the difference between the instantaneous breather frequency calculated at time @xmath117 , @xmath118 , and the frequency of the unperturbed breather @xmath119 . from this figure
, we can also conclude that there are two different mechanisms causing the shift of breather frequency .
firstly , there is an initial ( almost instantaneous ) rather large frequency shift , which can be interpreted as an adaption of the initially perturbed breather to the breather which is closest to the initial condition . as is shown below in sec .
[ subsec : higher - order ] , this time - independent frequency shift , which is observed to be always positive for the breathing mode , is a consequence of the overlap between the static part of the right - hand side of the second - order equation ( [ 2nd2 ] ) and the zero - frequency modes . secondly , there is the slow , continuous increase of the breather frequency which corresponds to the slow increase of @xmath120 in fig .
[ fig2 ] ( b ) , indicating a continuous transfer of norm from the internal mode to the breather .
it is described by the static part of the right - hand side of the fourth - order equation ( [ 4th ] ) ( see sec .
[ subsec : higher - order ] ) .
when the initial perturbation of the breather is taken in the direction of its internal pinning mode we observe , just as for the breathing mode , that the breather - internal mode interaction asymptotically _ always leads to breather growth_. an example is shown in fig .
[ fig3 ] , where the parameter values have been chosen so that the lowest harmonic that enters the phonon band is @xmath121 ( @xmath122 ) , we observe two qualitative differences compared to the case with breathing mode excitation .
firstly , since in this case the first phonon resonance occurs only in the third order equation ( [ 3rd ] ) , the decay of the internal mode amplitude will be slower , and a good fit is obtained by @xmath123 .
this agrees with the general result when @xmath124 is the lowest harmonic that enters the phonon band obtained for the continuum generalized nls equation in ref .
@xcite ; the derivation of the corresponding result for the discrete case ( see eq .
( [ ageneral ] ) ) is given in sec .
[ subsec : conservation ] . as a consequence of the slower decay of the internal mode amplitude
, the breather growth will also be slower when @xmath125 , as can be seen from fig .
[ fig3 ] ( c ) by comparing the time - scales with those of fig .
secondly , the initial shift of the breather frequency for a pinning mode excitation is always much smaller than for the breathing mode excitation ( also when @xmath100 is in the phonon band ) , and is when @xmath126 also observed to be negative .
the explanation for this is given in sec .
[ subsec : higher - order ] .
however , it is important to stress that also in the cases where the initial frequency shift is negative , we find that the continuous breather growth always will give an asymptotic frequency shift which is positive .      here , we will analyze the higher - order equations ( [ 2nd])-([4th ] ) by making use of the strategy of systematically removing the appearing divergent parts as outlined in sec .
[ subsec : strategy ] ( in analogy with the treatment of the continuous nls - type equations in e.g. refs .
first , we show how the dominating contribution to the time - independent frequency shift observed in the numerical simulations above can be calculated from the static part of the right - hand side of the second - order equation ( [ 2nd ] ) .
this frequency shift can be explicitly taken into account by replacing @xmath17 in eq .
( [ perturbation ] ) with @xmath127 , where @xmath119 is the unperturbed breather frequency and @xmath128 the second - order shift to be determined .
this implies that the additional term @xmath129 will be added to the right - hand side of eq .
( [ 2nd2 ] ) .
writing the response to the static part of ( [ 2nd2 ] ) as @xmath130 with real @xmath131 and @xmath132 then yields @xmath133 with @xmath46 as defined by eq .
( [ matrix ] ) .
if the expansion of the right - hand side of ( [ mstatic ] ) in the complete set of vectors consisting of the eigenvectors of @xmath46 ( including the phase mode ) and the growth mode contains some component on either of the two zero - frequency modes , the response @xmath134 will not remain bounded but diverge linearly with time .
thus , in order to remove this divergency , the frequency shift @xmath128 must be chosen so that both these components are identically zero .
the component on the growth mode is trivially zero , while the component on the phase mode is obtained by applying the pseudo - scalar product ( [ scalar ] ) with the vector corresponding to the growth mode and using ( [ overlap ] ) .
demanding this component to be zero yields @xmath135 this is typically positive for the breathing mode , since the dominating contribution to the sum in ( [ static ] ) comes from the central site @xmath136 , and @xmath137 is always positive . for the spatially antisymmetric pinning mode , there is no contribution to this sum from the central site , since @xmath138 and @xmath139 are zero . the change from a positive to a negative frequency shift when increasing @xmath17 in this case is related to a qualitative change of the nature of the growth mode @xmath140 : for @xmath141 we find that @xmath142 for all @xmath143 , so that all terms in the sum in eq .
( [ static ] ) are negative , while @xmath144 becomes positive also for sites in the neighborhood of @xmath136 for smaller @xmath17 .    for the rest of the analysis in this subsection , we assume for calculational simplicity that the internal mode frequency is such that @xmath105 is inside the phonon band ( and thus it is not applicable for the pinning mode excitation when @xmath145 ) .
then , the non - static part of the right - hand side of the second order equation ( [ 2nd2 ] ) will generally give rise to a non - localized response , which can be written in the form @xmath146 .
this response corresponds to the radiation field going out from the breather region , and since the right - hand side of ( [ 2nd2 ] ) is spatially localized and symmetric , this field should asymptotically correspond to two identical linear waves propagating to the left ( right ) for @xmath147 ( @xmath148 ) .
thus , the boundary conditions can be written as @xmath149 with @xmath150 according to ( [ adisp ] ) .
defining for general @xmath151 the matrix @xmath152 ( cf .
( [ matrix ] ) ) as @xmath153 the functions @xmath154 and @xmath155 are seen from ( [ 2nd2 ] ) to be determined by @xmath156 since for general @xmath151 every eigenvector of @xmath157 with eigenvalue @xmath158 is also an eigenvector of @xmath46 with eigenvalue @xmath159 , the right - hand side can be expanded on the basis of eigenvectors of @xmath46 ( including the zero - frequency modes ) .
the strength of the radiation field is then given by the expansion coefficient for the ( continuous spectrum ) eigenvector of @xmath46 with eigenvalue @xmath100 , since this corresponds to the eigenvalue @xmath160 of @xmath161 , and thus a spatially non - bounded response in ( [ m2omegap ] ) .
using the orthogonality relation ( [ orthogonal ] ) , this coefficient is simply the overlap between the right - hand side of ( [ m2omegap ] ) and the eigenvector of @xmath46 with eigenvalue @xmath100 calculated with the pseudo - scalar product ( [ scalar ] ) , which is generally nonzero .
next , we show how the dominating contribution to the decay of the internal mode amplitude as given by eq .
( [ decay ] ) is obtained from the condition that @xmath162 in the third - order equation ( [ 3rd ] ) should remain bounded . to this end , we assume a slow time - dependence of the internal mode amplitude of the form @xmath163 , and consider the response to the terms with frequency @xmath38 in the right - hand side of eq .
( [ 3rd ] ) corrected by the additional terms @xmath164 appearing as a consequence of including the time - dependence of @xmath51 and the second - order frequency shift @xmath128 from eq .
( [ static ] ) in the perturbation expansion ( [ perturbation ] ) .
writing the response to this part as @xmath165 yields @xmath166 a bounded response for @xmath167 exists only if the vector on the right - hand side of ( [ momegap ] ) has no component in the direction of the internal mode eigenvector @xmath39 , since this is the eigenvector corresponding to the eigenvalue zero of @xmath168 . using the orthogonality relation ( [ orthogonal ] ) , this component is obtained by application of the product ( [ scalar ] ) with the vector @xmath39 , and the condition that this component must be zero determines the time - evolution of @xmath51 . considering only the absolute value @xmath169 , the resulting equation has the form @xmath170 is given by @xmath171 } { \sum_n u_n w_n } = \frac{8 c \omega_p
|r_2|^2 \sin q_2 }   { \sum_n w_n { \cal l}_0 w_n } > 0 , \ ] ] where the second equality is obtained using eqs .
( [ w_n ] ) , ( [ bc ] ) and ( [ m2omegap ] ) , and the positivity of @xmath172 follows from the fact that , as mentioned in sec .  [ subsec : linear ] , the operator @xmath57 is positive definite for all @xmath173 @xcite .    finally , we show how the continuous increase of the breather frequency appears from the divergent response to the static part of the right - hand side of the fourth - order equation ( [ 4th ] ) . in its unmodified form ,
this part is given by @xmath174 now , it is clear from ( [ static ] ) and ( [ decay ] ) that the time - dependence of @xmath51 will induce a time - dependence of the second - order frequency shift @xmath128 of the form @xmath175 , so that we can express the total breather frequency up to order @xmath176 as @xmath177 , where a 4th order correction @xmath178 has also been included .
then , we must also take into account the time - dependence of the breather shape @xmath34 by writing @xmath179 . as a consequence
, the term @xmath180 will be added to the expression ( [ static4 ] ) for @xmath181 in the right - hand side of ( [ 4th ] ) , and writing the response to this total force as @xmath182 with real @xmath183 and @xmath184 yields @xmath185 the response @xmath186 will be bounded in time only if the right - hand side of ( [ m4 ] ) has no component either on the growth mode or on the phase mode , which gives two conditions for the determination of @xmath187 and @xmath178 . using ( [ scalar])-([overlap ] )
, we obtain by demanding the expansion coefficient for the growth mode to be zero : @xmath188 \right . \nonumber\\ \left .
+ \phi_n \left[2 u_n w_n   { \rm im}\left ( u_n^{(2)}\right)+\left(3w_n^2-u_n^2\right )   { \rm im}\left ( w_n^{(2)}\right ) \right ] \right\ } .\end{aligned}\ ] ] ( similarly , @xmath178 is obtained by demanding the component on the phase mode to be zero . ) thus , the dominating contribution to the frequency growth should be of order @xmath189 , so that with the approximate time - dependence ( [ decay ] ) of the internal - mode amplitude we obtain qualitatively @xmath190 which is in good agreement with the numerically observed time - dependence of the frequency shift as shown in fig .
[ fig2 ] ( c ) .
however , the positivity of @xmath187 is not easily seen from the expression ( [ lambda2dot ] ) , and therefore we will in the following subsection derive an alternative expression from which the positivity follows immediately , using the conservation laws for the norm resp .
hamiltonian .
we consider first the conservation law ( [ normcons ] ) for the total norm ( [ norm ] ) contained in any large but finite region around the breather .
averaging over a time - interval @xmath191 $ ] and using eqs .
( [ perturbation ] ) , ( [ subst ] ) , we can write the time - averaged norm to second order in @xmath192 ( putting @xmath193 ) as @xmath194 . \label{normav}\ ] ] ( note that there will be no contribution at order 2 from the static second - order correction @xmath131 , since the renormalization of the breather frequency @xmath17 according to ( [ static ] ) yields @xmath195 . ) in the case when @xmath105 belongs to the phonon band , we obtain the following balance equation , which is correct up to order @xmath196 : @xmath197 where we have used eqs .
( [ normcons ] ) , ( [ jndef ] ) and ( [ bc ] ) , and @xmath198 is as defined below eq .
( [ bc ] ) .
similarly , we can use the conservation law ( [ hamcons ] ) for the hamiltonian ( [ dnlsham ] ) together with the general expression for the hamiltonian flux density ( [ jhdef ] ) for a small - amplitude plane wave @xmath199 , @xmath200 to write the balance equation for the total time - averaged hamiltonian in the same region for the case of second order radiation : @xmath201 which is also correct to order @xmath196 . the lowest order contribution to the derivative @xmath202
can be obtained using the first equality in the equation of motion ( [ dnls ] ) and its complex conjugate as follows : @xmath203 which is always negative for an internal mode excitation since @xmath204 .
thus , we can combine the two balance equations ( [ dndt ] ) and ( [ dhdt ] ) , and using eqs .
( [ nh ] ) and ( [ dhda2 ] ) we obtain the expression ( [ decay ] ) for @xmath112 with @xmath172 as in ( [ gamma ] ) , together with the following expression for the frequency growth rate from which its positivity is immediately seen : @xmath205    this approach also has the advantage that it is easily generalized to the case where the lowest harmonic that enters the phonon band is @xmath107 with @xmath206 , i.e. , for the pinning mode excitation when @xmath207 .
then , we can write the boundary conditions at the infinities to lowest order in @xmath51 as @xmath208}$ ] , @xmath209 , where @xmath210 and @xmath211 according to ( [ adisp ] ) . consequently , we can proceed exactly as above , writing down the balance equations for the norm and hamiltonian to order @xmath212 just by modifying the right - hand sides of eqs .
( [ dndt ] ) resp .
( [ dhdt ] ) by replacing @xmath196 with @xmath212 , @xmath213 with @xmath214 , @xmath198 with @xmath215 , and @xmath105 with @xmath107 . combining the balance equations yields the following general expression for the time - dependence of the internal mode amplitude , @xmath216^ { 1/(2p-2)}}\;,\ ; \gamma_p = \frac{4pc|r_p|^2 \sin q_p } { \sum_n u_n w_n } > 0 , \label{ageneral}\ ] ] which is the analog to the expression obtained with similar arguments in @xcite for the continuum nls models .
and most importantly , we obtain a general expression for the breather frequency growth rate which is positive for all @xmath94 : @xmath217 thus , integrating eq .
( [ lambdadotgeneral ] ) using the time - dependence ( [ ageneral ] ) of the internal - mode amplitude , we obtain that the dominating contribution to the frequency growth generally can be written qualitatively as @xmath218    let us finally point out that the approach used in this section provides a simple physical interpretation of why the interaction between the breather and its internal mode in addition to generate radiation also should lead to breather growth . since the expression ( [ hamflowpw ] ) for the hamiltonian flux density of a small - amplitude plane wave is positive when @xmath219 and @xmath220 have the same sign , the hamiltonian energy for a plane wave always propagates in the same direction as the wave itself .
thus , the second ( or higher ) order radiation emitted from the breather region will always carry away a _ positive _ amount of the hamiltonian energy , or , equivalently , negative hamiltonian energy will flow into the breather region . moreover , from eq .
( [ dhda2 ] ) it is clear that the contribution to the hamiltonian from the internal mode always is negative and a monotonously _ decreasing _ function of its amplitude , and thus the decay of the internal mode would cause an _ increase _ of the hamiltonian in the breather region . consequently , since the hamiltonian of the pure breather is a monotonously decreasing function of the breather frequency , the breather should grow so that the total hamiltonian of ( breather+internal mode ) decreases .
a similar mechanism was recently found to cause soliton growth in the parametrically driven continuum nls equation in the regime of oscillatory instability @xcite , and this type of argument has also been used to explain the quasi - collapse of a broad excitation to a narrow localized state in the two - dimensional dnls equation @xcite .
we choose , as in the previous section , the initial perturbation @xmath77 to be an eigensolution of the linearized equations ( [ 1st ] ) , but now we consider the case of a spatially extended perturbation .
choosing without loss of generality a solution of the form ( [ subst ] ) with positive frequency @xmath38 yields the asymptotic behaviour @xmath221 where the wave vector @xmath220 @xmath222 is determined by the dispersion relation ( [ adisp ] ) , and @xmath223 is the phase shift across the breather .
thus , the excitation of an extended eigenmode corresponds to an interaction between the breather and a non - propagating , standing - wave phonon with small amplitude @xmath51 .
as was mentioned in the introduction , these standing - wave phonons are generally unstable , but since the instabilities become exponentially weak in the small - amplitude limit they are expected to have very little effect on the breather for the perturbation sizes and time - scales considered in this section .
we will return to the effects on the breather of these instabilities in sec .
[ sec : instabilities ] , where larger perturbations are considered .
in contrast to the case of excitation of a localized internal mode discussed in sec .  [
sec : internal ] , where higher order radiation always was emitted from the breather region , the condition ( [ pradiation ] ) yields that for the standing - wave perturbation the breather will radiate to higher order only if @xmath224 , so that both @xmath38 and @xmath225 are inside the phonon band . in terms of the phonon wave vector @xmath220 , this means that there is a critical value @xmath226 , @xmath227 such that for @xmath228 second order radiation will be emitted from the breather region , while for @xmath229 ( @xmath230 ) _ all _ multiples of @xmath38 are outside the phonon band , and no higher order radiation is emitted . as will be shown below , these two regions yield qualitatively different scenarios for the long - time evolution of the perturbed breather .
we also note that for @xmath231 we have @xmath232 , so that for the highly localized , high - frequency breathers _ no phonons can generate higher - order radiation _
( note also that there are no internal modes in this regime ) .
furthermore , we always have @xmath233 , so that the regime of higher - order radiation generation is a subset of the regime @xmath234 where modulational instability for _ travelling _ plane waves occurs @xcite .
let us first discuss the case @xmath235 .
a typical example of the long - time evolution of a breather interacting with a small - amplitude standing - wave phonon with @xmath235 is illustrated by fig .
[ fig4 ] . as is seen from fig .  [ fig4 ] ( a )
, the amplitude of the oscillations remains essentially constant in time , but a closer inspection reveals that the average value of @xmath236 asymptotically _ increases _ with an apparently constant rate ( see inset in fig .
[ fig4 ] ( b ) ) . similarly , fig .
[ fig4 ] ( b ) shows that also the total norm contained in any finite region around the breather asymptotically increases linearly with time .
we find that these results are generic for all cases when the phonon wave vector @xmath235 ( the spatial symmetry of the phonon is not important for the asymptotic behaviour ) , and thus we conclude that in this regime , the breather can pump energy from the phonon ( which is infinite for an infinite system ) , and thereby grow .    in the same spirit as for the internal mode excitation in sec .
[ subsec : conservation ] , we can give a simple argument based on the conservation laws to motivate why the generation of second order radiation should lead to breather growth . to this end , we assume that the initial standing - wave phonon is infinitely extended , and that far away from the breather a stationary regime will be reached corresponding to the following boundary conditions @xmath237 e ^{i\lambda t }   \ ; ; \ ; n \rightarrow \pm \infty .
% \label{scattering}\ ] ] thus , we have taken into account the second - order radiation with frequency @xmath100 generated at the breather region but neglected possible higher - order radiation ; moreover the resonance at the original phonon frequency @xmath38 in the third - order equation ( [ 3rd ] ) has been taken into account by allowing the incoming and outgoing complex amplitudes @xmath51 and @xmath238 to be different .
we can then , in analogy with eqs .
( [ dndt ] ) and ( [ dhdt ] ) , write the balance equations for the total norm and hamiltonian contained in a region around the breather averaged over a time - interval @xmath191 $ ] as @xmath239 ,   \label{dndtph}\ ] ] and @xmath240 ,   \label{dhdtph}\end{aligned}\ ] ] respectively . here , we have used the facts that the time - average of the norm current density ( [ jndef ] ) and the hamiltonian flux density ( [ jhdef ] ) are additive quantities for small - amplitude plane waves , and that in the stationary regime , the mode amplitudes @xmath51 , @xmath238 and @xmath213 are time - independent . combining eqs .
( [ dndtph ] ) and ( [ dhdtph ] ) and using ( [ nh ] ) , we obtain that the breather frequency grows with a constant rate given by @xmath241 the physical interpretation of this result is , similarly as for the case of internal mode excitation , that the generation of higher - order radiation results in a net flow of negative hamiltonian energy into the breather region , which is absorbed by the breather by increasing its frequency and maximum amplitude .
this process is similar to the one observed for the two - channel phonon scattering on breathers in a klein - gordon model with a morse potential in ref .
@xcite ; however in the latter case the second outgoing wave resulted from a resonance in the linearized equations and were therefore of the same order of magnitude as the incoming wave , and moreover the outcome in this case was breather decay since the energy of a klein - gordon breather is an increasing function of its amplitude .
considering now the regime @xmath242 where all multiples of @xmath38 are outside the phonon band , the most important conclusion from our extensive numerical investigations is that we _
never observe breather growth_. instead , we sometimes ( but not always ) observe a very slow decrease of @xmath120 , and an increase of the fluctuations around this mean - value .
this behaviour is illustrated by fig .
[ fig5 ] ( a ) and ( b ) .
a possible interpretation of these results is that , since the higher harmonics which are created by the breather - phonon interaction can not propagate , they stay trapped around the breather .
thus , this could lead to a transfer of energy from the pure breather , which acquires more and more internal frequencies and becomes a chaotic breather @xcite .
another possible interpretation is that the increase of the oscillation amplitude is connected with the oscillatory instabilities of the standing waves ; as we will show in the next section these instabilities provide a mechanism for breather decay . however in some cases , illustrated by fig .
[ fig5 ] ( c ) and ( d ) , the oscillation amplitude as well as its average value apparently approaches a constant limit value .
we have at present no explanation for this behaviour ( as will be discussed in the next section , there exist exact phonobreather solutions which could be candidates for such a final state , but they are unstable ) ; it is possible that the time ranges that we were able to study with sufficient numerical accuracy in these cases simply were too short to observe the scenario described by fig .
[ fig5 ] ( a ) and ( b ) .    to conclude this section , we repeat our main result that breather growth is observed if and only if @xmath235 , where @xmath226 is given by ( [ qc ] ) . the fact that @xmath232 for @xmath243 thus implies the existence of an upper limit beyond which the breather can not grow with the type of perturbations considered here .
we would also like to relate our results to recent numerical simulations of breathers interacting with _ propagating _ phonons in klein - gordon @xcite and fpu @xcite lattices .
for the klein - gordon lattice with a ( soft ) morse on - site potential , phonons with small wave vector @xmath220 were observed to yield breather growth , while phonons with large @xmath220 caused breather decay . for a fpu lattice with hard anharmonicity ,
the opposite situation was observed , i.e. , small-@xmath220 phonons caused breather decay and large-@xmath220 phonons breather growth .
the fact that the situation for the hard fpu lattice was opposite to that of the soft klein - gordon lattice could be expected , since in the former case the modulational instability occurs for large @xmath220 , whereas soft klein - gordon and dnls lattices with @xmath14 are modulationally unstable for small @xmath220 .
however , we stress that the relation between plane - wave modulational instability and breather growth is nontrivial , and at least for the case considered in this paper the critical value @xmath226 for breather growth from interaction with standing - wave phonons differs from the critical value @xmath244 for modulational instability of travelling waves .
in this section , with the aim at describing the interaction between a breather and a standing - wave phonon with non - negligible amplitude , we will take a slightly different point of view than in the preceding sections . instead of choosing as initial condition an exact breather solution and adding a perturbation corresponding to an eigenmode of the linearized equations
, we will here consider initial conditions which are exact _ phonobreather _ @xcite ( or _ nanopteron _ ) solutions . by definition ,
a phonobreather consists of a spatially localized breather on top of a spatially extended tail which is a nonlinear , standing - wave phonon .
( there are also solutions where the tail is a propagating wave @xcite , but we will not discuss them further here . )
phonobreathers exist generically for nonlinear lattice - equations ( see e.g. @xcite ) , but their existence normally requires an integer relationship between the breather and phonon frequencies .
however , for the dnls equation phonobreathers exist for _ any _ ( rational or irrational ) relation between the two frequencies @xcite , as a consequence of the additional invariance of the equation under global phase transformations .    since the phonobreathers are exact solutions consisting of a breather part and a standing - wave part , one could expect them to be attractors for the initial conditions considered in sec .
[ sec : phonons ] .
however , as was shown recently @xcite , generically for soft klein - gordon and dnls models with @xmath14 _ all _ phonobreathers with phonon wave vector @xmath245 will be _ linearly unstable _ if , for fixed phonon amplitude @xmath51 , the linear coupling @xmath246 is larger than some threshold value @xmath247 ( i.e. , away from the anticontinuous limit ) .
( for lattices with hard potentials and dnls with @xmath248 , the stable phonobreather has @xmath249 . )
these instabilities are caused by an oscillatory instability of the standing - wave phonon itself , which can be understood by considering the construction of a nonlinear standing wave with wave vector @xmath220 close to @xmath250 at the anticontinuous limit @xmath251 by introducing a periodic array of discommensurations or defects in the nonlinear phonon with wave vector @xmath250 and amplitude @xmath51 , @xmath252 , which is linearly stable for all @xmath51 and @xmath14 @xcite . in the anticontinuous limit @xmath251 ,
each defect consists when @xmath253 of one extra site with @xmath254 added to the @xmath250-phonon , which consequently suffers an additional phase shift of @xmath250 across each defect . for @xmath234 ,
each defect consists of several consecutive zero - amplitude sites with associated phase shifts ( a general method for generating the anticontinuous coding sequence for standing - wave phonons from a circle map is described in @xcite ) ; in this case it is also useful to consider the periodic repetition of sites with @xmath255 as defects of the zero - amplitude state .
the limit case of one isolated zero - amplitude defect , which is a discrete counterpart of the dark - soliton solution of the continuum nls equation , was investigated in @xcite .
the linear stability analysis of this mode showed that , although it is stable close to the anticontinuous limit , it suffers a bifurcation for @xmath256 where two pairs of eigenvalues of the eigenvalue problem ( [ matrix ] ) go out in the complex plane .
the resulting oscillatory instability occurs due to a resonance between a mode localized around the defect ( the defect pinning mode ) and linear radiation modes .
it was shown that for finite systems , the mode recovers its stability above some upper critical value of @xmath257 ( since the wavelength of the resonating linear modes becomes larger than the system size ) ; however this critical value increases with system size so that in the limit of an infinite system , the instability persists for all @xmath258 but with a growth - rate that decreases in an exponential - like fashion when approaching the continuum limit @xmath259 .
this instability was shown to result in the defect becoming mobile ( in nls terms , the stationary black soliton with zero minimum intensity transforms into a moving grey soliton with non - zero minimum intensity ) and radiation being emitted . in terms of the phase dynamics , this describes a moving , slowly spreading phase kink .
the instability scenario for the standing - wave phonons is basically the same as for the isolated defect , with the essential difference that the localized pinning modes associated with the individual , periodically repeated defects now will form a continuous defect band. in general , the krein signature @xcite of this defect band is opposite to that of the bands associated with the non - zero amplitude sites @xcite ( for @xmath234 there are generally several defect bands , but they will have the same krein signature ) , and as a consequence resonances between the bands will occur if the linear coupling @xmath246 is large enough , giving rise to similar oscillatory instabilities as described above ( details are given in @xcite ) .
we remark that earlier analysis @xcite of standing waves in nonlinear lattices , based on a quasi - continuum approximation , did not reveal these instabilities since their origin is the discrete nature of the lattice .
let us now return to the main objective of this section , namely to study the effect of the oscillatory standing - wave instabilities on the phonobreathers .
we find that the families of phonobreathers which are stable close to the anticontinuous limit and whose tails approach harmonic standing waves in the small - amplitude limit can be constructed from anticontinuous standing - wave solutions at @xmath251 , placing the breather at a zero - amplitude site of the phonon and adjusting it so that the resulting solution is either symmetric or antisymmetric around the breather site . denoting the anticontinuous breather amplitude by @xmath260 ( the phase of the breather site is unimportant when @xmath261 @xcite ) , this yields the following possibilities :    \(i ) for @xmath262 the antisymmetric anticontinuous solution ( here @xmath263 ) @xmath264 with the asymptotic behaviour @xmath265 in the continuum limit @xmath259 .
( note that this solution is antisymmetric only at @xmath251 . )
\(ii ) from ( [ i ] ) we can construct a symmetric solution for @xmath262 by introducing an additional phase shift of @xmath250 at one side of the breather site , giving for @xmath263 @xmath266 with the asymptotic behaviour @xmath267 in the continuum limit .
\(iii ) for @xmath268 , @xmath269 even , the number of consecutive zero - amplitude sites is odd and the antisymmetric anticontinuous solution is for e.g. @xmath270 @xmath271 behaving as @xmath272 in the continuum limit .
\(iv ) from ( [ iii ] ) the symmetric solution for @xmath268 , @xmath269 even , is constructed by a phase shift as above , giving for @xmath270 @xmath273 with the asymptotic behaviour @xmath274 in the continuum limit .
\(v ) for @xmath268 , @xmath269 odd , the number of consecutive zero - amplitude sites is even and we must add an extra site to obtain the antisymmetric anticontinuous solution , which for e.g. @xmath275 becomes @xmath276 behaving as @xmath277 in the continuum limit .
( a solution with similar properties is obtained by instead removing one zero - amplitude site ; however its symmetric counterpart is always unstable . )
\(vi ) the symmetric counterpart of ( [ v ] ) for @xmath268 , @xmath269 odd , is constructed by a phase shift as above , giving for @xmath275 @xmath278 with the asymptotic behaviour @xmath279 in the continuum limit .
a typical example on the time - evolution for an initially very weakly perturbed phonobreather with phonon wave vector @xmath262 and phonon amplitude small but non - negligible compared to the breather amplitude is illustrated in fig .
( the example in the figure belongs to type ( i ) , but similar dynamics is observed also for the spatially symmetric states of type ( ii ) . )
we can clearly distinguish two different steps leading to the final breather destruction .
the first step is the linear oscillatory instability described above , which leads to the generation of new internal frequencies of the breather , and to the movement of the defect sites in a similar way as for the case of an isolated defect . in the second step ,
the moving defects start interacting , and a close inspection of fig .
[ fig6 ] ( a ) shows that neighboring defects tend to merge and create regions of accumulated phase fluctuations travelling around in the lattice .
these will interact with the breather , and apparently cause its decay .
when the breather has decayed sufficiently to have an excitable pinning mode , it will start to move in the lattice but with rapidly decreasing amplitude , and it will finally be destroyed .
we have at present no complete understanding for the mechanism by which the interaction of the breather with the moving phase kinks cause its decay , but we remark that a similar scenario was observed when adding to the dnls equation an external , parametric white noise term @xcite . in that case , the white - noise approximation allowed a qualitative understanding of the breather decay as a consequence of phase fluctuations by using a collective coordinate approach .
however , to observe this scenario for breather destruction it is necessary ( at least for a finite system ) that the phonon amplitude is not too small compared to the breather amplitude .
if we increase the breather amplitude in fig .
[ fig6 ] ( or decrease the phonon amplitude ) sufficiently , we find that although the oscillatory instability develops , the fluctuations created in the second step will be too weak to cause the breather to decay , and it will live seemingly forever as a chaotic phonobreather. the absence of decay for small perturbations can be viewed as a consequence of the fact that the single - site dnls breather is nonlinearly ( lyapunov ) stable for norm - conserving perturbations , in the sense that @xmath280 remains arbitrarily close to the breather for all times if the initial perturbation is small enough @xcite .
thus , it is clear that for finite systems , the breather can not be destroyed unless the phonon amplitude exceeds some critical value , while nothing can be said about the infinite system since any infinitely extended phonon obviously has an infinite norm .    with @xmath234 ,
the first step resulting from the oscillatory instability occurs in a similar way as for @xmath262 : the breather acquires new frequencies and the small - amplitude sites of the phonon start moving . by instead making the interpretation that the sites with non - zero amplitude at the anticontinuous limit are defects in the zero - amplitude state ,
their movement can be seen as a consequence of the repulsive interaction between spatially separated , small - amplitude breathers with opposite phases observed e.g. in @xcite .
we also observe , similarly as for @xmath262 , the merging of neighboring defects , but in this case their interaction with the breather will not lead to breather decay , but rather to breather growth if the original amplitude of the breather is not too large .
a typical example is illustrated in fig .
it can be seen by a careful inspection of fig .
[ fig7 ] ( a ) how the merging of small - amplitude sites results in localized humps of larger amplitude reminiscent of small - amplitude moving breathers travelling around in the lattice .
the interaction of these humps with the original breather leads to growth of the latter in a similar way as observed in @xcite .
however , this growth stops when the breather amplitude has reached a critical value which is close to ( but apparently smaller than ) that corresponding to the limit value @xmath281 for small - amplitude perturbations found in sec .
[ sec : phonons ] ( the latter corresponds to @xmath282 ) .
the final state appears also here to be a chaotic phonobreather ; we have followed the time evolution of this kind of state for times up to @xmath283 without seeing any signs of decay .
also , if the initial breather frequency is chosen above the critical value @xmath281 , we typically do not observe breather growth ; instead the mean value of the chaotic amplitude oscillations resulting from the oscillatory instability remains close to the initial amplitude .
investigating the interaction between discrete nonlinear schrdinger breathers and small perturbations , we have found firstly that exciting an internal mode of the breather always leads to a slow energy transfer to the breather , i.e. , to breather growth . furthermore , we found that a dnls - breather can pump energy from a small - amplitude standing - wave phonon , provided that the phonon wave vector is smaller than the critical value @xmath226 given by eq .
( [ qc ] ) . in both cases
, the mechanism for breather growth involves the higher - order generation of radiating modes .
since this mechanism disappears at the threshold value @xmath281 of the breather frequency , it is impossible for a breather to grow beyond this value with the type of small - amplitude perturbations considered here . to analyze the interaction between breathers and standing - wave phonons of small but non - negligible amplitude , we considered the long - time evolution of weakly perturbed exact phonobreather solutions .
the instabilities of these , originating in oscillatory instabilities of the nonlinear phonons , where shown to lead to propagating inhomogeneities whose interaction with the breather provided a mechanism for breather decay and destruction ( when the phonon wave vector @xmath262 ) or growth ( @xmath284 and @xmath285 ) .    as was mentioned already in the introduction , the existence of the two conserved quantities ( [ dnlsham ] ) and ( [ norm ] ) makes the dnls equation non - generic among nonlinear lattice equations , and it is therefore necessary to investigate to what extent the results obtained in this paper apply also for klein - gordon and fpu lattices .
we plan to address these questions in a forthcoming publication , but let us stress already here that the perturbational approach used here for the dnls equation needs to be modified to account for the fact that generically , the dynamics of the breather involves also higher harmonics of its fundamental frequency . moreover , the approach used in secs .
[ subsec : conservation ] and [ sec : phonons ] based on the conservation laws can not be directly applied in the absence of a second conserved quantity . however , in view of the wide applicability of the dnls equation ( and in particular its appearance as a limit case of general lattice equations as mentioned in the introduction ) , we believe that the mechanisms for breather growth and destruction described in this paper are essential ingredients also for the corresponding processes in general lattice models .
we thank yu .
s. kivshar for giving us an early preprint of ref .
@xcite , i. v. barashenkov for directing our attention to ref .
@xcite , and a. m. morgante for discussions on phonobreathers and standing - wave instabilities .
m. j. acknowledges a marie curie research training grant from the european community .
a preliminary version of these results was presented at the conference nonlinearity  99 ( heraklion , may 10 - 14 , 1999 ) .
t. dauxois , m. peyrard , and c. r. willis , physica d * 57 * , 267 ( 1992 ) ; t. dauxois , m. peyrard , and a. r. bishop , phys . rev .
e * 47 * , 684 ( 1993 ) ; m. peyrard , in _ nonlinear cooperative phenomena in biological systems _ , edited by l. matsson ( world scientific , singapore , 1998 ) , p. 276 , and references therein .
a. m. morgante , oral communication at the conference _
nonlinearity  99 _ , heraklion , 10 - 14 may 1999 , organized by t. geisel and g. p. tsironis ; a. m. morgante , g. kopidakis , m. johansson , and s. aubry , in preparation .
b. gaididei , s. f. mingaleev , p. l. christiansen , and k.  .
rasmussen , phys .
e * 55 * , 6141 ( 1997 ) ; k.  .
rasmussen , p. l. christiansen , m. johansson , yu .
b. gaididei , and s. f. mingaleev , physica d * 113 * , 134 ( 1998 ) .                                , @xmath286 , and @xmath287 .
( a ) shows the time - evolution of the central - site intensity @xmath236 , ( b ) shows its time - average @xmath120 calculated using eq .
( [ average ] ) , while ( c ) shows the instantaneous shift of breather frequency @xmath288 .
the solid line in the main figure in ( c ) is a fit using eq .
( [ shift ] ) with @xmath289 , @xmath290 , @xmath291 , and @xmath292.,title="fig:",height=332 ] , @xmath286 , and @xmath287 .
( a ) shows the time - evolution of the central - site intensity @xmath236 , ( b ) shows its time - average @xmath120 calculated using eq .
( [ average ] ) , while ( c ) shows the instantaneous shift of breather frequency @xmath288 .
the solid line in the main figure in ( c ) is a fit using eq .
( [ shift ] ) with @xmath289 , @xmath290 , @xmath291 , and @xmath292.,title="fig:",height=332 ] , @xmath286 , and @xmath287 .
( a ) shows the time - evolution of the central - site intensity @xmath236 , ( b ) shows its time - average @xmath120 calculated using eq .
( [ average ] ) , while ( c ) shows the instantaneous shift of breather frequency @xmath288 .
the solid line in the main figure in ( c ) is a fit using eq .
( [ shift ] ) with @xmath289 , @xmath290 , @xmath291 , and @xmath292.,title="fig:",height=332 ]    , @xmath293 , and @xmath287 .
( a ) shows the time - evolution of @xmath294 , where @xmath136 is the central site of the breather , ( b ) ( main figure ) shows the instantaneous shift of breather frequency @xmath288 , and inset in ( b ) shows the time - average @xmath120.,title="fig:",height=521 ] , @xmath293 , and @xmath287 .
( a ) shows the time - evolution of @xmath294 , where @xmath136 is the central site of the breather , ( b ) ( main figure ) shows the instantaneous shift of breather frequency @xmath288 , and inset in ( b ) shows the time - average @xmath120.,title="fig:",height=521 ]    .
parameters are @xmath295 , @xmath296 ( @xmath297 ) , @xmath298 , and @xmath287 .
( a ) shows the time - evolution of the central - site intensity @xmath236 , ( b ) ( main figure ) shows the total norm contained in a region of 120 sites around the breather , and inset in ( b ) shows the time - average @xmath299.,title="fig:",height=498 ] .
parameters are @xmath295 , @xmath296 ( @xmath297 ) , @xmath298 , and @xmath287 .
( a ) shows the time - evolution of the central - site intensity @xmath236 , ( b ) ( main figure ) shows the total norm contained in a region of 120 sites around the breather , and inset in ( b ) shows the time - average @xmath299.,title="fig:",height=498 ]    .
parameters are @xmath295 , @xmath287 , and in ( a ) , ( b ) @xmath300 ( @xmath301 ) and @xmath302 , resp . in ( c ) , ( d ) @xmath303 ( @xmath304 ) and @xmath298 . ( a ) , ( c ) show the time - evolution of the central - site intensity @xmath236 , while ( b ) , ( d ) show its time - average.,title="fig:",height=302 ] .
parameters are @xmath295 , @xmath287 , and in ( a ) , ( b ) @xmath300 ( @xmath301 ) and @xmath302 , resp . in ( c ) , ( d ) @xmath303 ( @xmath304 ) and @xmath298 . ( a ) , ( c ) show the time - evolution of the central - site intensity @xmath236 , while ( b ) , ( d ) show its time - average.,title="fig:",height=302 ] .
parameters are @xmath295 , @xmath287 , and in ( a ) , ( b ) @xmath300 ( @xmath301 ) and @xmath302 , resp . in ( c ) , ( d ) @xmath303 ( @xmath304 ) and @xmath298 .
( a ) , ( c ) show the time - evolution of the central - site intensity @xmath236 , while ( b ) , ( d ) show its time - average.,title="fig:",height=302 ] .
parameters are @xmath295 , @xmath287 , and in ( a ) , ( b ) @xmath300 ( @xmath301 ) and @xmath302 , resp . in ( c ) , ( d ) @xmath303 ( @xmath304 ) and @xmath298 . ( a ) , ( c ) show the time - evolution of the central - site intensity @xmath236 , while ( b ) , ( d ) show its time - average.,title="fig:",height=302 ]    ) with @xmath263 for a system of 120 sites ( periodic boundary conditions ) , perturbed only by the numerical truncation errors .
the breather frequency is @xmath305 , phonon frequency @xmath306 ( phonon amplitude @xmath307 ) , and @xmath287 .
( a ) shows @xmath308 , while ( b ) ( main figure ) shows the intensity of the breather central site @xmath236 .
inset in ( b ) shows the inverse participation number @xmath309 , which gives a qualitative measure of the degree of localization.,title="fig:",height=476 ] ) with @xmath263 for a system of 120 sites ( periodic boundary conditions ) , perturbed only by the numerical truncation errors .
the breather frequency is @xmath305 , phonon frequency @xmath306 ( phonon amplitude @xmath307 ) , and @xmath287 .
( a ) shows @xmath308 , while ( b ) ( main figure ) shows the intensity of the breather central site @xmath236 .
inset in ( b ) shows the inverse participation number @xmath309 , which gives a qualitative measure of the degree of localization.,title="fig:",height=476 ]    ) with @xmath270 for a system of 120 sites ( periodic boundary conditions ) , perturbed only by the numerical truncation errors .
the breather frequency is @xmath310 , phonon frequency @xmath311 ( phonon amplitude @xmath312 ) , and @xmath287 .
( a ) shows @xmath308 , while ( b ) shows the intensity of the breather central site @xmath236.,title="fig:",height=519 ] ) with @xmath270 for a system of 120 sites ( periodic boundary conditions ) , perturbed only by the numerical truncation errors .
the breather frequency is @xmath310 , phonon frequency @xmath311 ( phonon amplitude @xmath312 ) , and @xmath287 .
( a ) shows @xmath308 , while ( b ) shows the intensity of the breather central site @xmath236.,title="fig:",height=519 ]