ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu

ATTRACTORS FOR DISSIPATIVE LATTICE DIFFERENTIAL EQUATIONS WITH LOCAL AND NONLOCAL

NONLINEARITIES

JARDEL MORAIS PEREIRA

Abstract. We study the dynamics of some dissipative lattice differential equations with local and nonlocal nonlinearities. Using a difference inequality due to Nakao [14] and suitable estimates, we prove the existence of global at- tractors. In addition, we briefly discuss the dynamics of two periodic lattice differential equations.

1. Introduction

Lattice differential equations (LDEs) are used to model various problems that occur in important areas of science and technology, for instance, in biophysics [13], electrical engineering [4], image processing [5], chemical reaction theory [10], etc.. They also naturally arise as spatial discretization of continuous models. The dynamics of nonlinear LDEs is a wide-ranging theme that includes as subjects the existence of solutions, existence and stability of traveling waves, asymptotic behavior, attractors and their properties, etc., see e.g. [1, 2, 3, 8, 20, 25] and the references therein. In particular, the existence of attractors for LDEs is a subject that attracts a great deal of attention. In this article, we study the dynamics of some dissipative LDEs with local and nonlocal nonlinearities in arbitrary spatial dimensions. Our main objective is to prove the existence of global attractors.

Firstly, we consider the following class of second order LDEs

¨

u_{n}(t) + (−1)^{p}∆^{p}_{d}u_{n}(t) +αu_{n}(t) +F(n, u_{n}(t),∇^{+}u_{n}(t)) +g(n,u˙_{n}(t)) =f_{n},
un(0) =u0,n, u˙n(0) =u1,n, (1.1)
where n ∈ Z^{d} and t ∈ R^{+}. In (1.1), α is a positive constant, p is any positive
integer, ∆^{p}_{d} = ∆d◦ · · · ◦∆d, p times, and ∆d denotes the d-dimensional discrete
Laplacian operator defined by ∆_{d}u_{n}=Pd

i=1(u_{n+e}_{i}+u_{n−e}_{i}−2un), where{ei}^{d}_{i=1}is
the canonical basis ofR^{n}. We assume that the nonlinear termF(n, u_{n}(t),∇^{+}u_{n}(t))
has the form

F(n, un(t),∇^{+}un(t)) =h0(n, un(t))−

d

X

i=1

∂^{−}_{i} hi(∂_{i}^{+}un(t)), (1.2)

2010Mathematics Subject Classification. 37L30, 37L60, 39A12.

Key words and phrases. Global attractors; dissipative lattice differential equations;

local and nonlocal nonlinearities; periodic lattice differential equations; difference inequality.

c

2021. This work is licensed under a CC BY 4.0 license.

Submitted June 15, 2020. Published July 19, 2021.

1

whereh0:Z^{d}×R→R,hi:R→R,i= 1, . . . , d, are functions satisfying appropriate
assumptions stated in Section 2 and∂_{i}^{+}un =un+e_{i}−un and ∂_{i}^{−}un =un−un−e_{i},
withi= 1, . . . , d. Moreover,g(n,u˙n(t)) is a nonlinear dissipation.

Many papers in the literature deal with the existence of attractors for LDEs.

Here, we just mention some works that are related to our model (1.1). The existence
of attractor and finite-dimensional approximations for the casep = 1,h_{i} ≡0 for
i = 1, . . . , d, was studied by Zhou [33]. The existence of a global attractor for
the class (1.1) when d= 1 and h_{i} ≡ 0 for i = 1, . . . , d, was first investigated by
Oliveira and Pereira [18]. Later, for this last class with a suitable delay term, the
existence of a pullback attractor was established by Wang and Bai [29]. Another
paper related to (1.1) is [32] in which the authors studied the existence and the
upper semi-continuity of attractors for second order lattices with retarded terms in
dimension one (d= 1) whenp= 1. In addition, the investigation of the existence
of attractors for second order LDEs in weighted spaces was considered in [7] by
Han and some contributions concerning non-autonomous and stochastic LDEs are
provided in [6] and [30], respectively.

An important special class of LDEs included in (1.1) is

¨

u_{n}(t) + (−1)^{p}∆^{p}_{d}u_{n}(t) +αu_{n}(t) +h_{0}(n, u_{n}(t)) +g_{0}( ˙u_{n}(t)) =f_{n}. (1.3)
Models of type (1.3) when p = 1 are known as discrete nonlinear Klein-Gordon
models and appear in different physical contexts, see e.g. [21, 23]. When p= 2,
equation (1.3) can be regarded as discrete versions of beam equations. An example
of dissipative term for which our results apply is g0 ∈ C^{1}(R;R), g0(0) = 0, and
g_{0}^{0}(s)≥c0 >0 for alls∈R. Observe that we do not demand any relation on the
parametersαandc0, andg^{0}_{0}(s) need not be bounded above.

The inclusion of the term−Pd

i=1∂_{i}^{−}hi(∂_{i}^{+}un) in (1.2) was motivated by some
studies on the dynamics of nonlinear LDEs in periodic spaces [16, 17, 19], and
continuous models of beam equations studied in the literature, see e.g. [12, 24, 31].

Note that if we choosep= 2,h_{i}(s) =|s|^{q−2}s,q≥2,i= 1, . . . , d, in (1.1), then we
obtain a discrete version of the beam equation

utt+ ∆^{2}u+αu+h0(x, u)−∆qu+g(x, ut) =f(x),
where x ∈ R^{d}, u = u(x, t) and ∆_{q}u = Pd

i=1

∂

∂x_{i} |_{∂x}^{∂u}

i|^{q−2}_{∂x}^{∂u}

i

is the usual q- Laplacian operator.

Secondly, we study the existence of global attractors for the following class of dissipative LDEs with a nonlocal nonlinearity

¨

u_{n}(t) + (−1)^{p}∆^{p}_{d}u_{n}(t) +αu_{n}(t) +F(u_{n}(t)) +g(n,u˙_{n}(t)) =f_{n},

un(0) =u0,n, u˙n(0) =u1,n, (1.4) where

F(un(t)) =h^{0}(un(t)) X

m∈Z^{d}

V(n−m)h(um(t)). (1.5)
In (1.4) and (1.5), g:Z^{d}×R→R, h: R→R^{+}, and V :Z^{d} →R^{+} are functions
satisfying suitable conditions stated in Section 2. To the best of our knowledge, the
existence of attractor for LDEs of type (1.4) with a nonlocal nonlinearity as in (1.5)
has not been considered before. The continuous convolution term corresponding
to (1.5) is known as Hartree-type nonlinearity and appears in Schr¨odinger and
Klein-Gordon equations, see [22, 27].

In addition to the problems described above we also briefly discuss the existence of global attractors for two periodic LDEs. We first consider a periodic problem for the LDE in (1.1) in arbitrary spatial dimension, then we consider the following one-dimensional periodic LDE withα= 0 and nonlinear dissipation

¨

un(t) + (−1)^{p}∆^{p}_{1}un(t)−∂^{−}_{1}h(∂_{1}^{+}un(t))−∂_{1}^{−}g(∂^{+}_{1}u˙n(t)) =fn, (1.6)
whereh:R→Randg:R→Rsatisfy appropriate conditions stated in Section 4.

Note that when g(s) =µ s, withµ >0, we obtain the strong dampingµ∆_{1}u˙_{n}(t).

Whenp= 2,g(s) =µs, andh(s) =α_{1}s^{5}−α_{2}s^{3}−α_{3}s, withα_{i}>0, i= 1,2,3, model
(1.6) can be regarded as a discrete version of the beam equation studied by Racke
and Shang in [24]. However, unlike the continuous model treated in [24], here, the
global attractor we obtain exists in the whole phase space where the dynamics is
considered.

Our main purpose in this paper is to prove the existence of global attractors
for the LDEs (1.1) and (1.4) by combining the use of a difference inequality due
to Nakao [14] with a method to derive “tail estimates of solutions” introduced by
Wang in [28]. As far as we know, in the context of the dynamics of discrete models,
this approach was only used for LDEs of type (1.1) in dimension one in [18] when
hi(s)≡0,i= 1, . . . , d, under conditions on the nonlinear termsh0(n, s) andg(n, s)
more restrictive than those used in this paper. We note that our assumptions do
not require any growth condition on the dissipative termg(n, s). In particular, our
results apply to cases with nonlinear dissipative terms such asg(n, s) =a_{n}(s^{3}+s)
org(n, s) =b_{n}sinhs, when the real constantsa_{n} andb_{n} are suitably chosen.

This paper is organized as follows. In Section 2, we state the assumptions on the
functionsh_{0},h_{i},i= 1, . . . , d, in (1.2),gin (1.1), andhandV in (1.5) that we need
to prove the existence of solutions and global attractors. Then, after introducing
some notation, we briefly discuss the global well-posedness of problems (1.1) and
(1.4). In Section 3, we establish the existence of global attractors for the semi-
groups generated by the solutions of (1.1) and (1.4). We first prove the existence of
absorbing sets, then we prove the asymptotic compactness of the semigroups. The
proofs are based on Nakao’s method [14] and suitable estimates. In Section 4, we
show how some arguments used in Section 3 can be adapted to prove the existence
of global attractors for the periodic problems described above. In the case of model
(1.6), we also used a Poincar´e inequality valid for the periodic space where the
problem is considered. Our results, in particular, generalize and complement the
studies of [16, 17, 18]. Finally, in appendix A, we present examples of functions
that satisfy some assumptions used in this paper.

2. Existence of solutions

In this section, we briefly discuss the existence of solutions for the initial value
problems (1.1) and (1.4). We begin establishing the appropriate assumptions on
the functions in (1.1), (1.2), and (1.5) and introducing some notation. We denote
by`^{p} the space of real sequencesu= (un)_{n∈Z}d such thatkuk`^{p}<∞, where

kuk_{`}p= X

n∈Z^{d}

|u_{n}|^{p}1/p

, if 1≤p <∞,
kuk`^{∞} = sup

n∈Z^{d}

|un|, ifp=∞.

Whenp= 2, `^{2} is a Hilbert space with the inner product
(u, v)_{`}2= X

n∈Z^{d}

unvn, u, v∈`^{2}.

In this case, we denote byk · kthe corresponding norm. Also, to simplify notation,
we denote a sequence (un)_{n∈Z}d by (un).

We recall that for the`^{p} spaces the following embedding relation holds:

`^{q} ⊂`^{p}, kuk`^{p}≤ kuk`^{q}, 1≤q≤p≤ ∞.

For the functionsh_{0}:Z^{d}×R→R,h_{i}:R→R, i= 1, . . . , d, andg:Z^{d}×R→R,
we assume that

(A1) For eachs0>0, there exist positive constantsLj=Lj(s0),j= 1,2,3, such that

(i) |h0(n, s1)−h0(n, s2)| ≤L1|s1−s2|,
(ii) |g(n, s1)−g(n, s2)| ≤L2|s1−s2|,
(iii) |hi(s_{1})−h_{i}(s_{2})| ≤L_{3}|s1−s_{2}|,

for all s_{1}, s_{2} in R, |s1| ≤s_{0}, |s2| ≤s_{0}, for all n∈Z^{d} andi = 1, . . . , d. In
addition,h_{0}(n,0) = 0,g(n,0) = 0,h_{i}(0) = 0, for alln∈Z^{d}andi= 1, . . . , d.

(A2) There exist sequences of nonnegative real numbers b_{1} = (b_{1,n}) ∈`^{1}, b_{2} =
(b_{2,n})∈`^{1}, and a positive constantk_{1} such that

sh0(n, s) +b1,n≥k1(˜h0(n, s) +b2,n)≥0, ∀s∈Randn∈Z^{d},
where ˜h0(n, s) =Rs

0 h0(n, σ)dσ.

(A3) There exist positive constantsk0,isuch that

sh_{i}(s)≥k_{0,i}˜h_{i}(s), ∀s∈Randi= 1, . . . , d, where ˜h_{i}(s) =
Z s

0

h_{i}(σ)dσ≥0.

(A4) There exist constantsk2>0,r≥0 and a positive integern0 satisfying
sg(n, s)≥k2|s|^{r+2}, if|n|0≤n0,

sg(n, s)≥k_{2}|s|^{2}, if|n|_{0}> n_{0},

for alls∈R, where|n|_{0}= max_{1≤i≤d}|n_{i}|, ifn= (n_{1}, . . . , n_{d}).

Regarding the LDE (1.1) with the nonlocal term (1.4), we assume thatgsatisfies
the same hypotheses above and thath:R→R^{+} andV :Z^{d} →R^{+} satisfy

(A5) h ∈ C^{2}(R;R^{+}), h(0) = h^{0}(0) = 0 and there exists a positive constant c_{1}
such thatsh^{0}(s)≥c_{1}h(s)≥0 for alls∈R.

(A6) V = (V(m))∈`^{2} andV(m) =V(−m) for allm∈Z^{d}.
The following notation will be useful.

∇^{+}un = (∂_{1}^{+}un, . . . , ∂^{+}_{d}un), ∇^{−}un= (∂_{1}^{−}un, . . . , ∂_{d}^{−}un),

∇^{+}u_{n}· ∇^{+}v_{n}=

d

X

i=1

∂_{i}^{+}u_{n}∂_{i}^{+}v_{n}, |∇^{+}u_{n}|^{2}=∇^{+}u_{n}· ∇^{+}u_{n},

D^{p}u_{n} =

(∆^{p/2}_{d} un, ifpis even

∇^{+} ∆

p−1 2

d un

, ifpis odd,

(2.1)
where ∆^{0}_{d}=I.

Lemma 2.1. If u= (u_{n})∈`^{2} thenP

n∈Z^{d}|D^{p}u_{n}|^{2}≤(4d)^{p}kuk^{2}.

Proof. Using the elementary inequality Pd
i=1a_{i}2

≤dPd

i=1a^{2}_{i},a_{i} ∈R, and argu-
ing by induction we find that

X

n∈Z^{d}

|∆^{k}_{d}un|^{2}≤(4d)^{2k} X

n∈Z^{d}

|un|^{2},∀k∈N. (2.2)
In (2.2) and hereafter N denotes the set of all positive integers. If pis odd then,
from (2.1) and (2.2), we have

X

n∈Z^{d}

|D^{p}u_{n}|^{2}= X

n∈Z^{d}

|∇^{+}(∆

p−1 2

d u_{n})|^{2}= X

n∈Z^{d}
d

X

i=1

|∂_{i}^{+}∆

p−1 2

d u_{n}|^{2}

≤4 X

n∈Z^{d}
d

X

i=1

|∆

p−1 2

d un|^{2}= 4d X

n∈Z^{d}

|∆

p−1 2

d un|^{2}

≤4d(4d)^{p−1}kuk^{2}= (4d)^{p}kuk^{2}.

Similarly, using (2.1) and (2.2) we can treat the case when p is even.

Lemma 2.2. For any u= (un)andv= (vn)in`^{2} we have
(−1)^{p} X

n∈Z^{d}

(∆^{p}_{d}un)vn =
(P

n∈Z^{d}∆^{p/2}_{d} un∆^{p/2}_{d} vn, if pis even
P

n∈Z^{d}∇^{+}(∆

p−1 2

d un)· ∇^{+}(∆

p−1 2

d vn), if pis odd.

Proof. Sinceu= (un) andv= (vn) belong to`^{2}, we have
X

n∈Z^{d}

(∆dun)vn=

d

X

i=1

X

n∈Z^{d}

(∂_{i}^{+}un)vn−

d

X

i=1

X

n∈Z^{d}

(∂_{i}^{−}un)vn

=

d

X

i=1

X

n∈Z^{d}

(∂_{i}^{+}u_{n})v_{n}−

d

X

i=1

X

n∈Z^{d}

(∂_{i}^{+}u_{n})v_{n+e}_{i}

=− X

n∈Z^{d}
d

X

i=1

∂_{i}^{+}un∂^{+}_{i} vn=−X

n∈Z^{d}

∇^{+}un· ∇^{+}vn.

This proves Lemma 2.2 ifp= 1. The general case follows by induction onp.

In what follows, given a sequenceu= (un), we will write
h0(u) = (h0(n, un)), g(u) = (g(n, un)),
Au= ((−1)^{p}∆^{p}_{d}un), B(u) =

−

d

X

i=1

∂_{i}^{−}hi(∂_{i}^{+}un)

. (2.3)

Also, we will use the Hilbert space H = `^{2}×`^{2} equipped with the usual inner
product and norm,

((u, v),(w, z))H = (u, w)_{`}2+ (v, z)_{`}2 and k(u, v)kH= kuk^{2}+kvk^{2}^{1/2}
,
for any (u, v) and (w, z) inH.

Lemma 2.3. Under assumption (A1), we have

(i) h0, B, andg are locally Lipschitz continuous maps from `^{2} into itself.

(ii) A:`^{2}→`^{2} is a bounded operator and kAuk ≤(4d)^{p/2}kuk for allu∈`^{2}.

Proof. (i) Givenu= (un) in`^{2}, using (A1), we have
kB(u)k^{2}= X

n∈Z^{d}

d

X

i=1

∂_{i}^{−}hi(∂_{i}^{+}un)

2

≤d X

n∈Z^{d}
d

X

i=1

∂_{i}^{−}hi(∂_{i}^{+}un)

2

≤d X

n∈Z^{d}
d

X

i=1

|hi(∂_{i}^{+}un)−hi(∂_{i}^{−}un)|^{2}

≤dL_{3}(2kuk)^{2} X

n∈Z^{d}
d

X

i=1

|∂_{i}^{+}u_{n}−∂_{i}^{−}u_{n}|^{2}

≤(4d)^{2}L_{3}(2kuk)^{2}kuk^{2}<∞.

If u= (un) and v = (vn) belong to `^{2}, with kuk ≤ R and kvk ≤ R, then, using
(A1) again, we have

kB(u)−B(v)k^{2}≤d X

n∈Z^{d}
d

X

i=1

|∂_{i}^{−}h_{i}(∂_{i}^{+}u_{n})−∂_{i}^{−}h_{i}(∂^{+}_{i} v_{n})|^{2}

= X

n∈Z^{d}
d

X

i=1

|hi(∂_{i}^{+}un)−hi(∂_{i}^{+}vn)−[hi(∂_{i}^{−}un)−hi(∂_{i}^{−}vn)]|^{2}

≤2dL3(2R)^{2} X

n∈Z^{d}
d

X

i=1

|∂_{i}^{+}(un−vn)|^{2}+|∂_{i}^{−}(un−vn)|^{2}

≤(4d)^{2}L3(2R)^{2}ku−vk^{2}.

This shows thatB is a locally Lipschitz continuous map from`^{2} into itself. Simi-
larly, we prove thath_{0}, g:`^{2}→`^{2}are locally Lipschitz continuous maps.

(ii) It follows immediately from (2.2).

Using Lemma 2.3 we can write the initial value problem (1.1) in`^{2} as

¨

u(t) +Au(t) +αu(t) +h0(u(t)) +B(u(t)) +g( ˙u(t)) =f, t >0,

u(0) =u_{0}, u(0) =˙ u_{1}, (2.4)
where u0 = (u0,n), u1 = (u1,n), f = (fn), u(t) = (un(t)), ˙u(t) = ( ˙un(t)) and

¨

u(t) = (¨un(t)).

Theorem 2.4. Assume that (A1)–(A3) hold and let u0, u1 and f belong to `^{2}.
Assume also that sg(n, s) ≥ 0 for all n ∈ Z^{d} and s ∈ R. Then the initial value
problem (2.4) has a unique solution u ∈ C^{2}(R^{+};`^{2}). Moreover, for each τ > 0,
the map I:H →C([0, τ];H), defined byI(u0, u1)(t) = (u(t),u(t)),˙ 0 ≤t ≤τ, is
continuous.

Proof. Introducing the change of variable ˙u= v we can rewrite problem (2.4) in the spaceH as

dw

dt(t) +Bw(t) = 0, t >0, w(0) =w0,

wherew= (u, v), ^{dw}_{dt} = ( ˙u,v),˙ w0= (u0, u1), and

B(w) = (−v, Au+αu+h_{0}(u) +B(u) +g(v)−f). (2.5)

Using Lemma 2.3 we can easily prove that the mapBdefined by (2.5) is a locally
Lipschitz continuous map from H into itself. Then, an application of the Theory
of Ordinary Differential Equations in Banach Spaces shows that the initial value
problem (2.4) has a unique solutionu∈C^{2}([0, τ_{max});`^{2}) such that eitherτ_{max}=∞
orτ_{max}<∞and lim_{t→τ}−

maxk(u(t),u(t))k˙ H =∞.

To extend the solution globally we proceed as follows. Taking the inner product
of equation (2.4) with ˙u(t) in`^{2} and using Lemma 2.2 we find

d

dtE(t) =−(g( ˙u(t)),u(t))˙ `^{2} ≤0, ∀0≤t < τmax, (2.6)
whereE(t) is the energy associated with the initial value problem (2.4) given by

E(t) = 1

2ku(t)k˙ ^{2}+1

2kD^{p}u(t)k^{2}+α

2ku(t)k^{2}+ X

n∈Z^{d}

˜h0(n, un(t))

+ X

n∈Z^{d}
d

X

i=1

˜hi(∂_{i}^{+}un)−(f, u(t))_{`}2.

(2.7)

In (2.7) and hereafter kD^{p}uk^{2} = P

n∈Z^{d}|D^{p}u_{n}|^{2}. Thus, E(t) ≤ E(0), for all
0 ≤t < τmax. Since |(f, u)_{`}2| ≤ ^{α}_{4}kuk^{2}+ _{α}^{4}kfk^{2}, using (A2), (A3) and (2.7) we
deduce that

k(u(t),u(t))k˙ ^{2}_{H} ≤α^{−1}_{0} E(t),˜ ∀0≤t < τmax, (2.8)
whereα_{0}= min{^{1}_{2},^{α}_{4}} and

E(t) =˜ E(t) + 4

αkfk^{2}+kb2k_{`}1, ∀0≤t < τmax. (2.9)
From (2.8) and (2.9) we conclude thatτmax =∞. Finally, under the assumptions
of Theorem 2.4, the continuity of I can be proved using (2.8) and the Gronwall
inequality. Since the arguments are well known, we omit the details here.

Now, let us consider the initial value problem (1.4) with the nonlocal term (1.5).

By assumptions (A5) and (A6) we can define the mapF:`^{2}→`^{2} by
F(u) =

h^{0}(u_{n}) X

m∈Z^{d}

V(n−m)h(u_{m})

, ∀u= (u_{n})∈`^{2}.
Then, using the above notation, we can write (1.4) in the space`^{2} as

¨

u(t) +Au(t) +αu(t) +F(u(t)) +g( ˙u(t)) =f, t >0,

u(0) =u0, u(0) =˙ u1, (2.10)
To see thatFis a locally Lipschitz continuous map from`^{2}into iself, letu= (un)
and v = (vn) in `^{2} such that kuk ≤ R and kvk ≤ R. Since h ∈ C^{2}(R;R^{+}),
h(0) = h^{0}(0) = 0 by (A5) and V = (V(m)) ∈ `^{2} by (A6) and |un| ≤ R and

|vn| ≤R, for alln∈Z^{d}, then
kF(u)−F(v)k^{2}≤2 X

n∈Z^{d}

|h^{0}(un)−h^{0}(vn)|^{2} X

m∈Z^{d}

V(n−m)|h(um)|^{2}
+ 2 X

n∈Z^{d}

|h^{0}(v_{n})|^{2} X

m∈Z^{d}

V(n−m)|h(u_{m})−h(v_{m})|2

≤2M_{1}^{2}M_{2}^{2}kVk^{2}kuk^{2}ku−vk^{2}+ 2M_{1}^{2}M_{2}^{2}kVk^{2}kvk^{2}ku−vk^{2}

≤4R^{2}M_{1}^{2}M_{2}^{2}kVk^{2}ku−vk^{2},

whereM1= max_{|s|≤2R}|h^{0}(s)|andM2= max_{|s|≤2R}|h^{00}(s)|.

Proceeding as we did for the initial value problem (1.1), we can prove the fol- lowing results.

Theorem 2.5. Assume (A1)-(ii), with g(n,0) = 0 for all n∈Z^{d},(A5), (A6)and
let u0, u1, and f belong to `^{2}. Assume also that sg(n, s) ≥0 for all n ∈Z^{d} and
s∈R. Then the initial value problem (2.10)has a unique solutionu∈C^{2}(R^{+};`^{2}).

Moreover, for each τ >0, the mapI:H →C([0, τ];H), defined by I(u0, u1)(t) = (u(t),u(t)),˙ 0≤t≤τ, is continuous.

In this case, as before, we obtain the identity (2.6) with the energy function E(t) =1

2ku(t)k˙ ^{2}+1

2kD^{p}u(t)k^{2}+α
2ku(t)k^{2}
+1

2 X

n∈Z^{d}

X

m∈Z^{d}

V(n−m)h(um(t))h(un(t))−(f, u(t))_{`}2.

(2.11)

The same inequality in (2.8) can be derived with E(t) =˜ E(t) + 4

αkfk^{2}, ∀0≤t <∞. (2.12)
3. Existence of global attractors

Our aim in this section is to prove the existence of global attractors for the semigroups generated by the solutions of the initial value problems (2.4) and (2.10).

Let us first consider the initial value problem (2.4). Using Theorem 2.4 we can define
a semigroup of continuous operators{S(t)}_{t≥0} onH as follows

S(t)(u_{0}, u_{1}) = (u(t),u(t)),˙ ∀(u_{0}, u_{1})∈H. (3.1)
To prove the existence of a global attractor for {S(t)}t≥0 in H it is sufficient
to prove that {S(t)}t≥0 has an absorbing set in H and that it is asymptotically
compact in H, see e.g. [26]. Our proofs are based on a difference inequality by
Nakao as stated in Lemma 3.1 below. This difference inequality was introduced
in [14] to study the existence of attractors for some nonlinear wave equations with
nonlinear dissipation. Some other applications to the study of the dynamics of
continuous models can be seen in [9, 11, 15]. In the context of LDEs, it was used
in [16, 17, 18].

Lemma 3.1(Nakao [14]). Letψ(t)be a nonnegative continuous function on[0, T), T >1, possibly T =∞, satisfying

sup

t≤s≤t+1

ψ(s)^{1+γ} ≤C[ψ(t)−ψ(t+ 1)] +K, ∀0≤t < T −1, (3.2)
with some C >0,K >0 andγ >0. Then

ψ(t)≤

C^{−1}γ(t−1)^{+}+ sup

0≤s≤1

ψ(s)−γ−_{γ}^{1}

+K^{γ+1}^{1} , 0≤t < T.

If (3.2)holds with γ= 0, then ψ(t)≤ sup

0≤s≤1

ψ(s) C C+ 1

^{[t]}

+K, 0≤t < T,

where[t]is the largest integer less than or equal to t andβ^{+}= max{β,0}.

Let us now introduce some notation that will be used in this section. We first observe that if (A4) holds then, in particular, we have

sg(n, s)≥0,∀n∈Z^{d} ands∈R.
Let

P(t)^{2}=E(t)−E(t+ 1) and Q(t) = X

n∈Z^{d}

˙

un(t)g(n,u˙n(t)), ∀t≥0.

According to (2.6),E(t) is a non-increasing function on [0,∞) and
P(t)^{2}=

Z t+1

t

Q(s)ds, ∀t≥0. (3.3)

Also, by the Mean Value Theorem for integrals, there exist real numbers t1 ∈
[t, t+^{1}_{4}] and t2∈[t+^{3}_{4}, t+ 1] such that

Z t+^{1}_{4}

t

Q(s)ds=1

4Q(t1) and

Z t+1

t+^{3}_{4}

Q(s)ds= 1

4Q(t2). (3.4)
To simplify notation we will write ˆh(u) =h_{0}(u) +B(u) for allu∈`^{2}.

Lemma 3.2. Assume that(A1)–(A4)hold and letu_{0},u_{1}, andf belong to`^{2}. Then
there exists a positive constant C_{0} such that

sup

t≤s≤t+1

E(s)˜ ≤C_{0} P(t)^{r+2}^{4} +P(t)^{4}+P(t)^{2}+kfk^{2}+kb1k`^{1}

,∀t≥0. (3.5)
Proof. Taking the inner product of 2.4 withu=u(t) in`^{2}, using Lemma 2.2, and
integrating the result over [t1, t2] we obtain

Z t_{2}

t1

kD^{p}u(s)k^{2}ds+α
Z t_{2}

t1

ku(s)k^{2}ds+
Z t_{2}

t1

(ˆh(u(s)), u(s))_{`}2ds

= ( ˙u(t1), u(t1))_{`}2−( ˙u(t2), u(t2))_{`}2+
Z t2

t1

ku(s)k˙ ^{2}ds

− Z t2

t_{1}

(g( ˙u(s)), u(s))_{`}2ds+
Z t2

t_{1}

(f, u(s))_{`}2ds.

(3.6)

Let us estimate the terms in the right hand side of (3.6). We initially write
( ˙u(tj), u(tj))_{`}2 = X

|n|0≤n0

˙

un(tj)un(tj) + X

|n|0>n0

˙

un(tj)un(tj), j = 1,2, (3.7) withn0 as in (A4). Using (2.8) we have

X

|n|0≤n0

|un(tj)|^{2}1/2

≤α^{−1/2}_{0} sup

t≤s≤t+1

E(s)˜ ^{1/2}. (3.8)
Using H¨older’s inequality, (A4), (3.3), and (3.4) we have

X

|n|0≤n0

|u˙n(tj)|^{2}≤(2n0+ 1)^{r+2}^{rd} X

|n|0≤n0

|u˙n(tj)|^{r+2}_{r+2}^{2}

≤(2n0+ 1)^{r+2}^{rd}

k_{2}^{−1} X

|n|0≤n0

g(n,u˙n(tj)) ˙un(tj)_{r+2}^{2}

≤(2n_{0}+ 1)^{r+2}^{rd} k^{−}

2 r+2

2 4^{r+2}^{2} P(t)^{r+2}^{4} .

(3.9)

Similarly, we can estimate the second term in (3.7) to find X

|n|0>n_{0}

|u˙_{n}(t_{j})u_{n}(t_{j})| ≤2α_{0}^{−1/2}k_{2}^{−1/2}P(t) sup

t≤s≤t+1

E(s)˜ ^{1/2}. (3.10)
From (3.7)-(3.10) we conclude that there is a positive constant C0,1 depending
only onn0,α,k2,r, anddsuch that

|(u(t_{1}),u(t˙ _{1}))_{`}2−(u(t_{2}),u(t˙ _{2}))_{`}2| ≤C_{0,1}

P(t)^{r+2}^{2} +P(t)
sup

t≤s≤t+1

E(s)˜ ^{1/2}. (3.11)
Proceeding as in (3.9), replacingtj bys∈[t1, t2], and noticing that

X

|n|0>n0

|u˙_{n}(s)|^{2}≤k_{2}^{−1}Q(s),
for anys∈[t_{1}, t_{2}] by (A4), we obtain

Z t_{2}

t_{1}

ku(s)k˙ ^{2}ds≤(2n0+ 1)^{r+2}^{rd} k^{−}

2 r+2

2

Z t_{2}

t_{1}

Q(s)^{r+2}^{2} ds+k_{2}^{−1}
Z t_{2}

t_{1}

Q(s)ds

≤(2n0+ 1)^{r+2}^{rd} k^{−}

2 r+2

2 P(t)^{r+2}^{4} +k^{−1}_{2} P(t)^{2}.

(3.12)

Next, for eacht≥0 fixed, we define the sets
I_{1}(t) =

n∈Z^{d}; |u˙_{n}(t)| ≤1 , I_{2}(t) =Z^{d}\I1(t).

Note that by (A1),|g(n,u˙_{n}(s))| ≤L_{2}(1)|u˙_{n}(s)|, whenevern∈I_{1}(s). Then
X

n∈I1(s)

un(s)g(n,u˙n(s))≤α

2ku(s)k^{2}+ 1
2α

X

n∈I1(s)

|g(n,u˙n(s))|^{2}

≤α

2ku(s)k^{2}+ 1

2αL_{2}(1) X

n∈I1(s)

|u˙_{n}(s)| |g(n,u˙_{n}(s))|

=α

2ku(s)k^{2}+ 1

2αL2(1) X

n∈I1(s)

˙

un(s)g(n,u˙n(s))

≤α

2ku(s)k^{2}+ 1

2αL_{2}(1)Q(s).

(3.13)

In addition, using (A4) and (2.8) we have X

n∈I2(s)

un(s)g(n,u˙n(s))≤ X

n∈I2(s)

|un(s)| |u˙n(s)| |g(n,u˙n(s))|

≤ ku(s)k X

n∈I_{2}(s)

|u˙_{n}(s)| |g(n,u˙_{n}(s))|

=ku(s)k X

n∈I2(s)

˙

un(s)g(n,u˙n(s))

≤ ku(s)kQ(s)≤α^{−1/2}_{0} Q(s) sup

t≤s≤t+1

E(s)˜ ^{1/2}.

(3.14)

It follows from (3.13), (3.14), and (3.3) that
Z t_{2}

t_{1}

(g( ˙u(s), u(s))`^{2}ds≤ α
2

Z t_{2}

t_{1}

ku(s)k^{2}ds+ 1

2αL2(1)P(t)^{2}
+α^{−1/2}_{0} P(t)^{2} sup

t≤s≤t+1

E(s)˜ ^{1/2}.

(3.15)

Finally, in view of (2.8), we easily see that

Z t_{2}

t_{1}

(f, u(s))_{`}2ds

≤α^{−1/2}_{0} kfk sup

t≤s≤t+1

E(s)˜ ^{1/2}. (3.16)
Substituting (3.11), (3.12), (3.15), and (3.16) into (3.6) we obtain the estimate

Z t_{2}

t1

kD^{p}u(s)k^{2}ds+α
2

Z t_{2}

t1

ku(s)k^{2}ds+
Z t_{2}

t1

J1(s)ds

≤C_{0,2}h

P(t)^{r+2}^{2} +P(t)^{2}+P(t) +kfk

× sup

t≤s≤t+1

E(s)˜ ^{1/2}+P(t)^{r+2}^{4} +P(t)^{2}i

+kb1k_{`}1,

(3.17)

whereC0,2 is a positive constant depending only on n0, α, k2, r, d, L2(1), and
J_{1}(t) = X

n∈Z^{d}

u_{n}(t)h_{0}(n, u_{n}(t)) +b_{1,n}

+ X

n∈Z^{d}
d

X

i=1

∂^{+}_{i} u_{n}(t)h_{i}(∂_{i}^{+}u_{n}(t))≥0,
for allt≥0, because of (A2) and (A3).

On the other hand, using hypotheses (A2) and (A3) we have X

n∈Z^{d}

˜h_{0}(n, u_{n}(t)) + X

n∈Z^{d}
d

X

i=1

h˜_{i}(∂_{i}^{+}u_{n}(t))

≤k^{−1}_{1} X

n∈Z^{d}

un(t)h0(n, un(t)) +b1,n

− kb2k`^{1}

+ X

n∈Z^{d}
d

X

i=1

k_{0,i}^{−1}

∂^{+}_{i} un(t)hi(∂_{i}^{+}un(t))

≤k0

h X

n∈Z^{d}

un(t)h0(n, un(t)) +b1,n

+ X

n∈Z^{d}
d

X

i=1

∂_{i}^{+}un(t)hi(∂_{i}^{+}un(t))i

− kb2k_{`}1,

(3.18)

wherek0= max{k^{−1}_{1} , k^{−1}_{0,i}, i= 1, . . . , d}.

Then, integrating (2.9) over [t_{1}, t_{2}] and using (2.7), (3.17), (3.12), and (3.18) we
deduce that

Z t_{2}

t_{1}

E(s)ds˜

=
Z t_{2}

t_{1}

E(s)ds+ 4 α

Z t_{2}

t_{1}

kfk^{2}ds+
Z t_{2}

t_{1}

kb2k_{`}1ds

≤max{k0,1}hZ t_{2}
t1

kD^{p}u(s)k^{2}ds+α
2

Z t_{2}

t1

ku(s)k^{2}ds+
Z t_{2}

t1

J1(s)dsi +1

2 Z t2

t1

ku(s)k˙ ^{2}ds+ 4
α

Z t2

t1

kfk^{2}ds

≤C_{0,3}

P(t)^{r+2}^{2} +P(t)^{2}+P(t) +kfk
sup

t≤s≤t+1

E(s)˜ ^{1/2}+P(t)^{r+2}^{4} +P(t)^{2}
+ 4

αkfk^{2}+ max{k_{0},1}kb_{1}k_{`}1, (3.19)

for some positive constantC0,3 depending only onn0, α, k0, k2, r, d, andL2(1).

Furthermore, by the Mean Value Theorem for integrals, there existst^{∗}∈[t1, t2]
such that

1 2

E(t˜ ^{∗})≤(t2−t1) ˜E(t^{∗}) =
Z t_{2}

t1

E(s)˜ ds. (3.20)

If ˜t∈[t, t^{∗}] then using (2.6) and (3.3) we obtain
E(˜˜ t) = ˜E(t^{∗}) +

Z t^{∗}

˜t

Q(s)ds≤E(t˜ ^{∗}) +P(t)^{2}.
Analogously, we obtain the same estimate if ˜t∈[t^{∗}, t+ 1]. Therefore,

E(s)˜ ≤E(t˜ ^{∗}) +P(t)^{2}, ∀s∈[t, t+ 1].

Using this fact, (3.19) and (3.20), we complete the proof.

In what follows we will writeK_{0}=kfk^{2}+kb1k`^{1} and assume thatK_{0}>0.

Lemma 3.3. Under the assumptions of Lemma 3.2 there exists ρ0>0 such that
B[0;ρ0] ={(w, z)∈H;k(w, z)kH ≤ρ0} is an absorbing set for{S(t)}t≥0 inH.
Proof. Let O be any bounded subset of H and let ρ = ρ(O) be a positive con-
stant such that k(w, z)kH ≤ ρ,∀(w, z) ∈ O. Assume that (u0, u1) ∈ O. Let us
first consider the case r > 0. Since ˜E(t) is a non-increasing function, then using
(A1)-(A3) and Lemma 2.1 we can find a positive constant µ0 depending only on
p, d, α, ρ, k1, k0,i,i= 1, . . . , d,kfk,kb1k`^{1}, andkb2k`^{1} such that

P(t)^{2}=E(t)−E(t+ 1) = ˜E(t)−E(t˜ + 1)≤2 ˜E(0)≤µ_{0}, ∀t≥0. (3.21)
It follows from (3.5) and (3.21) that

sup

t≤s≤t+1

E(s)˜ ^{1+}^{r}^{2} ≤C1,1[ ˜E(t)−E(t˜ + 1)] + 2C0K0

1+^{r}_{2}

, (3.22)

where C_{1,1} is a positive constant depending only on r, C_{0}, and µ_{0}. Applying the
first part of Lemma 3.1 to (3.22) with ψ(t) = ˜E(t), γ = ^{r}_{2}, C =C1,1, and K =

2C0K0

1+^{r}_{2}

we obtain E(t)˜ ≤h

C_{1,1}^{−1}r

2(t−1)^{+}+
sup

0≤s≤1

E(s)˜ −r/2i−2/r

+ 2C0K0. (3.23) From (3.23), using (3.21), we deduce that

E(t)˜ ≤C1(1 +t)^{−2/r}+ 2C0K0, (3.24)
whereC1 is a positive constant depending only onr, C0, andµ0. Combining (3.24)
with (2.8) we obtain

kS(t)(u0, u1)k^{2}_{H} ≤α^{−1}_{0} C1(1 +t)^{−}^{2}^{r} + 2α^{−1}_{0} C0K0 ∀t≥0. (3.25)
Consequently,kS(t)(u0, u1)kH ≤ρ0, for allt≥τ if we take

ρ_{0}= 2C0K0

α_{0}
1/2

, τ=τ(O) = maxn

0, C1

2C_{0}K_{0}
r/2

−1o .

This completes the proof of Lemma 3.3 if r >0. The proof for the case r = 0 is analogous. Indeed, using (3.5) and (3.21) we have

sup

t≤s≤t+1

E(s)˜ ≤C_{0}(2 +µ_{0})[ ˜E(t)−E(t˜ + 1)] +C_{0}K_{0},∀t≥0,

with C0 and µ0 as before. Then, applying the second part of Lemma 3.1 to this inequality and using (2.8) again, we deduce that

kS(t)(u0, u_{1})k^{2}_{H}≤α^{−1}_{0} µ_{0}
2

C_{2}
C2+ 1

^{[t]}

+α^{−1}_{0} C_{0}K_{0}

≤α^{−1}_{0} µ_{0}
2

C_{2}+ 1
C2

e^{−νt}+α^{−1}_{0} C_{0}K_{0}, ∀t≥0,

(3.26)

where C2 =C0(2 +µ0) andν = ln

C2+1
C_{2}

. This implies Lemma 3.3 in the case
r= 0 with ρ0= ^{2C}_{α}^{0}^{K}^{0}

0

^{1/2}

andτ =τ(O) = max

0,_{ν}^{1}ln ^{µ}_{2C}^{0}^{(C}^{2}^{+1)}

2C_{0}K_{0} .

Next, we prove that the semigroup {S(t)}_{t≥0} is asymptotically compact in H.
Here, the main step consists in using a method introduced by B. Wang in [28]

combined with Lemma 3.1 to derive an appropriate estimate for the “tail” of the solution of (2.4). More precisely, we will show that, for all >0, there existτ()>0 and a positive integerk() such that

X

|n|0≥k()

[( ˙u_{n}(t))^{2}+ (u_{n}(t))^{2}]< , for allt≥τ(),

whenever the initial data (u0, u1) belongs to the absorbing set B[0;ρ0]. To do this, we will need the following auxiliary lemma whose proof relies on the following elementary identities valid for any sequences w = (wn) and z = (zn) and i = 1, . . . , d.

∂_{i}^{+}(wnzn) = (∂_{i}^{+}wn)zn+e_{i}+wn∂^{+}_{i} zn, (3.27)

∂_{i}^{+}(wnzn) = (∂_{i}^{+}wn)zn+wn∂_{i}^{+}zn+∂_{i}^{+}wn∂_{i}^{+}zn, (3.28)
partial^{−}_{i} (wnzn) = (∂_{i}^{−}wn)zn+wn∂_{i}^{−}zn−∂_{i}^{−}wn∂_{i}^{−}zn, (3.29)

∆d(wnzn) = (∆dwn)zn+wn∆dzn+∇^{+}wn· ∇^{+}zn+∇^{−}wn· ∇^{−}zn. (3.30)
Note that (3.30) follows from (3.28) and (3.29).

Lemma 3.4. Let u= (un(t))belong toC^{1}(R^{+};`^{2}), and(θn) belong to`^{2}. Then
(−1)^{p} X

n∈Z^{d}

∆^{p}_{d}un(t) (θnu˙n(t)) = 1
2

d dt

X

n∈Z^{d}

θn|D^{p}un(t)|^{2}+ X

n∈Z^{d}
d

X

i=1

(∂_{i}^{+}θn)z_{p,n}^{(i)}(t),
where

X

n∈Z^{d}
d

X

i=1

|z^{(i)}_{p,n}(t)| ≤C(p, d)k(u(t),u(t))k˙ ^{2}_{H}, ∀t≥0,
for some positive constantC(p, d)depending onpandd.

Proof. A proof of this lemma whend= 1 was first presented in [18]. We will argue by induction. Consider first podd. Ifp= 1 then using Lemma 2.2 and (3.27) we have

− X

n∈Z^{d}

∆dun(θnu˙n) = X

n∈Z^{d}

∇^{+}un· ∇^{+}(θnu˙n)

= X

n∈Z^{d}
d

X

i=1

∂_{i}^{+}u_{n} θ_{n}∂^{+}_{i} u˙_{n}+∂_{i}^{+}θ_{n}u˙_{n+e}_{i}

= 1 2

d dt

X

n∈Z^{d}

θn|∇^{+}un|^{2}+ X

n∈Z^{d}
d

X

i=1

∂_{i}^{+}θnz_{1,n}^{(i)},

where z^{(i)}_{1,n} = (∂_{i}^{+}un) ˙un+ei satisfies P

n∈Z^{d}

Pd

i=1|z_{1,n}^{(i)}| ≤ 2dk(u,u)k˙ ^{2}_{H} by Lemma
2.1.

Assume now that Lemma 3.4 holds forp= 2k−1 with k ∈N. Using Lemma 2.2 again and (3.30) we have

(−1)^{2k+1} X

n∈Z^{d}

∆^{2k+1}_{d} u_{n}(θ_{n}u˙_{n})

=−X

n∈Z^{d}

∆^{2k}_{d} u_{n}∆_{d}(θ_{n}u˙_{n})

=−X

n∈Z^{d}

∆^{2k}_{d} un

(∆dθn) ˙un+θn∆du˙n+∇^{+}θn· ∇^{+}u˙n+∇^{−}θn· ∇^{−}u˙n
.

(3.31)

By Lemma 2.2 and (3.28) we also see that

− X

n∈Z^{d}

∆^{2k}_{d} un(∆dθn) ˙un

= X

n∈Z^{d}
d

X

i=1

∂_{i}^{+}θn(∆^{2k}_{d} un)∂_{i}^{+}u˙n+ X

n∈Z^{d}
d

X

i=1

∂_{i}^{+}θn(∂_{i}^{+}∆^{2k}_{d} un) ˙un

+ X

n∈Z^{d}
d

X

i=1

∂^{+}_{i} θ_{n}(∂_{i}^{+}∆^{2k}_{d} u_{n})∂_{i}^{+}u˙_{n}.

(3.32)

In addition,

− X

n∈Z^{d}

∆^{2k}_{d} un∇^{−}θn· ∇^{−}u˙n

=− X

n∈Z^{d}
d

X

i=1

∆^{2k}_{d} u_{n+e}_{i}∂_{i}^{+}θ_{n}∂_{i}^{+}u˙_{n}

=− X

n∈Z^{d}
d

X

i=1

∂_{i}^{+}∆^{2k}_{d} un∂_{i}^{+}θn∂^{+}_{i} u˙n− X

n∈Z^{d}
d

X

i=1

∆^{2k}_{d} un∂_{i}^{+}θn∂_{i}^{+}u˙n.

(3.33)

Substituting (3.32) and (3.33) into (3.31) results
(−1)^{2k+1} X

n∈Z^{d}

∆^{2k+1}_{d} un(θnu˙n)

=− X

n∈Z^{d}

∆^{2k}_{d} un(θn∆du˙n)

+ X

n∈Z^{d}
d

X

i=1

∂_{i}^{+}θ_{n}

(∂_{i}^{+}∆^{2k}_{d} u_{n}) ˙u_{n}−∆^{2k}_{d} u_{n}(∂_{i}^{+}u˙_{n})
.

(3.34)