Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: synchronization is a dynamical process where two or more interacting oscillatory systems end up with identical movement . in 1665 huygens experimented with maritime pendulum clocks and discovered the anti - phase synchronization of two pendulums clocks mounted on the common frame  @xcite . since then , synchronization has become a basic concept in nonlinear and complex systems @xcite . such systems include , but are not limited to , musical instruments , electric power systems , and lasers . there are numerous applications in mechanical  @xcite and electrical  @xcite engineering . new applications are found in mathematical biology such as synchronous variation of cell nuclei , firing of neurons , forms of cooperative behavior of animals and humans  @xcite . recently , huygen s experiment has been widely discussed and several experimental devices were built  @xcite . it was shown that two real mechanical clocks when mounted to a horizontally moving beam can synchronize both in - phase and anti - phase  @xcite . in all these experiments synchronization was achieved due to energy transfer via the oscillating beam , supporting huygen s intuition  @xcite . one of the rapidly developing areas in between physics and mathematics , is the topic of @xmath1-symmetry , which has started as a way to characterize non - hermitian hamiltonians in quantum mechanics  @xcite . the key idea is that a linear schrdinger operator with a complex - valued potential , which is symmetric with respect to combined parity ( @xmath2 ) and time - reversal(@xmath3 ) transformations , may have a real spectrum up to a certain critical value of the complex potential amplitude . in nonlinear systems , this distinctive feature may lead to existence of breathers ( time - periodic solutions localized in space ) as continuous families of their energy parameter . the most basic configuration having @xmath1 symmetry is a dimer , which represents a system of two coupled oscillators , one of which has damping losses and the other one gains some energy from external sources . this configuration was studied in numerous laboratory experiments involving electric circuits  @xcite , superconductivity  @xcite , optics  @xcite and microwave cavities  @xcite .    in the context of synchronization of coupled oscillators in a @xmath1-symmetric system , one of the recent experiments was performed by bender et al . these authors considered a @xmath1-symmetric hamiltonian system describing the motion of two coupled pendula whose bases were connected by a horizontal rope which moves periodically in resonance with the pendula . the phase transition phenomenon , which is typical for @xmath1-symmetric systems , happens when some of the real eigenvalues of the complex - valued hamiltonian become complex . the latter regime is said to have _ broken _ @xmath1 symmetry .    on the analytical side , dimer equations were found to be completely integrable  @xcite . integrability of dimers is obtained by using stokes variables and it is lost when more coupled nonlinear oscillators are added into a @xmath0-symmetric system . nevertheless , it was understood recently @xcite that there is a remarkable class of @xmath0-symmetric dimers with cross - gradient hamiltonian structure , where the real - valued hamiltonians exist both in finite and infinite chains of coupled nonlinear oscillators . analysis of synchronization in the infinite chains of coupled oscillators in such class of models is a subject of this work .    in the rest of this section , we describe how this paper is organized . we also describe the main findings obtained in this work . section [ model ] introduces the main model of coupled pendula driven by a resonant periodic movement of their common strings . see figure [ pendulapic ] for a schematic picture . by using an asymptotic multi - scale method , the oscillatory dynamics of coupled pendula is reduced to a @xmath0-symmetric discrete nonlinear schrdinger ( dnls ) equation with gains and losses . this equation generalizes the dimer equation derived in @xcite . section [ symmetry ] describes symmetries and conserved quantities for the @xmath0-symmetric dnls equation . in particular , we show that the cross - gradient hamiltonian structure obtained in @xcite naturally appears in the asymptotic reduction of the original hamiltonian structure of newton s equations of motion for the coupled pendula . section [ breathers ] is devoted to characterization of breathers , which are time - periodic solutions localized in the chain . we show that depending on parameters of the model ( such as detuning frequency , coupling constant , driving force amplitude ) , there are three possible types of breather solutions . for the first type , breathers of small and large amplitudes are connected to each other and do not extend to symmetric synchronized oscillations of coupled pendula . in the second and third types , large - amplitude and small - amplitude breathers are connected to the symmetric synchronized oscillations but are not connected to each other . see figure [ branches ] with branches ( a ) , ( b ) , and ( c ) , where the symmetric synchronized oscillations correspond to the value @xmath4 and the breather amplitude is given by parameter @xmath5 . section [ zero - equilibrium ] contains a routine analysis of linear stability of the zero equilibrium , where the phase transition threshold to the broken @xmath0-symmetry phase is explicitly found . breathers are only studied for the parameters where the zero equilibrium is linearly stable . section [ variational ] explores the hamiltonian structure of the @xmath0-symmetric dnls equation and characterizes breathers obtained in section [ breathers ] from their energetic point of view . we show that the breathers for large value of parameter @xmath6 appear to be saddle points of the hamiltonian function between continuous spectra of positive and negative energy , similar to the standing waves in the dirac models . therefore , it is not clear from the energetic point of view if such breathers are linearly or nonlinearly stable . on the other hand , we show that the breathers for smaller values of parameter @xmath6 appear to be saddle points of the hamiltonian function with a negative continuous spectrum and finitely many ( either three or one ) positive eigenvalues .    section [ stability ] is devoted to analysis of spectral and orbital stability of breathers . for spectral stability , we use the limit of small coupling constant between the oscillators ( the same limit is also used in sections [ breathers ] and [ variational ] ) and characterize eigenvalues of the linearized operator . the main analytical results are also confirmed numerically . see figure [ spectra ] for the three types of breathers . depending on the location of the continuous spectral bands relative to the location of the isolated eigenvalues , we are able to prove nonlinear orbital stability of breathers for branches ( b ) and ( c ) . we are also able to characterize instabilities of these types of breathers that emerge depending on parameters of the model . regarding branch ( a ) , nonlinear stability analysis is not available by using the energy method . our follow - up work @xcite develops a new method of analysis to prove the long - time stability of breathers for branch ( a ) . the summary of our findings is given in the concluding section [ conclusion ] , where the main results are shown in the form of table 1 . a simple yet universal model widely used to study coupled nonlinear oscillators is the frenkel - kontorova ( fk ) model @xcite . it describes a chain of classical particles coupled to their neighbors and subjected to a periodic on - site potential . in the continuum approximation , the fk model reduces to the sine - gordon equation , which is exactly integrable . the fk model is known to describe a rich variety of important nonlinear phenomena , which find applications in solid - state physics and nonlinear science  @xcite . we consider here a two - array system of coupled pendula , where each pendulum is connected to the nearest neighbors by linear couplings . figure [ pendulapic ] shows schematically that each array of pendula is connected in the longitudinal direction by the torsional springs , whereas each pair of pendula is connected in the transverse direction by a common string . newton s equations of motion are given by @xmath7 where @xmath8 correspond to the angles of two arrays of pendula , dots denote derivatives of angles with respect to time @xmath9 , and the positive parameters @xmath10 and @xmath11 describe couplings between the two arrays in the longitudinal and transverse directions , respectively . the type of coupling between the two pendula with the angles @xmath12 and @xmath13 is referred to as the _ direct coupling _ between nonlinear oscillators ( see section 8.2 in @xcite ) . we consider oscillatory dynamics of coupled pendula under the following assumptions .    * the coupling parameters @xmath10 and @xmath11 are small . therefore , we can introduce a small parameter @xmath14 such that both @xmath10 and @xmath11 are proportional to @xmath15 . * a resonant periodic force is applied to the common strings for each pair of coupled pendula . therefore , @xmath11 is considered to be proportional to @xmath16 , where @xmath17 is selected near the unit frequency of linear pendula indicating the @xmath18 parametric resonance between the force and the pendula . mathematically , we impose the following representation for parameters @xmath10 and @xmath19 : @xmath20 where @xmath21 are @xmath14-independent parameters , whereas @xmath14 is the formal small parameter to characterize the two assumptions ( a1 ) and ( a2 ) .    in the formal limit @xmath22 , the pendula are uncoupled , and their small - amplitude oscillations can be studied with the asymptotic multi - scale expansion @xmath23 + \mu^3 x_n(t;\mu ) , \\[3pt ] y_n(t ) = \mu \left [ b_n(\mu^2 t ) e^{i \omega t } + \bar{b}_n(\mu^2 t ) e^{-i \omega t } \right ] + \mu^3 y_n(t;\mu ) , \end{array } \right . \label{expansion}\ ] ] where @xmath24 are amplitudes for nearly harmonic oscillations and @xmath25 are remainder terms . in a similar context of single - array coupled nonlinear oscillators , it is shown in @xcite how the asymptotic expansions like ( [ expansion ] ) can be justified . from the conditions that the remainder terms @xmath25 remain bounded as the system evolves , the amplitudes @xmath24 are shown to satisfy the discrete nonlinear schrdinger ( dnls ) equations , which bring together all the phenomena affecting the nearly harmonic oscillations ( such as cubic nonlinear terms , the detuning frequency , the coupling between the oscillators , and the amplitude of the parametric driving force ) . a similar derivation for a single pair of coupled pendula is reported in @xcite .    using the algorithm in @xcite and restricting the scopes of this derivation to the formal level , we write the truncated system of equations for the remainder terms : @xmath26 where @xmath27 and @xmath28 are uniquely defined . bounded solutions to the linear inhomogeneous equations ( [ linear - equation ] ) exist if and only if @xmath29 for every @xmath30 . straightforward computations show that the conditions @xmath29 are equivalent to the following evolution equations for slowly varying amplitudes @xmath24 : @xmath31 the system ( [ dnls ] ) takes the form of coupled parametrically forced dnls equations ( see @xcite for references on parametrically forced nls equations ) . there exists an invariant reduction of system ( [ dnls ] ) given by @xmath32 to the scalar parametrically forced dnls equation . the reduction ( [ synchr-1 ] ) corresponds to the symmetric synchronized oscillations of coupled pendula of the model ( [ oscillators ] ) with @xmath33 in what follows , we consider a more general class of synchronized oscillations of coupled pendula . it turns out that the model ( [ dnls ] ) can be cast to the form of the parity  time reversal ( @xmath0 ) dnls equations @xcite . using the variables @xmath34 the system of coupled dnls equations ( [ dnls ] ) is rewritten in the equivalent form @xmath35 , \\[3pt ] 2 i \dot{v}_n = \epsilon \left ( u_{n+1 } - 2 u_n + u_{n-1 } \right )   + \omega u_n - i \gamma v_n + 2\left [ \left ( |u_n|^2 + 2 |v_n|^2 \right ) u_n + \bar{u}_n v_n^2 \right ] , \end{array } \right . \label{pt - dnls}\end{aligned}\ ] ] which is the starting point for our analytical and numerical work . the invariant reduction ( [ synchr-1 ] ) for system ( [ dnls ] ) becomes @xmath36    without loss of generality , one can scale parameters @xmath37 , @xmath38 , and @xmath39 by a factor of two in order to eliminate the numerical factors in the system ( [ pt - dnls ] ) . also in the context of hard nonlinear oscillators ( e.g. in the framework of the @xmath40 theory ) , the cubic nonlinearity may have the opposite sign compared to the one in the system ( [ pt - dnls ] ) . however , given the applied context of the system of coupled pendula , we will stick to the specific form ( [ pt - dnls ] ) in further analysis . the system of coupled dnls equations ( [ pt - dnls ] ) is referred to as the @xmath0-symmetric dnls equation because the solutions remain invariant with respect to the action of the parity @xmath2 and time - reversal @xmath3 operators given by @xmath41 = \left [ \begin{array}{c } v \\ u \end{array } \right ] , \qquad \mathcal{t } \left [ \begin{array}{c } u(t ) \\ v(t ) \end{array } \right ] = \left [ \begin{array}{c } \bar{u}(-t ) \\ \bar{v}(-t ) \end{array } \right].\ ] ] the parameter @xmath39 introduces the gain  loss coefficient in each pair of coupled oscillators due to the resonant periodic force . in the absence of all other effects , the @xmath39-term of the first equation of system ( [ pt - dnls ] ) induces the exponential growth of amplitude @xmath42 , whereas the @xmath39-term of the second equation induces the exponential decay of amplitude @xmath43 , if @xmath44 . the system ( [ pt - dnls ] ) truncated at a single site ( say @xmath45 ) is called the @xmath0-symmetric dimer . in the work of barashenkov _ et al . _ @xcite , it was shown that all @xmath0-symmetric dimers with physically relevant cubic nonlinearities represent hamiltonian systems in appropriately introduced canonical variables . however , the @xmath0-symmetric dnls equation on a lattice does not typically have a hamiltonian form if @xmath46 .    nevertheless , the particular nonlinear functions arising in the system ( [ pt - dnls ] ) correspond to the @xmath0-symmetric dimers with a cross  gradient hamiltonian structure @xcite , where variables @xmath47 are canonically conjugate . as a result , the system ( [ pt - dnls ] ) on the chain @xmath48 has additional conserved quantities . this fact looked like a mystery in the recent works @xcite . here we clarify the mystery in the context of the derivation of the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) from the original system ( [ oscillators ] ) . indeed , the system ( [ oscillators ] ) of classical newton particles has a standard hamiltonian structure with the energy function @xmath49 since the periodic movement of common strings for each pair of pendula result in the time - periodic coefficient @xmath19 , the energy @xmath50 is a periodic function of time @xmath9 . in addition , no other conserved quantities such as momenta exist typically in lattice differential systems such as the system ( [ oscillators ] ) due to broken continuous translational symmetry . after the system ( [ oscillators ] ) is reduced to the coupled dnls equations ( [ dnls ] ) with the asymptotic expansion ( [ expansion ] ) , we can write the evolution problem ( [ dnls ] ) in the hamiltonian form with the standard straight - gradient symplectic structure @xmath51 where the time variable @xmath9 stands now for the slow time @xmath52 and the energy function is @xmath53 the energy function @xmath54 is conserved in the time evolution of the hamiltonian system ( [ hamiltonian - dnls ] ) . in addition , there exists another conserved quantity @xmath55 which is related to the gauge symmetry @xmath56 with @xmath57 for solutions to the system ( [ dnls ] ) . when the transformation of variables ( [ variables ] ) is used , the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) is cast to the hamiltonian form with the cross - gradient symplectic structure @xmath58 where the energy function is @xmath59 the gauge - related function is written in the form @xmath60 the functions @xmath61 and @xmath62 are conserved in the time evolution of the system ( [ pt - dnls ] ) . these functions follow from ( [ energy - dnls ] ) and ( [ charge - dnls ] ) after the transformation ( [ variables ] ) is used . thus , the cross - gradient hamiltonian structure of the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) is inherited from the hamiltonian structure of the coupled oscillator model ( [ oscillators ] ) . we characterize the existence of breathers supported by the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) . in particular , breather solutions are continued for small values of coupling constant @xmath38 from solutions of the dimer equation arising at a single site , say the central site at @xmath45 . we shall work in a sequence space @xmath63 of square integrable complex - valued sequences . time - periodic solutions to the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) are given in the form @xcite : @xmath64 where the parameter @xmath6 is considered to be real , the factor @xmath65 is introduced for convenience , and the sequence @xmath66 is time - independent . the breather ( [ stationary ] ) is a localized mode if @xmath67 , which implies that @xmath68 as @xmath69 . the breather ( [ stationary ] ) is considered to be @xmath0-symmetric with respect to the operators in ( [ pt - symmetry ] ) if @xmath70 .    the reduction ( [ reduction ] ) for symmetric synchronized oscillations is satisfied if @xmath71 the time - periodic breathers ( [ stationary ] ) with @xmath72 generalize the class of symmetric synchronized oscillations ( [ reduction - e ] ) . the time - independent sequence @xmath67 can be found from the stationary @xmath0-symmetric dnls equation : @xmath73 , \\[3pt ] e v_n = \epsilon \left ( u_{n+1 } - 2 u_n + u_{n-1 } \right ) + \omega u_n - i \gamma v_n + 2\left [ \left ( |u_n|^2 + 2 |v_n|^2 \right ) u_n + \bar{u}_n v_n^2 \right ] . \end{array } \right . \label{eq : state}\end{aligned}\ ] ] the @xmath0-symmetric breathers with @xmath70 satisfy the following scalar difference equation @xmath74 note that the reduction ( [ reduction - e ] ) is compatible with equation ( [ eq : statpt ] ) in the sense that if @xmath4 and @xmath75 , then @xmath76 satisfies a real - valued difference equation . let us set @xmath77 for now and consider solutions to the dimer equation at the central site @xmath45 : @xmath78 the parameters @xmath39 and @xmath37 are considered to be fixed , and the breather parameter @xmath6 is thought to parameterize continuous branches of solutions to the nonlinear algebraic equation ( [ dimer - states ] ) . the solution branches depicted on figure [ branches ] are given in the following lemma .    ) . ] [ lemma - dimer ] assume @xmath46 . the algebraic equation ( [ dimer - states ] ) admits the following solutions depending on @xmath39 and @xmath37 :    * @xmath79 - two symmetric unbounded branches exist for @xmath80 , * @xmath81 - an unbounded branch exists for every @xmath82 , * @xmath83 - a bounded branch exists for @xmath84 ,    where @xmath85 .    substituting the decomposition @xmath86 with @xmath87 and @xmath88 into the algebraic equation ( [ dimer - states ] ) , we obtain @xmath89 excluding @xmath90 by using the fundamental trigonometric identity , we obtain the explicit parametrization of the solutions to the algebraic equation ( [ dimer - states ] ) by the amplitude parameter @xmath5 : @xmath91.\ ] ] the zero - amplitude limit @xmath92 is reached if @xmath93 , in which case @xmath94 , where @xmath85 . if @xmath95 , the solution branches ( if they exist ) are bounded away from the zero solution . now we analyze the three cases of parameters @xmath39 and @xmath37 formulated in the lemma .    * if @xmath79 , then the parametrization ( [ branch - e - a ] ) yields a monotonically increasing map @xmath96 because @xmath97 > 0.\ ] ] in the two asymptotic limits , we obtain from ( [ branch - e - a ] ) : @xmath98 see figure  [ branches](a ) . * if @xmath81 , the parametrization ( [ branch - e - a ] ) yields a monotonically increasing map @xmath99 , where @xmath100 indeed , we note that @xmath101 and @xmath102 so that the derivative in ( [ positivity - a ] ) is positive for every @xmath103 . we have @xmath104 see figure  [ branches](b ) . * if @xmath83 , then the parametrization ( [ branch - e - a ] ) yields a monotonically decreasing map @xmath105 , where @xmath106 in ( [ a - minus ] ) , the first choice is made if @xmath107 and the second choice is made if @xmath108 . both choices are the same if @xmath109 . we note that @xmath110 , therefore , the derivative ( [ positivity - a ] ) needs to be rewritten in the form @xmath111 < 0,\ ] ] where @xmath112 for both @xmath107 and @xmath113 . in the two asymptotic limits , we obtain from ( [ branch - e - a ] ) : @xmath114 see figure  [ branches](c ) . note that branches ( b ) and ( c ) coexist for @xmath83 . the reduction ( [ reduction - e ] ) corresponds to the choice : @xmath115 if @xmath44 , this choice corresponds to @xmath116 for @xmath117 , that is , the point @xmath4 on branch ( c ) . if @xmath118 , it corresponds to @xmath119 for any @xmath81 , that is , the point @xmath4 on branch ( b ) . every solution of lemma [ lemma - dimer ] can be extended to a breather on the chain @xmath48 which satisfies the spatial symmetry condition in addition to the @xmath0 symmetry : @xmath120 in order to prove the existence of the symmetric breather solution to the difference equation ( [ eq : statpt ] ) , we use the following implicit function theorem . + * implicit function theorem * ( theorem 4.e in @xcite ) . _ let @xmath121 and @xmath122 be banach spaces and let @xmath123 be a @xmath124 map on an open neighborhood of the point @xmath125 . assume that @xmath126 and that @xmath127 there are @xmath128 and @xmath129 such that for each @xmath130 with @xmath131 there exists a unique solution @xmath132 of the nonlinear equation @xmath133 with @xmath134 . moreover , the map @xmath135 is @xmath124 near @xmath136 . _    with two applications of the implicit function theorem , we prove the following main result of this section . [ theorem - existence - soliton ] fix @xmath46 , @xmath137 , and @xmath138 , where @xmath139 if @xmath93 . there exists @xmath140 sufficiently small and @xmath141 such that for every @xmath142 , there exists a unique solution @xmath143 to the difference equation ( [ eq : statpt ] ) satisfying the symmetry ( [ spatial - symmetry ] ) and the bound @xmath144 where @xmath5 and @xmath90 are defined in lemma [ lemma - dimer ] . moreover , the solution @xmath145 is smooth in @xmath38 .    in the first application of the implicit function theorem , we consider the following system of algebraic equations @xmath146 where @xmath147 is given , in addition to parameters @xmath39 , @xmath37 , and @xmath6 . let @xmath148 , @xmath149 , @xmath150 , @xmath151 , and @xmath152 . then , we have @xmath153 and the jacobian operator @xmath154 is given by identical copies of the matrix @xmath155,\ ] ] with the eigenvalues @xmath156 . by the assumption of the lemma , @xmath157 , so that the jacobian operator @xmath154 is one - to - one and onto . by the implicit function theorem , for every @xmath147 and every @xmath158 sufficiently small , there exists a unique small solution @xmath159 of the system ( [ fixed - point-1 ] ) such that @xmath160 where the positive constant @xmath161 is independent from @xmath38 and @xmath162 . thanks to the symmetry of the difference equation ( [ eq : statpt ] ) , we find that @xmath163 , @xmath164 satisfy the same system ( [ fixed - point-1 ] ) with @xmath165 , with the same unique solution .    in the second application of the implicit function theorem , we consider the following algebraic equation @xmath166 where @xmath167 depends on @xmath162 , @xmath39 , @xmath37 , and @xmath6 , satisfies the bound ( [ eq : u ] ) , and is uniquely defined by the previous result . let @xmath168 , @xmath169 , @xmath150 , @xmath151 , and @xmath170 . then , we have @xmath171 , where @xmath5 and @xmath90 are defined in lemma [ lemma - dimer ] . the jacobian operator @xmath172 is given by the matrix @xmath173 \biggr|_{u_0 = a e^{i \theta } } \\ \label{matrix - jacobian } & = & \left [ \begin{array}{cc } e - i \gamma - \frac{12 e a^2}{\omega + 8 a^2 } & -\omega - 12 a^2 \\ -\omega - 12 a^2 & e + i \gamma - \frac{12 e a^2}{\omega + 8 a^2 } \end{array } \right].\end{aligned}\ ] ] we show in lemma [ proposition - unbounded ] below that the matrix given by ( [ matrix - jacobian ] ) is invertible under the conditions @xmath46 and @xmath174 . by the implicit function theorem , for every @xmath158 sufficiently small , there exists a unique solution @xmath147 to the algebraic equation ( [ fixed - point-2 ] ) near @xmath175 such that @xmath176 where the positive constant @xmath177 is independent from @xmath38 . the bound ( [ bound - soliton ] ) holds thanks to the bounds ( [ eq : u ] ) and ( [ eq : u0 ] ) . since both equations ( [ fixed - point-1 ] ) and ( [ fixed - point-2 ] ) are smooth in @xmath38 , the solution @xmath145 is smooth in @xmath38 .    in the following result , we show that the matrix given by ( [ matrix - jacobian ] ) is invertible for every branch of lemma [ lemma - dimer ] with an exception of a single point @xmath4 on branch ( c ) for @xmath178 . [ proposition - unbounded ] with the exception of the point @xmath4 on branch ( c ) of lemma [ lemma - dimer ] for @xmath178 , the matrix given by ( [ matrix - jacobian ] ) is invertible for every @xmath46 . the matrix given by ( [ matrix - jacobian ] ) has zero eigenvalue if and only if its determinant is zero , which happens at @xmath179 eliminating @xmath180 by using parametrization ( [ branch - e - a ] ) and simplifying the algebraic equation for nonzero @xmath181 , we reduce it to the form @xmath182 we now check if this constraint can be satisfied for the three branches of lemma [ lemma - dimer ] .    * if @xmath79 , the constraint ( [ constraint - parameters ] ) is not satisfied because the left - hand side @xmath183 exceeds the right - hand side @xmath184 . * if @xmath81 and @xmath103 , where @xmath185 is given by ( [ a - plus ] ) , the constraint ( [ constraint - parameters ] ) is not satisfied because the left - hand side @xmath186 exceeds the left - hand side @xmath184 both for @xmath187 and for @xmath188 . * if @xmath83 and @xmath189 , where @xmath190 is given by ( [ a - minus ] ) , the constraint ( [ constraint - parameters ] ) is not satisfied because the left - hand side is estimated by @xmath191 in the first case , we have @xmath192 , so that the left - hand side is strictly smaller than @xmath193 . in the second case , we have @xmath194 , so that the left - hand side is also strictly smaller than @xmath193 . only if @xmath109 , the constraint ( [ constraint - parameters ] ) is satisfied at @xmath4 , when @xmath195 and @xmath196    hence , the matrix ( [ matrix - jacobian ] ) is invertible for all parameter values with one exceptional case .    in the asymptotic limit @xmath197 as @xmath198 , see lemma [ lemma - dimer ] , the matrix ( [ matrix - jacobian ] ) is expanded asymptotically as @xmath199 + \mathcal{o}(1 ) \quad \mbox{\rm as } \quad |e| \to \infty , \label{matrix - jacobian - asymptotics}\end{aligned}\ ] ] with the two eigenvalues @xmath200 and @xmath201 . thus , the matrix given by ( [ matrix - jacobian - asymptotics ] ) is invertible for every branch extending to sufficiently large values of @xmath6 . here we discuss the linear stability of the zero equilibrium in the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) . the following proposition yields a simple result . the zero equilibrium of the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) is linearly stable if @xmath202 , where @xmath203 . [ prop - zero - equilibrium ]    truncating the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) at the linear terms and using the fourier transform @xmath204 we obtain the linear homogeneous system @xmath205 = \left [ \begin{array}{c } 0 \\ 0 \end{array } \right ] , \quad \mbox{\rm where } \quad \hat{d}(k ) : = \left [ \begin{matrix } -2 \omega(k ) - i \gamma & -\omega + 4 \epsilon \sin^2(k/2 ) \\ -\omega + 4 \epsilon \sin^2(k/2 ) & -2 \omega(k ) + i \gamma \end{matrix } \right].\ ] ] the determinant of @xmath206 is zero if and only if @xmath207 is found from the quadratic equation @xmath208 for any @xmath202 , where @xmath209 is given by ( [ gamma-0 ] ) , the two branches @xmath210 found from the quadratic equation ( [ dispersion - relation ] ) are real - valued and non - degenerate for every @xmath211 $ ] . therefore , the zero equilibrium is linearly stable .    on the other hand , for any @xmath212 , the values of @xmath207 are purely imaginary either near @xmath213 if @xmath214 or near @xmath215 if @xmath188 . therefore , the zero equilibrium is linearly unstable .    the value @xmath209 given by ( [ gamma-0 ] ) represents the phase transition threshold and the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) is said to have broken @xmath0-symmetry for @xmath212 . if @xmath77 , the zero equilibrium is only linearly stable for @xmath216 . since the localized breathers can not be stable when the zero background is unstable , we shall study stability of breathers only for the case when @xmath216 , that is , in the regime of unbroken @xmath0-symmetry . it follows from theorem [ theorem - existence - soliton ] that each interior point on the solution branches shown on figure [ branches ] generates a fundamental breather of the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) . we shall now characterize these breathers as relative equilibria of the energy function . thanks to the cross - gradient symplectic structure ( [ hamiltonian - pt - dnls ] ) , the stationary @xmath0-symmetric dnls equation can be written in the gradient form @xmath217 keeping in mind the additional conserved quantity @xmath62 given by ( [ charge - pt - dnls ] ) , we conclude that the stationary solution @xmath66 is a critical point of the combined energy function given by @xmath218    if we want to apply the lyapunov method in order to study nonlinear stability of stationary solutions in hamiltonian systems , we shall investigate convexity of the second variation of the combined energy functional @xmath219 at @xmath66 . using the expansion @xmath220 , @xmath221 and introducing extended variables @xmath222 and @xmath223 with the blocks @xmath224 we can expand the smooth function @xmath219 up to the quadratic terms in @xmath223 : @xmath225 where @xmath226 is the self - adjoint ( hessian ) operator defined on @xmath63 and the scalar product was used in the following form : @xmath227    using ( [ energy - pt - dnls ] ) and ( [ charge - pt - dnls ] ) , the hessian operator can be computed explicitly as follows @xmath228 where blocks of @xmath229 at each lattice node @xmath30 are given by @xmath230\end{aligned}\]]and @xmath231 is the discrete laplacian operator applied to blocks of @xmath223 at each lattice node @xmath30 : @xmath232 in the expression for @xmath233 , we have used the @xmath0-symmetry condition @xmath70 for the given stationary solution @xmath66 . we study convexity of the combined energy functional @xmath219 at @xmath66 . since the zero equilibrium is linearly stable only for @xmath216 ( if @xmath77 ) , we only consider breathers of theorem [ theorem - existence - soliton ] for @xmath216 . in order to study eigenvalues of @xmath226 for small values of @xmath38 , we use the following perturbation theory . + * perturbation theory for linear operators * ( theorem vii.1.7 in @xcite ) . _ let @xmath234 be a family of bounded operators from banach space @xmath235 to itself , which depends analytically on the small parameter @xmath38 . if the spectrum of @xmath236 is separated into two parts , the subspaces of @xmath235 corresponding to the separated parts also depend analytically on @xmath38 . in particular , the spectrum of @xmath234 is separated into two parts for any @xmath158 sufficiently small . _    with an application of the perturbation theory for linear operators , we prove the following main result of this section .    [ theorem - critical - point ] fix @xmath46 , @xmath37 , and @xmath6 along branches of the @xmath0-symmetric breathers @xmath66 given by theorem [ theorem - existence - soliton ] such that @xmath93 and @xmath138 , where @xmath139 . for every @xmath237 sufficiently small , the operator @xmath226 admits a one - dimensional kernel in @xmath63 spanned by the eigenvector @xmath238 due to the gauge invariance , where the blocks of the eigenvector are given by @xmath239 in addition ,    * if @xmath240 , the spectrum of @xmath226 in @xmath63 includes infinite - dimensional positive and negative parts . * if @xmath241 and @xmath83 , the spectrum of @xmath242 in @xmath63 includes an infinite - dimensional negative part and either three or one simple positive eigenvalues for branches ( b ) and ( c ) of lemma [ lemma - dimer ] respectively .    if @xmath77 , the breather solution of theorem [ theorem - existence - soliton ] is given by @xmath243 for every @xmath244 and @xmath86 , where @xmath5 and @xmath90 are defined by lemma [ lemma - dimer ] . in this case , the linear operator @xmath245 decouples into @xmath246-by-@xmath246 blocks for each lattice node @xmath30 .    for @xmath45 , the @xmath246-by-@xmath246 block of the linear operator @xmath229 is given by @xmath247.\]]using relations ( [ parameterization - a ] ) and ( [ branch - e - a ] ) , as well as symbolic computations with maple , we found that the @xmath246-by-@xmath246 matrix block @xmath248 admits a simple zero eigenvalue and three nonzero eigenvalues @xmath249 , @xmath250 , and @xmath251 given by @xmath252 for each branch of lemma [ lemma - dimer ] with @xmath46 and @xmath138 , we have @xmath253 , so that @xmath254 . furthermore , either @xmath255 or @xmath256 if and only if @xmath257 expanding this equation for nonzero @xmath5 yields constraint ( [ constraint - parameters ] ) . with the exception of a single point @xmath4 at @xmath178 , we showed in lemma [ proposition - unbounded ] that the constraint ( [ constraint - parameters ] ) does not hold for any of the branches of lemma [ lemma - dimer ] . therefore , @xmath258 and @xmath259 along each branch of lemma [ lemma - dimer ] and the signs of @xmath249 , @xmath250 , and @xmath251 for each branch of lemma [ lemma - dimer ] can be obtained in the limit @xmath198 for branches ( a ) and ( b ) or @xmath260 for branch ( c ) . by means of these asymptotic computations as @xmath198 or @xmath260 , we obtain the following results for the three branches shown on figure [ branches ] :    * @xmath261 . * @xmath261 . * @xmath262 , @xmath263 , and @xmath264 .    for @xmath265 , the @xmath246-by-@xmath246 block of the linear operator @xmath229 is given by @xmath266.\ ] ] each block has two double eigenvalues @xmath267 and @xmath268 given by @xmath269 since there are infinitely many nodes with @xmath244 , the points @xmath267 and @xmath268 have infinite multiplicity in the spectrum of the linear operator @xmath229 . furthermore , we can sort up the signs of @xmath267 and @xmath268 for each point on the three branches shown on figure [ branches ] :    * if @xmath270 , then @xmath271 and @xmath272 . * if @xmath241 and @xmath83 , then @xmath273 .    by using the perturbation theory for linear operators , we argue as follows :    * since @xmath226 is hermitian on @xmath63 , its spectrum is a subset of the real line for every @xmath158 . * the zero eigenvalue persists with respect to @xmath158 at zero because the eigenvector ( [ kernel ] ) belongs to the kernel of @xmath226 due to the gauge invariance for every @xmath158 . * the other eigenvalues of @xmath229 are isolated away from zero . the spectrum of @xmath226 is continuous with respect to @xmath38 and includes infinite - dimensional parts near points @xmath267 and @xmath268 for small @xmath237 ( which may include continuous spectrum and isolated eigenvalues ) as well as simple eigenvalues near @xmath274 ( if @xmath274 are different from @xmath275 ) . the statement of the theorem follows from the perturbation theory and the count of signs of @xmath274 and @xmath275 above .    in the asymptotic limit @xmath197 as @xmath198 , we can sort out eigenvalues of @xmath226 asymptotically as : @xmath276 where the remainder terms are @xmath277 as @xmath278 . the values @xmath249 , @xmath251 , and @xmath267 are close to each other as @xmath278 . it follows from theorem [ theorem - critical - point ] that for @xmath240 , the breather @xmath66 is a saddle point of the energy functional @xmath219 with infinite - dimensional positive and negative invariant subspaces of the hessian operator @xmath226 . this is very similar to the hamiltonian systems of the dirac type , where stationary states are located in the gap between the positive and negative continuous spectrum . this property holds for points 1 and 3 on branches ( a ) and ( b ) shown on figure [ branches ] . no branches other than @xmath240 exist for @xmath79 . on the other hand , points 2 and 4 on branches ( b ) and ( c ) shown on figure [ branches ] satisfy @xmath241 and @xmath83 . the breather @xmath66 is a saddle point of @xmath219 for these points and it only has three ( one ) directions of positive energy in space @xmath63 for point 2 ( point 4 ) . spectral stability of breathers can be studied for small values of coupling constant @xmath38 by using the perturbation theory @xcite . first , we linearize the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) at the breather ( [ stationary ] ) by using the expansion @xmath279 , \quad v(t ) = e^{-\frac{1}{2 } i et } \left [ v + { \bf v}(t ) \right],\ ] ] where @xmath280 is a small perturbation satisfying the linearized equations @xmath281 , \\ 2 i \dot{{\bf v}}_n + e { \bf v}_n = \epsilon \left ( { \bf u}_{n+1 } - 2 { \bf u}_n + { \bf u}_{n-1 } \right )   + \omega { \bf u}_n - i \gamma { \bf v}_n \\ \phantom{text } + 2 \left [ 2 \left ( |u_n|^2 + |v_n|^2 \right ) { \bf u}_n + ( u_n^2 + v_n^2 ) \bar{{\bf u}}_n + 2 ( \bar{u}_n v_n + u_n \bar{v}_n ) { \bf v}_n + 2 u_n v_n \bar{{\bf v}}_n \right ] . \end{array } \right . \label{pt - dnls - lin}\end{aligned}\ ] ]    the spectral stability problem arises from the linearized equations ( [ pt - dnls - lin ] ) after the separation of variables : @xmath282 where @xmath283 is the eigenvector corresponding to the spectral parameter @xmath284 . note that @xmath285 and @xmath286 are no longer complex conjugate to each other if @xmath284 has a nonzero imaginary part . the spectral problem can be written in the explicit form @xmath287 , \\ ( e - i \lambda + i \gamma ) \psi_n - \omega \nu_n = \epsilon \left ( \nu_{n+1 } - 2 \nu_n + \nu_{n-1 } \right ) \\ \phantom{texttext } + 2 \left [ 2 |u_n|^2 ( \varphi_n + 2\nu_n ) + ( u_n^2 + \bar{u}_n^2 ) ( 2 \psi_n + \chi_n ) \right ] , \\ ( e + i \lambda + i \gamma ) \chi_n - \omega \varphi_n = \epsilon \left ( \varphi_{n+1 } - 2 \varphi_n + \varphi_{n-1 } \right ) \\ \phantom{texttext } + 2 \left [ 2 |u_n|^2 ( 2 \varphi_n + \nu_n ) + ( u_n^2 + \bar{u}_n^2 ) ( \psi_n + 2 \chi_n ) \right ] , \\ ( e - i \lambda - i \gamma ) \nu_n - \omega \psi_n = \epsilon \left ( \psi_{n+1 } - 2 \psi_n + \psi_{n-1 } \right ) \\ \phantom{texttext } + 2 \left [ 2 |u_n|^2 ( 2 \psi_n + \chi_n ) + ( u_n^2 + \bar{u}_n^2 ) ( \varphi_n + 2 \nu_n ) \right ] , \end{array } \right.\ ] ] where we have used the condition @xmath70 for the @xmath0-symmetric breathers . recalling definition of the hessian operator @xmath226 in ( [ pt - dnls - var ] ) , we can rewrite the spectral problem ( [ spectral - problem ] ) in the hamiltonian form : @xmath288 where @xmath289 is a symmetric matrix with the blocks at each lattice node @xmath30 given by @xmath290.\ ] ] we note the hamiltonian symmetry of the eigenvalues of the spectral problem ( [ eq : eigenval ] ) . [ prop - ham - sym ] eigenvalues of the spectral problem ( [ eq : eigenval ] ) occur either as real or imaginary pairs or as quadruplets in the complex plane . assume that @xmath291 is an eigenvalue of the spectral problem ( [ eq : eigenval ] ) with the eigenvector @xmath292 . then , @xmath293 is an eigenvalue of the same problem with the eigenvector @xmath294 , whereas @xmath295 is also an eigenvalue with the eigenvector @xmath296 .    if @xmath83 and @xmath297 ( points 2 and 4 shown on figure [ branches ] ) , theorem [ theorem - critical - point ] implies that the self - adjoint operator @xmath226 in @xmath63 is negative - definite with the exception of either three ( point 2 ) or one ( point 4 ) simple positive eigenvalues . in this case , we can apply the following hamilton  krein index theorem in order to characterize the spectrum of @xmath298 . + * hamilton  krein index theorem * ( theorem 3.3 in @xcite ) . _ let @xmath299 be a self - adjoint operator in @xmath300 with finitely many negative eigenvalues @xmath301 , a simple zero eigenvalue with eigenfunction @xmath302 , and the rest of its spectrum is bounded from below by a positive number . let @xmath303 be a bounded invertible skew - symmetric operator in @xmath300 . let @xmath304 be a number of positive real eigenvalues of @xmath305 , @xmath306 be a number of quadruplets @xmath307 that are neither in @xmath308 nor in @xmath309 , and @xmath310 be a number of purely imaginary pairs of eigenvalues of @xmath305 whose invariant subspaces lie in the negative subspace of @xmath299 . let @xmath311 be finite and nonzero . then , @xmath312 _    [ theorem - index ] fix @xmath46 , @xmath83 , and @xmath313 , where @xmath139 . for every @xmath237 sufficiently small , @xmath314 for branch ( b ) of lemma [ lemma - dimer ] and @xmath315 for branch ( c ) of lemma [ lemma - dimer ] with @xmath316 . for branch ( c ) with @xmath317 , there exists a value @xmath318 such that @xmath319 for @xmath320 and @xmath315 for @xmath321 . if @xmath46 , @xmath83 , @xmath241 , and @xmath237 is sufficiently small , theorem [ theorem - critical - point ] implies that the spectrum of @xmath226 in @xmath63 has finitely many positive eigenvalues and a simple zero eigenvalue with eigenvector @xmath238 . therefore , the hamilton  krein index theorem is applied in @xmath63 for @xmath322 , @xmath323 , and @xmath324 . we shall verify that @xmath325 where @xmath238 is given by ( [ kernel ] ) and @xmath326 denotes derivative of @xmath222 with respect to parameter @xmath6 . the first equation @xmath327 follows by theorem [ theorem - critical - point ] . by differentiating equations ( [ eq : state ] ) in @xmath6 , we obtain @xmath328 for every @xmath6 , for which the solution @xmath222 is differentiable in @xmath6 . for @xmath77 , the limiting solution of lemma [ lemma - dimer ] is differentiable in @xmath6 for every @xmath72 and @xmath138 . due to smoothness of the continuation in @xmath38 by theorem [ theorem - critical - point ] , this property holds for every @xmath237 sufficiently small .    by using ( [ generalized - kernel - lemma ] ) with @xmath329 , we obtain @xmath330 where we have used the definition of @xmath62 in ( [ charge - pt - dnls ] ) . we compute the slope condition at @xmath77 : @xmath331,\ ] ] where relations ( [ parameterization - a ] ) and ( [ branch - e - a ] ) have been used .    for branch ( b ) of lemma [ lemma - dimer ] with @xmath83 , we have @xmath332 and @xmath333 , so that @xmath334 . by continuity , @xmath335 remains strictly positive for small @xmath237 . thus , @xmath336 and @xmath314 by the hamilton  krein index theorem .    for branch ( c ) of lemma [ lemma - dimer ] with @xmath83 , we have @xmath337 and @xmath338 . therefore , we only need to inspect the sign of the expression @xmath339 . if @xmath316 , then for every @xmath340 , we have @xmath341 therefore , @xmath336 and @xmath315 by the hamilton  krein index theorem .    on the other hand , if @xmath342 , we have @xmath343 at @xmath92 ( @xmath344 ) and @xmath345 at @xmath116 ( @xmath4 ) . since the dependence of @xmath5 versus @xmath6 is monotonic , there exists a value @xmath318 such that @xmath319 for @xmath320 and @xmath315 for @xmath321 .    if @xmath315 and @xmath346 , orbital stability of a critical point of @xmath219 in space @xmath63 can be proved from the hamilton  krein theorem ( see @xcite and references therein ) . orbital stability of breathers is understood in the following sense . we say that the breather solution ( [ stationary ] ) is orbitally stable in @xmath63 if for every @xmath347 sufficiently small , there exists @xmath348 such that if @xmath349 satisfies @xmath350 , then the unique global solution @xmath351 , @xmath352 to the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) satisfies the bound @xmath353 [ def - stability ]    the definition of instability of breathers is given by negating definition [ def - stability ] . the following result gives orbital stability or instability for branch ( c ) shown on figure [ branches ] . [ theorem - bounded ] fix @xmath46 , @xmath83 , and @xmath313 . for every @xmath237 sufficiently small , the breather @xmath66 for branch ( c ) of lemma [ lemma - dimer ] is orbitally stable in @xmath63 if @xmath316 . for every @xmath317 , there exists a value @xmath318 such that the breather @xmath66 is orbitally stable in @xmath63 if @xmath321 and unstable if @xmath320 . the theorem is a corollary of lemma [ theorem - index ] for branch ( c ) of lemma [ lemma - dimer ] and the orbital stability theory from @xcite . orbital stability of breathers for branches ( a ) and ( b ) of lemma [ lemma - dimer ] does not follow from the standard theory because @xmath354 for @xmath240 and @xmath355 for branch ( b ) with @xmath241 . nevertheless , by using smallness of parameter @xmath38 and the construction of the breather @xmath66 in theorem [ theorem - existence - soliton ] , spectral stability of breathers can be considered directly . spectral stability and instability of breathers is understood in the following sense . [ def - spectral ] we say that the breather solution ( [ stationary ] ) is spectrally stable if @xmath356 for every bounded solution of the spectral problem ( [ eq : eigenval ] ) . on the other hand , if the spectral problem ( [ eq : eigenval ] ) admits an eigenvalue @xmath357 with an eigenvector in @xmath358 , we say that the breather solution ( [ stationary ] ) is spectrally unstable . the following theorem gives spectral stability of breathers for branches ( a ) and ( b ) shown on figure [ branches ] . [ theorem - stability ] fix @xmath46 , @xmath93 , and @xmath6 along branches ( a ) and ( b ) of lemma [ lemma - dimer ] with @xmath72 and @xmath138 . for every @xmath237 sufficiently small , the spectral problem ( [ eq : eigenval ] ) admits a double zero eigenvalue with the generalized eigenvectors @xmath359 where the eigenvector @xmath238 is given by ( [ kernel ] ) and the generalized eigenvector @xmath326 denotes derivative of @xmath222 with respect to parameter @xmath6 . for every @xmath6 such that the following non - degeneracy condition is satisfied , @xmath360 the breather @xmath66 is spectrally stable .    if @xmath77 , the breather solution of theorem [ theorem - existence - soliton ] is given by @xmath243 for every @xmath244 and @xmath86 , where @xmath5 and @xmath90 are defined by lemma [ lemma - dimer ] . in this case , the spectral problem ( [ spectral - problem ] ) decouples into @xmath246-by-@xmath246 blocks for each lattice node @xmath30 . recall that @xmath361 at @xmath77 . for @xmath45 , eigenvalues @xmath284 are determined by the @xmath246-by-@xmath246 matrix block @xmath362 . using relations ( [ parameterization - a ] ) and ( [ branch - e - a ] ) , as well as symbolic computations with maple , we found that the @xmath246-by-@xmath246 matrix block @xmath362 has a double zero eigenvalue and a pair of simple eigenvalues at @xmath363 , where @xmath364    for @xmath265 , eigenvalues @xmath284 are determined by the @xmath246-by-@xmath246 matrix block @xmath365 , where @xmath233 is given by ( [ block - l - n ] ) . if @xmath216 , @xmath72 , and @xmath138 , where @xmath85 , each block has four simple eigenvalues @xmath366 and @xmath367 , where @xmath368 so that @xmath369 . since there are infinitely many nodes with @xmath244 , the four eigenvalues are semi - simple and have infinite multiplicity . if @xmath237 is sufficiently small , we use perturbation theory for linear operators from section [ variational ] .    * the double zero eigenvalue persists with respect to @xmath158 at zero because of the gauge invariance of the breather @xmath66 ( with respect to rotation of the complex phase ) . indeed , @xmath327 follows from the result of theorem [ theorem - critical - point ] . the generalized eigenvector is defined by equation @xmath370 , which is equivalent to equation @xmath371 . differentiating equations ( [ eq : state ] ) in @xmath6 , we obtain @xmath372 . since @xmath373 = 1 $ ] and @xmath374 the second generalized eigenvector @xmath375 exists as a solution of equation @xmath376 if and only if @xmath377 . it follows from the explicit computation ( [ slope - explicit ] ) that if @xmath77 , then @xmath378 for every @xmath6 along branches ( a ) and ( b ) of lemma [ lemma - dimer ] . by continuity , @xmath378 for small @xmath237 . therefore , the zero eigenvalue of the operator @xmath379 is exactly double for small @xmath237 . * using the same computation ( [ slope - explicit ] ) , it is clear that @xmath380 for every @xmath6 along branches ( a ) and ( b ) of lemma [ lemma - dimer ] . assume that @xmath381 and @xmath382 , which is expressed by the non - degeneracy condition ( [ assumption - non - degeneracy ] ) . then , the pair @xmath383 is isolated from the rest of the spectrum of the operator @xmath379 at @xmath77 . since the eigenvalues @xmath363 are simple and purely imaginary , they persist on the imaginary axis for @xmath158 because they can not leave the imaginary axis by the hamiltonian symmetry of proposition [ prop - ham - sym ] . * if @xmath216 , @xmath72 , and @xmath138 , the semi - simple eigenvalues @xmath366 and @xmath367 of infinite multiplicity are nonzero and located at the imaginary axis at different points for @xmath77 . they persist on the imaginary axis for @xmath158 according to the following perturbation argument . first , for the central site @xmath45 , the spectral problem ( [ spectral - problem ] ) can be written in the following abstract form @xmath384 where @xmath385 denotes a continuation of @xmath248 in @xmath38 . thanks to the non - degeneracy condition ( [ assumption - non - degeneracy ] ) as well as the condition @xmath157 , the matrix @xmath386 is invertible . by continuity , the matrix @xmath387 is invertible for every @xmath38 and @xmath284 near @xmath77 and @xmath388 . therefore , there is a unique @xmath389 given by @xmath390 which satisfies @xmath391 near @xmath77 and @xmath388 , where @xmath10 is a positive @xmath38- and @xmath284-independent constant . next , for either @xmath164 or @xmath165 , the spectral problem ( [ spectral - problem ] ) can be represented in the form @xmath392 where @xmath393 denotes a continuation of @xmath233 given by ( [ block - l - n ] ) in @xmath38 , whereas the operator @xmath231 is applied with zero end - point condition at @xmath45 . we have @xmath394 and @xmath395 near @xmath77 and @xmath388 . therefore , up to the first order of the perturbation theory , the spectral parameter @xmath284 near @xmath396 is defined from the truncated eigenvalue problem @xmath397 which is solved with the discrete fourier transform ( [ discrete - fourier ] ) . in order to satisfy the dirichlet end - point condition at @xmath45 , the sine  fourier transform must be used , which does not affect the characteristic equation for the purely continuous spectrum of the spectral problem ( [ truncated - problem ] ) . by means of routine computations , we obtain the characteristic equation in the form , see also equation ( [ dispersion - relation ] ) , @xmath398 where @xmath211 $ ] is the parameter of the discrete fourier transform ( [ discrete - fourier ] ) . solving the characteristic equation ( [ continuous - spectrum ] ) , we obtain four branches of the continuous spectrum @xmath399 where the two sign choices are independent from each other . if @xmath93 is fixed and @xmath237 is small , the four branches of the continuous spectrum are located on the imaginary axis near the points @xmath366 and @xmath367 given by ( [ lambda - plus - minus - eig ] ) . + in addition to the continuous spectrum given by ( [ continuous - spectrum ] ) , there may exist isolated eigenvalues near @xmath366 and @xmath367 , which are found from the second - order perturbation theory @xcite . under the condition @xmath72 and @xmath138 , these eigenvalues are purely imaginary . therefore , the infinite - dimensional part of the spectrum of the operator @xmath400 persists on the imaginary axis for @xmath158 near the points @xmath366 and @xmath367 of infinite algebraic multiplicity . the statement of the lemma follows from the perturbation theory and the fact that all isolated eigenvalues and the continuous spectrum of @xmath401 are purely imaginary .    in the asymptotic limit @xmath197 as @xmath198 , the eigenvalues @xmath402 and @xmath396 defined by ( [ lambda-0-eig ] ) and ( [ lambda - plus - minus - eig ] ) are given asymptotically by @xmath403 where the remainder terms are @xmath277 as @xmath278 . the values @xmath402 , @xmath404 , and @xmath405 are close to each other as @xmath406 . computations in the proof of theorem [ theorem - stability ] can be extended to the branch ( c ) of lemma [ lemma - dimer ] . indeed , @xmath380 for branch ( c ) with either @xmath316 or @xmath317 , and @xmath6 near @xmath407 . on the other hand , @xmath408 if @xmath317 and @xmath6 near @xmath409 . as a result , branch ( c ) is spectrally stable in the former case and is spectrally unstable in the latter case , in agreement with theorem [ theorem - bounded ] .    observe in the proof of theorem [ theorem - stability ] that @xmath410 if @xmath95 . in this case , branch ( b ) of lemma [ lemma - dimer ] is spectrally unstable . this instability corresponds to the instability of the zero equilibrium for @xmath95 , in agreement with the result of proposition [ prop - zero - equilibrium ] .    before presenting numerical approximations of eigenvalues of the spectral problem ( [ eq : eigenval ] ) , we compute the krein signature of wave continuum . this helps to interpret instabilities and resonances that arise when isolated eigenvalues @xmath383 cross the continuous bands near points @xmath366 and @xmath367 . the krein signature of simple isolated eigenvalues is defined as follows . [ def - krein ] let @xmath411 be an eigenvector of the spectral problem ( [ eq : eigenval ] ) for an isolated simple eigenvalue @xmath380 . then , the energy quadratic form @xmath412 is nonzero and its sign is called the krein signature of the eigenvalue @xmath402 .    definition [ def - krein ] is used to simplify the presentation . similarly , one can define the krein signature of isolated multiple eigenvalues and the krein signature of the continuous spectral bands in the spectral problem ( [ eq : eigenval ] ) @xcite . the following lemma characterizes krein signatures of the spectral points arising in the proof of theorem [ theorem - stability ] . [ lemma - krein ] fix @xmath46 , @xmath93 , and @xmath413 with @xmath138 . assume the non - degeneracy condition ( [ assumption - non - degeneracy ] ) . for every @xmath237 sufficiently small , we have the following for the corresponding branches of lemma [ lemma - dimer ] :    * the subspaces of @xmath379 in @xmath63 near @xmath366 , @xmath367 , and @xmath383 have positive , negative , and positive krein signature , respectively ; * the subspaces of @xmath379 in @xmath63 near @xmath366 , @xmath367 , and @xmath383 have negative , positive ( if @xmath414 ) or negative ( if @xmath415 ) , and positive krein signature , respectively ; * all subspaces of @xmath379 in @xmath63 near @xmath366 , @xmath367 , and @xmath383 ( if @xmath380 ) have negative krein signature . we proceed by the perturbation arguments from the limit @xmath77 , where @xmath416 is a block - diagonal operator consisting of @xmath417 blocks . in particular , we consider the blocks for @xmath265 , where @xmath233 is given by ( [ block - l - n ] ) . solving ( [ spectral - problem ] ) at @xmath77 and @xmath388 , we obtain the eigenvector @xmath418 as a result , we obtain for the eigenvector @xmath419 , @xmath420    for branch ( a ) , @xmath79 and @xmath414 . therefore , @xmath421 for @xmath422 and @xmath423 for @xmath424 . for branch ( b ) , @xmath83 and either @xmath414 or @xmath425 . in either case , @xmath423 for @xmath422 . on the other hand , for @xmath424 , @xmath421 if @xmath414 and @xmath423 if @xmath425 .    for branch ( c ) , @xmath83 and @xmath425 . in this case , @xmath423 for either @xmath422 or @xmath424 .    finally , the krein signature for the eigenvalue @xmath426 denoted by @xmath427 follows from the computations of eigenvalues @xmath274 in the proof of theorem [ theorem - critical - point ] . we have @xmath428 for branches ( a ) and ( b ) because @xmath429 and we have @xmath430 for branch ( c ) because @xmath431 , whereas the eigenvalue @xmath263 is controlled by the result of lemma [ theorem - index ] . the signs of all eigenvalues are nonzero and continuous with respect to parameter @xmath38 . therefore , the count above extends to the case of small nonzero @xmath38 . the spectrum of @xmath401 is shown at figure  [ spectra ] . panels ( a ) , ( b ) and ( c ) correspond to branches shown at figure  [ branches ] .     for different branches of breathers . ]    * we can see on panel ( a ) of figure [ spectra ] that @xmath432 do not intersect for every @xmath414 and are located within fixed distance @xmath277 , as @xmath433 . note that the upper - most @xmath402 and @xmath404 have positive krein signature , whereas the lowest @xmath405 has negative krein signature , as is given by lemma [ lemma - krein ] . * we observe on panel ( b ) of figure [ spectra ] that @xmath434 intersects @xmath402 , creating a small bubble of instability in the spectrum . the insert shows that the bubble shrinks as @xmath435 , in agreement with theorem [ theorem - stability ] . there is also an intersection between @xmath436 and @xmath402 , which does not create instability . these results are explained by the krein signature computations in lemma [ lemma - krein ] . instability is induced by opposite krein signatures between @xmath404 and @xmath402 , whereas crossing of @xmath405 and @xmath402 with the same krein signatures is safe of instabilities . note that for small @xmath6 , the isolated eigenvalue @xmath402 is located above both the spectral bands near @xmath404 and @xmath405 . the gap in the numerical data near @xmath344 indicates failure to continue the breather solution numerically in @xmath38 , in agreement with the proof of theorem [ theorem - existence - soliton ] . * we observe from panel ( c ) of figure [ spectra ] that @xmath402 and @xmath437 intersect but do not create instabilities , since all parts of the spectrum have the same signature , as is given by lemma [ lemma - krein ] . in fact , the branch is both spectrally and orbitally stable as long as @xmath380 , in agreement with theorem [ theorem - bounded ] . on the other hand , there is @xmath318 , if @xmath317 , such that @xmath408 for @xmath438 , which indicates instability of branch ( c ) , again , in agreement with theorem [ theorem - bounded ] .    as we see on panel ( b ) of figure [ spectra ] , @xmath402 intersects @xmath404 for some @xmath439 . in the remainder of this section , we study whether this crossing point is always located on the right of @xmath440 . in fact , the answer to this question is negative . we shall prove for branch ( b ) that the intersection of @xmath402 with either @xmath404 or @xmath441 occur either for @xmath442 or for @xmath443 , depending on parameters @xmath39 and @xmath37 . fix @xmath46 , @xmath83 , and @xmath413 along branch ( b ) of lemma  [ lemma - dimer ] . there exists a resonance @xmath444 at @xmath445 with @xmath446 if @xmath447 and @xmath448 if @xmath449 , where @xmath450 moreover , if @xmath451 , there exists a resonance @xmath452 at @xmath445 with @xmath448 . let us first assume that there exists a resonance @xmath444 at @xmath453 and find the condition on @xmath39 and @xmath37 , when this is possible . from the definitions ( [ lambda-0-eig ] ) and ( [ lambda - plus - minus - eig ] ) , we obtain the constraint on @xmath181 : @xmath454 after canceling @xmath455 since @xmath456 with @xmath185 given by ( [ a - plus ] ) , we obtain @xmath457 which has two roots @xmath458 since @xmath103 , the lower sign is impossible because this leads to a contradiction @xmath459 the upper sign is possible if @xmath460 . using the parametrization  , we substitute the root for @xmath181 to the equation @xmath461 and simplify it : @xmath462 \\ & = & \frac{2 \sqrt{\omega^2 - 4 \gamma^2 } \left(2 |\omega| + \sqrt{\omega^2 - 4 \gamma^2}\right)^2}{\left(|\omega| + \sqrt{\omega^2 - 4 \gamma^2}\right)}.\end{aligned}\ ] ] this equation further simplifies to the form : @xmath463 squaring it up , we obtain @xmath464 which has only one positive root for @xmath465 given by @xmath466 this root yields a formula for @xmath467 in ( [ omega - star ] ) . since there is a unique value for @xmath468 , for which the case @xmath469 is possible , we shall now consider whether @xmath446 or @xmath470 for @xmath471 or @xmath447 .    to inspect the range @xmath470 , we consider a particular case , for which the intersection @xmath472 happens at @xmath4 . in this case , @xmath473 given by ( [ a - plus ] ) , so that the condition @xmath474 can be rewritten as @xmath475 there is only one negative root for @xmath37 and it is given by @xmath476 . by continuity , we conclude that @xmath444 for @xmath477 and @xmath452 for @xmath478 , both cases correspond to @xmath448 . finally , we verify that the case @xmath444 occurs for @xmath446 if @xmath447 . indeed , @xmath479 and @xmath480 as @xmath198 , so that @xmath481 as @xmath482 . on the other hand , the previous estimates suggest that @xmath483 for every @xmath425 if @xmath447 . therefore , there exists at least one intersection @xmath444 for @xmath446 if @xmath447 . we have reduced newton s equation of motion for coupled pendula shown on figure [ pendulapic ] under a resonant periodic force to the @xmath0-symmetric dnls equation ( [ pt - dnls ] ) . we have shown that this system is hamiltonian with conserved energy ( [ energy - pt - dnls ] ) and an additional constant of motion ( [ charge - pt - dnls ] ) . we have studied breather solutions of this model , which generalize symmetric synchronized oscillations of coupled pendula that arise if @xmath4 . we showed existence of three branches of breathers shown on figure [ branches ] . we also investigated their spectral stability analytically and numerically . the spectral information on each branch of solutions is shown on figure [ spectra ] . for branch ( c ) , we were also able to prove orbital stability and instability from the energy method . the technical results of this paper are summarized in table 1 and described as follows .              for branch ( a ) , we found that it is disconnected from the symmetric synchronized oscillations at @xmath4 . along this branch , breathers of small amplitudes @xmath5 are connected to breathers of large amplitudes @xmath5 . every point on the branch corresponds to the saddle point of the energy function between two wave continua of positive and negative energies . every breather along the branch is spectrally stable and is free of resonance between isolated eigenvalues and continuous spectrum . in the follow - up work @xcite , we will prove long - time orbital stability of breathers along this branch . for branch ( b ) , we found that the large - amplitude breathers as @xmath482 are connected to the symmetric synchronized oscillations at @xmath4 , which have the smallest ( but nonzero ) amplitude @xmath119 . breathers along the branch are spectrally stable except for a narrow instability bubble , where the isolated eigenvalue @xmath402 is in resonance with the continuous spectrum . the instability bubble can occur either for @xmath414 , where the breather is a saddle point of the energy function between two wave continua of opposite energies or for @xmath415 , where the breather is a saddle point between the two negative - definite wave continua and directions of positive energy . when the isolated eigenvalue of positive energy @xmath402 is above the continuous spectrum near @xmath404 and @xmath367 , orbital stability of breathers can be proved by using the technique in @xcite , which was developed for the dnls equation .    finally , for branch ( c ) , we found that the small - amplitude breathers at @xmath484 are connected to the symmetric synchronized oscillations at @xmath4 , which have the largest amplitude @xmath116 . breathers are either spectrally stable near @xmath344 or unstable near @xmath4 , depending on the detuning frequency @xmath37 and the amplitude of the periodic resonant force @xmath39 . when breathers are spectrally stable , they are also orbitally stable for infinitely long times .    * acknowledgements . * the authors thank i.v . barashenkov and p.g . kevrekidis for useful discussions regarding this project . the work of a.c . is supported by the graduate scholarship at mcmaster university . the work of d.p . is supported by the ministry of education and science of russian federation ( the base part of the state task no . 2014/133 , project no . 2839 ) . c. huygens , _ the pendulum clock or geometrical demonstrations concerning the motion of pendula as applied to clocks _ , iowa state university press series in the history of technology and science ( iowa state university press , ames , ia , 1986 ) . mohamed , m.a . rahman and r. seethapathy , _ robust line - voltage sensorless control and synchronization of lcl - filtered distributed generation inverters for high power quality grid connection _ , power electronics , ieee transactions * 27 * , 87 ( 2012 ) .
Abstract:

Answer: we derive a hamiltonian version of the @xmath0-symmetric discrete nonlinear schrdinger equation that describes synchronized dynamics of coupled pendula driven by a periodic movement of their common strings . in the limit of weak coupling between the pendula , we classify the existence and spectral stability of breathers ( time - periodic solutions localized in the lattice ) supported near one pair of coupled pendula . orbital stability or instability of breathers is proved in a subset of the existence region .