Ka Ïs Ammari
email: kais.ammari@fsm.rnu.tn

AND Fathi Hassine
email: fathi.hassine@fsm.rnu.tn

Luc Robbiano
email: luc.robbiano@uvsq.fr

STABILIZATION OF FRACTIONAL EVOLUTION SYSTEMS WITH MEMORY

Keywords: Mathematics Subject Classification. 35A01, 35A02, 35M33, 93D20 stabilization, abstract-wave equation, fractional-evolution

This paper is devoted to the analysis of the problem of stabilization of fractional (in time) partial differential equations. We consider the following equation

with the initial data u(0) = u 0 , where A is a unbounded operator in Hilbert space and ∂ α,η t stands for the fractional derivative. We provide two main results concerning the behavior of the solutions when t -→ +∞. We look first to the case η > 0 where we prove that the solution of this problem is exponential stable then we consider the case η = 0 when we prove under some consideration on the resolvent that the energy of the solution goes to 0 as t goes to the infinity as 1/t α .

For three centuries the theory of fractional derivatives developed mainly as a pure theoretical field of mathematics useful only for mathematicians. However, in the last few decades many authors pointed out that derivatives and integrals of non-integer order are very suitable for the description of properties of various real materials, e.g. polymers. It has been shown that new fractional-order models are more adequate than previously used integer-order models.

Fractional derivatives provide an excellent instrument for the description of memory and hereditary properties of various materials and processes. This is the main advantage of fractional

Introduction

Integer-order derivatives and integrals have clear physical interpretation and are used for describing different concepts in classical physics. For example, the position of a moving object can be represented as a function of time, the object velocity is then the first derivative of the function, the acceleration is the second derivative and so on. Fractional derivatives and integrals, being generalization of the classical derivative and integrals are expected to have even broader meaning. unfortunately, there is no such result in the literature until now. derivatives in comparison with classical integer-order models, in which such effects are in fact neglected. The advantages of fractional derivatives become apparent in modeling mechanical and electrical properties of real materials, as well as in the description of rheological properties of rocks, and in many other fields.

Fractional integrals and derivatives also appear in the theory of control of dynamical systems, when the controlled system or/and the controller is described by a fractional differential equation [START_REF] Ammari | Fractional-feedback stabilization for a class of evolution systems[END_REF][START_REF] Cuesta | Asymptotic behaviour of the solutions of fractional integro-differential equations and some time discretizations[END_REF][START_REF] Qi | On the lack of controllability of fractional in time ODE and PDE[END_REF].

Fractional calculus includes various extensions of the usual definition of derivative from integer to real order, including the Riemann-Liouville derivative, the Caputo derivative, the Riesz derivative, the Weyl derivative, etc. In this paper,we only consider the Caputo derivative that leads to an initial condition which is physically meaningful [START_REF] Samko | Fractional integrals and derivatives, Theory and applications[END_REF].

These models are relevant, in particular, in the context of spatially disordered systems, porous media, fractal media, turbulent fluids and plasmas, biological media with traps, binding sites or macro-molecular crowding, stock price movements, etc. We refer the readers to [START_REF] Ben-Avraham | Diffusion and reactions in fractals and disordered systems[END_REF][START_REF] Bouchaud | Anomalous diffusion in disordered media: statistical mechanisms, models and physical applications[END_REF][START_REF] Metzler | The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics[END_REF] and the rich references therein for the motivation and description of the model. On the other hand, we refer to [START_REF] Caffarelli | Regularity theory for parabolic nonlinear integral operators[END_REF][START_REF] Caffarelli | Drift diffusion equations with fractional diffusion and the quasigeostrophic equation[END_REF][START_REF] Caffarelli | Regularity estimates for the solution and the free boundary of the obstacle problem for the fractional Laplacian[END_REF][START_REF] Caffarelli | An extension problem related to the fractional Laplacian[END_REF][START_REF] Constantin | Diffusion and mixing in fluid flow[END_REF][START_REF] Kiselev | Global well-posedness for the critical 2D dissipative quasigeostrophic equation[END_REF] and the rich references therein for mathematical analysis of these models.

Let H be a Hilbert space equipped with the norm . H , and let A : D(A) ⊂ H -→ H be a closed and densely defined operator on H. We consider the following Cauchy problem described by the mean of the fractional derivative as follow:

(1.1)    ∂ α,η t u(t) = Au(t) - η Γ(1 -α) t 0 (t -s) -α e -η(t-s) u(s) ds, t > 0, u(0) = u 0 ,
where ∂ α,η t denoted the fractional derivative defined by

(1.2) ∂ α,η t v(t) = 1 Γ(1 -α) t 0 (t -s) -α e -η(t-s) v (s) ds, 0 < α < 1, η ≥ 0.
The main result of this paper concerns the precise asymptotic behavior of the solutions of (1.1). In [START_REF] Cuesta | Asymptotic behaviour of the solutions of fractional integro-differential equations and some time discretizations[END_REF], the author consider a fractional integro-differential equation equivalent in some sense to system (1.1) but for the case 1 < α < 2 and η = 0 where he shows how the asymptotic behavior of the continuous solution depends on some parameter ω where it's assumed that A is a sectoral operator with a sector depends on ω. Precisely, if ω ≥ 0, then the continuous solutions are bounded by an exponential of type e ω 1 α t and if ω < 0, then the solutions show a merely algebraic decay of order o 1 ωt α . In this work we prove under some consideration on the resolvent behavior (weaker then the one considered in [START_REF] Cuesta | Asymptotic behaviour of the solutions of fractional integro-differential equations and some time discretizations[END_REF]) that the second kind of behavior still true if η = 0 and one shows also that the energy is exponentially stable if η > 0. This paper is organized as follows: In section 2 we prove the well-posedness of system (1.1) and gives an exponential stability result in the case η > 0. In section 3 we prove that the energy of system (1.1) is polynomially stable for the case η = 0 while in section 4 we consider an integro-differential equation where we prove the well-posedness of the equation and a polynomial decay rate of the energy.

Well-posedness and exponential stabilization

We define the convolution product of a and u by

a * u(t) = t 0 a(t -s)u(s) ds, ∀ t ∈ R + , a ∈ L 1 loc (R + ), u ∈ L p loc (R + , H), p ∈ [1, +∞[,
and for β > 0 the functions g β are given by

g β (t) =    1 Γ(β) t β-1 t > 0, 0 t ≤ 0.
Noting that these functions satisfy the semigroup property, namely (2.1)

g β * g γ (t) = g β+γ (t), ∀ t > 0, β, γ > 0.
The Riemann-Liouville fractional integral of order 0 < α < 1 is defined as follow

J α u(t) = g α * u(t), ∀ u ∈ L 1 loc (R + ; H), t > 0.
Using (2.1) we follow that J α verifying the semigroup property, (2.2)

J β J γ u(t) = J β+γ u(t), ∀ t > 0, β, γ > 0.
For every u ∈ L 1 loc (R + ; H) such that g 1-α * u ∈ W 1,1 loc (R + ; H), the Riemann-Liouville fractional derivative of order α is defined by

D α t u(t) = (g 1-α * u) (t) = (J 1-α u) (t), ∀ t > 0. The operator D α
t is the left inverse and right invertible of J α (see [3, Theorem 1.5]) more precisely, we have

(2.3) D α t J α u(t) = u(t), ∀ u ∈ L 1 loc (R + ; H), t > 0, and for every u ∈ L 1 loc (R + ; H) such that g 1-α * u ∈ W 1,1 loc (R + ; H), (2.4) J α D α t u(t) = u(t) ∀ t > 0.
We recall that the Caputo fractional derivative of order α > 0 when η = 0 is defined by

D α t u(t) = J 1-α u (t), ∀ u ∈ W 1,1 loc (R + ;
H), t > 0, then by integration by parts we follow that when u ∈ W 1,1 loc (R + ; H), we have

D α t u(t) = D α t (u -u(0))(t), ∀ t > 0.
Hence, the Caputo derivative D α t is a left inverse of J α but in general it is not a right inverse, namely using (2.3) and (2.4) we have

(2.5) D α t J α u(t) = u(t), ∀ u ∈ L 1 loc (R + ), t > 0, and (2.6) J α D α t u(t) = u(t) -u(0), ∀ u ∈ C(R + ; H), g 1-α * (u -u(0)) ∈ W 1,1 loc (R + ; H), t > 0. By setting v(t) = u(t)e ηt , (1.1) is equivalent to the following problem (2.7) D α t v(t) = Av(t), t > 0, v(0) = u 0 .
Applying J α in both sides of the first line of (2.7) and using (2.5) and (2.6), we conclude that when

u ∈ C(R + ; H) satisfying g 1-α * (v -u 0 ) ∈ W 1,1 loc (R + ; H) (2.7) is equivalent to the following integral differential equation (2.8) u(t) = e -ηt u 0 + e -ηt (g α * (e η. Au))(t), ∀ t ≥ 0.
The well-posedness of a system such as (1.1) is related to the notion of what is called solution operator defined as follow:

Definition 2.1. A family (S α,η (t)) t≥0 ∈ L(H) (denoted simply by (S α (t)) t≥0 if η = 0)
S α,η (t)x = e -ηt x + e -ηt (g α * (e η. AS α,η (.)x))(t), ∀ x ∈ D(A), t ≥ 0. Definition 2.2. A function u ∈ C(R + , H) is called a strong solution of (2.8) if u ∈ C(R + , D(A))
and (2.8) holds on R + .

Definition 2.3. The problem (2.8) is called well-posed if for any u 0 ∈ D(A), there is a unique strong solution u(t, u 0 ) of (2.8), and when (2.10)

u 0 n ∈ D(A) such that u 0 n -→ 0 as n +∞ then u(t, u 0 n ) -→ 0 as n +∞ in H, uniformly on compact intervals. Definition 2.4. A function u ∈ C(R + , H) is called a strong solution of (1.1) if u ∈ C(R + , D(A)), g 1-α * ((e η. u) -u 0 ) ∈ C 1 (R + , H)
v(t) = u 0 + g α * Av(t), t ≥ 0,
where we denoted by v(t) = e ηt u(t). Following to [22, Proposition 1.1], (2.10) is wellposed if and only if it admits a solution operator Sα (t). The result follow easily by observing that S α,η (t) = Sα (t)e -ηt is a solution operator of (2.8) and then (2.9) holds. 2/ The first implication follows easily from the definitions. Now suppose that (2.8) is wellposed and let u its solution and let v as above solution of (2.10). To prove the result we only have to prove that g 1-α * (v -u 0 ) ∈ C 1 (R + , H). Since v ∈ C(R + , D(A)) then the convolution of (2.10) with g 1-α gives

g 1-α * (v -u 0 )(t) = A t 0 v(s) ds, t ≥ 0,
where we have used (2.1). Then it is easy to show that g

1-α * (v -u 0 ) ∈ C 1 (R + , H).
This completes the proof.

Definition 2.6. The solution operator S α (t) (η = 0) is called exponentially bounded if there are M ≥ 1 and ω ≥ 0 such that

(2.11) S α (t) L(H) ≤ M e ωt .
The operator A is said to belong to C α (M, ω) if the problem (1.1) has a solution operator S α (t)

satisfying (2.11). Denote C α (ω) = ∪ {C α (M, ω) : M ≥ 1} and C α = ∪ {C α (ω) : ω ≥ 0}.
For θ ∈ [0, π) we denote by

Σ θ = {z ∈ C * : | arg(z)| < θ} .
Definition 2.7. A solution operator S α (t) of (2.7) is called analytic if S α (t) admits an analytic extension to a sector Σ θ0 for some

θ 0 ∈ (0, π 2
]. An analytic solution operator is said to be of analytic type (θ 0 , ω 0 ) if for each θ < θ 0 and ω > ω 0 there is M = M (θ, ω) such that

S α (t) L(H) ≤ M e ωRe(t) , ∀ t ∈ Σ θ .
Where Re(λ) stands for the real part of λ. where Re(λ) stands for the real part of λ. The set of all operators A ∈ C α , generating analytic solution operator S α,η (t) of type (θ 0 , ω 0 ) is denoted by A α (θ 0 , ω 0 ). In addition, denote A α (θ 0 ) = ∪ {A α (θ 0 , ω 0 ) :

ω 0 ∈ R + } and A α = ∪ A α (θ 0 ) : θ 0 ∈]0, π 2 ] .
Proposition 2.2. [3, Corollary 2.17] Suppose that {λ : Re(λ) > 0} ⊂ ρ(A) and for some C > 0 we have

(2.12) (λI -A) -1 L(H) ≤ C Re(λ) , ∀ λ ∈ ρ(A), Re(λ) > 0.
Then for any α ∈ (0, 1),

A ∈ A α min 1 α -1 , 1 π 2 , 0 .
Theorem 2.1. Suppose that A is a m-dissipative operator on the Hilbert space H, then A generates a solution operator S α,η (t) and system (1.1) is well-posed. In particular, when η > 0 and α ∈ (0, 1) system (1.1) is exponentially stable and for some M > 0 we have

(2.13) S α,η (t)u 0 H ≤ M e -ηt u 0 H , ∀ u 0 ∈ D(A), t ≥ 0.
Proof. Since the operator A is m-dissipative then by [25, Proposition 3.1.9] property (2.12) holds.

According to Proposition 2.2, A is a generator of solution operator S α (t) of (2.7), therefore Proposition 2.1 leads to the well-posedness of the problem (2.7) and consequently the wellposedness of (1.1) since v(t) = u(t)e ηt . Besides, for some constant M > 0 we have v(t) = S α (t)u 0 satisfies S α (t)u 0 H ≤ M u 0 H , ∀ t ≥ 0. Therefore, the solution operator since u(t) = v(t)e -ηt and (2.13) holds. This completes the proof.

Examples 2.1. As examples we consider the following systems:

(2.14)        ∂ α,η t u(t, x) + (i∆ + a(x))u(t, x) + η Γ(1 -α) t 0 (t -s) -α e -η(t-s) u(s, x) ds = 0, t > 0, x ∈ Ω, u = 0, (0, +∞) × ∂Ω, u(0, x) = u 0 (x), x ∈ Ω, and
(2.15)        ∂ α,η t u(t, x) + (i∆ + i)u(t, x) + η Γ(1 -α) t 0 (t -s) -α e -η(t-s) u(s, x) ds = 0, t > 0, x ∈ Ω, ∂ ν u = ib(x) u, (0, +∞) × ∂Ω, u(0, x) = u 0 (x), x ∈ Ω,
where Ω is a smooth bounded open domain of R n , ∂ ν = ν.∇ is the derivative along ν, the unit normal vector pointing outward of Ω and a ∈ L ∞ (Ω), b ∈ L ∞ (∂Ω) are non-identically zero and non-negative functions.

By a direct implication of Theorem 2.1 we obtain for η > 0 exponential stability results for (2.14) and (2.15) without any geometric conditions (see [START_REF] Bardos | Sharp sufficient conditions for the observation, control, and stabilization of waves from the boundary[END_REF] for example) on the supports of a and b.

Polynomial stabilization

The aim of this section is to establish a polynomial stabilization result of the system (1.1) only for the case η = 0. For this purpose we introduce first some properties of the Mittag-Leffler function ([13, Chapter XVIII] and [20, chapter 1]) E α,β defined by

E α,β (z) = +∞ n=0 z n Γ(αn + β) = 1 2iπ C µ α-β e µ µ α -z dµ, ∀ z ∈ C, α, β > 0,
where C is a contour which starts and ends at -∞ and encircles the disc

D = {µ ∈ C : |µ| ≤ |z| 1 α } counter-clockwise.
For short, we denote E α (z) = E α,1 (z). The first property claims (see [START_REF] Podlubny | Fractional Differential Equations[END_REF]Theorem 1.6]) that for every β > 0 and 0 < α < 2, there exists a constant c > 0, such that

(3.1) |E α,β (-t)| ≤ c 1 + t , ∀ t > 0.
Consider also the function of Wright type φ γ (see [START_REF] Gorenflo | Analytical properties and applications of the Wright function[END_REF][START_REF] Mainardi | Fractional relaxation-oscillation and fractional diffusion-wave phenomena[END_REF][START_REF] Wright | The generalized Bessel function of order greater than one[END_REF]) given by

Φ γ (z) = +∞ n=0 (-z) n n!Γ(1 -γ(n + 1)) = 1 2iπ C µ γ-1 e µ-zµ γ dµ, 0 < γ < 1,
where C is a contour which starts and ends at -∞ and encircles the origin once counterclockwise. The relationship between the Mittag-Leffler function E γ and the function of Wright type Φ γ is given by

(3.2) E γ (z) = +∞ 0 Φ γ (t) e zt dt, ∀ z ∈ C, 0 < γ < 1.
That is, E γ (-z) is the Laplace transform of Φ γ in the whole complex plane. Therefore, Φ γ is a probability density function,

(3.3) Φ γ (t) ≥ 0, ∀ t > 0; and +∞ 0 Φ γ (t) dt = 1.
One of the main ingredients of this section is the following proposition.

Proposition 3.1. [3, Theorem 3.1] Let 0 < α < β ≤ 2, γ = α β and ω ≥ 0. If A ∈ C β (ω) then A ∈ C α (ω 1 γ
) and the following representation holds

(3.4) S α (t) = +∞ 0 ϕ t,γ (s)S β (s) ds, ∀ t > 0,
where ϕ t,γ (s) = t -γ Φ γ (st -γ ). The identity (3.4) holds in the strong sense.

The main result of this section is given by the following theorem.

Theorem 3.1. We suppose that A generates a C 0 -semigroup on the Hilbert space H such that the following properties hold,

(3.5) Re(λ) < 0, ∀ λ ∈ σ(A), and sup
Re(λ)≥0 (λI -A) -1 L(H) < +∞.
Then (1.1) admits a solution operator S α (t) such that there exists c > 0

(3.6) S α (t) L(H) ≤ c 1 + t α , ∀ t ≥ 0.
Proof. With the assumption made on the theorem we can apply Proposition 3.1 with β = 1 and of course γ = α then (1.1) admits a solution operator S α (t) given by the following formula, where we denoted S 1 (t) simply by S(t) which is the C 0 -semigroup generated by A. Thanks to the assumptions (3.5), then according to [START_REF] Huang | Characteristic conditions for exponential stability of linear dynamical systems in Hilbert space[END_REF][START_REF] Prüss | On the spectrum of C 0 -semigroups[END_REF] the uniform stabilization holds, that is there exist ω 0 > 0 and K > 0 such that

(3.8) S(t) L(H) ≤ Ke -ω0t ∀ t > 0.
Performing a change of variable on (3.7) and using (3.3) and (3.8) we obtain

(3.9) S α (t) L(H) = +∞ 0 Φ α (s)S(st α ) ds L(H) ≤ +∞ 0 Φ α (s)e -ω0st α ds ∀ t > 0.
Following to (3.2) and (3.9) we find

S α (t) L(H) ≤ E α (-ω 0 t α ) ∀ t > 0.
Estimate (3.6) follows now from (3.1) and this completes the proof.

Examples 3.1. As examples we consider here the same systems as above but with η = 0:

(3.10)

   ∂ α,η t u(t, x) + (i∆ + a(x))u(t, x) = 0, t > 0, x ∈ Ω, u = 0, (0, +∞) × ∂Ω, u(0, x) = u 0 (x), x ∈ Ω, and (3.11)
   ∂ α,η t u(t, x) + (i∆ + i)u(t, x) = 0, t > 0, x ∈ Ω, ∂ ν u = ib(x) u, (0, +∞) × ∂Ω, u(0, x) = u 0 (x), x ∈ Ω,
where Ω is a smooth bounded open domain of R n , ∂ ν = ν.∇ is the derivative along ν, the unit normal vector pointing outward of Ω and a ∈ L ∞ (Ω), b ∈ L ∞ (∂Ω) are non-identically zero and non-negative functions.

By a direct implication of Theorem 3.1 we obtain polynomial stability results for (3.10) and (3.11) under geometric conditions G.C.C. (see, respectively, [START_REF] Lebeau | Équation des ondes amorties, Algebraic and geometric methods in mathematical physics[END_REF] and [START_REF] Bardos | Sharp sufficient conditions for the observation, control, and stabilization of waves from the boundary[END_REF], for example) on the supports of a and b.

Extension to some integro-differential equation

Let X be a Hilbert space equipped with the norm . X , and let A : D(A) ⊂ X → X be a closed, self-adjoint and strictly positive operator on X with dense domain. We introduce the scale of Hilbert spaces X β , β ∈ R, as follows: for every β ≥ 0, X β = D(A β ), with the norm z β = A β z X . The space X -β is defined by duality with respect to the pivot space H as follows: X -β = X * β for β > 0. The operator A can be extended (or restricted) to each X β , such that it becomes a bounded operator

A : X β → X β-1 , ∀ β ∈ R.
Let a bounded linear operator B : U → X - 1 2 , where U is another Hilbert space which will be identified with its dual.

We consider the following integro-differential equation (4.1)

D α t u(t) + g α * Au(t) + BB * u(t) = 0, t > 0, u(0) = u 0 .
We set the Hilbert space H = X 1 2 ×X and we consider the unbounded operator A : D(A) -→ H defined by

A = 0 I -A -BB * , where D(A) = {(v, u) ∈ H : u ∈ X 1 2 , Av + BB * u ∈ X}. It is well known (see [1, 2]) that A is a generator of a C 0 -semigroup of contractions on H. Definition 4.1. A function u ∈ C(R + , X) such that g α * u ∈ C(R + , X 1
2 ) is called a strong solution of (4.1) if the couple

g α * u u ∈ C(R + , D(A)), g 1 * u g 1-α * (u -u 0 ) ∈ C 1 (R + , H)
and (4.1) holds on R + with u 0 ∈ X.

Definition 4.2. The problem (4.1) is called well-posed if for any u 0 ∈ X 1 2 such that BB * u 0 ∈ X 1 2 , there is a unique strong solution u(t, u 0 ) of (4.1), and when u 0 n ∈ D(A) such that u 0 n -→ 0 as n +∞ then u(t, u 0 n ) -→ 0 in X and g α * u(., u 0 n )(t) -→ 0 in X 1 2 as n +∞ in H, uniformly on compact intervals. Theorem 4.1. Under the above assumptions made on the operator A, system (4.1) is wellposed in such away if u 0 ∈ X 1 2 such that BB * u 0 ∈ X and we have the following regularity of the solution

u ∈ C(R + , X 1
2 ), g α * u ∈ C(R + , X 1 2 ), g 1-α * (u -u 0 ) ∈ C 1 (R + , X) If in addition, the following properties hold, Then for some constant C > 0 and for any data u 0 ∈ X 1 2 , the solution u(t) of (4.1) satisfies the following asymptotic estimates (4.3) u(t) X ≤ C 1 + t α u 0 X , ∀ t ≥ 0, and (4.4)

g α * u(t) X 1 2 ≤ C 1 + t α u 0 X , ∀ t ≥ 0.
Remark 4.1. We note that in the case where B ∈ L(U, X) and if u 0 ∈ X 1 2 we have immediately that BB * u 0 ∈ X.

Proof. Let's consider the following equation (4.5) D α t U (t) = AU (t), t ≥ 0, U (0) = U 0 , where we have denoted by

U (t) = v(t) u(t) and U 0 = v 0 u 0 .
Since the operator A is m-dissipative (see [START_REF] Ammari | Fractional-feedback stabilization for a class of evolution systems[END_REF]) then according to Theorem 2.1, system (4.5) is well-posed as given by Definition 2.5. In the other hand, according to section 2 system (4.5) is equivalent to the following integral equation (4.6) U (t) = U 0 + g α * AU (t).

Equation (4.6) can be also writing as follow:

v(t) u(t) = v 0 u 0 + g α * 0 I -A -BB * v(t) u(t)
.

  is called a solution operator (or a resolvent) for (1.1) or for (2.8) if the following conditions are satisfied (a) S α,η (t) is strongly continuous for t ≥ 0 and S α,η (0) = I. (b) S α,η (t)(D(A)) ⊂ D(A) and AS α,η (t)x = S α,η (t)Ax for all x ∈ D(A) and t ≥ 0. (c) The resolvent equation holds

(3. 7 )

 7 S α (t) = +∞ 0 ϕ t,α (s)S(s) ds, ∀ t > 0,

( 4 . 2 )

 42 iR ⊂ ρ(A), and lim sup µ∈R,|µ|→+∞(iµI -A) -1 L(H) < +∞.

  and (1.1) holds on R + . Equation (2.8) is well-posed if and only if it admits a resolvent S α,η (t). If this is the case we have in addition (g α * (e η. S α,η (.)x))(t) ∈ D(A) for all x ∈ H and t ≥ 0 and we have

	Definition 2.5. The problem (1.1) is called well-posed if for any u 0 ∈ D(A), there is a unique
	strong solution u(t, u 0 ) of (1.1), and when u 0 n ∈ D(A) such that u 0 n -→ 0 as n u(t, u 0 n ) -→ 0 as n +∞ in H, uniformly on compact intervals.	+∞ then
	Proposition 2.1. 1/ (2.9) S

α,η (t)x = e -ηt x + e -ηt A(g α * (e η. S α,η (.

)x))(t), ∀ x ∈ H, t ≥ 0. 2/ System (1.1) is well-posed if and only if (2.8) is well-posed.

Proof.

1/ Problem (2.8) can be written as follows:

Equivalently, we have (4.7) u(t) = u 0 -g 2α * Au(t) -g α * Av 0 -g α * BB * u(t),

where we have used the semigroup property (2.1). By taking v 0 = 0 and u 0 ∈ X 1 2 such that BB * u 0 ∈ X in (4.7), we obtain

Since in this case U (t) is given by the couple g α * u(t) u(t) where u(t) is the solution of the system (4.8) and the problem (4.5) is well-posed then by Definition 2.5, U ∈ C(R, D(A)) and

2 ) and g 1-α * (uu 0 ) ∈ C 1 (R + , X) then by applying the operator D α t on both sides of(4.8) we find that (4.8) is equivalent to

where we have used here the semigroup property (2.1) and (2.5). Now we have proved that system (4.5) with U 0 = 0 u 0 with u 0 ∈ X 1 2 such that BB * u 0 ∈ X is equivalent to the equation (4.8). Since (4.5) is well-posed in the sens of Definition 2.5 then problem (4.1) is also well-posed in the sens of Definition 4.2 and the regularities of the solution u(t) of (4.1) given by the theorem hold. Now if assumptions (4.2) hold, then according to Theorem 3.1 the solution U of (4.5) satisfies the following estimation

for some constant C > 0, then we follow

This implies in particular the estimates (4.3) and (4.4) and completes the proof.

Example 4.1. We consider the following integro-differential equation:

where Ω is a smooth bounded open domain of R n and a ∈ L ∞ (Ω) is non-identically zero and positive function.

By a direct implication of Theorem 4.1 we obtain a polynomial stability result for (4.9) under a geometric condition G.C.C (see [START_REF] Lebeau | Équation des ondes amorties, Algebraic and geometric methods in mathematical physics[END_REF] for more details) on the support of a.