Global Existence and Compact Attractors for the Discrete Nonlinear Schrödinger equation

Nikos I. Karachalios and Athanasios N. Yannacopoulos Department of Statistics and Actuarial Science University of the Aegean Karlovassi, 83200, Samos, Greece

November 27, 2006

Abstract
We study the asymptotic behavior of solutions of discrete nonlinear Schrödinger-type (DNLS) equations. For a conservative system, we consider the global in time solvability and the question of existence of standing wave solutions. Similarities and differences with the continuous counterpart (NLS-partial differential equation) are pointed out.
For a dissipative system we prove existence of a global attractor and its stability under finite dimensional approximations.
Similar questions are treated in a weighted phase space. Finally, we propose possible extensions for various types of DNLS equations.

1 Introduction

The discrete nonlinear Schrödinger equation (DNLS), is a very popular model with a great variety of applications, ranging from physics to biology. To name but a few, the DNLS has been successfully applied to the modelling of localized pulse propagation in optical fibers and wave guides, to the study of energy relaxation in solids, the behaviour of amorphous material, to the modelling of self-trapping of vibrational energy in proteins or studies related to the denaturation of the DNA double strand (for an account of possible applications see e.g.
the recent review [23). The DNLS serves both as a model in its own right (modelling cases where the nature of the problem is inherently discrete), or as an approximation of the continuous nonlinear Schrödinger through a numerical scheme.
The study of the dynamics of the discrete nonlinear Schrödinger equation, has been an active theme of research in the past decade. The major part of the activity, has focused on the study of localized excitations and nonlinear waves (solitons and breathers), in the conservative DNLS. This activity has started with the ingenious construction introduced in [26, of localized in space time periodic or time-quasiperiodic solutions of dicrete lattice systems (which include the DNLS as special case), using the so called anti-integrable limit or anti-continuous limit. Such solutions are essentially found for the weak coupling limit between the lattice sites.
Other constructions are available using dynamical systems methods, such as the reduction to a symplectic map and the homoclinic trajectory method. The stability of such solutions has been studied using both analytic and numerical techniques and they have been shown in many cases to be robust under perturbations. They can be appplied to the study of diverse physical and biological phenomena such as lattice dynamics of solids, selective bond excitations, localization of electromagnetic waves in photonic crystals with nonlinear response etc. For detailed reviews on this theme we refer to [4, [14, [18. There have also been recent interesting works on more complicated solutions arising in higher dimensions than one, such as discrete vortex solutions (see eg. [24or [27).
In this work we mainly consider the asymptotic behavior of solutions for the following complex lattice dynamical system (LDS)
i u ˙ n + 1 ε ( u n 1 2 u n + u n + 1 ) + i δ u n + F ( u n ) = g n , n Z , ε > 0 , (1)
u n ( 0 ) = u n , 0 , n Z (2)
The real parameter δ > 0   introduces weak damping, while the case δ = 0   and g n = 0 ,   for all n Z   , corresponds to a conservative (NLS) lattice system.
Using techniques from the theory of infinite dimensional systems, we are able to show the global in time existence of solutions to both the conservative and the dissipative DNLS equation. This ensures the well posedeness of the model. This result is interesting as it provides a rigorous proof of the fact that blow up is not expected to occur in the discrete models as it may happen in the continuous model, a fact that has been observed numerically by a number of authors [7. Futhermore, for a finite lattice (assuming Dirichlet boundary conditions), and by using a discrete version of the mountain pass Theorem, we prove the existence of nontrivial standing waves, in the conservative case. For the infinite dimensional case, we show non-existence of nontrivial standing waves for certain values of energy, by using a fixed point argument. For the dissipative case we may further prove a more strong result; the existence of a global attractor, that attracts all bounded sets of the infinite dimensional phase space. It can further be shown that the dynamics on this global attractor may be approximated by the dynamics of finite dimensional systems. Results of this kind, for lattice dynamical systems of first and second order are presented in the recent works [8, 38, 39.
Many of the physically interesting solutions of the DNLS equation present strong spatial localization properties.
This is true for instance for soliton solutions or breathers. To cover such situations, we study the DNLS equation in weighted spaces, with properly chosen weight functions. For the dissipative case we may prove global existence of solutions in weighted spaces as well as the existence of a global attractor. Finally, we comment on possible generalizations of the above results in higher spatial dimensions and for more general models.

2 Formulation of the problem

Problem ( 1 )-( 2 ) is an infinite system of ordinary differential equations (infinite dimensional dynamical system). We will formulate the problem in an appropriate functional setting. In what follows, a complex Hilbert space X   with the sesquilinear form B X ( , )   and the norm | | | | X   , will be considered as a real Hilbert space with the scalar product ( , ) X = R e B X ( , )   .Let T : D ( T ) X X   a C   -linear, self-adjoint 0   operator with dense domain D ( T )   on the Hilbert space X   , equipped with the scalar product ( , ) X   . The space X T   is the completion of D ( T )   in the norm | | u | | T 2 = | | u | | X 2 ( T u , u ) X   , for u X T   , and we denote by X T *   its dual and by T *   the extension of T   to the dual of D ( T )   , denoted by D ( T ) *   . A function F : X X   is Lipschitz continuous on bounded subsets of X   , if for all M > 0   there exists a constant L ( M )   such that | | F ( u ) F ( v ) | | X L ( R ) | | u v | | X   , for all u , v B M   the closed ball of X   , of center 0   and of radius R   .
For the treatment of ( 1 )-( 2 ), we shall use complexifications of the usual real sequence spaces denoted by
p : = { u = ( u n ) n Z C : | | u | | p : = ( n Z | u n | p ) 1 p < } . (3)
Let us recall that between p   , spaces the following elementary embedding relation [19,pg.145holds,
q p , | | u | | p | | u | | q 1 q p . (4)
We will see in the sequel that relation ( 4 ) has important effects concerning the global in time solvability of DNLS equation, especially in the conservative case. For p = 2   we get the usual Hilbert space of square-summable (complex) sequences endowed with the real scalar product
( u , v ) 2 = R e n Z u n v n ¯ , u , v 2 . (5)
We describe next the discrete differential operators employed. For any u , v 2   we consider the linear operators A , B , B * : 2 2   ,
( B u ) n Z = u n + 1 u n , ( B * u ) n Z = u n 1 u n , (6)
( A u ) n Z = ( u n 1 2 u n + u n + 1 ) . (7)
Moreover it holds that
( B u , v ) 2 = ( u , B * v ) 2 , (8)
( A u , v ) 2 = ( B u , B v ) 2 , u , v 2 . (9)
The functional set up will also include (see [38) the following bilinear form and induced norm
( u , v ) 1 2 : = ( B u , B v ) 2 + ( u , v ) 2 , (10)
| | u | | 1 2 2 : = | | B u | | 2 2 + | | u | | 2 2 . (11)
We denote by 1 2   the Hilbert space with scalar product ( 10 ) and norm ( 11 ). The usual norm of 2   and ( 11 ) are equivalent, since
| | u | | 2 2 | | u | | 1 2 2 c | | u | | 2 2 , c > 0 . (12)
For local existence, we examine the following examples of nonlinearities F : C C   :
( N 1 )   F ( 0 ) = 0   and there exist constants c > 0   , β 0   such that | F ( z 1 ) F ( z 2 ) | c ( | z 1 | β + | z 2 | β ) | z 1 z 2 |   , or alternatively ( N 2 )   F ( z ) = f ( | z | 2 ) z   where f , f : R R   , are continuous.
Both cases include the power-law nonlinearity F ( z ) = | z | 2 σ z   (for ( N 2 )   we must have 1 < σ <   ). Since we intend to apply the general theory of abstract Schrödinger equations [11, 12, the following Lemma will be useful.
Lemma 2.1 Let F : C C   satisfy ( N 1 )   or ( N 2 )   . Then the function F   defines an operator (still denoted by F   ) F : 2 2 , ( F ( u ) ) n Z : = F ( u n ) ,   which is Lipschitz continuous on bounded sets of 2   .
Proof: We consider first the case ( N 1 )   . Let u B R   a closed ball in 2   of center 0   and radius R   . We have from ( 3 ) that
| | F ( u ) | | 2 2 c 2 n Z | u n | 2 β + 2 = c 2 | | u | | 2 β + 2 2 β + 2 c 2 | | u | | 2 2 β + 2 (13)
hence F : 2 2   , is bounded on bounded sets of 2   .
For u , v B R   , we observe that
| | F ( u ) F ( v ) | | 2 2 c 2 n Z ( | u n | β + | v n | β ) 2 | u n v n | 2 (14)
c 2 sup n Z [ ( | u n | β + | v n | β ) 2 ] n Z | u n v n | 2 c 2 4 R 2 β | | u v | | 2 2 , (15)
justifying that the map F : 2 2   is Lipschitz continuous on bounded sets of 2   with Lipschitz constant L ( R ) = c 2 R β   .
For the case ( N 2 )   we have
| | F ( u ) | | 2 2 = n Z | f ( | u n | 2 ) | 2 | u n | 2 . (16)
Since f : R R   is continuous, there exists a monotone increasing C 1   -function g : R R   such that
| f ( ρ ) | g ( | ρ | ) , for all ρ R . (17)
(see e.g [37,p.g796). Note also that | u n | 2 | | u n | | 2 2   for all n Z   . Now ( 16 ) and ( 17 ) imply the inequality
| | F ( u ) | | 2 2 n Z g ( | u n | 2 ) 2 | u n | 2 n Z g ( | | u | | 2 2 ) 2 | u n | 2
{ m a x ρ [ 0 , R 2 ] g ( ρ ) } 2 n Z | u n | 2
K ( R ) 2 | | u | | 2 2 . (18)
with K ( R ) = g ( R 2 )   .
To check the Lipschitz property, we may argue as for the proof of ( 18 ). For some θ ( 0 , 1 )   and an appropriate C 1   -function g 1 : R R   , we get the inequality
| | F ( u ) F ( v ) | | 2 2 2 n Z | f ( | u n | 2 ) | 2 | u n v n | 2 + 2 n Z | f ( | u n | 2 ) f ( | v n | 2 ) | 2 | v n | 2 , (19)
n Z | f ( | u n | 2 ) f ( | v n | 2 ) | 2 | v n | 2 = n Z | f ( θ | u n | 2 + ( 1 θ ) | v n | 2 ) | 2 ( | u n | + | v n | ) 2 ( | u n | | v n | ) 2 | v n | 2
{ m a x ρ [ 0 , 2 R 2 ] g 1 ( ρ ) } 2 sup n Z [ ( | u n | + | v n | ) 2 | v n | 2 ] n Z | u n v n | 2 . (20)
A combination of ( 19 ) and ( 20 ), implies that there exists a constant L ( R ) = ( 2 g 2 ( R 2 ) + 8 R 4 g 1 ( 2 R 2 ) ) 1 / 2   such that
| | F ( u ) F ( v ) | | 2 L ( R ) | | u v | | 2 . (21)
The Lemma is proved.   We observe by ( 9 ), that the operator A   satisfies the relations
( A u , u ) 2 = | | B u | | 2 2 0 , (22)
( A u , v ) 2 = ( u , A v ) 2 , (23)
therefore defines a self-adjoint operator on D ( A ) = X = 2   and A 0   . Note that the graph norm | | u | | D ( A ) = | | A u | | 2 2 + | | u | | 2 2 ,   is also an equivalent norm with the 2   -norm since | | u | | 2 2 n Z | u n + 1 2 u n + u n 1 | 2 + n Z | u n | 2 c | | u | | 2 2 .   In our case, as it is indicated by ( 10 )-( 22 )-( 23 ), we may choose X A = 1 2   equipped with the norm | | u | | A 2 = | | u | | X 2 ( A u , u ) X | | u | | 1 2   , for u 2   . Moreover, D ( A ) = X = 2 = D ( A ) *   . Obviously A * = A   and the operator i A : 2 2   defined by ( i A ) u = i A u   for u 2   , is C   -linear and skew-adjoint and i A   generates a group ( T ( t ) ) t R   , of isometries on 2   . It is easy to check that the same properties hold for the operator ε 1 A   with ε > 0   (i.e., we may consider the equivalent norm ε 1 | | B u | | 2 2 + | | u | | 2   ).
Thus, for fixed T > 0   and ( u n , 0 ) n Z : = u 0 2   , a function u C ( [ 0 , T ] , 2 )   is a solution of ( 1 )-( 2 ) if and only if
u ( t ) = T ( t ) u 0 + i 0 t T ( t s ) F 1 ( u ( s ) ) d s , ( F 1 ( u ) ) n Z : = F 1 ( u n ) = i δ u n + F ( u n ) g n . (24)
Our local existence result can be stated as follows:
Theorem 2.1 Let assumptions ( N 1 )   or ( N 2 )   be satisfied, assume that g : = ( g n ) n Z 2   . Then there exists a function T * : 2 ( 0 , ]   with the following properties:
(a) For all u 0 2   , there exists u C ( [ 0 , T * ( u 0 ) ) , 2 )   such that for all 0 < T < T * ( u 0 )   , u   is the unique solution of ( 1 )-( 2 ) in C ( [ 0 , T ] , 2 )   (well posedeness).
(b) For all t [ 0 , T * ( u 0 ) )   ,
T * ( u 0 ) t 1 2 ( L 1 ( R ) + 1 ) : = T R , R = | | g | | 2 + 2 | | u ( t ) | | 2 , L 1 ( R ) = δ + L ( R ) , (25)
where L ( R )   is the Lipschitz constant for the map F : 2 2   . Moreover the following alternative holds: (i) T * ( u 0 ) =   , or (ii) T * ( u 0 ) <   and lim t T * ( u 0 ) | | u ( t ) | | 2 =   (maximality).
(c) T * : 2 ( 0 , ]   is lower semicontinuous. In addition, if { u n 0 } n N   is a sequence in 2   such that u n 0 u 0   and if T < T * ( u 0 )   , then S ( t ) u 0 n S ( t ) u 0   in C ( [ 0 , T ] , 2 )   , where S ( t ) u 0 = u ( t )   , t [ 0 , T * ( u 0 ) )   , denotes the solution operator (continuous dependence on initial data).
Proof: Since the result is an application of the results in [11,p.g.56-59, we present only an outline of the proof, for the sake of completeness. It follows from Lemma  2.1 that F 1 : 2 2   is Lipschitz continuous on bounded sets of 2   with Lipschitz constant L 1 ( R ) = δ + L ( R )   . For all u 0 2   with | | u 0 | | 2 M   we set R = 2 M + | | F 1 ( 0 ) | | 2 = 2 M + | | g | | 2   . Following closely the lines of [11,Lemma4.3.2-Proposition4.3.3and using ( 24 ) we may show that the map Φ u ( t ) : = T ( t ) u 0 + i 0 t T ( t s ) F 1 ( u ( s ) ) d s ,   is a contraction on the complete metric space ( X R , d )   where X R : = { u C ( [ 0 , T R ] , 2 ) : | | u ( t ) | | 2 R , t [ 0 , T R ] } , T R = 1 2 L 1 ( R ) + 1 , d = max t [ 0 , T R ] | | u ( t ) v ( t ) | | 2 .   For u 0 2   we define  24  T * ( u 0 ) = sup { T > 0 : S ( t ) u 0 = u ( t ) C ( [ 0 , T ] , 2 ) solution of ( ) }   . Inequality ( 25 ) follows by the contradiction argument of [11,Theorem4.3.4: We assume instead of ( 25 ), that for R = | | u ( t ) | | 2   , T * ( u 0 ) t < T R   and we consider the solution v C ( [ 0 , T R ] , 2 )   of v ( s ) = T ( s ) u ( t ) + i 0 s T ( s σ ) F 1 ( v ( σ ) ) d σ , s [ 0 , T R ] .   Then we can define a solution w C ( [ 0 , t + T R ] , 2 )   of ( 24 ) at T = t + T R   , defined by
w ( s ) = { u ( s ) , s [ 0 , t ] , v ( s t ) , s [ t , t + T R ] .
which is in contradiction with the assumption T * ( u 0 ) < t + T R   . To prove (c), as in [11,Proposition4.3.7, we set R = 2 sup t [ 0 , T ] | | u ( t ) | | 2   and define τ n = sup { t [ 0 , T * ( u 0 n ) ) : | | S ( s ) u 0 n | | 2 2 R , s [ 0 , t ] }   . Since u 0 n u 0   in 2   , it follows that for sufficiently large n   , | | u 0 n | | 2 < R   . Therefore τ n > T R > 0   . Moreover, for all t T   and t τ n   , we get from ( 24 ) and Gronwall's Lemma the inequality
| | S ( t ) u 0 n S ( t ) u 0 | | 2 e L 1 ( 2 R ) T | | u 0 n u 0 | | 2 .
Letting n   , we get that S ( t ) u 0 n S ( t ) u 0   in 2   and τ n > T   , which implies that T * ( u 0 n ) > T   .  

3 Global existence for the conservative case

In this section we discuss the global in time solvability of ( 1 )-( 2 ) for the case δ = 0   and g n = 0   , for all n Z   (conservative case). We consider as a model problem, the following DNLS equation with arbitrary power nonlinearity
i u ˙ n + 1 ε ( u n 1 2 u n + u n + 1 ) + | u n | 2 σ u n = 0 , n Z , 0 σ < , ε > 0 , (26)
u n ( 0 ) = u n , 0 , n Z (27)
We also consider the (NLS) partial differential equation
i t u + u x x + | u | 2 σ u = 0 , x R , t > 0 , (28)
u ( x , 0 ) = u 0 ( x ) .
We will show that there is a vast difference concerning the global solvability of the discrete (NLS) equation ( 26 ) with ε > 0   and its continuum limit as ε 0   given by ( 28 ). Concerning ( 26 ), we have the following Theorem.
Theorem 3.1 Let u 0 2   . For any ε > 0   , and 0 σ <   , there exists a unique solution of ( 26 ) such that u ( t ) C 1 ( [ 0 , ) , 2 )   .
Proof: The result is a consequence of discrete conservation laws satisfied by the solutions of ( 26 )-( 27 ). We take the scalar product of ( 26 ) with i u   . By using ( 22 ) we obtain
d d t | | u ( t ) | | 2 2 = 0 , or | | u ( t ) | | 2 2 = | | u 0 | | 2 2 , for every t [ 0 , T * ( u 0 ) ) , (29)
i.e. | | u ( t ) | | 2   is uniformly bounded on the maximal interval of existence. It follows then by Theorem  2.1 (b)-(c) that T * ( u 0 ) =   and s u p { | | u ( t ) | | 2 , t [ 0 , ) } <   .  
Remark 3.1 Theorem  3.1 can be also established by a discrete version of the conservation of energy, satisfied by solutions of ( 26 ). This approach elucidates the role of the nonlinearity exponent, in the global in time solvability of DNLS: It can be easily checked that
E ( u ( t ) ) = E ( u 0 ) , E ( u ( t ) ) : = 1 ε | | u ( t ) | | 1 2 2 1 σ + 1 | | u ( t ) | 2 σ + 2 2 σ + 2 , t [ 0 , T * ( u 0 ) ) . (30)
Then, by using ( 4 ) and ( 30 ) we may derive the estimate
| | u ( t ) | | 1 2 2 | | u 0 | | 1 2 2 + 2 ε σ + 1 | | u 0 | | 2 2 σ + 2 . (31)
As a consequence of ( 31 ) we obtain once again that T * ( u 0 ) =   and s u p { | | u ( t ) | | 1 2 , t [ 0 , ) } <   .
Clearly, Theorem  3.1 , covers the case of DNLS with a nonlinearity satisfying ( N 2 )   . The case of a globally Lipschitz continuous function F : [ 0 , ) R   with F ( 0 ) = 0   is also included. Following [11,Section7.2, [12,p.g53, such a function can be extended to the complex function F : C \ { 0 } C   by setting F ( z ) = z | z | F ( | z | )   . It can be shown as in Lemma  2.1 that F : 2 2   is locally Lipschitz and ( F ( u ) , i u ) 2 = 0   which implies ( 29 ). Setting F * ( z ) = 0 | z | F ( s ) d s   it can be shown as in Section 4, Lemma  4.2 that the functional ( V ( u ) ) n : = n Z F * ( u n )   is a C 1   -functional on 2   and V ( u ) = F ( u )   . In this case, we observe conservation of energy defined by E ( u ( t ) ) : = 1 ε | | u ( t ) | | 1 2 2 + V ( u )   .  
For a comparison between the discrete (NLS) equation ( 26 ) and its continuous counterpart ( 28 ), we refer to the main results concerning ( 28 ) (see [11, 12, 28).
For u 0 H 1 ( R )   and 0 σ <   there exists a unique maximal solution of ( 28 ), u ( t ) C ( [ 0 , T m a x ) , H 1 ( R ) ) C 1 ( [ 0 , T m a x ) , L 2 ( R ) )   . In addition: If σ < 2   then T max =   and u   is bounded in H 1 ( R )   . Let σ 2   . Assume that u 0 H 1 ( R )   such that R | x | 2 | u 0 | 2 d x <   (initial data with finite variance) and E ( u 0 ) <   . Then T max <   .
On the other hand if | | u 0 | | H 1   is sufficiently small, T m a x =   and u   is bounded in H 1 ( R )   .
It follows from Theorem  3.1 , that solutions of ( 26 ) with ε > 0   , exist globally in 1 2   -norm unconditionally with respect to the degree of the nonlinearity the size of the initial data and the sign of the initial energy E ( u 0 )   . This is in agreement with the numerical observations in [7: According to [7, for a discrete system of the form ( 26 ), numerical simulations provide evidence that the solution of the discrete equation still exists after localization, while the solution of the continuum system blows-up in finite time. Note that the discrete version of initial data with finite variance reads as n Z n 2 | u n , 0 | 2 <   and such data belong to 2   .

4 Existence of standing wave solutions for the conservative DNLS

We conclude the discussion on the conservative DNLS, with a discussion on the existence of standing wave solutions. The standing wave solution for fixed ω 2   is given by the ansatz
u n ( t ) = e i ω 2 t φ n , n Z , φ n independent of t . (32)
We do not restrict φ n   to be real. We are focused on DNLS ( 26 ). It can be easily seen, that any standing wave solution of ( 26 ), satisfies the following equation
1 ε ( φ n 1 2 φ n + φ n + 1 ) + ω 2 φ n = | φ n | 2 σ φ n , n Z . (33)
We will study in this section two related problems. The first problem we will study is given by equation ( 33 ), with n   taking finite values. We consider the finite dimensional subspace of 2   C 2 m + 1 : = { ψ n 2 : ψ ( m + 1 ) = ψ m + 1 = 0 } .   We consider the following nonlinear system in C 2 m + 1   ,
1 ε ( φ n 1 2 φ n + φ n + 1 ) + ω 2 φ n = | φ n | 2 σ φ n , | n | m , (34)
φ ( m + 1 ) = φ m + 1 = 0 . (35)
This is a finite dimensional problem, and is related to the problem of existence of standing wave solutions for the DNLS equation in a finite lattice with Dirichlet boundary conditions. Although this problem is finite dimensional, it is still of interest as any numerical approximation of the continuous space NLS, or even of the infinite dimensional DNLS will necessarily lead to a finite dimensional problem of the above type.
The set C 2 m + 1   endowed with the discrete inner product and induced norm
( φ , ψ ) 2 : = R e n = m n = m φ n ψ n ¯ , | | ψ | | 2 : = n = m n = m | ψ n | 2 , φ , ψ C 2 m + 1 , (36)
is a (finite dimensional) Hilbert space. We consider now the operators ( A 1 ψ ) | n | m : = ψ n 1 2 ψ n + ψ n + 1 , ( B 1 ψ ) | n | m = ψ n + 1 ψ n .   It can be easily checked (see also [3,pg.117) that
( A 1 ψ , ψ ) 2 = n = m n = m | ψ n + 1 ψ n | 2 , ( A 1 φ , ψ ) 2 = ( B 1 φ , B 1 ψ ) 2 . (37)
Hence, we may also consider the inner product and the corresponding norm in C 2 m + 1  
( φ , ψ ) 1 , 2 : = ( B 1 φ , B 1 ψ ) 2 + ( φ , ψ ) 2 , | | ψ | | 1 , 2 : = n = m n = m ( | ψ n + 1 ψ n | 2 + | ψ n | 2 ) . (38)
The norm in ( 38 ) is equivalent with ( 36 ), and the constants in the equivalence inequality, can be chosen in order to be independent of m   .
The second is the problem given by equation ( 33 ), which is an infinite dimensional system. This is related with the problem of existence of standing wave solutions of the conservative DNLS equation in the infinite lattice.
In what follows we will use a variational principle to show existence of nontrivial standing wave solutions for ( 34 )-( 35 ) and a fixed point argument to show non-existence of non-trivial standing waves, for certain parameter values for the infinite dimensional problem ( 33 ).

4.1 Existence of non trivial standing wave solutions for the finite dimensional problem

We will treat this problem using a variational approach. We start with some observations and results which are common to both problems. Unless the opposite is explicitly stated, everything that follows will hold for both the finite dimensional and the infinite-dimensional system.
Solutions of ( 33 ) are critical points of the functional
E ( φ ) = 1 2 ε n Z | ( B φ ) n | 2 + ω 2 2 n Z | φ n | 2 1 2 σ + 2 n Z | φ n | 2 σ + 2 . (39)
To establish differentiability of the functional E : 2 R   , we shall use the following discrete version of the dominated convergence Theorem, provided by [9.
Lemma 4.1 Let { ψ i , k }   be a double sequence of summable functions (i.e i Z | ψ i , k | <   ) and lim k ψ i , k = ψ i   , for all i Z   . If there exists a summable sequence { g i }   such that | ψ i , k | g i   for all i , k   's, we have that lim k i Z ψ i , k = i Z ψ i   .
We then have the following Lemma.
Lemma 4.2 Let ( φ n ) n Z = φ 2 σ + 2   for some 0 < σ <   . Then the functional V ( φ ) = n Z | φ n | 2 σ + 2 ,   is a C 1 ( 2 σ + 2 , R )   functional and
< V ( φ ) , ψ > = ( 2 σ + 2 ) R e n Z | φ n | 2 σ φ n ψ n ¯ , ψ = ( ψ n ) n Z 2 σ + 2 . (40)
Proof: We assume that φ , ψ 2 σ + 2   . Then for any n Z   , 0 < s < 1   , we get from the mean value Theorem
V ( φ + s ψ ) V ( ψ ) s = 1 s R e n Z 0 1 d d θ | φ n + θ s ψ n | 2 σ + 2 d θ
= ( 2 σ + 2 ) R e n Z 0 1 | φ n + s θ ψ n | 2 σ ( φ n + s θ ψ n ) ψ n ¯ d θ . (41)
For the rhs of ( 41 ), we have the estimate
( 2 σ + 2 ) n Z | φ n + θ s ψ n | 2 σ + 1 | ψ n | ( 2 σ + 2 ) n Z ( | φ n | + | ψ n | ) 2 σ + 1 | ψ n | = n Z z n
The sequence ( z n ) n Z   is summable since
n Z | z n | c ( n Z ( | φ n | + | ψ n | ) 2 σ + 2 ) 2 σ + 1 2 σ + 2 ( n Z | ψ n | 2 σ + 2 ) 1 2 σ + 2 .
Letting s 0   , the existence of the Gateaux derivative ( 40 ) of the functional V : 2 σ + 2 R   , follows from Lemma  4.1 (discrete dominated convergence).
We show next that the functional V : 2 σ + 2 2 σ + 2 2 σ + 1   is continuous. For ( φ n ) n Z = φ 2 σ + 2   , we set F 1 ( φ n ) = | φ n | 2 σ φ n   . We consider a sequence φ m 2 σ + 2   such that φ m φ   in 2 σ + 2   . We get then the inequality
| V ( φ m ) V ( φ ) , ψ | c | | F 1 ( φ m ) F 1 ( φ ) | | q | | ψ | | p , q = 2 σ + 2 2 σ + 1 , p = 2 σ + 2 . (42)
We denote by ( φ m ) n   the n   -th coordinate of the sequence φ m 2   . Since F 1   satisfies condition ( N 1 )   with β = 2 σ   , by setting Φ n = | ( φ m ) n | 2 σ + | φ n | 2 σ   , we get from Hölder's inequality that
n Z | F 1 ( ( φ m ) n ) F 1 ( φ n ) | q c n Z ( Φ n ) q | ( φ m ) n φ n | q
c ( n Z | ( φ m ) n φ n | 2 σ + 2 ) 1 2 σ + 1 ( n Z ( Φ n ) σ + 1 σ ) 2 σ 2 σ + 1 0 , as m .
By ( 22 ),( 23 ) and Lemma ( 4.2 ), it follows that the functional E   defined by ( 39 ) is C 1 ( 2 , R )   . For convenience, we recall [13,Definition4.1,pg.130(PS-condition) and [13,Theorem6.1,pg.140or [31,Theorem6.1,pg. 109(Mountain Pass Theorem of Ambrosseti-Rabinowitz [1).
Definition 4.1 Let X   be a Banach space and E : X R   be C 1   . We say that E   satisfies condition ( P S )   if, for any sequence { u n } X   such that | E ( u n ) |   is bounded and E ( u n ) 0   as n   , there exists a convergent subsequence. If this condition is only satisfied in the region where E α > 0   (resp E α < 0   ) for all α > 0   , we say E   satisfies condition ( P S + )   (resp.
( P S )   ).
Theorem 4.1 Let E : X R   be C 1   and satisfy (a) E ( 0 ) = 0   , (b) ρ > 0   , α > 0 : | | u | | X = ρ   implies E ( u ) α   , (c) u 1 X : | | u 1 | | X ρ   and E ( u 1 ) < α   . Define Γ = { γ C 0 ( [ 0 , 1 ] , X ) : γ ( 0 ) = 0 , γ ( 1 ) = u 1 } .   Let F γ = { γ ( t ) X : 0 t 1 }   and = { F γ : γ Γ }   . If E   satisfies condition ( P S )   , then β : = inf F γ sup { E ( v ) : v F γ } α   is a critical point of the functional E   .
We shall verify definition  4.1 and the assumptions of Theorem  4.1 to show that problem ( 34 )-( 35 ), has a nontrivial solution.
Theorem 4.2 Consider the finite dimensional problem ( 34 )-( 35 ). Let 0 < σ <   and ε > 0   . Then for any ω 0   , there exists a non-trivial solution of ( 34 )-( 35 ).
Proof (Condition ( P S )   ) We observe that the norm
| | φ | | ε 2 = ( 1 ε | | B φ | | 2 2 + ω 2 | φ | | | 2 2 ) 1 2 , (43)
is equivalent with the norm of 1 2   defined by ( 11 ) since
min { 1 ε , ω 2 } | | φ | | 1 2 2 | | φ | | ε 2 2 max { 1 ε , ω 2 } | | φ | | 1 2 2 , (44)
and is also equivalent with the usual norm of 2   , as it follows from ( 12 ). Now we let a sequence φ m   of 2   be such that | E ( φ m ) | < M   for some M > 0   and E ( φ m ) 0   as m   . By using ( 39 ) and Lemma  4.2 , we observe that for m   sufficiently large
M E ( φ m ) 1 2 σ + 2 E ( φ m ) , φ m = ( 1 2 1 2 σ + 2 ) | | φ m | | ε 2 2 . (45)
Therefore the sequence φ m   is bounded. Thus, we may extract a subsequence, still denoted by φ m   , such that
φ m φ i n 2 , as m . (46)
Focusing now on ( 34 )-( 35 ), we observe that since in the finite dimensional space C 2 m + 1   the weak convergence coincides with the strong, φ m   converges strongly to φ   in C 2 m + 1   . Therefore the functional E   associated with the boundary value problem ( 34 )-( 35 ), satisfies condition ( P S )   .
(Mountain Pass assumptions) We now check the conditions for the validity of the Mountain Pass Theorem for the functional E   . We use the same notations for the equivalent norms in the finite dimensional space C 2 m + 1   , with those of the norms of the infinite dimensional spaces, having in mind that we deal with finite dimensional sums.
For every φ C 2 m + 1   it holds that
| | φ | | 2 σ + 2 κ 1 | | φ | | ε 2 , κ 1 = 1 min { 1 / ε , ω } .
The following inequality holds
E ( φ ) = 1 2 | | φ | | ε 2 2 1 2 σ + 2 | | φ | | 2 σ + 2 2 σ + 2
1 2 | | φ | | ε 2 2 κ 1 2 σ + 2 2 σ + 2 | | φ | | ε 2 2 σ + 2 .
Choosing φ C 2 m + 1   such that | | φ | | ε 2 = r   , we observe that if 0 < r < ( σ + 1 κ 1 2 σ + 2 ) 1 2 σ   , then
E ( φ ) α > 0 , α = r 2 ( 1 2 κ 1 2 σ + 2 2 σ + 2 r 2 σ ) .
This establishes the first condition for the validity of the Mountain Pass Theorem, for the functional E   . To check the behaviour of E   at infinity, let us consider φ = t e   where t R   and e C 2 m + 1   , | | e | | ε 2 = 1   . We have that
E ( t e ) = t 2 2 1 2 σ + 2 t 2 σ + 2 | | e | | 2 σ + 2 2 σ + 2 ,
as t +   . This ensures the validity of the second condition of the Mountain Pass Theorem. For fixed φ 0   and choosing t   sufficiently large, we may set u 1 = t φ   to obtain the existence of a non-trivial solution for ( 34 )-( 35 ).
 
Remark 4.1 In theorem  4.2 we have to state explicitly that we deal with the finite dimensional problem, since otherwise we may not prove the validity of the PS condition. In the infinite dimensional case, due to the lack of the Schur property for the space 2   (in contrast with the space 1   which posses this property-weak convergence coincides with strong convergence), we may not conclude the strong convergence of the sequence, from its weak convergence. The generalization of the mountain pass argument in the case of the infinite dimensional conservative system is a subject of future research, and will be considered elsewhere.  

4.2 Non-existence of non trivial standing waves

We now provide some non-existence results for nontrivial standing waves. The results that follow are valid for both problems ( 34 )-( 35 ) and ( 33 ).
We recall [37,Theorem18.E,pg.68(Theorem of Lax and Milgram), which will be used to establish existence of solutions for an auxiliary linear, non-homogeneous problem, related to ( 33 ).
Theorem 4.3 Let X   be a Hilbert space and A : X X   be a linear continuous operator. Suppose that there exists c * > 0   such that
R e ( A u , u ) X c * | | u | | X 2 , for all u X . (47)
Then for given f X   , the operator equation A u = f , u X   , has a unique solution
Theorem 4.4 There exist no nontrivial standing waves of energy less than E c ( ω , σ ) : = ( ω 4 4 ) 1 / 4 σ   .
Proof: Let ω R   , ω 0   , be fixed. For any ε > 0   , the operator A ω : 2 2   , defined by
( A ω φ ) n Z = 1 ε ( φ n 1 2 φ n + φ n + 1 ) + ω 2 φ n , (48)
is linear and continuous and satisfies assumption ( 47 ) of Theorem  4.3 : Using ( 44 ) and ( 11 ) we get that
( A ω φ , φ ) 2 = 1 ε | | B φ | | 2 2 + ω 2 | | φ | | 2 ω 2 | | φ | | 2 2 for all φ 2 . (49)
For given z 2   , we consider the linear operator equation
( A ω φ ) n Z = | z n | 2 σ z n , φ 2 . (50)
The function F 1 : C C   , F 1 ( z ) = | z n | 2 σ z n   satisfies Lemma  2.1 , with β = 2 σ   . Therefore assumptions of Theorem  4.3 are satisfied, and ( 50 ) has a unique solution φ 2   . For some R > 0   , we consider the closed ball of 2   , B R : = { z 2 : | | z | | 2 R }   , and we define the map P : 2 2   , by P ( z ) : = φ   where φ   is the unique solution of the operator equation ( 50 ). Clearly the map P   is well defined.
Let z   , ξ B R   such that φ = P ( z )   , ψ = P ( ξ )   . The difference χ : = φ ψ   satisfies the equation
( A ω χ ) n Z = F 1 ( z ) F 1 ( ξ ) . (51)
Taking now the scalar product of ( 51 ) with χ   and using ( 14 ), we derive
1 ε | | B χ | | 2 2 + ω 2 | | χ | | 2 2 | | F 1 ( z ) F 1 ( ξ ) | | 2 | | χ | | 2
2 R 2 σ | | z ξ | | 2 | | χ | | 2
ω 2 2 | | χ | | 2 2 + 2 ω 2 R 4 σ | | z ξ | | 2 2 . (52)
From ( 52 ), ( 11 ), ( 44 ), we obtain the inequality
| | χ | | 2 2 = | | P ( z ) P ( ξ ) | | 2 2 4 ω 4 R 4 σ | | z ξ | | 2 2 . (53)
Since P ( 0 ) = 0   , from inequality ( 53 ), we derive that for R < ( ω 4 / 4 ) 1 / 4 σ : = E c ( ω , σ )   , the map P : B R B R   and is a contraction. Therefore P   , satisfies the assumptions of Banach Fixed Point Theorem and has a unique fixed point, the trivial one. Hence, for R < E c ( ω , σ )   the only standing wave is the trivial.  
Remark 4.2 The above theorem is interesting when viewed in the following way: The 2   norm corresponds to the energy of the excitations of the lattice. If the energy of the excitation is less that E c ( ω , σ )   the lattice may not support a standing wave of frequency ω   . This relation may be seen as some kind of dispersion relation of frequency vs energy for the standing wave solutions. It contains information on the type of nonlinearity, through its dependence on the nonlinearity exponent σ   . Let us note that a similar computation to ( 52 ) shows that P   maps also the ball B R   with R ( ω 4 ) 1 / 4 σ : = R 1   , to itself and obviously E c < R 1   , but wthin this range, the map is not essentially a contraction.  
Remark 4.3 Arguments of the same type as above, may be used for the proof of existence of unique steady states for the case of nonlinearities that do not vanish when φ = 0   (see the end of Remark  3.1 ).  

5 Existence of global attractor in 2   for the dissipative case

In this section we study the asymptotic behavior of solutions of ( 1 )-( 2 ) in the dissipative case, δ > 0   and g n 0   , i.e. taking into account the effect of a weak dissipation and of an external excitation. This time, as a model problem we shall employ the dissipative case of ( 26 )-( 27 )
i u ˙ n + 1 ε ( u n 1 2 u n + u n + 1 ) + i δ u n + | u n | 2 σ u n = g n , n Z , 0 σ < , ε > 0 , (54)
u n ( 0 ) = u n , 0 , n Z . (55)
We also consider the weakly damped driven NLS partial differential equation
i t u + u x x + i δ u + | u | 2 σ u = g , x R , t > 0 , (56)
u ( x , 0 ) = u 0 ( x ) .
It is well known that the dissipation and forcing terms have important effects in the long-time behavior of ( 56 ). As an example, we refer to the case σ = 2   , for which has been observed numerically the existence of chaotic attractors [10, 29, in constrast with the conservative case which is completely integrable by the inverse scattering theory, for sufficiently smooth initial data. We refer to [15, 16, 21, 22, 34for results on the existence, finite dimensionality, and regularity of global attractors for ( 56 ).

5.1 Existence of absorbing set

This subsection is devoted to the existence of an absorbing set in 2   . We have the following Lemma.
Lemma 5.1 Let ( u 0 , n ) n Z = u 0 , ( g n ) n Z = g 2   . For any ε > 0   and 0 σ <   , the dynamical system defined by ( 54 )-( 55 )
S ( t ) : u 0 2 u ( t ) 2 , (57)
possesses a bounded absorbing set 0   in 2   : For every bounded set   of 2   , there exists t 0 ( , 0 )   such that for all t t 0 ( , 0 )   , it holds S ( t ) 0   .
Proof: Ẇe take the scalar product of ( 54 ) with i u ( t )   and we get
1 2 d d t | | u | | 2 2 + δ | | u | | 2 2 = I m n Z u ¯ n g n δ 2 | | u | | 2 2 + 1 2 δ | | g | | 2 2 . (58)
and therefore
d d t | | u | | 2 2 + δ | | u | | 2 2 1 δ | | g | | 2 2 . (59)
From ( 59 ) we derive that u L ( [ 0 , ) , 2 )   : Gronwall's Lemma implies that
| | u ( t ) | | 2 2 | | u 0 | | 2 2 exp ( δ t ) + 1 δ 2 | | g | | 2 2 { 1 exp ( δ t ) } . (60)
Letting t   we infer that
limsup t | | u ( t ) | | 2 2 1 δ 2 | | g | | 2 2 .
Setting ρ 2 = | | g | | 2 2 / δ 2   , it follows that for any number ρ 1 > ρ   the ball 0   of 2   centered at 0   of radius ρ 1   is an absorbing set for the semigroup S ( t )   : If   is a bounded set of 2   included in a ball of 2   centered at 0   of radius R   , then for t t 0 ( , 0 )   where
t 0 = 1 δ log R 2 ρ 1 2 ρ 2 , (61)
it holds | | u ( t ) | | 2 2 ρ 1 2   , i.e. S ( t ) 0   . Note that in the absence of external excitation, the dynamical system exhibits trivial dynamics, in the sense that limsup t | | u ( t ) | | 2 2 = 0   , as ( 60 ) clearly shows.  
Remark 5.1 The existence of the absorbing set, can be shown again by working in a different way, in the equivalent norm of 1 2   : By taking the scalar product of ( 54 ) with u ˙ δ u   we may derive the energy equation
1 2 d d t J ( u ( t ) ) + δ J ( u ( t ) ) = Λ ( u ( t ) ) , J ( u ) : = | | u | | 1 2 2 ε { 1 σ + 1 | | u | | 2 σ + 2 2 σ + 2 2 ( g , u ) } , (62)
Λ ( u ) : = ε δ { σ σ + 1 | | u | | 2 σ + 2 2 σ + 2 + ( g , u ) 2 } + I m n Z u ¯ n g n . (63)
From ( 4 ) and ( 61 ) we obtain the following estimate for Λ ( u )  
Λ ( u ) ε δ { σ σ + 1 ρ 1 2 σ + 2 + | | g | | 2 ρ 1 } + | | g | | 2 ρ 1 = λ 1 . (64)
Gronwall's inequality applied in ( 62 ) implies that
J ( u ( t ) ) J ( u 0 ) exp ( 2 δ t ) + λ 1 δ { 1 exp ( 2 δ t ) } (65)
Letting t   we get limsup t J ( u ( t ) ) λ 1 / δ   . Hence, for t t 0  
| | u ( t ) | | 1 2 2 ε ρ 1 2 σ + 2 + 3 ε | | g | | 2 ρ 1 + 1 δ | | g | | 2 ρ 1 : = ρ 2 2 (66)
Relation ( 66 ), gives an estimate of the radius of the absorbing ball, involving the degree of the nonlinearity and lattice spacing. Note for example, that for sufficiently small 0 < ε < ρ 1 | | g | | 2 δ 1 ρ 1 2 σ + 1 + 3 | | g | | 2 ,   we get that ρ 2 < ρ 1   .  
Throughout the rest of the paper, for simplicity reasons we set ε = 1   .

5.2 Asymptotic compactness of the semigroup and existence of global attractor

We now prove the asymptotic compactness of the semigroup generated by the DNLS equation ( 54 )-( 55 ).
To establish this property, we follow the approach introduced in [8and applied also in [38for a second order lattice dynamical system. This approach is based on the derivation of estimates on the tail ends of solutions of ( 54 )-( 55 ).
Lemma 5.2 Let ( u 0 , n ) n Z = u 0   where   is a bounded set of 2   and ( g n ) n Z = g 2   . For any 0 σ <   , and η > 0   , there exist T ( η )   and K ( η )   such that the solution ( u n ) n Z = u   of ( 54 )-( 55 ) satisfies for all t T ( η )   , the estimate
| n | > 2 M | u n ( t ) | 2 2 η δ , for any M > K ( η ) . (67)
Proof: Choose a smooth function θ C 1 ( R + , R )   satisfying the following properties
{ θ ( s ) = 0 , 0 s 1 0 θ ( s ) 1 , 1 s 2 θ ( s ) = 1 , s 2 .
and
| θ ( s ) | C 0 , s R + , (68)
for some C 0 R   . We shall use the shorthand notation θ n = θ ( | n | M )   and R e ( u n ) = u 1 , n   , I m ( u n ) = u 2 , n   We now multiply the DNLS equation with the function θ n u ¯ n   , n Z   , and we sum over all sites and keep the imaginary part. The effect of the function θ   on u n   is to cut-off all sites within a ball of radius comparable to M   and to take into account only the remote sites. This results to
1 2 d d t n Z θ n | u n | 2 + δ n Z θ n | u n | 2 n Z { ( B u 2 ) n ( B θ u 1 ) n ( B u 1 ) n ( B θ u 2 ) n } = I m n Z g n θ n u ¯ n . (69)
The nonlinear term contributes a purely real term so it drops out of the above equality. We must now estimate the remaining terms. Using the Cauchy-Schwartz inequality we find for the term containing the forcing that
I m n Z g n θ n u ¯ n δ 2 n Z θ n | u n | 2 + 1 2 δ | n | M | g n | 2 .
We now estimate the last term on the lhs of equation ( 69 ). After some algebra we obtain
I 1 = n Z { ( B u 2 ) n ( B θ u 1 ) n ( B u 1 ) n ( B θ u 2 ) n } = n Z { ( θ n + 1 θ n ) ( u 1 , n u 2 , n + 1 u 2 , n u 1 , n + 1 ) } .
We may further estimate I 1   as follows,
| I 1 | sup n Z | θ n + 1 θ n | n Z | ( u 1 , n u 2 , n + 1 u 2 , n u 1 , n + 1 ) |
C 1 M | | u | | 2 2 C 1 M ρ 1 2 ,
where C 1   depends only on C 0   . For the last inequality we applied the mean value theorem and ( 68 ) and the restriction of the solution to the absorbing ball in 2   ( ρ 1   is the radius of the absorbing ball of Lemma  5.1 ).
Therefore we obtain the following differential inequality
1 2 d d t n Z θ n | u n | 2 + δ 2 n Z θ n | u n | 2 C 1 M ρ 1 2 + 1 2 δ | n | > M | g n | 2 .
Using the Gronwall inequality we obtain the following estimate
n Z θ n | u n | 2 e δ ( t t 0 ) n Z θ n | u n ( t 0 ) | 2 + 1 δ ( 2 C 1 M ρ 1 2 + 1 δ | n | > M | g n | 2 ) ,
for t > t 0   where t 0   is the time of entry of initial data bounded in 2   , to the absorbing ball of radius ρ 1   in 2   .
Since g 2   , then for all η > 0   , there exists K ( η )   such that
2 C 1 M ρ 1 2 + 1 δ | n | > M | g n | 2 η , M > K ( η ) .
Therefore, for all η   and for t > t 0   and M > K ( η )   we obtain that
n Z θ n | u n | 2 e δ ( t t 0 ) ρ 1 2 + 1 δ η
Choosing t   large enough we may then obtain
| n | > 2 M | u n | 2 n Z θ n | u n | 2 2 η δ .
This estimate holds as long as t T ( η )   where T ( η ) = t 0 + 1 δ ln ( δ ρ 1 2 η )   and M > K ( η )   . This concludes the proof of the Lemma.   Remark The above Lemma has an additional implication, since it allows us to obtain an estimate for the time scales required for the localization of the energy to the part of the lattice with | n | M   .  
Proposition 5.1 Let 0 σ <   . The semigroup S ( t )   is asymptotically compact in 2   , that is, if the sequence φ n   is bounded in 2   and t n   , then S ( t n ) φ n   is precompact in 2   .
Proof: Let us consider a bounded sequence φ n 2   , such that | | φ n | | 2 r   , r > 0   , n = 1 , 2 , . . .   By Lemma  5.1 , there exists T ( r ) > 0   such that S ( t ) φ n 0   , t T ( r )   .
As t n   there exists an integer N 1 ( r )   such that t n T ( r )   for n N 1 ( r )   and
S ( t n ) φ n 0 , n N 1 ( r ) . (70)
Therefore S ( t n ) φ n   is weakly relatively compact. Thus, there exists φ 0 0   and a subsequence of S ( t n ) φ n   (not relabelled) such that S ( t n ) φ n φ 0   weakly in 0   . To ensure precompactness we need to show that the above convergence is strong, that is for all η > 0   , there exists N ( η ) N   such that | | S ( t n ) φ n φ 0 | | 2 η , n N ( η ) .   From ( 70 ) and the tail estimates provided by Lemma  5.2 , we may conclude the existence of some K 1 ( η )   and T 1 ( η )   such that
| i | K 1 ( η ) | ( S ( t ) S ( T r ) φ n ) i | 2 η 2 8 , t T 1 ( η ) ,
where by ( S ( t ) φ n ) i   we denote the i   -th coordinate of the infinite sequence S ( t ) φ n 2   . Since t n   we may find N 2 ( r , η ) N   such that t n T r + T ( η )   if n N 2 ( r , η )   . Hence,
| i | K 1 ( η ) | ( S ( t n ) φ n ) i | 2 = | i | K 1 ( η ) | ( S ( t n T r ) S ( T r ) φ n ) i | 2 η 2 8
where we used the semigroup property for S ( t )   .
Since φ 0 2   we have the following tail estimate | i | K 2 ( η ) | ( φ 0 ) i | 2 η 2 8 .   Let us now choose K ( η ) = m a x ( K 1 ( η ) , K 2 ( η ) )   . Since S ( t n ) φ n φ 0   in 2   , we may see that ( S ( t n ) φ n ) i ( φ 0 ) i   strongly in C 2 K ( η ) + 1   . To prove the strong convergence we will break | | | | 2   into two distinct parts, the finite dimensional part (corresponding to C 2 K ( η ) + 1   ) and the infinite dimensional part corresponding to the tails. Thus,
| | S ( t n ) φ n φ 0 | | 2 2 = | i | K ( η ) | ( S ( t n ) φ n φ 0 ) i | 2 + | i | > K ( η ) | ( S ( t n ) φ n φ 0 ) i | 2
η 2 2 + 2 | i | > K ( η ) ( | ( S ( t n ) φ n ) i | 2 + | ( φ 0 ) i | 2 ) η 2
where the first estimate comes from the strong convergence in the finite dimensional space C 2 K ( η ) + 1   .   Using Proposition  5.1 and Theorem 1.1.1 of [33, we may arrive at the main result of this section which can be stated as follows:
Theorem 5.1 Let 0 σ <   . The semigroup S ( t )   associated to ( 54 )-( 55 ) possesses a global attractor A = ω ( 0 ) 0 2   which is compact, connected and maximal among the functional invariant sets in 2   .

6 Finite dimensional approximation of the global attractor

This section is devoted to a result of stability of the global attractor of the dynamical system generated by ( 54 )-( 55 ) under its approximation by a global attractor of an appropriate finite dimensional dynamical system.
Let us note first that problem ( 54 )-( 55 ) could be viewed as a boundary value problem satisfying the boundary condition that the solutions are bounded at infinity. Since u ( t ) C 1 ( R + , 2 )   clearly lim n u n ( t ) = 0 , t 0 .   We seek for an approximation of ( 54 )-( 55 ) by a finite dimesional system of ordinary differential equations.
We consider the boundary value problem in C 2 m + 1   ,
i ψ ˙ n + ( ψ n 1 2 ψ n + ψ n + 1 ) + i δ ψ n + | ψ n | 2 σ ψ n = g n , | n | m , 0 σ < , (71)
ψ ( m + 1 ) ( ) = ψ ( m + 1 ) ( ) = 0 , (72)
ψ n ( 0 ) = ψ n , 0 , | n | m . (73)
In a similar manner to the infinite dimensional system we may show that for the well posedness and asymptotic behavior of solutions of ( 71 ) the following result holds.
Proposition 6.1 Let ψ 0 : = ( ψ n , 0 ) | n | m C 2 m + 1   . For 0 σ <   , there exists a unique solution of ( 71 )-( 73 ) such that ψ C 1 ( [ 0 , ) , C 2 m + 1 )   . The dynamical system defined by ( 71 )-( 73 )
S m ( t ) : ψ 0 C 2 m + 1 ψ ( t ) C 2 m + 1 , (74)
possesses a bounded absorbing set O 0   in C 2 m + 1   and a global attractor A m O 0 C 2 m + 1   : For every bounded set O   of 2   , there exists t 1 ( O , O 0 )   such that for all t t 1 ( O , O 0 )   , it holds that S m ( t ) O O 0   and for every t 0   S m ( t ) A m = A m   .
Note that the estimates used for the proof of Proposition  6.1 (similar to those used in the proof Lemma  5.1  are independent of m   .
Following [38, 39(see also [6, 20for a similar idea applied to pdes considered in all of R N   ), we observe that the C 2 m + 1   -solution of ( 71 )-( 73 ) can be extended naturally in the infinite dimensional space 2   , as
( φ m ( t ) ) m Z : = { ψ ( t ) : = ( ψ n ( t ) ) | n | m , | n | m , 0 , | n | > m . (75)
The extension ( 75 ) will be used in order to verify that the global attractor A   of the semigroup S ( t )   associated with ( 54 )-( 55 ), can be approximated by the global attractor A m   of S m ( t )   associated to ( 71 )-( 73 ) as m   .
We recall that the semidistance of two nonempty compact subsets of a metric space X   endowed with the metric d X ( , )   is defined as
d ( 1 , 2 ) = sup x 1 inf y 2 d X ( x , y ) .
Theorem 6.1 For any 0 σ <   , the global attractor A m   converges to A   in the sense of the semidistance related to 2   : lim m d ( A m , A ) = 0   .
Proof: We denote by U   an open-neighborhood of the absorbing ball 0   of S ( t )   . Obviously A   attracts U   . For arbitrary m Z   , we consider the semigroup S m ( t )   defined by Proposition  6.1 and its global attractor A m   . Exactly as in Lemma  5.1 , it can be shown that 0 C 2 m + 1   is also an absorbing set for S m ( t )   . Therefore A m 0 C 2 m + 1 U C 2 m + 1 ,   and A m   attracts U C 2 m + 1   . In the light of Proposition  6.1 and [33,TheoremI1.2,pg.28, it remains to verify that for every compact interval I   of R +   ,
δ m ( I ) : = sup ψ 0 U C 2 m + 1 sup t I d ( S m ( t ) ψ 0 , S ( t ) ψ 0 ) 0 , as m . (76)
We consider the corresponding solution ψ ( t ) = S m ( t ) ψ 0   , ψ ( 0 ) = ψ 0   , in C 2 m + 1   through ( 71 )-( 73 ). Then by Proposition  6.1 , it follows that ψ ( t ) A m   for any t R +   . Therefore, if ρ > 0   is the m   -independent radius of the absorbing ball O m   in C 2 m + 1   , then for every t R +   , | | ψ ( t ) | | 2 2 ρ 2   . Using ( 75 ), we clearly observe that φ m ( t )   also satisfies the estimates
| | φ m ( t ) | | 2 2 ρ 2 , | | φ ˙ m ( t ) | | 2 2 C ( ρ , | | g | | ) , (77)
the latter derived by ( 71 ). According to [33,Theorem10.1pg.331-332or [39,Lemma4,pg.60, for the justification of ( 76 ) it suffices to show that φ m ( t )   converges to a solution φ ( t )   of ( 54 )-( 55 ) in an arbitrary compact interval of R +   , and φ 0 = φ ( 0 )   in a bounded set of 2   . Let I   be an arbitrary compact interval of R +   .
From estimates ( 77 ), we may extract a subsequence φ μ   of φ m   , such that
φ μ ( t ) φ ( t ) , in 2 , as μ , for every t D , (78)
where D   denotes a countable dense subset of I   .
For any t I   we consider the sequence γ m ( t ) : = ( φ m ( t ) , z ) 2 , z 2 ,   which is diferrentiable by ( 77 ) as a function of t   and γ m ( t ) = ( φ ˙ m ( t ) , z ) 2   . By the mean-value Theorem, there exists ξ I   such that, for fixed t , s I  
| ( φ m ( t ) φ m ( s ) , z ) 2 | = | γ m ( t ) γ m ( s ) | = | ( ψ ˙ m ( t ) , z ) 2 | | t s |
sup ξ I | | φ ˙ m ( ξ ) | | 2 | | z | | 2 | t s | C | t s | .
It follows then, that there exists an m   -independent constant C 1   such that
| | φ m ( t ) φ m ( s ) | | 2 C 1 | t s | , (79)
i.e the sequence γ m   is equicontinuous. Hence by Ascoli's Theorem, it follows that the convergence ( 78 ), holds uniformly on I   as μ   . Summarizing, we obtain that for the subsequence φ μ   hold the convergence relations
φ μ φ in C ( I , 2 ) ,
φ μ * φ in L ( I , 2 ) , (80)
φ ˙ μ * φ ˙ in L ( I , 2 ) .
For every z 2   and ω ( t ) C 0 ( I )   we consider the formula (see [12,p.g.59)
I ( i φ ˙ μ ( t ) , z ) 2 ω ( t ) d t I ( B φ μ ( t ) , B z ) 2 ω ( t ) d t + δ I ( i φ μ ( t ) , z ) 2 ω ( t ) d t
+ I ( f ( φ μ ( t ) ) , z ) 2 ω ( t ) d t = I ( g , z ) 2 ω ( t ) d t , (81)
where f ( φ μ ) : = | φ μ | 2 σ φ μ   . Needless to say that any solution of ( 81 ), is in our case, a solution of ( 54 )-( 55 ) and vice versa. Using ( 80 ) we may pass to the limit in ( 81 ) . Note for example, that since f : 2 2   is Lipschitz continuous on bounded sets of 2   , there exists from ( 77 ), a constant c ( ρ )   such that | | f ( φ μ ) f ( φ ) | | 2 c ( ρ ) | | φ μ φ | | 2   . Then from ( 80 ) we infer
| I ( f ( φ μ ( t ) ) f ( φ ( t ) ) , z ) ω ( t ) d t | I | | f ( φ μ ( t ) ) f ( φ ( t ) ) | | 2 | | z | | 2 ω ( t ) d t
c I | | φ μ ( t ) φ ( t ) | | 2 | | z | | 2 ω ( t )
c sup t I | | φ μ ( t ) φ ( t ) | | 2 | | z | | 2 I | ω ( t ) | d t 0 , as μ .
Now we conclude by the same arguments as in [39,Lemma4,pg.61: Since I   is arbitrary, ( 81 ) is satisfied for all t R +   , i.e. φ ( t )   solves ( 54 )-( 55 ). Moreover by ( 80 ) we get that φ ( t )   is bounded in 2   for all t R +   . Therefore φ ( t ) A   , which implies that φ μ ( 0 ) φ ( 0 )   and φ 0   , is at least, in a bounded set of 2   . Since the convergence holds for any other subsequence having the above formulated properties, by a contradiction argument using uniqueness, we may deduce that the convergence holds for the original sequence φ m   . Condition ( 76 ) is proved.
  .

7 The case of weighted spaces: Existence of spatially exponentially localized solutions

We now turn to the properties of the solutions of the DNLS ( 54 )-( 55 ), in weighted spaces. We consider weight function w n   which is an increasing function of | n |   , satisfying for all n Z   , the following condition:
( W ) { 1 w n | w n + 1 w n | d 1 w n d ̲ 2 w n w n + 1 .
the space w 2   , w 2 = { u n C : | | u | | w 2 2 : = n Z w n | u n | 2 < } .   It can easily be seen that the space w 2   is a Hilbert space with the norm | | | | w 2   . Such spaces are the discrete analogue of weighted L 2   spaces. A choice for such a function may be the exponential function w n = e x p ( λ | n | )   for λ > 0   . Existence of solutions in such spaces will provide us with the existence of (exponentially) localized solutions for the DNLS equation ( 54 )-( 55 ). The use of such spaces is important in the study of existence of soliton solutions or breathers. An instance where such spaces have been used is in [5where the existence of exponentially localized solutions has been studied in conservative lattices using a continuation argument, related to the anti-integrable limit.
Since the operator A   is not symmetric in the space w 2   we cannot apply Theorem  2.1 for the local existence in such spaces. We thus have to resort to general existence Theorems in Banach spaces. To this end we need to show that the operators involved are Lipschitz. This time as a model case we consider ( 1 )-( 2 ) with a nonlinearity satisfying ( N 2 )   .
We have the following,
Lemma 7.1 Let condition ( N 2 )   be fulfilled. The operator T : w 2 w 2   , defined by T ( z ) = f ( | z | 2 ) z   satisfies the following properties:
  • (i) T   is bounded on bounded sets of w 2   and
  • (ii) T   is locally Lipschitz continuous.
Proof: (i) Let u B R   , a closed ball of w 2   of radius R   . It follows from condition ( W )   that | z n | 2 | | z | | w 2 2   for all n Z   . Since f   is continuous, we may argue as for the proof of Lemma  2.1 , in order to get the inequality
| | T ( z ) | | w 2 2 = n Z w n | f ( | z n | 2 ) | 2 | z n | 2 n Z w n g ( | z n | 2 ) 2 | z n | 2
n Z w n g ( | | z | | w 2 2 ) 2 | z n | 2 { max ρ [ 0 , R 2 ] g ( ρ ) } 2 n Z w n | z n | 2 c ( R ) | | z | | w 2 .
for some positive constant c ( R )   . Thus we conclude that the operator T   is bounded on bounded sets of w 2   .
(ii) Since ( N 2 )   holds, for some θ ( 0 , 1 )   we have
| | T ( z ) T ( z ) | | w 2 2 n Z w n | f ( | z n | 2 ) | 2 | z n z n | 2
+ 2 n Z w n | f ( θ | z n | 2 + ( 1 θ ) | z n | 2 ) | 2 ( | z n | + | z n | ) 2 | z n | 2 | z n z n | 2 .
A similar inequality to ( 19 ) should be obtained: We may see that
| | T ( z ) T ( z ) | | l w 2 2 2 n Z w n | f ( | | z | | w 2 2 ) | 2 | z n z n | 2
+ 2 { max ρ [ 0 , 2 R 2 ] g 1 ( ρ ) } 2 c ( R ) n Z w n | z n z n | 2
c 1 ( R ) ( | | z z | | w 2 2 .
This concludes the proof of the Lemma.  
Lemma 7.2 The operator A : w 2 w 2   , defined by ( A u ) n = u n + 1 2 u n + u n 1   is globally Lipschitz on w 2   .
Proof: Let u , v B R   . Then ( A u ) n Z ( A v ) n Z = ( u n + 1 v n + 1 ) 2 ( u n v n ) + ( u n 1 v n 1 )   and it follows that | | A u A v | | w 2 4 | | u v | | w 2   .   Using the above two Lemmas we may restate the local existence result, in the case of weighted spaces.
Theorem 7.1 Let assumption ( N 2 )   be satisfied, assume that u 0 w 2   . There exists T * ( u 0 ) > 0   , such that for all 0 < T < T * ( u 0 )   , there exists a unique solution of the problem ( 1 )-( 2 ), u ( t ) C 1 ( [ 0 , T ] , w 2 )   .
Proof: This time, we write ( 1 )-( 2 ), as an ordinary differential equation in w 2  
u ˙ ( t ) + Φ ( u ( t ) ) = 0 ,
u ( 0 ) = u 0 ,
where Φ ( u ) = T ( u ) + A ( u ) + i δ u g   and u ( t )   lies in w 2   . Lemmas  7.1 - 7.2 suffice for the application of standard existence and uniqueness Theorems for ordinary differential equations in Banach spaces [37,pg.78-82.   In order to prove the existence of global solutions we need to obtain some a priori estimates for the solution. In the dissipative case this is achieved by proving the existence of a globally attracting ball of finite radius in w 2   .
Lemma 7.3 Assume condition ( W )   on the weight function and that the damping coefficient satisfies
δ 2 2 d 1 d ̲ 2 1 / 2 0 . (82)
Let ( u 0 , n ) n Z = u 0 , ( g n ) n Z = g w 2   and condition ( N 2 )   be satisfied. A dynamical system can be defined by ( 1 )-( 2 ),
S ( t ) : u 0 w 2 u ( t ) w 2 , (83)
possessing a bounded absorbing set 0   in w 2   : For every bounded set   of w 2   , there exists t 0 ( , 0 )   such that for all t t 0 ( , 0 )   , it holds S ( t ) 0   .
Proof: We multiply ( 1 ) with w n u ¯ n   , n Z   add over all lattice sites and keep the imaginary part. Working similarly as in Lemma  5.2 , we obtain
1 2 d d t n Z w n | u n | 2 n Z ( w n + 1 w n ) ( u 1 , n u 2 , n + 1 u 2 , n u 1 , n + 1 ) + δ n Z w n | u n | 2
δ 2 n Z w n | u n | 2 + 1 2 δ n Z w n | g n | 2
Using the assumptions on the weight function ( W )   , we find that
n Z | ( w n + 1 w n ) ( u 1 , n u 2 , n + 1 u 2 , n u 1 , n + 1 ) | d 1 n Z w n | ( u 1 , n u 2 , n + 1 u 2 , n u 1 , n + 1 ) |
d 1 ( n Z w n | u 1 , n u 2 , n + 1 | + n Z w n | u 2 , n u 1 , n + 1 | ) (84)
We now estimate the sums in the above inequality as follows:
n Z w n | u 1 , n u 2 , n + 1 | = n Z w n 1 / 2 w n 1 / 2 | u 1 , n u 2 , n + 1 | ( n Z w n | u 1 , n | 2 ) 1 / 2 ( n Z w n | u 2 , n + 1 | 2 ) 1 / 2
d ̲ 2 1 / 2 ( n Z w n | u 1 , n | 2 ) 1 / 2 ( n Z w n + 1 | u 2 , n + 1 | 2 ) 1 / 2 d ̲ 2 1 / 2 | | u | | w 2 2 ,
and similarly,
n Z w n | u 2 , n u 1 , n + 1 | d ̲ 2 1 / 2 | | u | | w 2 2 .
Therefore,
n Z | ( w n + 1 w n ) ( u 1 , n u 2 , n + 1 u 2 , n u 1 , n + 1 ) | 2 d 1 d ̲ 2 1 / 2 | | u | | w 2 2 .
Thus, we obtain the differential inequality
1 2 d d t | | u | | w 2 2 + ( δ 2 2 d 1 d ̲ 2 1 / 2 ) | | u | | w 2 2 1 2 δ | | g | | w 2 2
Using the Gronwall inequality, if δ 2 2 d 1 d ̲ 2 1 / 2 0   we obtain the existence of the attracting ball in w 2   .
This concludes the proof of the Proposition.  
Remark 7.1 In the case of exponential weight w n = e x p ( λ n )   , it appears that d 1 = e x p ( λ ) 1   and d ̲ 2 = e x p ( λ )   . The condition for existence of attracting ball becomes 8 s i n h ( λ / 2 ) δ   .
We may further prove the existence of a global attractor attracting all bounded sets of w 2   . The next Lemma provides us with tail estimates in the weighted space.
Lemma 7.4 Let ( u 0 , n ) n Z = u 0   where   is a bounded set of w 2   and ( g n ) n Z = g w 2   . Under assumptions ( N 2 )   on the nonlinearity, ( W )   on the weight function and assumption ( 82 ) on the dissipation, for any η > 0   , there exist T ( η )   and K ( η )   such that the solution ( u n ) n Z = u   of ( 1 )-( 2 ) satisfies for all t T ( η )   , the estimate
| n | > 2 M w n | u n | 2 2 η δ , for any M > K ( η ) . (85)
Proof: The proof follows closely the proof of Lemma  5.2 for the corresponding tail estimates in the space 2   , only that we now multiply the equation by w n θ n u ¯ n   , instead of multiplying simply by θ n u ¯ n   . We obtain the following inequality
1 2 d d t n Z w n θ n | u n | 2 I 2 + δ n Z w n θ n | u n | 2 δ 2 n Z w n θ n | u n | 2 + 1 2 δ n Z w n θ n | g n | 2 , (86)
where I 2 = n Z { ( θ n + 1 w n + 1 θ n w n ) ( u 1 , n u 2 , n + 1 u 2 , n u 1 , n + 1 ) } .   We rewrite ( θ n + 1 w n + 1 θ n w n ) = ( θ n + 1 θ n ) w n + 1 + θ n ( w n + 1 w n ) .   Using arguments similar to those used in the proof of Lemma  5.2 and Lemma  7.3 we see that
n Z | θ n + 1 θ n | | w n + 1 | | u 1 , n u 2 , n + 1 u 2 , n u 1 , n + 1 | C 2 ρ 2 2 M , (87)
where C 2   depends only on d 1 , d ̲ 2   and ρ 2   denotes the absorbing ball in w 2   . We also have the inequality
n Z θ n | w n + 1 w n | | u 1 , n u 2 , n + 1 u 2 , n u 1 , n + 1 | d 1 n Z θ n w n | u 1 , n u 2 , n + 1 | + d 1 n Z θ n w n | u 2 , n u 1 , n + 1 | . (88)
For the first term of the rhs of ( 88 ), holds the estimate
d 1 n Z θ n w n | u 1 , n u 2 , n + 1 | d 1 ( n Z θ n w n | u 1 , n | 2 ) 1 / 2 ( n Z θ n w n | u 2 , n + 1 | 2 ) 1 / 2
d 1 d ̲ 2 1 / 2 ( n Z θ n w n | u 1 , n | 2 ) 1 / 2 ( n Z θ n w n + 1 | u 2 , n + 1 | 2 ) 1 / 2 . (89)
An application of Young's inequality to the rhs of ( 89 ), implies that
d 1 n Z θ n w n | u 1 , n u 2 , n + 1 | 1 2 d 1 d ̲ 2 1 / 2 n Z θ n w n | u 1 , n | 2 + d 1 2 d ̲ 2 1 2 d 1 d ̲ 2 1 / 2 n Z θ n w n + 1 | u 2 , n + 1 | 2
1 2 d 1 d ̲ 2 1 / 2 n Z θ n w n | u 1 , n | 2 + 1 2 d 1 d ̲ 2 1 / 2 n Z θ n + 1 w n + 1 | u 2 , n + 1 | 2
+ 1 2 d 1 d ̲ 2 1 / 2 n Z ( θ n θ n + 1 ) w n + 1 | u 2 , n + 1 | 2
d 1 d ̲ 2 1 / 2 n Z θ n w n | u 1 , n | 2 + C 3 ρ 2 2 M , (90)
and similarly for the second term of the rhs of ( 88 )
d 1 n Z θ n w n | u 2 , n u 1 , n + 1 | d 1 d ̲ 2 1 / 2 n Z θ n w n | u 2 , n | 2 + C 4 ρ 2 2 M . (91)
Using ( 86 )-( 91 ) we derive that
1 2 d d t n Z w n θ n | u n | 2 + ( δ 2 2 d 1 d ̲ 2 1 / 2 ) n Z w n θ n | u n | 2 + C ρ 2 2 M + 1 2 δ n Z w n θ n | g n | 2 .
with C   depending only on d 1 , d ̲ 2   . The rest follows by direct generalization of the arguments of Lemma  5.2 .   In complete analogy as before we may prove the asymptotic compactness of the semigroup in the weighted spaces,
Proposition 7.1 The semigroup S ( t )   is asymptotically compact in w 2   , that is if the sequence φ n   is bounded in w 2   and t n   then S ( t n ) φ n   is precompact in w 2   .
This result leads to the existence of a global attractor in the weighted spaces.
Theorem 7.2 Let conditions ( N 2 ) , ( W )   and ( 82 ) be fulfilled. The semigroup S ( t )   associated to ( 1 )-( 2 ) possesses a global attractor A = ω ( 0 ) 0 w 2   which is compact, connected and maximal among the functional invariant sets in w 2   .

8 Remarks on possible extensions

We conclude, by mentioning some other examples of DNLS-type equations for which, extensions of the results on the global solvability and the existence of global attractors, could be investigated.
A first example is provided by the DNLS with potential [35, 36, i.e ( 1 )-( 2 ) with a nonlinearity of the form F ( u ) n Z = F ( u n ) + V ( n ) u   , where V : R R   is a real valued potential | V ( s ) | c   for every s R   . The potential function expresses the inhomogeneity properties of the medium.
Another interesting DNLS model, may be given by the spatial discretization of a modified NLS equation with viscocity
i u t + 1 2 u x x | u | 2 u = i δ u x x , δ > 0 (92)
which describes the light propagation in an array of optical fibers in a weakly lossy medium. Taking into account some external excitation and considering a general nonlinearity, the discretized counterpart of ( 92 ) reads as (see [30)
i u ˙ n + ( 1 i δ ) ( u n + 1 2 u n + u n 1 ) + f ( | u n | 2 ) u n = g n . (93)
Typical DNLS-models include equations of the form ([17, 25)
i U ˙ n + ( α n + 1 U n + 1 + β n + 1 U n + 1 ) + i δ n U n + γ n f ( | U n | 2 ) U n = g n .
For the coupling and nonlinear strength and dissipation parameter α n , β n , γ n , δ n R   we assume that there exist some constants c 1 , c 2 , c 3   such that
| α n | , | β n | , | γ n | c 1 , c 2 δ n c 3 , n Z .
The validity of the results could be examined for DNLS-type equations considered in Z N   , N > 1   ,
i u ˙ n + ( A u ) n + i δ n u n + f ( | u n | 2 ) u n = g n , (94)
u n ( 0 ) = u n , 0 , n = ( n 1 , n 2 , . . . n N ) Z N . (95)
where dissipation satisfies ( 94 ) and the linear operator A   has a decomposition introduced in ([39): It is assumed that
A = A 1 + A 2 + . . . + A N (96)
and
A j = B j * B j = B j * B j , | | B j | | ( 2 ) M , j = 1 , 2 , . . . , N , (97)
for some bounded linear operators B j : 2 2   , and its adjoint B j *   , defined as
( B j u ) n = l = m 0 l = m 0 C j , l u n j l , for all u = ( u n ) n Z N , j = 1 , 2 , . . . N (98)
( B j * u ) n = l = m 0 l = m 0 C j , l u n j l , n j l = ( n 1 , n 2 , . . . , n j 1 , n j + l , n j + 1 , . . . , n N ) Z N . (99)
Under assumptions ( 96 )-( 98 ) similar relations to ( 8 )-( 9 ) are satisfied,
( B u , v ) 2 = ( u , B * v ) 2 , ( A u , v ) 2 = ( B u , B v ) 2 , u , v 2 .
Especially for the conservative case, δ n = 0   , n Z N   , and the power-law nonlinearity, ( 96 )-( 98 ) are sufficient for the extension of Theorem  3.1 to ( 94 )-( 95 ). Although this is a simple observation, it is important, since it is a rigorous verification that spatial dimension does not play a particular role on the global existence of the solutions as it happens in the continuous model [7, [11, 12. For a discussion on the existence of global attractors for ( 94 )-( 95 ), one could be based on the framework provided by [39.
Acknowledgements We would like also to thank Professors D. Frantzeskakis, K. Hizanidis, and N. Stavrakakis for helpful discussions.
References

  1. A. Ambrosetti and P. H. Rabinowitz, Dual Variational Methods in Critical Point Theory and Applications, J. Funct. Anal. 14, (1973), 349-381.
  2. V. S. Afraimovich, S.N. Chow, J. K. Hale, Synchronization in lattices of coupled oscillators, Phys. D 103 (1997), no. 1-4, 442–451.
  3. G. Akrivis, Finite Difference Discretization of the cubic Schrödinger equation, IMA J. Numer. Anal., 13, (1993), 115-124.
  4. S. Aubry, Breathers in nonlinear lattices: Existence, linear stability and quantization, Physica D 103 (1997), 201-250.
  5. C. Baesens and R. S. MacKay, Exponential localization of linear response in networks with exponentially decaying coupling, Nonlinearity 10, (1997) 931-940
  6. A. V. Babin and M. I. Vishik, Attractors for Partial Differential Evolution Equations in an Unbounded Domain, Proc. Roy. Soc. Edinb., 116A, (1990), 221-243.
  7. O. Bang, J. Rasmussen, P. Christiansen, Subcritical localization in the discrete nonlinear Schrödinger equation with arbitrary nonlinearity, Nonlinearity 7, (1994), 205-218.
  8. P. W. Bates, K. Lu, B. Wang, Attractors for lattice dynamical systems, Internat. J. Bifur. Chaos Appl. Sci. Engrg. 11 (2001), no. 1, 143–153.
  9. P. W. Bates, A. Chmaj, A Discrete Convolution Model for Phase Transitions, Arch. Rational Mech. Anal. 150 (1999) 281-305.
  10. K. J. Blow, N. J. Doran, Global and Local Chaos in the Pumped Nonlinear Schrödinger Equation, Physical Review Letters, Vol. 52, No. 7, (1984), 526-539.
  11. T. Cazenave, A. Haraux, Introduction to Semilinear Evolution Equations, Oxford Lecture Series in Mathematics and its Applications 13, 1998.
  12. T. Cazenave, An introduction to Nonlinear Schrödinger equations, Textos des Métodos Matemáticos 26, I.M.U.F.R.J., Rio de Janeiro, 1996.
  13. S. N. Chow, J. K. Hale Methods of Bifurcation Theory, Grundlehren der mathematischen Wissenschaften-A series of Comprehensive Studies in Mathematics 251, Springer-Verlag, New-York, 1982.
  14. S. Flach and C. R. Willis Discrete Breathers Phys. Reports 295 (1998) 181-264
  15. J M Ghidaglia, Finite Dimensional Behavior for Weakly Damped Driven Schrödinger Equation, Ann. Inst. Henri Poincaré, 5, (1988), 365-405.
  16. O. Goubet, Regularity of the attractor for a weakly damped nonlinear Schrodinger equation, Appl. Anal. 60 (1996), no. 1-2, 99–119.
  17. B. C. Gupta, K. Kundu, Localised states in a 1-D nonlinear chain, Phys. Lett. A 235, (1997), 176-182.
  18. D. Hennig and G. P. Tsironis, Wave transmission in nonlinear lattices, Phys. Rep. 307 (1999) 333-432
  19. F. Hirsch, G. Lacombe, Elements of Functional Analysis, Graduate Texts in Mathematics 192, Springer-Verlag, New-York, 1999.
  20. N. Karachalios and N. Stavrakakis, Existence of Global Attractors for Semilinear Dissipative Wave Equations on R N   , J. Diff. Eqns, 157 (1999), 183-205.
  21. N. Karachalios and N. Stavrakakis, Global Attractor for the Weakly Damped Driven Schrödinger Equation in H 2 ( R )   , NoDEA Nonlinear Differential Equations Appl., Birkhäuser, 9 (2002), 347-360.
  22. P. Laurençot, Long Time Behavior for Weakly Damped Driven Nonlinear Schrödinger Equation in R N , N 3   , NoDEA Nonlinear Differential Equations Appl., Birkhäuser, 2 (1995), 357-369.
  23. P. G. Kevrekidis, K. O. Rasmussen and A. R. Bishop, The discrete nonlinear Schrödinger equation: A survey of recent results, Int. Journal of Modern Physics B, 15 (2001), 2833-2900
  24. P. G. Kevrekidis, B. A. Malomed, A. R. Bishop and D. J. Frantzeskakis, Localized vortices with a semi-integer charge in nonlinear dynamical lattices, Phys. Rev. E 65 (2002) 016605
  25. S. W. Kim, S. Kim, The structure of eigenmodes and phonon scattering by discrete breathers in the discrete nonlinear Schrödinger chain, Phys. D 141, (2000), 91-103.
  26. R. S. Mackay and S. Aubry, Proof of existence of breathers for time reversible or hamiltonian networks of weakly coupled oscillators Nonlinearity 7 (1994) 1623
  27. B. A. Malomed, P. G. Kevrekidis, D. J. Frantzeskakis, H. E. Nistazakis and A. N. Yannacopoulos, 1 and 2-D solitons in second harmonic generating lattices, Phys. Rev. E 65 (2002) 056606
  28. Y. Martel, Blow-up for the nonlinear equation in nonisotropic spaces, Nonlinear Analysis, TMA 28, no. 12, (1997), 1903-1908.
  29. K. Nozaki, N. Bekki, Low Dimensional Chaos in a Driven Damped Nonlinear Schrödinger Equation, Physica D 21, (1986), 381-393.
  30. M. Salerno, B. M. Malomed, V. Konotop, Shock wave dynamics in a discrete nonlinear Schrödinger equation with internal losses, Phys. Review E 62, No.6, (2000), 8651-8656.
  31. M. Struwe, Variational Methods-Applications to Nonlinear Partial Differential Equations and Hamiltonian Systems (2nd edition), A Series of Modern Surveys in Mathematics, Vol. 34, Springer-Verlag, (1996).
  32. C. Sulem, P. L. Sulem, The Nonlinear Schrödinger Equation: Self-Focusing and Wave Collapse, Applied Mathematical Sciences 139, Springer-Verlag, New-York, 1999.
  33. R. Temam, Infinite-Dimensional Dynamical Systems in Mechanics and Physics, 2nd edition, Springer-Verlag, New York, 1997.
  34. X. Wang, An Energy Equation for the Weakly Damped Driven Nonlinear Schrödinger equation and its application to their Attractors, Physica D 88 (1995), 167-175.
  35. A. N. Yannacopoulos, D. J. Frantzeskakis, C. Polymilis and K. Hizanidis, Conditions for soliton trapping in random potentials using Lyapunov exponents of stochastic ODEs, Phys. Lett. A 271 (2000) pp. 334-340
  36. A. N. Yannacopoulos, D. J. Frantzeskakis, C. Polymilis and K. Hizanidis, Motion of Schrödinger solitary waves in the presence of random external potentials, Physica Scr. 65 (2002) pp. 363-368
  37. E. Zeidler, Nonlinear Functional Analysis and its Applications, Vols I, II, (Fixed Point Theorems, Monotone Operators), Springer-Verlag, Berlin, 1990.
  38. S. Zhou, Attractors for Second order Dynamical Systems, J. Diff. Eqns 179, (2002), 605-624.
  39. S. Zhou, Attractors for first order dissipative lattice dynamical systems, Physica D 178 (2003), 51-61.