Romain Joly
email: romain.joly@univ-grenoble-alpes.fr

Camille Laurent
email: laurent@ann.jussieu.fr

Decay of semilinear damped wave equations: cases without geometric control condition

Keywords: damped wave equations, stabilization, semi-uniform decay, unique continuation property, small trapped sets, weak attractors

We consider the semilinear damped wave equation

In this article, we obtain the first results concerning the stabilization of this semilinear equation in cases where γ does not satisfy the geometric control condition. When some of the geodesic rays are trapped, the stabilization of the linear semigroup is semi-uniform in the sense that ∥e At A -1 ∥ ≤ h(t) for some function h with h(t) → 0 when t → +∞. We provide general tools to deal with the semilinear stabilization problem in the case where h(t) has a sufficiently fast decay.

Introduction

We consider the semilinear damped wave equation

   ∂ 2 tt u(x, t) + γ(x)∂ t u(x, t) = ∆u(x, t) -αu(x, t) -f (x, u(x, t)) (x, t) ∈ Ω × (0, +∞) u |∂Ω (x, t) = 0 (x, t) ∈ ∂Ω × (0, +∞) (u(•, t = 0), ∂ t u(•, t = 0)) = U 0 = (u 0 , u 1 ) ∈ H 1 0 (Ω) × L 2 (Ω) (1.
1) in the following general framework:

(i) the domain Ω is a two-dimensional smooth compact and connected manifold with or without smooth boundary. If Ω is not flat, ∆ has to be taken as Beltrami Laplacian operator.

(ii) the constant α ≥ 0 is a non-negative constant. We require that α > 0 in the case without boundary to ensure that ∆ -αId is a negative definite self-adjoint operator.

(iii) the damping γ ∈ L ∞ (Ω, R + ) is a bounded function with non-negative values. Since we want to consider a damped equation, we will assume that γ does not vanish everywhere.

(iv) the non-linearity f ∈ C 1 (Ω × R, R) is of polynomial type in the sense that there exists a constant C and a power p ≥ 1 such that for all (x, u)

∈ Ω × R, |f (x, u)| + |∇ x f (x, u)| ≤ C(1 + |u|) p and |f ′ u (x, u)| ≤ C(1 + |u|) p-1 . (1.2)
Moreover, in most of this paper, we will be interested in the stabilization problem and we will also assume that

∀(x, u) ∈ Ω × R , f (x, u)u ≥ 0 . (1.3)
We introduce the space X = H 1 0 (Ω) × L 2 (Ω) and the operator A defined by

D(A) = (H 2 (Ω) ∩ H 1 0 (Ω)) × H 1 0 (Ω) A = ( 0 Id ∆ -αId -γ(x)
)

.

In this paper, we are interested in the cases where the linear semigroup e At has no uniform decay, that is that ∥e At ∥ L(X) does no converge to zero. We only assume a semi-uniform decay, but sufficiently fast in the following sense.

(v) There exist a function h(t) such that

∀U 0 ∈ D(A) , ∥e At U 0 ∥ X ≤ h(t) ∥U 0 ∥ D(A) (1.4)
and there is σ h ∈ (0, 1] such that lim t→∞ h(t) = 0 and ∀σ ∈ [0, σ h ] ,

∫ ∞ 0 h(t) 1-σ dt < ∞ . (1.5)
Condition (1.5) requires a decay rate fast enough to be integrable. Roughly speaking, this article shows that this condition, together with a suitable unique continuation property, are sufficient to obtain a stabilization of the semilinear equation. The relevant unique continuation property is explained in Proposition 3.5 below. We present two general results where it can be obtained.

Our first result concerns analytic nonlinearities and smooth dampings. b) the damping γ is of class C 1 or at least that there exists γ ∈ C 1 (Ω, R + ) such that (v) holds with γ replaced by γ and such that the support of γ is contained in the support of γ.

c) the power p of f in (1.2) and the decay rate h(t) of the semigroup in (1.4) satisfy h(t) = O(t -β ) with β > 2p.

Then, any solution u of (1.1) satisfies

∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 --------→ t-→+∞ 0 .
Moreover, for any R and σ > 0, there exists h R,σ (t) which goes to zero when t goes to +∞ such that the following stabilization hold. For any U 0 ∈ H 1+σ 0 (Ω) × H σ (Ω), if u is the solution of (1.1), then

∥(u 0 , u 1 )∥ H 1+σ ×H σ ≤ R =⇒ ∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 ≤ h R,σ (t) --------→ t-→+∞ 0 .
Our assumptions (v) and c) on the decay of the linear semigroup may seem strong. They are satisfied in the cases where the set of trapped geodesics, the ones which do not meet the support of the damping, is small and hyperbolic in some sense. Several geometries satisfying (v) and c) have been studied in the literature, see the concrete examples of Figure 1 and the references therein. Notice in particular that the example of domain with holes is particularly relevant for applications where we want to stabilize a nonlinear material with holes by adding a damping or a control in the external part. There is a huge literature about the damped wave equation and the purpose of the examples presented here is mainly to illustrate our theorem to non specialists. Moreover, the subject is growing fastly, giving more and more examples of geometries where we understand the effect of the damping and where we may be able to apply our results. We do not pretend to exhaustivity and refer to the bibliography of the more recent [START_REF]Long Jin Damped wave equations on compact hyperbolic surfaces[END_REF] for instance.

In some cases, the unique continuation property required in Proposition 3.5 can be obtained without considering analytic nonlinearities or conditions on the growth of f as Hypothesis c) of Theorem 1.1. Instead, we require a particular geometry, which will be introduced more precisely in Section 6.

Theorem 1.2. Consider the damped wave equation (1.1) in the framework of Assumptions (i)- (v). Assume in addition that:

a) the function (x, u) → f (x, u) is of class C 1 (Ω × R, R).
b) there exists a pseudo-convex foliation of Ω in the sense of Definitions 6.4 or 6.6.

Then, the conclusions of Theorem 1.1 hold. This result can be applied in several situations of Figure 1: the "disk with two holes", the "peanut of rotation" and the "open book". In these cases, the stabilization holds for any natural nonlinearity.

We expect that the decay rate h R,σ (t) is related to the linear decay rate h(t) of Assumption (v). We are able to obtain this link for the typical decays of the examples of Figure 1. • if the decay rate of Assumption (v) satisfies h(t) = O(t -α ) with α > 1, then the nonlinear equation admits a decay of the type h R,σ (t) = O(t -σα ).

• if the decay rate of Assumption (v) satisfies h(t) = O(e -at 1/β ) with a > 0 and β > 0, then the nonlinear equation admits a decay of the type h R,σ (t) = O(e -bσt 1/(β+1) ) for some b > 0.

Notice that this result is purely local in the sense that the decay rate is obtained when the solution is close enough to 0. Our proofs do not provide an explicit estimate of the time needed to enter this small neighborhood of 0. Also notice that the loss in the power of the second case of Proposition 1.3 is due to an abstract setting: in the concrete examples, we may avoid this loss, see the remark below Lemma 4.2 and the concrete applications to the examples of Figure 1.

To our knowledge, Theorems 1.1 and 1.2 are the first stabilization results for the semilinear damped wave equation when the geometric control condition fails. This famous geometric control condition has been introduced in the works of Bardos, Lebeau, Rauch and Taylor (see [START_REF] Bardos | Sharp sufficient conditions for the observation, control and stabilization of waves from the boundary[END_REF]) and roughly requires that any geodesic of the manifold Ω meets the support of the damping γ. This condition implies that the linear semigroup of the damped wave equation satisfies a uniform decay ∥e At ∥ L(H 1 0 ×L 2 ) ≤ M e -λt . In this context, the stabilization of the semilinear damped wave equation has been studied since a long time, see for example [START_REF] Haraux | Stabilization of trajectories for some weakly damped hyperbolic equations[END_REF][START_REF] Zuazua | Exponential decay for semilinear wave equations with localized damping[END_REF][START_REF] Dehman | Stabilisation pour l'équation des ondes semi-linéaire[END_REF][START_REF] Belhassen Dehman | Stabilization and control for the subcritical semilinear wave equation[END_REF][START_REF] Joly | Stabilization for the semilinear wave equation with geometric control condition[END_REF]. Under this condition and for large data, the proof often divides into a part dealing with high frequencies with linear arguments and another

Disk with holes

We set Ω to be a convex flat surface with a damping γ efficient near the boundary. Typically, we may take the flat disk B(0, 1) of R 2 and assume that there exist r ∈ (0, 1) and γ > 0 such that γ(x) ≥ γ > 0 for |x| > r. Inside the interior zone without damping, assume that at least two holes exist; to simplify, we also assume that these holes are disks and are small in a sense specified later. Notice that there exist some periodic geodesics which do not meet the support of the damping. This example has been studied in [START_REF] Burq | Contrôle de l'équation des plaques en présence d'obstacles strictement convexes[END_REF][START_REF] Burq | Geometric control in the presence of a black box[END_REF].

' & $ %

The peanut of rotation

We consider a compact twodimensional manifold without boundary. We assume that the damping γ is effective, that is uniformly positive, everywhere except in the central part of the manifold. This part is a manifold of negative curvature and invariant by rotation along the y-axis. More precisely, let us set this part to be equivalent to the cylinder endowed with the metric g(y, θ) = dy 2 + cosh 2 (y)dθ 2 . This central part admits a unique (up to change of orientation) periodic geodesic which is unstable; any other geodesic meets the support of the damping. This example has been studied for example in [START_REF] Christianson | From resolvent estimates to damped waves[END_REF], [START_REF] Schenck | Exponential stabilization without geometric control[END_REF] and the references therein.

' & $ % γ(x 1 , x 2 ) = |x 1 | β

The open book

We consider the torus T 2 with flat geometry. The damping γ is assumed to depend only on the first coordinate and to be of the type γ(x) = |x 1 | β with β > 0 to be chosen small enough. In this case, there is a unique (up to change of orientation) geodesic which does not meet the support of the damping. This example has been studied in [START_REF] Léautaud | Sharp polynomial energy decay for locally undamped waves, Séminaire Laurent Schwartz -EDP et applications[END_REF].

' & $ %

Hyperbolic surfaces

We consider a compact connected hyperbolic surface with constant negative curvature -1 (for example a surface of genus 2 cut out from Poincaré disk). The damping is any non zero function γ(x) ≥ 0. This example has been studied in [START_REF]Long Jin Damped wave equations on compact hyperbolic surfaces[END_REF] following the fractal uncertainty principle of [START_REF] Bourgain | Dyatlov Spectral gaps without the pressure condition[END_REF]. We also refer to other results in any dimension with pressure conditions [START_REF] Schenck | Energy Decay for the Damped Wave Equation Under a Pressure Condition[END_REF] following ideas of [START_REF] Anantharaman | Entropy and the localization of eigenfunctions[END_REF].

Figure 1: the main applications of Theorems 1.1 and 1.2 presented in this paper. The gray parts show the localization of the damping (white=no damping). The more geometrically constrained Theorem 1.2 apply to the "disk with two holes", the "peanut of rotation" and the "open book". one dealing with low frequencies that often requires a unique continuation argument. The high frequency problem was solved by Dehman [START_REF] Dehman | Stabilisation pour l'équation des ondes semi-linéaire[END_REF] with important extension by Dehman-Lebeau-Zuazua [START_REF] Belhassen Dehman | Stabilization and control for the subcritical semilinear wave equation[END_REF] using microlocal defect measure. Yet, the unique continuation was proved by classical Carleman estimates (see Section 6.2 below) which restricted the generality of the geometry. Using techniques from dynamical systems applied to PDEs, the authors of the present article proved in [START_REF] Joly | Stabilization for the semilinear wave equation with geometric control condition[END_REF] a general stabilization result under Geometric Control Condition, at the cost of an assumptions of analyticity of the nonlinearity. Theorem 1.1 is in the same spirit as [START_REF] Joly | Stabilization for the semilinear wave equation with geometric control condition[END_REF] and intends to prove that related techniques can be extended to a weaker damping. Theorem 1.2 is more in the spirit of the other references, taking advantage of particular geometries, but avoiding analyticity.

Notice that the "disk with one hole" satisfies this geometric control condition and thus it is not considered in this paper.

In the cases where the geometric control condition fails, the decay of the linear semigroup is not uniform. At least, if γ does not vanish everywhere, it is proved in [START_REF] Dafermos | Asymptotic behavior of solutions of evolution equations[END_REF] (see also [START_REF] Haraux | Stabilization of trajectories for some weakly damped hyperbolic equations[END_REF]) that the trajectories of the linear semigroup goes to zero (see Theorem 2.1 below). In fact, the decay can be estimated with a loss of derivative as

∥e At U ∥ H 1 ×L 2 ≤ h(t)∥U ∥ H 2 ×H 1 with h(t) -----→ t→+∞ 0 . (1.6)
In the general case, as soon as γ ̸ ≡ 0, the decay rate can be taken as h(t) = O(ln(ln(t))/ ln t) as shown in [START_REF] Lebeau | Equation des ondes amorties[END_REF][START_REF] Lebeau | Stabilisation de l'équation des ondes par le bord[END_REF]. In some particular situations, γ misses the geometric control condition but very closely: typically there is only one (up to symmetries) geodesic which does not meet the support of the damping and this geodesic is unstable. In this case, we may hope a better decay than the O(ln(ln(t))/ ln t) one, see for example [START_REF] Christianson | From resolvent estimates to damped waves[END_REF][START_REF] Léautaud | Sharp polynomial energy decay for locally undamped waves, Séminaire Laurent Schwartz -EDP et applications[END_REF] and the other references of Figure 1.

To our knowledge, until now, there was no result concerning the semilinear damped wave equation (1.1) when the geometric control condition fails. Thus Theorems 1.1 and 1.2 provides the first examples of semi-uniform stabilization for the semilinear damped wave equation. Notice that our results deeply rely on the fact that the decay rate of (1.6) is integrable. Typically, for the situations of Figure 1, it is of the type h(t) = O(e -λt α ) or h(t) = O(1/t β ) with sufficiently large β > 0.

Theorems 1.1 and 1.2 concern the stabilization of the solutions of (1.1) in the sense that their H 1 × L 2 -norm goes to zero. Notice that, since the energy of the damped wave equation is non-increasing (see Section 2), we knew that this H 1 × L 2 -norm is at least bounded. Such a uniform bound is not clear a priori for the H 2 × H 1 -norm. However, basic arguments provide this bound as a corollary of Theorems 1.1 and 1.2 if the decay is fast enough, which is the case of the "disk with holes", the "peanut of rotation" and the "hyperbolic surfaces" of Figure 1.

Theorem 1.4. Consider the damped wave equation (1.1) in the framework of Theorems

1.1 or 1.2. Assume that for all R > 0, the decay rate h R,1 (t) is faster than polynomial, i.e. h R,1 (t) = o(t -k ) for any k ∈ N. Also assume that γ is of class C 1 and f is of class C 2 (Ω × R, R).
Then the H 2 (Ω) × H 1 (Ω)-norm of the solutions are bounded in the following sense. For any R > 0, there exists C(R) > 0 such that, for any

U 0 ∈ (H 2 (Ω) ∩ H 1 0 (Ω)) × H 1 0 (Ω) such that ∥U 0 ∥ H 2 ×H 1 ≤ R, the solution u of (1.1) satisfies sup t≥0 ∥(u, ∂ t u)(t)∥ H 2 ×H 1 ≤ C(R) .
Note that, in the case without damping, this result is sometimes expected to be false. It is related to the weak turbulence, described as a transport from low frequencies to high frequencies.

The main purpose of this paper is to obtain new examples of stabilization for the semilinear damped wave equation and to introduce the corresponding methods and tools. We do not pretend to be exhaustive and the method may be easily used to obtain further or more precise results. For example:

• the boundary condition may be modified, typically in the case of the disk with holes, Neumann boundary condition may be chosen at the exterior boundary.

• for simplicity, the examples of Figure 1 and the main results of this article concern two-dimensional manifolds. However, the arguments of this paper can be used to deal with higher-dimensional manifolds. There are some technical complications, mainly due to the Sobolev embeddings. For example, in dimension d = 3, the degree of f in (1.2) should satisfy p < 3 (p < 5 if we use Strichartz estimates as done in [START_REF] Joly | Stabilization for the semilinear wave equation with geometric control condition[END_REF] using [START_REF] Belhassen Dehman | Stabilization and control for the subcritical semilinear wave equation[END_REF]) and the order β of the vanishing of γ in the example of the open book should not be too large. To simplify, we choose to state our results in dimension d = 2. However, several intermediate results in this article are stated for dimensions d = 2 or d = 3.

• It is also possible to combine the strategy of this paper with other tricks and technical arguments. For example, we may consider unbounded manifolds or manifolds of dimension d = 3 with nonlinearity of degree p ∈ [START_REF] Anatolii | Regular attractors of semigroups and evolution equations[END_REF][START_REF] Charles | Non-uniform stability for bounded semigroups on Banach spaces[END_REF], which are supercritical in the Sobolev sense. This would requires to use Strichartz estimates in addition to Sobolev embeddings as done in [START_REF] Belhassen Dehman | Stabilization and control for the subcritical semilinear wave equation[END_REF] or [START_REF] Joly | Stabilization for the semilinear wave equation with geometric control condition[END_REF].

• assume that we replace the sign condition (1.3) by an asymptotic sign condition

∃R > 0 , ∀(x, u) ∈ Ω × R , |u| ≥ R ⇒ f (x, u)u ≥ 0 .
Then they may exist several equilibrium points and the stabilization to zero cannot be expected. However, the arguments of this paper show that the energy E introduced in Section 2 is a strict Lyapounov functional and that any solution converges to the set of equilibrium points. We can also show the existence of a weak compact attractor in the sense that there is an invariant compact set A ⊂ H 1 0 (Ω)×L 2 (Ω), which consists of all the bounded trajectories and such that any regular set B bounded in (H 2 (Ω)∩H 1 0 (Ω))×H 1 0 (Ω) is attracted by A in the topology of H 1 0 (Ω)×L 2 (Ω). Notice that this concept of weak attractor is the one of Babin and Vishik in [START_REF] Anatolii | Regular attractors of semigroups and evolution equations[END_REF]. At this time, the asymptotic compactness property of the semilinear damped wave equation was not discovered and people thought that a strong attractor (attracting bounded sets of H 1 0 (Ω) × L 2 (Ω)) was impossible due to the lack of regularization property for the damped wave equation. Few years later, Hale [START_REF] Hale | Asymptotic behavior of dissipative systems[END_REF] and Haraux [START_REF] Haraux | Stabilization of trajectories for some weakly damped hyperbolic equations[END_REF] obtain this asymptotic compactness property and the existence of a strong attractor. Thus this notion of weak attractor has been forgotten. It is noteworthy that it appears again here. Notice that we cannot hope a better attraction property since even in the linear case, {0} is not an attractor in the strong sense.

The organization of this paper follows the proof of stabilization of the examples of Figure 1. We add step by step the techniques required to deal with our guiding examples, from the simplest to the most complicated one.

Sections 2 and 3 contain the basic notations and properties. The asymptotic compactness of the semilinear dynamics is proved and the problem is reduced to a unique continuation property. In Section 4, we show the estimations of Proposition 1.3. Section 5 then proves the nonlinear stabilization in the "open book" case, where the unique continuation property is trivial. Section 6 stated several unique continuation results, enabling to prove Theorem 1.2. We obtain as a consequence the stabilization in the case of the "peanut of rotation" in Section 7. Section 8 studies the linear semigroup for the case of the "disk with holes" before we apply Theorem 1.2 in the case of the "disk with two holes" in Section 9. Theorem 1.1 is proved in Section 10 by showing an asymptotic analytic regularization. It is applied to the "disk with three or more holes" and hyperbolic surfaces, assuming f analytic in u, in Sections 11 and 12. In Section 13, we show how to obtain Theorem 1.4 as a corollary of Theorems 1.1 or 1.2. This article finishes with three appendices on the links between the decay of the semigroup e At and the resolvent (A -iµId) -1 .

Acknowledgements:

The authors deeply thank Matthieu Léautaud for his contributions to the appendices. They are also grateful to Nicolas Burq for several discussions and the suggestion of Theorem 1.4. Part of this work has been made in the fruitful atmosphere of the Science Center of Benasque Pedro Pascual and has been supported by the project ISDEEC ANR-16-CE40-0013.

Notations and basic facts

We use the notations of Equation (1.1), of Assumptions (i)-(v) and of the introduction. In particular, we recall that X = H 1 0 (Ω) × L 2 (Ω) and

D(A) = (H 2 (Ω) ∩ H 1 0 (Ω)) × H 1 0 (Ω) A = ( 0 Id ∆ -αId -γ(x)
)

.

The operator A is the classical linear damped wave operator corresponding to the linear part of (1.1). Due to Lumer-Phillips theorem, we know that this operator generates a linear semigroup of contractions e At on X and on D(A) and that ∀t ≥ 0 , ∥e At ∥ L(X) ≤ 1 and ∥e At ∥ L(D(A)) ≤ 1 .

Notice that the second estimate is a direct consequence of the commutation of A and e At and does not require any regularity on γ.

For any σ ∈ [0, 1], we set

X σ = (H 1+σ (Ω) ∩ H 1 0 (Ω)) × H σ 0 (Ω) .
Thus X 0 = X and X 1 = D(A) and X σ is an interpolation space between X 0 and X 1 . In particular, by interpolation, e At is defined in X σ and we have

∀σ ∈ [0, 1] , ∀t ≥ 0 , ∥e At ∥ L(X σ ) ≤ 1 . (2.1)
We set F ∈ C 0 (X) to be the function

F : U = ( u v ) ∈ X -→ ( 0 -f (•, u) ) ∈ X . (2.2) Notice that, if Ω is two-dimensional, H 1 0 (Ω) → L p (Ω) for any p ∈ [1, +∞) and if Ω is three- dimensional H 1 0 (Ω) → L 6 (Ω). Thus, f (u) is well defined in L 2 (Ω) due to Assumption (1.2)
if dim(Ω) = 2 or if dim(Ω) = 3 and p ≤ 3. Moreover, for any R, u and v with ∥u∥

H 1 ≤ R and ∥v∥ H 1 ≤ R, we have ∥f (•, u) -f (•, v)∥ L 2 = (u -v) ∫ 1 0 f ′ u (•, u + s(u -v))ds L 2 ≤ C(R)∥u -v∥ H 1
and so F is lipschitzian on the bounded sets of X. As a consequence, the damped wave equation (1.1) is well posed in X and admits local solutions if dim(Ω) = 2 or if dim(Ω) = 3 and p ≤ 3.

With the above notation, our main equation writes

∂ t U = AU + F (U ) U (t = 0) = U 0 ∈ X . (2.3)
In particular, Duhamel's formula yields

U (t) = e At U 0 + ∫ t 0 e A(t-s) F (U (s)) ds .
We introduce the potential

V (x, u) = ∫ u 0 f (x, s)ds .
Due to (1.2) and the above arguments, V (•, u) defines a Lipschitz function from the bounded sets of H 1 0 (Ω) into L 1 (Ω). The classical energy associated to (1.1) is defined along a trajectory U = (u, ∂ t u) as

E(U ) = ∫ Ω 1 2 (|∇u| 2 + α|u| 2 + |∂ t u| 2 ) + V (x, u) .
The damping effect appears by the computation

∂ t E(U (t)) = - ∫ Ω γ(x)|∂ t u| 2 .
(2.4)

In particular, the energy E is non-increasing along the trajectories. Moreover, the sign assumption (1.3) yields that V (x, u) ≥ 0. Thus, we have that E(U ) ≥ C∥U ∥ 2 X and that E(t), t ≥ 0, is bounded on the bounded sets of X. All together, the above properties show that for any U 0 ∈ X, the solution U = (u, ∂ t u) of (1.1) is defined for all non-negative times and remains in a bounded set of X, which only depends on ∥U 0 ∥ X .

A fundamental question of this paper concerns the solution for which the energy is constant: are they equilibrium points or may they be moving trajectories? At least, the answer is known for the linear equation, see [START_REF] Dafermos | Asymptotic behavior of solutions of evolution equations[END_REF] and also [START_REF] Haraux | Stabilization of trajectories for some weakly damped hyperbolic equations[END_REF].

Theorem 2.1. Dafermos (1978).

Assume that the damping γ ≥ 0 does not vanish everywhere. Then, for any U 0 ∈ X, we have e At U 0 --------→ t-→+∞ 0 in X .

Asymptotic compactness and reduction to a unique continuation problem

In this section, we assume a fast enough semi-uniform linear decay as described by (1.4) and (1.5). We first notice that, by linear interpolation, we have the following result.

Proposition 3.1. For any σ 1 , σ 2 such that 0 ≤ σ 2 < σ 1 ≤ 1, the linear semigroup is well defined from X σ 1 in X σ 2 and we have

∀U 0 ∈ X σ 1 , ∥e At U 0 ∥ X σ 2 ≤ h(t) σ 1 -σ 2 ∥U 0 ∥ X σ 1 .
Proof: We interpolate the estimates (2.1) for σ = 1 and (1.4) with respective weights (σ 2 /σ 1 , 1 -σ 2 /σ 1 ) and obtain

∀U 0 ∈ D(A) , ∥e At U 0 ∥ X σ 2 /σ 1 ≤ h(t) 1-σ 2 /σ 1 ∥U 0 ∥ D(A) .
It remains to interpolated the above estimate and (2.1) for σ = 0 with respective weights (σ 1 , 1 -σ 1 ). □

We also need some regularity properties for F . The following properties depend on Sobolev embeddings and so of the dimension d of Ω. For d = 2, which is the case in our examples, the properties are general. For d = 3, they are more restrictive but they are shown in the same way. We choose to also consider this case in our paper for possible later uses. Proposition 3.2. Assume that dim(Ω) = 2. Then for any σ ∈ [0, 1), the function F maps any bounded set

B of X = H 1 0 (Ω) × L 2 (Ω) in a bounded set F (B) contained in X σ = (H 1+σ (Ω) ∩ H 1 0 (Ω)) × H σ 0 (Ω). Moreover, F (B) has compact enclosure in X σ . If dim(Ω) = 3 and if (1.
2) holds for some p ∈ [0, 3), then the same properties hold for σ ∈ [0, (3 -p)/2).

Proof: Assume that dim(Ω) = 2. First notice that we only have to show that f (•, u) is compactly bounded in H σ 0 (Ω) since the first component of F is zero. Also notice that f (x, 0) = 0 due to the sign assumption (1.3), thus the Dirichlet boundary condition possibly contained in H σ 0 (Ω) will be fulfilled by f (x, u) if u ∈ H 1 0 (Ω). Due to the Sobolev embeddings, and since Ω is compact, it is sufficient to show that F (B) is bounded in W 1,q (Ω) for all q ∈ [1, 2) to obtain compactness in H σ (Ω) for any σ ∈ [0, 1). Since f is of polynomial type due to (1.2), we know that f (x, u) is bounded in L q (Ω) for any q ∈ [1, 2). On the other hand, using (1.2), we have

∥∇(f (x, u))∥ L q ≤ ∥∇ x f (x, u)∥ L q + ∥f ′ u (x, u)∇u∥ L q ≤ C∥(1 + |u|) p ∥ L q + C∥(1 + |u|) p-1 ∇u∥ L q ≤ C(1 + ∥u∥ p L pq + ∥∇u∥ L 2 ∥u∥ p-1 L r )
with r = (p -1) 2q 2-q defined as soon as q < 2. This shows that f (•, u) belongs to W 1,q (Ω) for any q ∈ [1, 2) and concludes the proof for dim(Ω) = 2.

The case dim(Ω) = 3 is similar once we use the suitable Sobolev embeddings. □

The main results of this section are the following asymptotic compactness properties.

Proposition 3.3. If dim(Ω) = 2, set σ * = σ h . If dim(Ω) = 3, assume that p ∈ [0, 3) in (1.2) and that σ h ∈ ((p -1)/2, 1) in (1.5) and set σ * = σ h -(p -1)/2.
Let U (t) = (u, ∂ t u) where u solves (1.1) and let (t n ) be a sequence of times such that t n → +∞. Then, there exist a subsequence (t φ(n) ) and a solution W (t) = (w, ∂ t w)(t) of (1.1) defined for all t ∈ R, such that

∀t ∈ R , U (t φ(n) + t) --------→ n-→+∞ W (t) in X = H 1 0 (Ω) × L 2 (Ω) .
Moreover, the solution W is globally bounded in X σ for all σ ∈ [0, σ * ) and the energy E(W (t)) is constant.

Proof: Assume first that dim(Ω) = 2. We have

U (t n ) = e Atn U 0 + ∫ tn 0 e A(tn-s) F (U (s)) ds . (3.1)
Due to Theorem 2.1, the term e Atn U 0 goes to zero in X. Thus, it remains to show that ∫ tn 0 e A(tn-s) F (U (s)) ds is a compact term in X. First notice that U (s) is uniformly bounded for s ≥ 0 due to the non-increasing energy E (see Section 2). Due to Proposition 3.2, F (U (s)) thus belongs to a bounded set of X σ 1 for all σ 1 ∈ [0, 1). By Assumption (1.5) and Proposition 3.1, e A(tn-s) F (U (s)) has an integral in [0, t n ] bounded in X σ 2 uniformly with respect to n, for any σ 2 ∈ (0, σ h ). Thus ∫ tn 0 e A(tn-s) F (U (s)) ds is a compact sequence in X σ for any σ ∈ [0, σ 2 ). As a consequence, for any σ as close as wanted to σ h , we may extract a subsequence (t φ(n) ) such that ∫ t φ(n) 0 e A(t φ(n) -s) F (U (s)) ds converges to some limit W (0) in X σ . Since the linear term of (3.1) goes to zero in X for t φ(n) → +∞, U (t φ(n) ) converges to W (0) ∈ X σ for the norm of X.

Let W (t) = (w, ∂ t w)(t) be the maximal solution of the damped wave equation (1.1) corresponding to the initial data W (0). Let t ∈ R, for n large enough t φ(n) + t ≥ 0 and thus U (t φ(n) + t) is well defined and uniformly bounded in X. Since our equation is well posed, the solution is continuous with respect to the initial data. Thus, since U (t φ(n) ) converges to W (0) in X, we have that U (t φ(n) + t) converges to W (t) for all t such that W (t) is well defined. But due to the uniform bound on U (t φ(n) + t), W (t) is uniformly bounded and thus the solution may be extended to a global solution W (t), t ∈ R. In addition, the X σ -bound obtained above for W (0) only depends on the X-bound on U (s) which is uniform due to non-increase of the energy of U (t). Thus, the same arguments applied to the convergence U (t φ(n) + t) → W (t) give the same X σ -bound for W (t) for all t ∈ R. Finally, since the energy of U (t) is non-increasing and non-negative, for any t ∈ R, we must have E(U (t n + t)) -E(U (t n )) → 0 when n → 0 (since t n goes to +∞). This shows that E(W (t)) is constant and finishes the proof.

The case dim(Ω) = 3 is similar once we take into account the constraints given by Proposition 3.2.

□ Proposition 3.4. If dim(Ω) = 2, set σ * = σ h . If dim(Ω) = 3, assume that p ∈ [0, 3) in (1.2) and that σ h ∈ ((p -1)/2, 1) in (1.5) and set σ * = σ h -(p -1)/2. Let σ > 0 and R > 0. Let U n (t) = (u n , ∂ t u n ) a sequence of solutions u n of (1.1) such that (U n (0)) ⊂ X σ and ∥U n (0)∥ X σ ≤ R. Let (t n ) be a sequence of times such that t n → +∞ and let σ ′ ∈ [0, min(σ, σ * )). Then, there exist subsequences (t φ(n) ) and (U φ(n) ) and a solution W (t) = (w, ∂ t w)(t) of (1.1) defined for all t ∈ R, such that ∀t ∈ R , U φ(n) (t φ(n) + t) --------→ n-→+∞ W (t) in X σ ′ .
Moreover, the solution W is globally bounded in X σ ′ and the energy E(W (t)) is constant.

Proof:

The arguments are similar as the ones of the above proof of Proposition 3.3. The term e Atn U n (0) goes to zero in X σ ′ due to Proposition 3.1 because σ ′ < σ. We bound the integral ∫ tn 0 e A(tn-s) F (U n (s)) ds as in the proof of Proposition 3.3: U n (s) is uniformly bounded in X, so F (U n (s)) is uniformly bounded in X η with η < 1 in dimension 2 or η < (3 -p)/2 in dimension 3. Proposition 3.1 together with (1.5) implies that the integral is uniformly bounded in X σ ′ with σ ′ < σ * . The compactness follows by leaving any small amount of regularity in the process. To obtain the convergence to W (t) for all t, we use the same argument as the one of the proof of Proposition 3.3 to first show the convergence in X. Then, the above arguments also show the compactness of U (t φ(n) + t) in X σ ′ and thus the convergence to W (t) also holds in X σ ′ . The last property is the same as the ones of Proposition 3.3. □

The conclusions of Theorem 1.1 then follow from Propositions 3.3 and 3.4 as soon as we can prove that W ≡ 0 for any subsequences of any sequences (t n ) and U n . To this end, notice that E(W (t)) is constant and its derivative (2.4) implies that ∫ Ω γ(x)|∂ t w| 2 = 0 for all time. To formulate this property as a unique continuation property, we set as usual z = ∂ t w and notice that z solves

   ∂ 2 tt z(x, t) = ∆z(x, t) -αz(x, t) -f ′ u (x, w(x, t))z (x, t) ∈ Ω × R z |∂Ω (x, t) = 0 (x, t) ∈ ∂Ω × R z(x, t) ≡ 0 (x, t) ∈ support(γ) × R (3.2)
If this implies z ≡ 0 everywhere, this means that w(x, t) = w(x) is constant in time and

solves ∆w(x) -αw(x) = f (x, w(x)) .
Multiplying by w and integrating, we obtain

∥∇w∥ 2 + α∥w∥ 2 = - ∫ Ω f (x, w(x))w(x) dx .
By the sign Assumption (1.3), this yields w ≡ 0. Thus, it only remains to study this unique continuation property.

Rate of the nonlinear decay: proof of Proposition 1.3

The purpose of this section is to prove Proposition 1.3. When estimating the decay rate of the nonlinear system, we will not exactly need the decay of the linear semigroup but more precisely the decay rate of the linearization at u = 0. This is not difficult since an estimate as (1.4) is a high-frequency result: the behavior of the high frequencies is the difficult part and we only need that the low frequencies do not lie on the imaginary axes.

The polynomial case

The case of polynomial decay is obtained as follows.

Lemma 4.1. Assume the sign hypothesis (1.3) and assume that (1.4) holds with

h(t) = O(t -α ), with α > 1. Set à = A + ( 0 -f ′ u (x, 0) ) = ( 0 Id ∆ -αId -f ′ u (x, 0) -γ(x)
) .

Then, there exists C > 0 such that

∀t ≥ 0 , ∀U 0 ∈ D(A) , ∥e Ãt U 0 ∥ X ≤ C (1 + t) α ∥U 0 ∥ D(A) .
Proof: Due to (1.3), we have that f (x, 0) = 0 and that f ′ u (x, 0) ≥ 0. Since à is a compact perturbation of A, we do not expect that the behavior for high frequencies should be modified. For low frequencies, the sign of f ′ u (x, 0) is sufficient to avoid eigenvalues on the imaginary axes.

To prove rigorously these facts as quickly as possible, we use the results stated in Appendix with

H = L 2 (Ω), L = -∆ + α, B = γ and V = f ′ u (x, 0). Due to Theorem A.4, there exist µ 0 and C > 0 such that ∀µ ∈ R with |µ| ≥ µ 0 , ∥(A -iµ) -1 ∥ L(L 2 ) ≤ C|µ| 1/α .
We now use Proposition B.4 to obtain that the same estimate holds for the resolvents ( Ã -iµ) -1 for large µ. Moreover, for µ in a compact interval, Proposition B.1 ensures that the resolvent ( Ã -iµ) -1 is well defined. Applying Theorem A.4 in the converse sense concludes the proof. □

Proof of the first case of Proposition 1.3:

we assume that the conclusions of Theorem 1.1 hold. In particular, the trajectory of a ball of X σ of radius R is attracted by {0} in X σ ′ for a small enough σ ′ > 0. Thus, it is sufficient to prove that the decay has the same rate as the linear one, as soon as we start from a small ball of X σ ′ of radius ρ > 0 and stay in it.

We consider the linearization of our equation near the stable state u = 0. We set à be as in Lemma 4.1 and

F (u) = (0, f (x, u) -f ′ x (x, u)u). Since H 1+σ ′ (Ω) is embedded in L ∞ (Ω)
and is an algebra, we may bound the derivatives of f and by linearization, for any small δ > 0, we may work with

U (t) in a ball of X σ ′ of radius ρ, which is such that ∥ F (U )∥ X 1 ≤ δ∥U ∥ X .
Let U (t) be a trajectory in the small ball of X σ ′ with ∥U 0 ∥ X σ ≤ ρ. We have

(1 + t) σα U (t) = (1 + t) σα e Ãt U 0 + (1 + t) σα ∫ t 0 e Ã(t-s) F (U (s)) ds .
The term (1 + t) σα e Ãt U 0 is bounded by the linear decay (see Lemma 4.1 and Proposition 3.1). By using the above estimate on F (U (s)), we get that

∥(1 + t) σα U (t)∥ X ≤ C + (1 + t) σα ∫ t 0 C (1 + (t -s)) α δ∥U (s)∥ X ds . Thus, max t∈[0,T ] (1 + t) σα ∥U (t)∥ X ≤ C + δC ( max s∈[0,T ] (1 + s) σα ∥U (s)∥ X ) × max t∈[0,T ] ∫ t 0 ( 1 + t 1 + s ) σα ds (1 + (t -s)) α
where C is a constant independent on T when T goes to +∞ and on the radius ρ of the starting ball when ρ goes to 0. The limit T → +∞ will prove our theorem as soon as we can show that the integral term is bounded uniformly in t. Indeed, up to work with ρ small enough, we may assume that δ is such that the whole last term is less than 1/2 max s∈[0,T ] (1 + s) σα ∥U (s)∥ X and may be absorbed by the left hand side.

To estimate the integral, we use the change of variable τ = (1 + s)/(t + 2), for which 1 + (t -s) = (t + 2)(1 -τ ). We obtain that

I(t) = ∫ t 0 ( 1 + t 1 + s ) σα ds (1 + (t -s)) α = ( 1 + t 2 + t ) σα 1 (t + 2) α-1 ∫ 1-1/(t+2) 1/(t+2) dτ τ σα (1 -τ ) α .
Recall that α > 1 and that σ ≤ 1. The integral ∫ 1 0 dτ τ σα (1-τ ) α does not converge at least close to 1 and if it also diverges close to 0, the blow up is slower or equal to the one occurring close to 1. Thus,

∫ 1-1/(t+2) 1/(t+2) dτ τ σα (1-τ ) α is of order O(t α-1
) when t goes to +∞. This shows that the whole integral I(t) is bounded uniformly in t. □

The exponential case

The following lemma is similar to Lemma 4.1, except that we cannot use the result of Borichev and Tomilov recalled in Theorem A.4. If the decay is not polynomial, then we must accept a logarithmic loss and use the results of Batty and Duyckaerts, Theorems A.2 and A.3.

Lemma 4.2. Assume the sign hypothesis (1.3) and assume that (1.4) holds with h(t) =

O(e -at 1/β ), with a > 0 and β > 0. Set

à = A + ( 0 -f ′ u (x, 0) ) = ( 0 Id ∆ -αId -f ′ u (x, 0) -γ(x)
) .

Then, there exists C > 0 and b > 0 such that

∀t ≥ 0 , ∀U 0 ∈ D(A) , ∥e Ãt U 0 ∥ X ≤ Ce -bt 1/(β+1) ∥U 0 ∥ D(A) .
Proof: As in the proof of Lemma 4.1, we use the results stated in Appendix with

H = L 2 (Ω), L = -∆ + α, B = γ and V = f ′ u (x, 0). Using Theorem A.2, we obtain that ∀µ ∈ R with |µ| ≥ µ 0 , ∥(A -iµ) -1 ∥ L(L 2 ) ≤ C| ln µ| β .
As in Lemma 4.1, we use Proposition B.4 to obtain that the same estimate holds for the resolvents ( Ã-iµ) -1 for large µ and Proposition B.1 to deal with the low frequencies. The difference is that Theorem A.3 yields a logarithmic loss when going back to the estimate of the semigroup (see the definition of M log ), leading to the exponent t 1/(β+1) . □

Remark:

We have seen that there is a logarithmic loss in our estimate. However, in the applications, we will obtain a better result. Indeed, this loss was already present in the original estimate for the linear semigroup because of the additional log in M log of Theorem A.3. In some sense, the above abstract result makes an additional use of the back and forth Theorems A.2 and A.3. We can improve our estimate by a shortcut: we go back to the estimate of the resolvent in the original proof of the linear decay, before the authors apply Theorem A.3, and we directly apply the above arguments to estimate ( Ã -iµ) -1 and then apply Theorem A.3. With this trick, we do not add a second logarithmic loss to the one of the original proof dealing with the linear semigroup. However, we can do this only in the concrete situations and not in an abstract result as Proposition 1.3.

Proof of the second case of Proposition 1.3:

the method is exactly the same as in the first case. The only difference is that, instead of bounding

∫ t 0 ( 1+t 1+s
) σα ds (1+(t-s)) α , we must here bound an integral of the type

I(t) = ∫ t 0 e σc(t γ -s γ ) e -c(t-s) γ ds
for some c > 0 and γ ∈ (0, 1). We set τ = s/t and obtain

I(t) = t ∫ 1 0 e ct γ (σ-στ γ -(1-τ ) γ ) dτ ≤ t ∫ 1 0 e cσt γ (1-τ γ -(1-τ ) γ ) dτ
and by symmetry

I(t) ≤ 2t ∫ 1/2 0 e σct γ (1-τ γ -(1-τ ) γ ) dτ . We notice that τ → 1 -τ γ -(1 -τ ) γ is decreasing for τ ∈ [0, 1/2] since its derivative is γ((1 -τ ) γ-1 -τ γ-1 ) with γ -1 < 0. Moreover, 1 -τ γ -(1 -τ ) γ ∼ -τ γ for small τ . Thus, there exists ν > 0 small enough such that 1 -τ γ -(1 -τ ) γ ≤ -ντ γ for τ ∈ [0, 1/2].
We get

I(t) ≤ 2t ∫ 1/2 0 e σct γ (-ντ γ ) dτ = 2 ∫ t/2 0 e -σcνs γ ds .
The integrand of these last bound is integrable on R + , thus I(t) is bounded uniformly with respect to t. Arguing as in the proof of the polynomial case, this proves the second part of Proposition 1.3. □

Application 1: the open book

In this section, we consider the third example of Figure 1. Let Ω = T 2 be the twodimensional torus and let α > 0 (there is no boundary and so no Dirichlet boundary condition). Assume that

γ(x 1 , x 2 ) = |x 1 | β .
We have the following decay estimate proved in [30,

U 0 ∥ X ≤ C (1 + t) 1+2/β ∥U 0 ∥ D(A) .
Of course, this estimate implies the decay assumptions (1.4) and (1.5) for σ h < 2/(2 + β). Since the support of γ is T 2 , the unique continuation property is trivial and Proposition 3.5 implies that the conclusions of Theorem 1.1 holds in this case. Moreover, Proposition 1.3 provides an explicit decay rate, which is optimal (since it is the same as the linear one). Due to the trivial unique continuation property, this case is far simpler than the general results Theorem 1.1 and 1.2. Nevertheless, it seems the first non-linear stabilization and decay estimate in a case where the linear semigroup has only a polynomial decay.

Theorem 5.2. Consider the damped wave equation

(1.1) in Ω = T 2 and with γ(x 1 , x 2 ) = |x 1 | β (β > 0), α > 0 and f satisfying (1.2) and (1.3). Then, any solution u of (1.1) satisfies ∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 --------→ t-→+∞ 0 .
Moreover, for any R and σ ∈ (0, 1], there exists C R,σ such that, for any solution u with

U 0 ∈ X σ , ∥(u 0 , u 1 )∥ H 1+σ ×H σ ≤ R =⇒ ∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 ≤ C R,σ (1 + t) σ(1+2/β) .

Unique continuation theorems

As proved in Section 3, the last step to prove stabilization is the unique continuation property: if z is a global solution of

   ∂ 2 tt z(x, t) = ∆z(x, t) -αz(x, t) -f ′ u (x, w(x, t))z (x, t) ∈ Ω × R z |∂Ω (x, t) = 0 (x, t) ∈ ∂Ω × R z(x, t) ≡ 0 (x, t) ∈ support(γ) × R (6.1)
then z ≡ 0. Except for the example of Section 5, the property is often difficult to obtain. The purpose of this section is to gather several results yielding this property. The first known result has been proved by Ruiz in [START_REF] Ruiz | Unique continuation for weak solutions of the wave equation plus a potential[END_REF]. It stated the unique continuation property in a bounded domain Ω ⊂ R d as soon as the support of γ contains a neighborhood of the boundary ∂Ω. This result has been generalized in [START_REF] Kazemi | Stability estimates for ill-posed Cauchy problems involving hyperbolic equations and inequalities[END_REF] (see also [START_REF] Lasiecka | global uniqueness and observability in one shot[END_REF] for Neumann boundary conditions). However, this kind of results is not relevant in this paper. Indeed, their geometric settings implies the uniform decay of the semigroup e At and we are interested here in cases where it is not satisfied. We need sharper results.

Unique continuation with coefficients analytic in time

A very general unique continuation property holds if the coefficients of a linear wave equation as (6.1) are analytic in time. This is a consequence of local continuation results proved by by Hörmander in [START_REF] Hörmander | On the uniqueness of the Cauchy problem under partial analyticity assumptions, Geometrical optics and related topics[END_REF] and generalized by Tataru in [START_REF] Tataru | Unique continuation for operators with partially analytic coefficients[END_REF] and also independently proved by Robbiano and Zuily in [START_REF] Robbiano | Uniqueness in the Cauchy problem for operators with partially holomorphic coefficients[END_REF]. These results concern in fact a very general setting but we restrict here the statements at the case of the wave equation. The application to the wave equation and the proof that the local results yield a global one are classical and straightforward, see for example [START_REF] Joly | Stabilization for the semilinear wave equation with geometric control condition[END_REF]Corollary 3.2] for the details.

{ ∂ 2 tt z = ∆z + b(x, t)∂ t z + c(x, t).∇z + d(x, t)z (x, t) ∈ Ω × (-T, T ) z |∂Ω (x, t) = 0 (x, t) ∈ ∂Ω × R . (6.2)
Let O be a non-empty open subset of Ω and assume that z( T,T ) and O T = Ω, then z ≡ 0 everywhere.

x, t) = 0 in O × (-T, T ). Then z(x, 0) ≡ 0 in O T = {x 0 ∈ Ω , d(x 0 , O) < T }. As consequences if z ≡ 0 in O × (-

Unique continuation through pseudo-convex surfaces without boundary

If the coefficients of (6.2) are not analytic in time, the geometry of the problem is more constrained. However, it could still include cases where the geometric control condition of [START_REF] Bardos | Sharp sufficient conditions for the observation, control and stabilization of waves from the boundary[END_REF] does not hold and thus where the semigroup e At is not uniformly stable, see the examples below.

We consider here Hörmander framework (see [START_REF] Hörmander | The analysis of linear partial differential operators. IV: Fourier integral operators[END_REF] for example). The principal symbol of the differential operator of (6.1) is of order two and writes locally

p(x, t, ξ, τ ) = ξ ⊺ .A(x).ξ -|τ | 2
where A(x) is a smooth family of positive definite symmetric matrices coding the Beltrami Laplacian operator in a local chart. Let ϕ(x, t, ξ, τ ) be a locally C 1 -function, we introduce the Poisson bracket

H p (ψ) = {p, ψ} = ∇ ξ p∇ x ψ + ∂ τ p∂ t ψ -∇ x p∇ ξ ψ -∂ t p∂ τ ψ .
Let ψ be a smooth function defined in a neighborhood O ⊂ R d+1 of (x 0 , t 0 ). Assume that (∇ x ψ, ∂ t ψ)(x 0 , t 0 ) ̸ = 0 so that Σ = {(x, t), ψ(x, t) = 0} defines a smooth hypersurface near (x 0 , t 0 ). Definition 6.2. The local hypersurface Σ is said to be non-characteristic at (x 0 , t 0 ) if

p(x 0 , t 0 , ∇ψ(x 0 , t 0 ), ∂ t ψ(x 0 , t 0 )) ̸ = 0 .
Moreover, Σ is said to be strongly pseudo-convex at (x 0 , t 0 ) if for any (ξ, τ ) ̸ = 0 such that p(x 0 , t 0 , ξ, τ ) = 0 and H p (ψ)(x 0 , t 0 ) = 0, we have

H 2 p (ψ)(x 0 , t 0 ) > 0 .
Notice that the above definition of strongly pseudo-convexity is adapted to the case of a real differential operator of order two. Thus it is perfectly adapted to the situation of this paper where the wave operator is p = ξ ⊺ .A(x).ξ -|τ | 2 . However, we emphasize that, in the general case, the assumption of pseudo-convexity is more complex, see [START_REF] Hörmander | The analysis of linear partial differential operators. IV: Fourier integral operators[END_REF].

The geometrical interpretation of Definition 6.2 is as follows. First, the fact that the surface is non-characteristic says that |∂ t ψ| 2 ̸ = ∇ψ ⊺ .A(x).∇ψ. This means that the surface is not moving at the exact same speed as the sound waves.

The pseudo-convexity is slightly more involved. Consider the total Hamiltonian flow σ → φ σ defined by

φ 0 (x, t, ξ, τ ) = (x, t, ξ, τ ) ∂ σ φ σ (x, t, ξ, τ ) = (∇ ξ p, ∂ τ p, -∇ x p, -∂ t p)(φ σ ) .
Since p is independent of t, τ is constant and thus t(σ) = t -2τ σ, meaning that σ is a simple new parametrization of time. Moreover, ∇ ξ p and ∇ x p are independent of t and τ . Thus, (x, ξ)(σ) follows the geodesic flow

∂ σ (x, ξ) = (∇ ξ g, -∇ x g)(x, ξ)
where g(x, ξ) = ξ ⊺ .A(x).ξ is the symbol of the local metric. Assume that p(x, t, ξ, τ ) = 0 at σ = 0. The Hamiltonian being conserved, we always have p(x, t, ξ, τ ) = 0 and |τ | 2 = ξ ⊺ A(x)ξ is constant in σ: the point x(σ) is moving along a geodesic of the metric at a speed which is of constant norm |τ | with respect to the metric. Let h be a function of (x, t, ξ, τ ), then the Poisson bracket {p, h} is the derivative at

σ = 0 of h(φ σ (x, t, ξ, τ )). Thus {p, h} = {g, h} -2τ ∂ t h = {g, h} + ∂ σ h
where {g, h} is the derivative along the geodesic (x, ξ)(σ) of the metric starting at x with speed ξ. Thus, the strongly pseudo-convexity condition H p ψ = 0 ⇒ H 2 p ψ > 0 means that if a geodesic of the surface is tangent to Σ in the space-time sense, then it must be contained in a non-degenerated sense in the half-space ψ(x, t) > 0 for t ̸ = 0. Finally notice that if ψ does not depend on time t (as in Definitions 6.4 and 6.6), then the strongly pseudo-convexity is a classical strong convexity: if a classical geodesic of the metric g is tangent to the surface ψ(x) = 0, it must be contained in the half-space ψ(x) > 0.

Theorem 28.4.3 of [START_REF] Hörmander | The analysis of linear partial differential operators. IV: Fourier integral operators[END_REF] mainly comes from [START_REF] Lerner | Unicité de Cauchy pour des opérateurs de type principal[END_REF] and is stated as follows. Let O be a small open neighborhood of a point (x 0 , t 0 ) in R d × R and let A(x) be a smooth family of positive definite symmetric matrices. Let b, c and d be bounded coefficients. Assume that z is a mild solution of

∂ 2 tt z = div A(x)∇z + b(x, t)∂ t z + c(x, t).∇z + d(x, t)z (x, t) ∈ O . (6.3)
Let Σ = {(x, t), ψ(x, t) = 0} be a smooth surface containing (x 0 , t 0 ) which is non-characteristic and strongly pseudo-convex in the sense of Definition 6.2. Then, if u(x, t) = 0 for all (x, t) ∈ O such that ψ(x, t) ≥ 0, we have u(x, t) ≡ 0 in a neighborhood of (x 0 , t 0 ). (i) the family of surfaces is smooth in the sense that it is locally described as level sets {x, ψ λ (x) = 0} where (x, λ) → ψ λ (x) is a local smooth function with ∇ x ψ λ ̸ = 0.

(ii) each surface is globally oriented in the sense that there exist disjoint sets

Σ ± λ such that locally {x ∈ Ω, ±ψ λ (x) > 0} ⊂ Σ ± λ and such that Ω = Σ - λ ∪ Σ λ ∪ Σ + λ .
(iii) for each λ, (x, t) → ψ λ (x) is pseudo-convex in the sense of Definition 6.2 as a function independent of t. Equivalently, Σ - λ is locally strictly convex in a neighborhood of its boundary Σ λ for the metric g: for each x ∈ Σ λ , a geodesic through x which is tangent at Σ λ is locally included in Σ + λ , x excepted.

(iv) the surfaces Σ λ are compact and have no boundary or equivalently do not meet ∂Ω.

A typical example of such oriented pseudo-convex foliation without boundary is given in Figure 2.

By a classical argument, we may state a global version of Theorem 6.3 as follows.

Theorem 6.5. Let Ω be a smooth compact manifold (with or without boundaries) and let ω ⊂ Ω be an open set. Assume that there exists an oriented pseudo-convex foliation without boundary (Σ λ ) λ∈[0, [START_REF] Anantharaman | Entropy and the localization of eigenfunctions[END_REF] of Ω in the sense of Definition 6.4. Also assume that

Σ + 0 ⊂ ω and ω ∪ ( ∪ λ∈[0,1) Σ + λ )
covers Ω up to a set of zero measure. Let b, c and d be bounded coefficients. Assume that z is a global mild solution of

{ ∂ 2 tt z = ∆z + b(x, t)∂ t z + c(x, t).∇z + d(x, t)z (x, t) ∈ Ω × R z |ω (x, t) = 0 (x, t) ∈ ω × R
with any suitable boundary conditions on ∂Ω such that the wave equation is well-posed and where ∆ is the Laplace-Beltrami operator related to Ω. Then z ≡ 0 everywhere.

Proof: We will show that z(•, t = 0) vanishes in ω ∪Σ + λ , for all λ ∈ [0, 1), which shows that z(t = 0) ≡ 0 and thus that z ≡ 0 due to the uniqueness properties of the wave equation. By assumption, z ≡ 0 in Σ + 0 ⊂ ω. Let λ 0 ∈ (0, 1) and let h α,T (t) = α(1 -t 2 /T 2 ). We consider the family of surfaces t ∈ [-T, T ] → Σ h α,T (t) which is locally parametrized by functions (x, t) → ψ h α,T (t) (x). Notice that it is a smooth family of smooth surfaces since due to Assumption (i) of Definition 6.4. Also notice that the larger is T , the smaller are the derivatives of these functions with respect to t. By assumption, each function x → ψ λ (x) is non-characteristic and strongly pseudo-convex as a function independent of t. By compactness, there exists T large enough such that (x, t) → ψ h α,T (t) (x) defines local surfaces which are non-characteristic and pseudo-convex for all α ∈ [0, λ 0 ] and t ∈ [-T, T ]. The parameter T is fixed in the remaining part of the proof and we may omit it in the notations.

Notice that, for any α, the family of set t ∈ [-T, T ] → Σ + h α (t) starts inside ω at t = -T and finishes inside ω at t = T . Moreover, for any small α, these sets always stay inside ω where z vanishes. Assume that there exist (x, t) and α ∈ [0, λ 0 ] such that x ∈ Σ + h α (t) and z(x, t) ̸ = 0. We set

α 0 = min{α ∈ (0, λ 0 ] , ∃t ∈ [-T, T ], ∃x ∈ Σ + h α (t) such that z(x, t) ̸ = 0} . (6.4)
By continuity, we know that z(x, t) = 0 for all x ∈ Σ + h α 0 (t) . Moreover, there exists t 0 ∈ (-T, T ) and x 0 ∈ Σ hα 0 (t 0 ) such that z is not identically zero in any neighborhood of (x 0 , t 0 ). Indeed, otherwise, by compactness, we may extend the set where z vanishes and contradict (6.4).

To conclude, it remains to use the local unique continuation property of Theorem 6.3 at (x 0 , t 0 ) with the time-space surface defined by (x, t) → ψ h α 0 (t) (x). The continuation implies that z vanishes near (x 0 , t 0 ) which contradicts the construction. Thus, z(x, t) = 0 for all t ∈ [-T, T ] and x ∈ ∪ α Σ + hα(t) . In particular z(•, t = 0) ≡ 0 in

∪ α Σ + hα(0) = ∪ λ≤λ 0 Σ + λ .
Since these arguments hold for all λ 0 < 1 and since ω ∪ λ∈[0,1) Σ + λ is Ω up to a set of measure zero, we have that z(•, t = 0) ≡ 0 in Ω. Well-posedness of the linear wave equation concludes that z ≡ 0 everywhere. □

Notice that, as it is stated, this unique continuation result needs an infinite time to be efficient, where Theorem 6.1 only need a finite explicit time. In fact, a careful look to the proof shows that a finite time is sufficient once we know that the family of surface is pseudo-convex and non-characteristic in a uniform way. However, such a bound of convexity is difficult to obtain in general cases and may be even impossible as for the example studied in Section 7.

A typical example of application is given in Figure 2: if Ω is a sphere and ω covers more than an hemisphere, then if z is a global solution of a linear wave equation which vanishes in ω for all times, then z ≡ 0. Notice that, in this case, the family of surfaces is uniformly pseudo-convex and the unique continuation holds in fact in finite time even if z is not a global in time solution.

Unique continuation through pseudo-convex surfaces with boundary

The case where the pseudo-convex surfaces Σ λ meet the boundary ∂Ω is more involved. Theorem 6.3 has been generalized to this case by Tataru (see [START_REF] Tataru | Unique continuation for solutions to PDE's; between Hörmander's theorem and Holmgren's theorem[END_REF][START_REF] Tataru | Carleman estimates and unique continuation for solutions to boundary value problems[END_REF][START_REF] Tataru | Unique continuation for operators with partially analytic coefficients[END_REF]). The boundary conditions are more difficult to describe geometrically, so we will only deal here with the case of flat geometry, that is g(x, ξ) = |ξ| 2 , and the case of Dirichlet boundary condition. Definition 6.6. A family of surfaces (Σ λ ) λ∈[0,1) is an oriented pseudo-convex foliation with boundary in a flat manifold Ω if:

(i) the family of surfaces is smooth in the sense that it is locally described as level sets {x, ψ λ (x) = 0} where (x, λ) → ψ λ (x) is a local smooth function with ∇ x ψ λ ̸ = 0.

(ii) each surface is globally oriented in the sense that there exist disjoint sets

Σ ± λ such that locally {x ∈ Ω, ±ψ λ (x) > 0} ⊂ Σ ± λ and such that Ω = Σ - λ ∪ Σ λ ∪ Σ + λ . (iii) for each λ, (x, t) → ψ λ (x) is pseudo-convex in
the sense of definition 6.2 as a function independent of t. Equivalently, Σ - λ is locally strictly convex in a neighborhood of its boundary: the tangent space to Σ λ at x 0 is locally included in Σ + λ , x 0 excepted. (iv) if a surface Σ λ meet ∂Ω at x, then ∂ ν ψ λ (x) < 0. Equivalently, the angle formed by Σ λ and ∂Ω in the region Σ - λ is strictly less than π/2. A typical example of such oriented pseudo-convex foliation with boundary is given in Figure 3. Notice that the condition at the boundary is consistent with the one inside the domain. Indeed, the geodesics are straight lines which bounce at the boundary according to Newton's laws. Geometrically, we ask that any geodesic either crosses Σ λ in a transversal way, or stay locally inside Σ + λ . By the same arguments as the ones in the proof of Theorem 6.5 and using the result of Tataru, we obtain a global unique continuation result. Theorem 6.7. Let Ω ⊂ R d be a compact domain and let ω ⊂ Ω be an open set. Assume that there exists an oriented pseudo-convex foliation (Σ λ ) λ∈[0, [START_REF] Anantharaman | Entropy and the localization of eigenfunctions[END_REF] of Ω in the sense of Definition 6.6. Also assume that Σ + 0 ⊂ ω and ω ∪

( ∪ λ∈[0,1) Σ + λ )
covers Ω up to a set of zero measure.

Let b, c and d be bounded coefficients. Assume that z is a global mild solution of

   ∂ 2 tt z = ∆z + b(x, t)∂ t z + c(x, t).∇z + d(x, t)z (x, t) ∈ Ω × R z |∂Ω (x, t) = 0 (x, t) ∈ ∂Ω × R z |ω (x, t) = 0 (x, t) ∈ ω × R
Then z ≡ 0 everywhere.

A typical example of application is given in Figure 3: if Ω is a disk and ω covers more than half of the boundary, then if z is a global solution of a linear wave equation which vanishes in ω for all times, then z ≡ 0.

Proof of Theorem 1.2

Theorem 1.2 is then a direct consequence of the unique continuation results stated in this Section: Proposition 3.5 and Theorems 6.5 and 6.7 imply Theorem 1.2.

Application 2: the peanut of rotation

We consider in this section the example of the peanut of rotation: a two-dimensional manifold where a central part is equivalent to the cylinder {x = (y, θ) ∈ (-1, 1) × S} endowed with the metric g(y, θ) = dy 2 + cosh 2 (y)dθ 2 (see Figure 1). The damping γ is assumed to be positive, except in a part x ∈ (-ℓ, ℓ) of the central part (ℓ ∈ (0, 1)). The decay of the linear damped wave semigroup has been established in [START_REF] Christianson | From resolvent estimates to damped waves[END_REF] and [START_REF] Schenck | Exponential stabilization without geometric control[END_REF]. In the setting of the peanut of rotation, there exist two positive constants C and λ such that the semigroup e At satisfies

∥e At U 0 ∥ X ≤ Ce -λ √ t ∥U 0 ∥ D(A) .
The decay rate of Theorem 7.1 obviously satisfies (1.4) and (1.5). Thus, once the unique continuation property is obtained, Proposition 3.5 yields the conclusion of Theorem 1.1 for the framework of the peanut of rotation. To obtain the unique continuation property, we will apply Theorem 6.5 with the family of pseudo-convex surfaces shown in Figure 4. Applying Theorem 1.2 and the ideas of Proposition 1.3, we obtain the following result. λ . Notice that the set Σ 0 consists in both exterior disk and the disks Σ λ get closer to x = 0 when λ get closer to 1. Thus, the central circle x = 0 is not included in ∪ λ∈[0,1) Σ + λ but this is not important since this set is of zero measure.

∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 --------→ t-→+∞ 0 . 21 Σ λ Σ λ Σ + λ Σ + λ Σ - λ y = 0
Moreover, for any R and σ ∈ (0, 1], there exists C R,σ such that, for any solution u with

U 0 ∈ X σ , ∥(u 0 , u 1 )∥ H 1+σ ×H σ ≤ R =⇒ ∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 ≤ C R,σ e -σ λ√ t
where λ is the linearized rate given in Lemma 4.2.

Proof: Let us first formally check that the family of disks introduced in Figure 4 is a suitable pseudo-convex foliation without boundary. We use the cylindrical coordinates x = (y, θ) ∈ (-1, 1)×S with associated tangent variables ξ = (ζ, Θ). By symmetry, we only consider the right-hand-side circles which are defined by ψ λ (x) = 0 with ψ λ (x) = y-(1-λ).

The circle Σ 0 corresponding to y = 1 is in the interior of the region ω where the damping is positive. When λ get closer to 1, the circle Σ λ get closer to y = 0. We are obviously in the setting of Definition 6.4 and thus of Theorem 1.2, except maybe for the assumption of strong pseudo-convexity. We already give a geometrical insight of this assumption, but let us check it formally. The local metric is given by g(y, θ) = dy 2 + cosh 2 (y)dθ 2 . The Laplace-Beltrami operator is thus given by

∆ g = ∂ 2 yy + 2 tanh(y)∂ y + 1 cosh 2 (y) ∂ 2 θθ .
The principal part of the wave operator is then

p(y, θ, t, ζ, Θ, τ ) = |ζ| 2 + 1 cosh 2 (y) |Θ| 2 -|τ | 2 .
Thus H p (ψ λ ) = 2ζ and

H 2 p (ψ λ ) = 4 sinh(y) cosh(y) 3 |Θ| 2 .
The pseudo-convexity condition is then checked. Indeed, if H p (ψ λ ) = 0 then ζ = 0 and since ξ = (ζ, Θ) must be non-zero, we must have Θ ̸ = 0. As ψ λ (y, θ) = 0 with λ ∈ [0, 1), we have y > 0 and thus H 2 p (ψ λ ) > 0. Looking carefully to the computations, one notes that, in fact, we only need that the radius cosh(y) of the cylindrical part is increasing for y > 0 and decreasing for y < 0 to obtain the unique continuation property.

The above arguments show the stabilization of the semilinear damped wave equation. Moreover, Proposition 3.4 shows the uniform convergence to 0 in X σ ′ for initial data in a more regular space X σ (σ > σ ′ ).

To obtain the decay estimate of Theorem 7.2, we argue as in the proof of the second case of Proposition 1.3. However, we claim that we can avoid the loss in the power by following the remark below Lemma 4.2. Indeed, Theorem 5.1 of [START_REF] Christianson | From resolvent estimates to damped waves[END_REF] implies in our framework that, for large µ,

∥(-∆ + µ 2 -iµγ) -1 ∥ L(L 2 ) ≤ C ln |µ| |µ| .
We argue as in the proof of Lemma 4.2. Applying Propositions B.1, B.2 and B.4, we obtain

∥( Ã -iµ) -1 ∥ L(L 2 ) ≤ C | ln µ| .
and by Theorem A.3, we get the conclusion of Lemma 4.2 in the form

∥e Ãt U 0 ∥ X ≤ Ce -b √ t ∥U 0 ∥ D(A) .
In other words, we avoid an additional logarithmic loss by directly dealing with the estimates of [START_REF] Christianson | From resolvent estimates to damped waves[END_REF] instead of using the back and forth implications of Theorems A.2 and A.3. It is then sufficient to follow the proof of the exponential case of Proposition 1.3 with the exponent γ = 1/2. □

Decay estimate in the disk with holes

In the previous examples of application, the decay of the semigroups was explicitly written in previous papers. In the case of a disk with several holes, we are not aware of a paper where an explicit decay is written. The corresponding scattering problem has been studied by Ikawa in [START_REF] Ikawa | Decay of solutions of the wave equation in the exterior of two convex obstacles[END_REF][START_REF] Ikawa | Decay of solutions of the wave equation in the exterior of several convex bodies[END_REF]. Many further studies have been published. In this article, we will use an estimate and a "black box argument" introduced by Burq and Zworski in [START_REF] Burq | Geometric control in the presence of a black box[END_REF].

Combining them with the results in Appendix, we obtain the following decay.

O i ∩ O j = ∅ for i ̸ = j, (b) the convex hull convhull(∪ i O i ) of the obstacles is contained in O, (c) no obstacle is in the convex hull of two others, that is that convhull(O i ∪ O j ) ∩ O k = ∅ for i, j,
At U 0 ∥ X ≤ Ce -λt 1/3 ∥U 0 ∥ D(A)
with C and λ two positive constants.

A typical domain consists in a smooth domain with several small holes as in Figure 5. Typically, if O i are balls of center c i and radius r i and if there is no triplet of aligned centers, then Assumption (d) holds for r i small enough since κ becomes large whereas L stay bounded.

Remarks:

• We do not claim that the decay rate e -t 1/3 is optimal. In fact, our proof uses rough arguments leading to logarithmic losses. We strongly believe that the right decay rate is e -t 1/2 . A strategy of proof may be to follow the arguments of Datchev and Vasy [START_REF] Datchev | Propagation through trapped sets and semiclassical resolvent estimates[END_REF], adding the presence of a boundary. However, this improvement is not central and would use techniques too far from the spirit of this paper.

• It is also certainly possible to relax the assumption about γ to the following weaker assumption. There exists a neighborhood K of O = ∪ i O i such that any geodesic ray starting in Ω \ K reach a point where γ ≥ ε > 0 before meeting K, for at least the backward or the forward flow.

Based on [START_REF] Ikawa | Decay of solutions of the wave equation in the exterior of several convex bodies[END_REF], the following estimate appears in [

(∆ R d \O + µ 2 ) -1 from Im(µ) > 0, where ∆ R d \O is the Laplacian operator on the exterior domain R d \ O = R d \ (∪ i O i ) with Dirichlet boundary condition. Then, for any cut-off function χ ∈ C ∞ c (R d ), we have ∥χR bb (µ)χ∥ L(L 2 (R d \O)) ≤ C ln(1 + |µ|) 1 + |µ| .
Using this black box in the same spirit of [START_REF] Burq | Geometric control in the presence of a black box[END_REF], we obtain the following observation estimate.

Lemma 8.3. Assume that the assumptions of Theorem 8.1 hold. Then there exists

C > 0 such that ∀µ ∈ R , ∀u ∈ D(∆) , ∥u∥ L 2 (Ω) ≤ C ln |µ| |µ| ∥(∆ + µ 2 )u∥ L 2 (Ω) + Cln |µ|∥ √ γu∥ L 2 (Ω)
where ∆ is the Laplacian operator on the bounded domain with holes Ω = O \ O with Dirichlet boundary conditions.

Proof: By compactness, there exists ε > 0 such that γ ≥ ε > 0 in a neighborhood of the exterior boundary ∂O. Let χ ∈ C ∞ 0 (R d ) be a smooth cut-off function equal to 1 in a neighborhood of convhull(∪ i O i ), equal to 0 outside O and such that 1 -χ is supported where γ ≥ ε > 0. This is possible by Assumption (b). We have immediately

∀u ∈ L 2 (Ω) , ∥(1 -χ)u∥ L 2 (Ω) ≤ 1 √ ε ∥ √ γu∥ L 2 (Ω) . (8.1)
Let χ 0 be a smooth cut-off function supported in Ω so that χ 0 ≡ 1 in a neighborhood of the support of χ and let χ 1 be another smooth cut-off function supported in Ω so that

χ 1 ≡ 1 in a neighborhood of the support of χ 0 . For all u ∈ H 2 (Ω) ∩ H 1 0 (Ω), we extend χ i u as a function in H 2 (R d \ O) ∩ H 1 0 (Ω).
Regarding χ i u, applying ∆ R d \O or ∆ gives the same result. We can thus apply the "black box" estimate of Proposition 8.2 as follows.

∥χu∥ L 2 (Ω) = ∥χχ 2 0 u∥ L 2 (Ω) = ∥χχ 0 R bb (µ)(∆ + µ 2 )χ 0 u∥ L 2 (Ω) ≤ ∥χχ 0 R bb (µ)χ 0 (∆ + µ 2 )u∥ L 2 (Ω) + ∥χχ 0 R bb (µ)[∆, χ 0 ]u∥ L 2 (Ω) ≤ ∥χ 0 R bb (µ)χ 0 ∥ L(L 2 (R d \O)) ∥(∆ + µ 2 )u∥ + ∥χχ 0 χ 1 R bb (µ)χ 1 [(∆ + µ 2 ), χ 0 ]u∥ ≤ ∥χ 0 R bb (µ)χ 0 ∥ L(L 2 (R d \O)) ∥(∆ + µ 2 )u∥ + ∥χ 1 R bb (µ)χ 1 ∥ L(L 2 (R d \O)) ∥[∆, χ 0 ]u∥ ≤ C ln(1 + |µ|) 1 + |µ| ( ∥(∆ + µ 2 )u∥ L 2 (Ω) + ∥u∥ H 1 (supp(∇χ 0 )) ) .
By interpolation and elliptic regularity, we have

∥u∥ 2 H 1 (supp(∇χ 0 )) ≤ C∥u∥ L 2 (supp(∇χ 0 )) ∥u∥ H 2 (supp(∇χ 0 )) ≤ C∥u∥ L 2 (supp(∇χ 0 )) ( ∥∆u∥ L 2 (supp(∇χ 0 )) + ∥u∥ L 2 (supp(∇χ 0 )) ) ≤ C∥u∥ L 2 (supp(∇χ 0 )) ( ∥(∆ + µ 2 )u∥ L 2 (supp(∇χ 0 )) + (1 + µ 2 )∥u∥ L 2 (supp(∇χ 0 )) ) ≤ C ( 1 2(1 + µ 2 ) ∥(∆ + µ 2 )u∥ L 2 (supp(∇χ 0 )) + 3 2 (1 + µ 2 )∥u∥ L 2 (supp(∇χ 0 ))
)

Since the support of ∇χ 0 is included in the place where γ ≥ ε > 0, both previous estimates yield

∥χu∥ L 2 (Ω) ≤ C ln(1 + |µ|) 1 + |µ| ( ∥(∆ + µ 2 )u∥ L 2 (Ω) + (1 + |µ|)∥ √ γu∥ L 2 (Ω)
) .

With (8.1)

Application 3: the disk with two holes

In the previous section, we have obtained a sufficiently fast decay rate for the semigroup of the damped wave equation in a disk with several holes as in Figure 1. If we prove the unique continuation property of Proposition 3.5 in this situation, then we would obtain the desired stabilization. To obtain the unique continuation property, we would like to use Theorem 6.7, that is to exhibit an oriented pseudo-convex foliation (Σ λ ) λ∈[0,1) with Σ + 0 included in a neighborhood of the boundary and ∪ λ Σ + ∞ covering almost all Ω. This is possible in the case where there is at most two holes in the disk and impossible if there are more holes, as shown in Figure 6. Left, the disk with two holes may be covered by a pseudo-convex foliation starting in a neighborhood of the exterior boundary, except for a single line, which is of measure zero. Right, an attempt to cover a disk with three holes with a pseudo-convex foliation. We easily notice that the central zone cannot be covered and Theorem 6.7 cannot be applied in this case.

0000
Thus, in the case where there is only two holes, Theorem 6.7 enables to use Proposition 3.5 and to obtain the conclusions of Theorem 1.1. Moreover, notice that in this case, there is no technical assumptions, neither (c) nor (d), in Theorem 8.1. We thus obtain the following result, as an application of Theorem 1.2, Proposition 1.3 and mutatis mutandis the same use of the remark below Lemma 4.2 as in the proof Theorem 7.2. The proof is left to the reader.

∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 --------→ t-→+∞ 0 .
Moreover, there exists λ such that, for any R and σ ∈ (0, 1], there exists C R,σ such that, for any solution u with U 0 ∈ X σ ,

∥(u 0 , u 1 )∥ H 1+σ ×H σ ≤ R =⇒ ∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 ≤ C R,σ e -σ λt 1/3 .
10 Analytic regularization and proof of Theorem 1.1

In Figure 6, we have seen that, in some situations, the unique continuation property of Lerner-Robbiano-Hörmander stated in Theorem 6.7 is not useful. To deal with these situations, we need another unique continuation property: Theorem 6.1 of Robbiano-Zuily-Hörmander to (3.2). This is only possible if the coefficients of (3.2) are analytic with respect to the time t. Thus, we need to choose f ′ u analytic with respect to u and to prove that the global solution w appearing in (3.2) is analytic in time. This is the basic idea leading to Theorem 1.1.

However, even if f is analytic, the damped wave equation does not regularize its solutions. To overcome this problem, we use the following fact known since the work of Hale and Raugel in [START_REF] Hale | Regularity, determining modes and Galerkin methods[END_REF]: the globally bounded solutions of the damped wave equation are as smooth as the non-linearity f . This asymptotic regularization property is linked to the asymptotic smoothness or compactness property (see Section 3). The idea that this asymptotic smoothing of the damped wave equation may be used to apply analytic unique continuation theorems originates from the work of Hale and Raugel, even if they did not publish this idea. The first published occurrence appears in [START_REF] Joly | New examples of damped wave equations with gradient-like structure[END_REF] (see also [START_REF] Joly | Stabilization for the semilinear wave equation with geometric control condition[END_REF]).

The article [START_REF] Hale | Regularity, determining modes and Galerkin methods[END_REF] contains several abstract theorems. They apply for linear semigroup with uniform decay, that is ∥e At ∥ L(X) ≤ M e -λt . The purpose of the present article is to study cases where this uniform decay fails, so [START_REF] Hale | Regularity, determining modes and Galerkin methods[END_REF] does not directly apply: we need to extend its results in our case where the semigroup has a weaker decay. Extending these results in the most general framework will lead to heavy notations and assumptions. That is why, we only consider here damped wave equations in a simple setting and in low dimension.

Analytic regularization of global bounded solutions

Let d = 2 or 3 and let Ω be a smooth manifold of dimension d with or without boundary and such that Ω is compact. Let γ ∈ C 1 (Ω, R) be a nonnegative damping, let ∆ be the Laplacian operator with Dirichlet boundary condition and let f ∈ C ∞ (Ω × R, R) be a smooth nonlinearity. We assume that f is of polynomial type in the sense that there exist p > 0 and C > 0 such that (1.2) holds. We consider global solutions of the damped wave equation

{ ∂ 2 tt u(x, t) + γ(x)∂ t u(x, t) = ∆u(x, t) -αu(x, t) -f (x, u(x, t)) (x, t) ∈ Ω × R u |∂Ω (x, t) = 0 (x, t) ∈ ∂Ω × R (10.1)
We use the notations of Section 2. In particular, for σ ∈ [0, 1], we set

X σ = (H 1+σ (Ω) ∩ H 1 0 (Ω)) × H σ 0 (Ω), X 0 = X = H 1 0 (Ω) × L 2 (Ω) and X 1 = D(A). Also notice that H 1+σ (Ω) is a subspace of C 0 (Ω) for d = 2 and σ > 0 or for d = 3 and σ > 1/2.
The purpose of this section is to prove the following result Theorem 10.1. Assume that the above setting holds, in particular assume that γ ∈ C 1 (Ω, R). Let U (t) ∈ C 0 (R, X) be a mild solution of (10.1) and assume moreover that (i) There exists σ 0 ∈ (0, 1) such that U (t) is defined for all t ∈ R and uniformly bounded in X σ 0 , that is that there exists C > 0 such that

∀t ∈ R , ∥U (t)∥ X σ 0 ≤ C . If d = 3, assume in addition that σ 0 > 1/2.
(ii) The linear semigroup satisfies the decay estimate

∀U 0 ∈ D(A) , ∥e At U 0 ∥ X ≤ M t β ∥U 0 ∥ D(A) (10.2)
with β > max(2p, 2/(1 -σ 0 )), where p is the polynomial growth of f in (1.2).

(iii) The function u ∈ R → f (x, u) ∈ R is analytic with respect to u,
Then the mapping t ∈ R → U (t) ∈ X σ 0 is analytic with respect to t.

We expect that the condition β > 2p in Assumption (ii) is not optimal: at least β > p should be sufficient if we get rid of the losses in the too general proofs of the auxiliary results in appendix and β > 2/(1 -σ 0 ) could be omitted by assuming more regularity on γ. Since our concrete applications have a linear decay of the type O(e -t β ), we let these probable improvements for later study.

The remaining part of this section is devoted to the proof of Theorem 10.1.

• Step 1: a trick to satisfy the boundary condition

We would like that u → f (•, u) maps H 1+σ (Ω) ∩ H 1 0 (Ω) into H σ 0 (Ω). For d = 2 or d = 3 and σ > 1/2, the regularity part is trivial since H 1+σ (Ω) ⊂ C 0 (Ω) is a normed algebra. However, the boundary condition is not necessarily fulfilled since f (x, 0) may be different from 0. Let us describe here a trick introduced in [START_REF] Hale | Regularity, determining modes and Galerkin methods[END_REF] to deal with this problem. Let e be the (time-independent) solution of ∆e

-αe = f (x, 0) in H 2 (Ω) ∩ H 1 0 (Ω). We notice that ũ = u -e solves ∂ 2 tt ũ + γ(x)∂ t ũ = ∆ũ -αũ -f (x, ũ + e(x)) + f (x, 0) . If we set f (x, ũ) = f (x, ũ + e(x)) -f (x, 0)
we obtain a function f as smooth as f which is also analytic with respect to ũ. Moreover, f (x, 0) = 0 for x ∈ ∂Ω, which shows that ũ → f (•, ũ) maps H 1+σ (Ω) ∩ H 1 0 (Ω) into H σ 0 (Ω), including the boundary condition. If we prove Theorem 10.1 for ũ and f , it clearly yields Theorem 10.1 for u and f . Thus, we may assume that f (x, 0) = 0 at the boundary and we will forget the tilde sign to lighten the notations in what follows. In particular, F = (0, f ) maps X σ into X 1+σ .

• Step 2: the decay of the high-frequencies semigroups

From now on, we also use the notations of the appendices, Sections B and C. In particular, we set L = -∆ + α. Since L is self-adjoint, positive and with compact resolvent, there exists an orthonormal basis (ϕ k ) k≥0 of eigenfunctions of L corresponding to the eigenvalues (λ k ) k≥0 . As in Section C, we introduce the high-frequencies truncations Q n , that are the projectors on the space span{ϕ k , k ≥ n} and we set

Q n = (Q n , Q n ) on X.
The proof of Theorem 10.1 is based on the following splitting. We introduce P n = Id -Q n and P n = Id -Q n , which are low-frequencies projections with finite rank. We consider the splitting

U = P n U + Q n U := V + W . ( 10.3)
Then (10.1) writes

{ ∂ t V = (P n AP n )V + (P n AQ n )W + P n F (V + W ) ∂ t W = (Q n AQ n )W + (Q n AP n )V + Q n F (V + W ) (10.4)
As a consequence of the results in the appendices, we have the following decay estimates.

Proposition 10.2. Assume that Hypothesis (ii) of Theorem 10.1 holds. Then, for all σ ∈ [0, 1), ν ∈ (0, 1 -σ] and ε > 0, there exists C > 0 such that, for all n ∈ N,

∀t ≥ 0 , ∀U ∈ Q n X σ+ν e QnAQnt U X σ ≤ C t νβ/2-ε ∥U ∥ X σ+ν .
(10.5)

Proof: Using the arguments of Proposition 3.1, it is sufficient to obtain the decay for σ = 0 and ν = 1, that is the decay estimate from D(A) into X.

In the appendices, we recall the result of Borichev and Tomilov in [START_REF] Borichev | Optimal polynomial decay of functions and operator semigroups[END_REF] (see Theorem A.4). In the context of the damped wave equation, we may also consider Proposition 2.4 of [START_REF] Anantharaman | Sharp polynomial decay rates for the damped wave equation on the torus[END_REF] by Anantharaman and Léautaud. We obtain that Hypothesis (ii) of Theorem 10.1 implies the estimate

∥u∥ L 2 ≤ C ( |µ| 1/β-1 ∥P B (µ)u∥ L 2 + |µ| 1/2β ∥ √ Bu∥ ) (10.6)
with B is the multiplication by γ and P B (µ) = -L -iµB + µ 2 Id. Note that the unique continuation assumed in [2, Proposition 2.4] is satisfied because Hypothesis (ii) is satisfied. Thus, we obtain that

∥u∥ L 2 ≤ C ( |µ| 1/β-1 ∥P (µ)u∥ L 2 + ( |µ| 1/β ∥ √ B∥ L(L 2 ) + |µ| 1/2β ) ∥ √ Bu∥ ) (10.7)
with

P (µ) = -L + µ 2 Id. Then, Propositions C.2 and C.1 show that ∀n ∈ N , ∀µ ∈ R , (Q n AQ n -iµ) -1 L(QnX) ≤ K |µ| b
with b = 2/β. At this point, we may use Theorem A.4 in appendix to obtain the decay of the linear semigroup. However, this result of [START_REF] Borichev | Optimal polynomial decay of functions and operator semigroups[END_REF] (as the one of [START_REF] Anantharaman | Sharp polynomial decay rates for the damped wave equation on the torus[END_REF]) is not stated with explicit constants and we need to be sure that these constants are uniform in n (even if this is surely the case). Thus, we accept here a small loss (harmless for the proof of Theorem 10.1) and we use the explicit statement of Theorem A.3 to obtain that, for any ε > 0, there exists a constant C such that

e QnAQnt U X ≤ C t β/2-ε ∥U ∥ D(A)
which concludes the proof. □

• Step 3: the finite determining modes

In this step, we follow the arguments of [START_REF] Hale | Regularity, determining modes and Galerkin methods[END_REF] with the main modifications coming from the weaker decay of the linear semigroup. We consider the complex setting, that is that the functions in X are complex valued. We recall that a function Ψ between two complex Banach spaces Y and Z is said to be holomorphic if its Fréchet derivative exists for any y ∈ Y . We introduce the notation

B M,δ (Y ) = {U (t) ∈ C 0 (R, Y ) | ∀t ∈ R, ∥ℜ(U (t))∥ Y ≤ M and ∥ℑ(U (t)∥ Y ≤ δ} .
The space B M,δ (Y ) is naturally endowed with the L ∞ (Y )-norm. We assumed that σ 0 > 0 for d = 2 and σ 0 > 1/2 for d = 3 as in Hypothesis (iii) of Theorem 10.1, so that H 1+σ 0 (Ω) ⊂ C 0 (Ω). We will use the holomorphic extension of f in a technical setting stated in the following lemma. Except this particular setting, the result is a straightforward consequence of the analyticity of f and we omit the proof.

Lemma 10.3. Assume u ∈ R → f (x, u) ∈ R is analytic with respect to u. Denote κ the injection constant ∥ • ∥ L ∞ ≤ κ∥ • ∥ H 1+σ 0 .
Let M be given and M ′ 0 ≥ 0. Then, there exists M ′ ≥ M ′ 0 , as well as two small positive constants δ and δ ′ such that the following holds. The function

z ∈ R → f (•, z) ∈ C 1 (Ω, R) has a holomorphic extension in {z ∈ C, |ℜ(z)| ≤ κ(M + M ′ + δ) and |ℑ(z)| ≤ κ(δ + δ ′ )}.
We apply the above lemma to obtain the following result. Proposition 10.4. Let M and M ′ 0 ≥ 0 be given and let n 0 ∈ N. Let M ′ , δ and δ ′ the constants given by the previous lemma. Then, there exist n ≥ n 0 so that for any function V (t) in the complex set B M +δ,δ (P n X σ 0 ), there exists a unique bounded solution W in

B M ′ ,δ ′ (Q n X σ 0 ) of ∂ t W = (Q n AQ n )W + (Q n AP n )V + Q n F (V + W ) .
(10.8)

In addition, the mapping V → W (V ) is lipschitzian and holomorphic.

Proof: Assume that W solution of (10.8) exists and is bounded in

C 0 (R, Q n X σ 0 ). Then, W (t) =e QnAQn(t-t 0 ) W (t 0 ) + ∫ t-t 0 0 e QnAQns ((Q n AP n )V (t -s) + Q n F (V + W )(t -s)) ds.
Using Proposition 10.2, when t 0 goes to -∞, we get

W (t) = ∫ ∞ 0 e QnAQns ((Q n AP n )V (t -s) + Q n F (V + W )(t -s)) ds .
Conversely, it is easy to see that a solution of the previous integral equation is a solution of (10.8). To prove Proposition 10.4, we set up a fixed point theorem for contracting maps. We introduce the map Φ V , defined for W bounded in

C 0 (R, Q n X σ 0 ) by Φ V (W )(t) = ∫ ∞ 0 e QnAQns ((Q n AP n )V (t -s) + Q n F (V + W )(t -s)) ds .
During the first part of the proof, we consider real valued functions, so that the terms including the function f are well defined. Let ν ∈ (0, 1 -σ 0 ] to be fixed later. Let

V = (v 1 , v 2 ), we have Q n AP n V = (0, Q n γ(x)P n v 2 ). Thus, using that γ is of class C 1 , ∥Q n AP n V ∥ X ν+σ 0 ≤ C∥P n v 2 ∥ H ν+σ 0 ≤ C|λ n | ν/2 ∥P n v 2 ∥ H σ 0
where (-λ k ) denotes the eigenvalues of the Laplacian operator L. Proposition 10.2 and the bound on V show that

∫ ∞ 0 e QnAQns (Q n AP n )V (t -s) ds L ∞ (R,X σ 0 ) ≤ C 1 (ν)|λ n | ν ∥V ∥ L ∞ (R,X σ 0 ) (10.9)
as soon as there is ε > 0 such that 1/t νβ 2 -ε is integrable in a neighborhood of +∞, that is for νβ > 2.

For any function

Z(t) = (z 1 , z 2 ) bounded in X σ 0 , we have Q n F (Z) = (0, Q n f (z 1 )
). Since H 1+σ 0 (Ω) is a normed algebra, using the polynomial growth of f stated in (1.2), we get

∥Q n F (Z)∥ X ν+σ 0 = ∥Q n f (z 1 )∥ H ν+σ 0 ≤ C |λ n | 1-ν ∥Q n f (z 1 )∥ H 1+σ 0 ≤ C(1 + ∥Z∥ p X σ 0 ) |λ n | 1-ν .
Once again, Proposition 10.2 shows that

∫ ∞ 0 e QnAQns Q n F (Z)(t -s) ds L ∞ (R,X σ 0 ) ≤ C 2 (ν)(1 + ∥Z∥ p L ∞ (R,X σ 0 ) ) |λ n | 1-ν (10.10)
as soon as νβ > 2.

In the same way, using the control of f ′ u stated in (1.2), we prove that, if

νβ > 2, ∫ ∞ 0 e QnAQns Q n ( F (Z) -F (Z ′ ) ) (t -s) ds L ∞ (R,X σ 0 ) ≤ C 3 (ν) ( 1 + ∥Z∥ p-1 L ∞ (R,X σ 0 ) + ∥Z ′ ∥ p-1 L ∞ (R,X σ 0 )
)

|λ n | 1-ν ∥Z -Z ′ ∥ L ∞ (R,X σ 0 ) . (10.11)
Gathering (10.9), (10.10) and (10.11), we obtain that, for any real functions V and W with ∥V ∥ L ∞ (R,X σ 0 ) ≤ M and ∥W ∥ L ∞ (R,X σ 0 ) ≤ M ′ , Φ is well defined, bounded and locally lipschitzian. To apply the fixed point theorem for contracting maps, we need that the Lipschitz constant is smaller than 1, which is implied by

C 3 (ν) ( 1 + M p-1 + M ′ p-1
)

|λ n | 1-ν ≤ 1 2 (10.12)
and that Φ V maps B M ′ ,0 (Q n X σ 0 ) into itself, which is implied by

C 1 (ν)|λ n | ν M + C 2 (ν)(1 + M p + M ′ p )) |λ n | 1-ν ≤ M ′ 2 . ( 10.13)
To this end, we choose ν > 0 such that (p -1)ν < 1 -ν and we fix M ′ to be equal to 4C 1 (ν)|λ n | ν M , so that the first term of (10.13) is smaller but satisfies the same growth than to the bound M ′ /2. Then, since λ n goes to +∞ when n → +∞, one can find n large enough such that (10.12) and (10.13) hold, since M = o(M ′ ) and

M ′ p-1 = o(|λ n | 1-ν )
when n goes to +∞. Taking n larger if needed, the bounds n ≥ n 0 and M ′ ≥ M ′ 0 are easily fulfilled.

It remains to check that we can find ν satisfying all the required conditions. The bound (p -1)ν < 1 -ν is equivalent to ν < 1/p. It is compatible with νβ > 2 since we assumed 2 β < 1 p . Moreover, we also need that ν ∈ (0, 1 -σ 0 ], which is possible since 2 β < 1 -σ 0 . Now, we extend our functions in a complex strip. By the previous bounds, if it is real, (v 1 +w 1 )(x, t) always stays smaller than κ(M +M ′ ) where κ is the is the injection constant ∥ • ∥ L ∞ ≤ κ∥ • ∥ H 1+σ 0 . Since f is analytic, it has a holomorphic extension in a complex neighborhood of the real interval [-κ(M + M ′ ); κ(M + M ′ )]. Thus, one can also consider functions V and W with small imaginary parts in X σ 0 . All the above estimates extend by continuity in this complex strip and, since (10.12) and (10.13) contain some margin, for δ, δ ′ > 0 small enough, Φ V can be extended as a contraction map from B M ′ ,δ ′ (Q n X σ 0 ) into itself for all V ∈ B M +δ,δ (P n X σ 0 ). Proposition 10.4 then follows from the fact that Φ V has a unique fixed point W (V ), which corresponds to the unique solution of (10.8).

To conclude, we notice that the above estimates also show that Φ V (W ) is Lipschitz continuous with respect to V and thus that the fixed point W (V ) is Lipschitz continuous with respect to V . To obtain that the fixed point W (V ) depends holomorphically of V , we have to show that the map Φ is holomorphic with respect to V and W . Then, one can conclude by using the implicit function theorem. To show that Φ is holomorphic, we have to show that it has Fréchet derivatives. This can be obtained by using the fact that f is holomorphic and arguments similar to the above ones.

□

The proof of Proposition 10.4 also yields the property of finite determining modes. Proposition 10.5. Assume that Hypothesis (ii) of Theorem 10.1 holds and let M be given. Let σ 0 > 0 for d = 2 or σ 0 > 1/2 for d = 3. Then there exists n ∈ N such that the following holds. Let U 1 (t) and U 2 (t) be two global solutions of (10.1) such that ∥U i (t)∥ X σ 0 ≤ M for all t ∈ R. If P n U 1 (t) = P n U 2 (t) for all times t ∈ R, then U 1 (t) ≡ U 2 (t) for all t.

Proof: We consider the projections V i = P n U i and W i = Q n U i . For any n, we have that ∥V i (t)∥ X σ 0 ≤ M and ∥W i (t)∥ X σ 0 ≤ M for all t ∈ R. We argue as in the proof of Proposition 10.4 with M = M ′ 0 and n 0 = 0. The fixed point argument of Proposition 10.4 can be applied for n large enough. It shows that, if the low frequencies V 1 = V 2 are known, then there is only one possible high frequencies part W (t), solution of (10.8). Since (10.3) and (10.4) are equivalent, this unique function W must be equal to W 1 as well as W 2 . So W 1 = W 2 and thus U 1 = U 2 . □

• Step 4: End of the proof of Theorem 10.1.

Let U (t) be the mild solution of Theorem 10.1 and let U = V + W be the splitting of (10.3) for some n ∈ N. We have that U (t) is uniformly bounded in X σ 0 by some constant M , thus, for all n ∈ N, V (t) and W (t) are also bounded by M , independent of n. From now on, we fix n, δ, M ′ and δ ′ as prescribed by Proposition 10.4 for such M , for M ′ 0 = M and n 0 = 0. We have the existence of a lipschitzian and holomorphic map Ṽ → W ( Ṽ ) defined in a neighborhood of V . We consider the ordinary differential equation

{ ∂ s Ṽ (s) = (P n AP n ) Ṽ (s) + P n AQ n W ( Ṽ ) + P n F ( Ṽ + W ( Ṽ )) Ṽ (0)(t) = V (t) ∈ C 0 (R, P n X σ 0 ) (10.14)
defined in the Banach space C 0 (R, P n X σ 0 ) of finite dimension. Notice that (10.14) is really an ODE since P n has finite rank and thus (P n A) is a bounded operator. Proposition 10.4 shows that Ṽ → W ( Ṽ ) is lipschitzian and holomorphic in a neighborhood of the initial data V (t). It also ensures that f , and thus F , are holomorphic in the complex set where Ṽ + W ( Ṽ ) takes values. As a consequence, by the classical theory of ODE's in Banach spaces, (10.14) admits a unique solution Ṽ (s) for small s ∈ C, |s| ≤ ϵ, and this solution is holomorphic with respect to s. Notice that the construction is just made such that, for any t ∈ R and for s ∈ [-ϵ, ϵ] real, s → Ũ (s)(t) = ( Ṽ (s) + W ( Ṽ (s)))(t) is a solution of (10.1) with P n Ũ (s) = Ṽ (s). By uniqueness of the solution of (10.14), using the translation invariance, we have Ũ (0)(t + s) = Ũ (s)(t). Thus t → Ũ (0)(t) is a mild solution of (10.1) with P n ( Ũ (0)(t)) = V (t). Due to Proposition 10.5, Ũ (0)(t) = U (t) for all t ∈ R. Since, for small s ∈ [-ϵ, ϵ] real, we have U (t + s) = Ũ (0)(t + s) = Ũ (s)(t), we get that s → U (t + s) is an analytic function and thus t → U (t) ∈ X σ 0 is an analytic function.

Proof of Theorem 1.1

Due to Proposition 3.5, we only need to obtain a unique continuation property. The difference with Theorem 1.2 is that the geometric background is quite general and we cannot use Theorems 6.5 or 6.7. Our goal is thus to use the analytic unique continuation of [START_REF] Robbiano | Uniqueness in the Cauchy problem for operators with partially holomorphic coefficients[END_REF] stated in Theorem 6.1.

Let W = (w, w t ) be a globally bounded solution of (1.1) as in Propositions 3.3 and 3.4. We have that W (t) is globally bounded in X σ 0 for t ∈ R and σ 0 > 0. We want to apply Theorem 10.1. The assumption β > 2p ≥ 2 is already contained in the assumptions of Theorem 1.1. Since we work here with Ω of dimension d = 2, we may choose σ 0 as small as wanted and in particular we can fulfill the condition β > 2/(1 -σ 0 ). If γ is of class C 1 , we may directly apply Theorem 10.1 and obtain that W (t) is analytic in time. If γ is not smooth, Assumption b) of Theorem 1.1 provides a C 1 damping γ with a smaller support but with the same decay properties for the corresponding damped wave semigroup. Since the energy of W (t) is constant, W (t) vanishes in the support of the damping γ and satisfies ∂ tt w(x, t) = ∆w(x, t) -αw(x, t) -f (x, w(x, t)) . (10.15) We may thus replace γ by the regular damping γ and W (t) is still a solution of the corresponding damped (or free) wave equation. We apply Theorem 10.1 in this setting and still get that W (t) is analytic in time.

Since w is smooth with respect to t and f is smooth, (10.15) yields that (∆ -α)w is in L 2 (Ω) and thus w is in H 2 (Ω). We differentiate the above equation to obtain

(∆ -α) 2 w = (∆ -α)(∂ 2 tt w + f (x, w)) = ∂ 2 tt (∆ -α)w + (∆ -α)f (x, w) = ∂ 4 tttt w + f ′ u (x, w)∂ 2 tt w + f ′′ uu (x, w)|∂ t w| 2 + (∆ -α)f (x, w)
showing that w belongs to H 4 (Ω). The process can be used as many times as wanted, showing that w(x, t) is also smooth with respect to x. Thus, the coefficients of (3.2) are smooth in x and analytic in t and the unique continuation property of Theorem 6.1 applies. Then Theorem 1.1 is a direct consequence of Proposition 3.5.

Application 4: the disk with many holes

In Section 9, we have proved the semi-uniform stabilization for the semilinear damped wave equation in the case of the disk with two holes. In Figure 6, we have seen that if the disk has three holes or more, Theorem 1.2 does not apply. In this case we assume that u → f (x, u) is analytic and apply Theorem 1.1. Notice that we may take γ ∈ L ∞ , since the second part of Assumption b) of Theorem 1.1 is satisfied (see the geometric conditions in Theorem 8.1). Once again, we obtain an estimation of the decay which is better than the one given by Proposition 1. Then, any solution u of (1.1) satisfies

∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 --------→ t-→+∞ 0 .
Moreover, there exists λ such that, for any R and σ ∈ (0, 1], there exists C R,σ such that, for any solution u with U 0 ∈ X σ ,

∥(u 0 , u 1 )∥ H 1+σ ×H σ ≤ R =⇒ ∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 ≤ C R,σ e -σ λt 1/3 .
12 Application 5: Hyperbolic surfaces

In the case of hyperbolic surfaces, some recent results in Jin [START_REF]Long Jin Damped wave equations on compact hyperbolic surfaces[END_REF] following the fractal uncertainty principle ideas of Bourgain-Dyatlov [START_REF] Bourgain | Dyatlov Spectral gaps without the pressure condition[END_REF] give a very good decay for any non trivial damping. In our nonlinear setting, the application of our previous results give the following result. Then, any solution u of (1.1) satisfies

∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 --------→ t-→+∞ 0 .
Moreover, there exists λ such that, for any R and σ ∈ (0, 1], there exists C R,σ such that, for any solution u with U 0 ∈ X σ ,

∥(u 0 , u 1 )∥ H 1+σ ×H σ ≤ R =⇒ ∥(u, ∂ t u)(t)∥ H 1 0 ×L 2 ≤ C R,σ e -σ λ√ t .
The result follows from an application of Theorem 1.1. The decay of the semigroup can be found in [START_REF]Long Jin Damped wave equations on compact hyperbolic surfaces[END_REF]. Again, it is quite sure that we could avoid the loss and obtain a better decay e -σ λt by following the arguments inside of the proof in [START_REF]Long Jin Damped wave equations on compact hyperbolic surfaces[END_REF] for a slightly modified operator.

Note that the results in the references involve the case α = 0, but similar result can be obtained for α ≥ 0 for the linear semigroup as in Lemma 4.2. We also refer to other results with pressure conditions Schenck [START_REF] Schenck | Energy Decay for the Damped Wave Equation Under a Pressure Condition[END_REF] following ideas of Anantharaman [START_REF] Anantharaman | Entropy and the localization of eigenfunctions[END_REF]. We also want to stress that the result of [START_REF]Long Jin Damped wave equations on compact hyperbolic surfaces[END_REF] follows several deep progress in the subject for ergodic flow and with various assumptions on the damping, but it would be impossible to make a complete bibliography. We refer to the bibliography in [START_REF]Long Jin Damped wave equations on compact hyperbolic surfaces[END_REF] for instance, or the survey [START_REF] Nonennmacher | Resonances in hyperbolic dynamics[END_REF] for a history of resolvent estimates that can lead to such result of damping.

A uniform bound for the H 2 × H 1 -norm

This section is devoted to the proof of Theorem 1.4. We will assume in the whole section that the conclusions of Theorem 1.1 hold.

The following lemma is very general and does not depend on the geometric setting.

= (u 0 , u 1 ) belongs to D(A) = (H 2 (Ω) × H 1 0 (Ω)) × H 1 0 (Ω), then U (t) = (u, ∂ t u)(t)
also belongs to D(A). Moreover, there exists β > 0 such that, for all R > 0, there exists C(R) > 0 such that,

∥(u 0 , u 1 )∥ H 2 ×H 1 ≤ R =⇒ ∥(u, ∂ t u)(t)∥ H 2 ×H 1 ≤ C(R)(1 + t) β .
For d = 2, the exponent β is as close to 1 as wanted.

Proof: The proof of this result is classical. Assume that U 0 = (u 0 , u 1 ) belongs to

D(A) = (H 2 (Ω) × H 1 0 (Ω)) × H 1 0 (
Ω), we deal with the Cauchy problem by classical arguments since the linear semigroup e At is well defined on D(A), H 2 (Ω) ⊂ L ∞ (Ω) and f (x, 0) = 0 due to (1.3), preserving the Dirichlet boundary conditions. Thus, the solution U (t) is locally well defined in D(A).

We recall that the damped wave equation admits the physical energy as a Lyapounov function

E(U ) = ∫ Ω 1 2 (|∇u| 2 + α|u| 2 + |∂ t u| 2 ) + V (x, u) .
As noticed in Section 2, this energy is non-increasing in time. Using the Sobolev embeddings and Assumptions (1.2) and (1.3), we obtain that ∥U (t)∥ H 1 ×L 2 is uniformly bounded if U 0 belongs to a bounded set of H 1 0 (Ω)×L 2 (Ω) and in particular is bounded by a constant

C 1 (R) if ∥(u 0 , u 1 )∥ H 2 ×H 1 ≤ R.
We introduce a energy of higher order

F(U ) = 1 2 ∫ Ω (|∆u| 2 + α|∇u| 2 + |∇∂ t u| 2 ) - ∫ Ω f (x, u)∆u .
Notice that this energy is well defined for U (t) ∈ D(A). To be more precise, we have

∫ f (x, u)∆u ≤ 1 4 ∫ |∆u| 2 + ∫ |f (x, u)| 2 and since H 1 (Ω) ⊂ L 2p (Ω), ∫ |f (x, u)| 2 is controlled by ∥u∥ 2p H 1 and thus by C 1 (R) 2p . In particular, there exists C 2 (R) such that 1 4 ∥U (t)∥ 2 H 2 ×H 1 -C 2 (R) ≤ F(U (t) ≤ 3 4 ∥U (t)∥ 2 H 2 ×H 1 + C 2 (R) .
We have

∂ t F(U ) = ∫ ∆u∆∂ t u + α∇u∇∂ t u + ∇∂ t u.∇∂ 2 tt u -f ′ u (x, u)∂ t u∆u -f (x, u)∆∂ t u = ∫ (∆u -αu -f (x, u) -∂ tt u)∆∂ t u -f ′ u (x, u)∂ t u∆u = ∫ γ(x)∂ t u∆∂ t u -f ′ u (x, u)∂ t u∆u = ∫ -γ(x)|∇∂ t u| 2 -∂ t u∇γ∇∂ t u -f ′ u (x, u)∂ t u∆u ≤ ∫ |∂ t u||∇γ||∇∂ t u| + ∫ |f ′ u (x, u)||∂ t u||∆u| (13.1)
The first term is bounded by ∥∇γ∥

L ∞ ∥U ∥ H 1 ×L 2 ∥U ∥ H 2 ×H 1 and so by C 3 (R)∥U ∥ H 2 ×H 1 .
The second term of (13.1) is bounded by C

1 (R)∥f ′ u (•, u)∥ L ∞ ∥U ∥ H 2 ×H 1 . We bound ∥f ′ u (•, u)∥ L ∞ as follows for d = 3: ∥f ′ u (•, u)∥ L ∞ ≤ (1 + ∥u∥ L ∞ ) p-1 ≤ (1 + ∥u∥ H 3/2+ε ) p-1 ≤ (1 + ∥u∥ 1/2-ε H 1 ∥u∥ 1/2+ε H 2 ) p-1
where ε ∈ (0, 1/2] can be chosen small enough so that θ = (p -1)(1/2 + ε) < 1 since p < 3. Thus the second term of (13.1) is bounded by

C 4 (R)(1 + ∥U ∥ 1+θ H 2 ×H 1 ). We finally obtain that ∂ t F(U ) ≤ C 5 (R)(C 6 (R) + F(U )) η with η = (1 + θ)/2 < 1. This show that F(U ) ≤ C(R)(1 + t) δ with δ = 1/(1 -θ).
In the case d = 2, the bound of the second term of (13.1) is of the type (1 + ∥u∥ Remark: in the case of the open book, the decay of the semigroup is polynomial and Theorem 1.4 does not apply. To adapt the above arguments and still get integrability where it is needed, we should assume that the vanishing order β is small enough. However, this constraint seems not compatible with γ of class C 1 . It could be possible to find sharper arguments but this is not the purpose of this paper.

1-ε H 1 ∥u∥ ε H 2 ) p-1 with ε ∈ (0,
∈ (0, 1/2), F (U (t)) = (0, f (•, u(t))) is uniformly bounded in (H 2+ε (Ω) ∩ H 1 0 (Ω) × H 1+ε 0 (

A Estimates of the resolvent and decay of the semigroup

The decay rate of a linear semigroup e At is closely related to the control of the resolvent (A -iµ) with µ ∈ R, that is the resolvent along the imaginary axis. A famous result of [START_REF] Gearhart | Spectral theory for contraction semigroups on Hilbert spaces[END_REF], [START_REF] Prüss | On the spectrum of C 0 -semigroups[END_REF] and [START_REF] Huang | Characteristic conditions for exponential stability of linear dynamical systems in Hilbert spaces[END_REF] is as follows.

Theorem A.1. Gearhart-Prüss-Huang Let e At be a C 0 -semigroup in a Hilbert space X and assume that there exists a positive constant M > 0 such that ∥e At ∥ L(X) ≤ M for all t ≥ 0. Then there exist C > 0 and λ > 0 such that ∀U ∈ X , ∥e At U ∥ X ≤ Ce -λt ∥U ∥ X if and only if iR ⊂ ρ(A) and sup µ∈R ∥(A -iµId) -1 ∥ L(X) < +∞ .

In the case of the weak stabilization, the resolvent (A -iµId) -1 is no more uniformly bounded for µ ∈ R. The rate of blow-up of this resolvent when µ → ±∞ is related to the decay of e At in L(D(A), X). A general relation has been obtained by Batty and Duyckaerts in [START_REF] Charles | Non-uniform stability for bounded semigroups on Banach spaces[END_REF]. The first implication is the following. The second implication is more useful but is not optimal in general due to the logarithmic loss in M log which is not expected. log is the inverse of M log which maps (M log (0), +∞) onto (0, +∞).

If X is a Hilbert space and M polynomial, we can get rid of the logarithmic term in M log as proved in [START_REF] Borichev | Optimal polynomial decay of functions and operator semigroups[END_REF] (see also [START_REF] Anantharaman | Sharp polynomial decay rates for the damped wave equation on the torus[END_REF] in the framework of a damped wave system). Let e At be a bounded C 0 -semigroup on a Hilbert space X with generator A such that iR ∩ σ(A) = ∅. Then, for a fixed α > 0, the following conditions are equivalent: (i) for large µ ∈ R, ∥(A -iµ) -1 ∥ L(X) = O(|µ| 1/α ), (ii) for large t ≥ 0, ∥e At A -1 ∥ L(X) = O(1/t α ), (iii) for all x ∈ H and for large t ≥ 0, ∥e At A -1 x∥ X = o(1/t α ).

B Estimates of the resolvent of abstract damped wave equations

In this section, we consider an abstract damped wave equation. Let H be Hilbert spaces and let L : D(L) → H be a positive self-adjoint operator with compact resolvent. Let B ∈ L(H) be a damping operator which is bounded, self-adjoint and non-negative. We set X = D(L 1/2 ) × H and

A = ( 0 Id -L -B ) D(A) = D(L) × D(L 1/2 ) .
For µ ∈ R, we also introduce the operator P B (µ) : D(L) → H defined by

P B (µ) = -L -iµB + µ 2 Id .
In Theorems A.3 and A.4, we have seen the importance of estimating the resolvent (Aiµ) -1 . In this section, we recall the equivalence with estimating P B (µ) -1 , which is often more convenient. This type of arguments is very classical and may be found in many articles dealing with the stabilization of damped wave equations. We present them here for sake of completeness and because we will need to generalize most of them in the next appendix.

We begin with the resolvent (A -iµ) -1 for fixed µ ∈ R, that is that we consider the low frequencies. Since L has compact resolvent, the P B (µ) is invertible if and only if its kernel is reduced to {0} and if it does, P B (µ) -1 is bounded from H into D(L). Thus (B.1) yields that (i)⇔(ii). Moreover, if P B (µ)u = 0 with u ̸ = 0 then taking the scalar product with u and considering the imaginary part, we get that ⟨Bu|u⟩ = 0 and so B 1/2 u = 0. Thus, Bu = 0 and Lu = µ 2 u showing that (iii) fails. In the converse way, if we assume that (iii) fails, the corresponding solution u also solves P B (µ)u = 0 showing that (ii) fails. □

To study the high-frequencies, we have to estimate the behaviors for large µ. We have for large µ PB (µ) = P B (µ)(Id -P B (µ) -1 V ) ,

showing that PB (µ) is invertible for large µ since P B (µ) -1 goes to 0. In addition, it shows that the estimates for ∥ PB (µ) -1 ∥ and ∥P B (µ) -1 ∥ are equivalent. Then the reverse implication of Proposition B.2 finishes the proof. □

C Estimates for the high-frequencies projections

In Section 10, we need to estimate the decay of the semigroup projected into the eigenspaces corresponding to the high frequencies of the Laplacian operator. This estimation is not direct in the cases where the projectors on high-frequencies do not commute with A. The purpose of this Section is to prove results yielding quickly to estimates of the decay of the high-frequencies by using the above results Theorem A.3 and A.4. In particular, we generalized some results of Appendix B by showing that they hold uniformly with respect to cut-off frequency of the projection. Notice that this type of decay estimates for the high-frequency part of the solutions of the damped wave equation is related to Theorem 10 of [START_REF] Burq | Injections de Sobolev probabilistes et applications[END_REF], which shows that the eigenspaces corresponding to the high frequencies of the Laplacian operator are mainly preserved by the flow of the damped wave equation. We use the notations of Appendix B. Since L is self-adjoint, positive and with compact resolvent, there exists an orthonormal basis (ϕ k ) k≥0 of eigenfunctions of L. We introduce the high-frequencies truncations Q n , that are the projectors on the space Span{ϕ k , k ≥ n}

Q n u = ∑ k≥n ⟨u|ϕ k ⟩ϕ k .
We also introduce the sequence of high-frequencies projections

Q n = (Q n , Q n ) on X.
We consider in Q n H the operators

P QnBQn (µ) = Q n P B (µ)Q n = -L -iµQ n BQ n + µ 2 Id .
and the projection of A on the high frequencies: that is, for any V ∈ Q n X,

Q n AQ n = ( 0 Id -L -Q n BQ n ) .
We prove a generalization of the classical implication of Proposition B.2, which is uniform with respect to the high-frequencies projections. Then, there exists K > 0 such that, for all n ∈ N and all |µ| ≥ µ 0 ,

(Q n AQ n -iµ) -1 L(QnX)
≤ KM (|µ|) .

It remains to estimate ∥u 1 ∥ D(L 1/2 ) . To this end, we take the scalar product of second line of (C.2) with u 1 and consider the real part to obtain

∥u 1 ∥ 2 D(L 1/2 ) -µ 2 ∥u 1 ∥ 2 H ≤ O(µ)∥V ∥ X ∥u 1 ∥ H
and thus, due to the above estimates, ∥u 1 ∥ 2 D(L 1/2 ) ≤ O(M (|µ|) 2 )∥V ∥ 2 X . □

To obtain estimates as (C.1), it is convenient to generalize Proposition B.3 to the case where high-frequencies projections appear. In this way, we can use the classical observability estimate without projections to study the decay of the high-frequencies semigroup.

Proposition C.2. We set P (µ) = -L + µ 2 Id := P 0 (µ) .

Assume that there exist two positive functions f and g and µ 0 ≥ 0 such that, for any µ with |µ| ≥ µ 0 and any u ∈ D(L), We bound the sum on the right by three times the largest of the three terms. Depending on which one is the largest one, we get three different bounds, which can be gathered in (C.5) and (C.6). □

∥u∥

Theorem 1 . 1 .

 11 Consider the damped wave equation (1.1) in the framework of Assumptions (i)-(v). Assume in addition that: a) the function (x, u) → f (x, u) is smooth and analytic with respect to u.

Proposition 1 . 3 .

 13 Consider a situation where the stabilization stated in Theorems 1.1 or 1.2 holds. Then,

Proposition 3 . 5 .

 35 Assume that z ≡ 0 is the only global solution of (3.2). Then the decay assumptions (1.4) and (1.5) imply the conclusions of Theorem 1.1.

Theorem 6 . 1 .

 61 Robbiano-Zuily, Hörmander (1998) Let T > 0 (or T = +∞) and let b, c and d be smooth coefficients. Assume moreover that b, c and d are analytic in time and that z is a strong solution of

Theorem 6 . 3 .

 63 Lerner and Robbiano (1985), Hörmander.

Theorem 6 .

 6 3 states a local unique continuation property through pseudo-convex surfaces. To use it, it is more convenient to have a global version. This kind of global foliation has already been introduced in [40] by Stefanov and Uhlmann. Definition 6.4. A family of surfaces (Σ λ ) λ∈[0,1) is an oriented pseudo-convex foliation without boundary in a compact manifold Ω if:

Figure 2 :

 2 Figure2: A oriented pseudo-convex foliation without boundary in the sphere S 2 . The surfaces Σ λ forms a smooth family of vertical circles inside a hemisphere and Σ + λ and Σ - λ are respectively the large and the small spherical caps. The geodesics beings the great circles, the ones which are tangent to a surface Σ λ stay inside Σ + λ .

Figure 3 :

 3 Figure 3: A oriented pseudo-convex foliation with boundary in the disk. The surfaces Σ λ forms a smooth family of curves inside a semidisk. The surfaces Σ - λ are strictly convex and the angle formed by Σ λ and the boundary of the disk is less than π/2 on Σ - λ side. The geodesics are straight lines bouncing at the boundary according to Newton's laws. The ones which are tangent to a surface Σ λ stay inside Σ + λ .

Theorem 7 . 1 .

 71 Christianson, Schenck, Vasy & Wunsch, 2014.

Theorem 7 . 2 .

 72 Consider the damped wave equation (1.1) in the framework of the peanut of rotation introduced above. Let α > 0 and f satisfying (1.2) and (1.3). Then, any solution u of (1.1) satisfies

Figure 4 :

 4 Figure 4: An oriented pseudo-convex foliation without boundary of the central part of the peanut. Each surface Σ λ consists in two vertical circles, the set Σ - λ being the interior part surrounded by these circles. Since the central part of the peanut is negatively curved, the geodesics tangent to the vertical circles stay in the exterior domain Σ +λ . Notice that the set Σ 0 consists in both exterior disk and the disks Σ λ get closer to x = 0 when λ get closer to 1. Thus, the central circle x = 0 is not included in ∪ λ∈[0,1) Σ + λ but this is not important since this set is of zero measure.

Theorem 8 . 1 .

 81 Let O ⊂ R d be a smooth bounded open set. For i = 1 . . . p, let O i ⊂ O be smooth strictly convex obstacles satisfying: (a) the obstacles are disjoint:

Figure 5 :

 5 Figure 5: A domain satisfying Hypotheses (a)-(d) of Theorem 8.1. Notice the presence of periodic geodesics bouncing on the obstacles and never meeting the support of the damping γ. The geometric control condition thus fails and we cannot hope a uniform decay of the semigroup in X.

Figure 6 :

 6 Figure 6: Two examples of disks with holes and associated attempts to draw a suitable pseudo-convex foliation covering the whole domain. Left, the disk with two holes may be covered by a pseudo-convex foliation starting in a neighborhood of the exterior boundary, except for a single line, which is of measure zero. Right, an attempt to cover a disk with three holes with a pseudo-convex foliation. We easily notice that the central zone cannot be covered and Theorem 6.7 cannot be applied in this case.

Theorem 9 . 1 .

 91 Consider the damped wave equation(1.1) in the framework of a disk with two convex holes and assume that the damping γ is strictly positive in a neighborhood of the exterior boundary. Let α > 0 and f satisfying (1.2) and (1.3). Then, any solution u of (1.1) satisfies

Theorem 12 . 1 .

 121 Let M be a compact connected hyperbolic surface with constant negative curvature -1. Assume that (a) the damping γ ∈ L ∞ (Ω, R + ) is non zero and α > 0, (b) the nonlinearity f is smooth, satisfies (1.2) and (1.3) and u → f (x, u) is analytic.

0 e

 0 Ω) by a constant C(R). We use a last time the formula of variation of the constantU (t) = e At U 0 + ∫ t At F (U (t -s))dsand the weak decay estimates ∥e At ∥ L(H 2+ε (Ω)∩H 1 0 (Ω)×H 1+ε 0 (Ω),D(A)) ≤ h(t) with h(t) integrable on [0, +∞) to obtain that ∥U (t)∥ D(A) is uniformly bounded. □

Theorem A. 2 .

 2 Batty-Duyckaerts (2008) [5, Proposition 1.3]. Let e At be a semigroup of operators on a space X and assume thatm(t) = sup s≥t ∥e At (A -1) -1 ∥ L(X)goes to 0 as t → +∞. Then iR belongs to the resolvent set of A and there exist µ 0 andC > 0 such that ∀µ ∈ R with |µ| ≥ µ 0 , ∥(A -iµ) -1 ∥ L(X) ≤ 1 + Cm -1 r ( 1 2(|µ| + 1) ) where m -1r is a right inverse of m, which maps (0, m(0)] onto [0, +∞).

Theorem A. 3 .

 3 Batty-Duyckaerts (2008) [5, Theorem 1.5].Let e At be a semigroup of operators on a space X such that ∥e At ∥ L(X) ≤ C for all t ≥ 0, and such that iR ∩ σ(A) = ∅. We setM (µ) = sup |τ |≤µ ∥(A -iµ) -1 ∥ L(X) and M log (µ) = M (µ)[ln(1 + M (µ)) + ln(1 + |µ|)] .Then, for any k ∈ N \ {0}, there exist C k and T k , depending only on C, k and M , such that ∀t ≥ T k , e At (A -1) -k ≤ C k (

Proposition B. 1 .

 1 Let µ ∈ R, the three following propositions are equivalent(i) (A -iµId) is invertible in L(X), (ii) P B (µ) is invertible in L(H),(iii) for any u ̸ = 0 solution of Lu = µ 2 u, we have ⟨Bu|u⟩ ̸ = 0.Proof: Let U = (u 1 , u 2 ) and V = (v 1 , v 2 ) be vectors of X = D(L 1/2 ) × H such that (A -iµ)U = V . We have equivalently { u 2 -iµu 1 = v 1 in D(L 1/2 ) P B (µ)u 1 = v 2 + Bv 1 + iµv 1 in H . (B.1)

Proposition B. 2 .

 2 With the above notations, both estimations are equivalent(i) for large µ ∈ R, ∥(A -iµ) -1 ∥ L(X) = O(M (|µ|)), (ii) for large µ ∈ R, ∥(P B (µ)) -1 ∥ L(H) = O( M (|µ|)|µ| ).

Proposition C. 1 .

 1 Assume that there exist a function M (µ), uniformly positive, and µ 0 ≥ 0 such that, for any µ with |µ| ≥ µ 0 , n ∈ N and u ∈ Q n D(L), we have ∥u∥ H ≤ M (|µ|) |µ| ∥P QnBQn (µ)u∥ H . (C.1)

10] Proposition 8.2. Ikawa's black box (Section 6.2 of [10]) Let

O i be obstacles in R d satisfying the Assumptions (a),(c) and (d) of Theorem 8.1. Let R bb (µ) be the outgoing resolvent of the Laplacian operator outside the obstacles, that is the meromorphic continuation of

  , this concludes the proof.Then, we apply Theorem A.3 of Batty and Duyckaerts stated in Appendix to obtain the decay with rate e -λt 1/3 . Notice that Proposition B.3 and Theorem A.3 contain some losses transforming the rate ln |µ| of Lemma 8.3 into first ln 2 |µ| and then ln 3 |µ|. This is responsible of the power 1/3 in the decay rate. □

□ Proof of Theorem 8.1: Applying Propositions B.3 and B.2 in Appendix, the observability estimate of Lemma 8.3 implies that there exists C > 0 such that ∀µ ∈ R , ∥(A -iµ)∥ L(X) ≤ C ln 2 (2 + |µ|) .

  3 because we follow the idea of the remark below Lemma 4.2. The details are left to the reader.

	(b) no obstacle is in the convex hull of two others, that is that convhull(O i ∪ O j ) ∩ O k = ∅
	for i, j, k different,
	(c) if κ denotes the infimum of the curvatures of the boundaries ∂O i of the obstacles and
	L the minimal distance between two obstacles, and assume that κL > p.

Theorem 11.1. Let O ⊂ R 2 be a smooth convex bounded open set. For i = 1 . . . p, let O i ⊂ O be smooth strictly convex obstacles satisfying: (a) the obstacles are disjoint:

O i ∩ O j = ∅ for i ̸ = j, Let Ω = O \ (∪ i O i ) be

the convex domain with holes and let α > 0. Assume moreover that (d) the damping γ ∈ L ∞ (Ω, R + ) is strictly positive in a neighborhood of the exterior boundary ∂O, (e) the nonlinearity f is smooth, satisfies (1.2) and (1.3) and u → f (x, u) is analytic.

  Lemma 13.1. Let u(t) be a solution of the damped wave equation (1.1) with γ of class C 1 and f of class C 2 satisfying (1.2) and (1.3). Assume moreover that Ω is of dimension d = 2 or of dimension d = 3 and in this last case assume that p < 3. Then, if U 0

  1] as small as needed. Thus the growth of F(U ) is of type (1 + t) δ with δ as close to 2 as wanted. Since F(U ) is equivalent to ∥U ∥ 2 H 2 ×H 1 , we obtain the polynomial growth of Lemma 13.1 with β as close to 1 as wanted.

□

Proof of Theorem 1.4: By the previous lemma, we know that the H 2 × H 1 -norm of U (t) has a at most polynomial growth. By assumption, we know that the H 1 × L 2 -norm of U (t) goes to zero faster than any polynomial decay. By interpolation, this shows that the norm ∥U (t)∥ H 1+ε ×H ε is bounded for ε ∈ (0, 1). Since H 1+ε (Ω) is an algebra, we have that f (•, u) is uniformly bounded in H 1+ε (Ω). Thus, for any initial data satisfying ∥U 0 ∥ H 2 ×H 1 ≤ R, and for any ε

  Then, for any µ with |µ| ≥ µ 0 , any n ∈ N and any u ∈ Q n D(L),In particular, the imaginary part of ⟨P QnBQn (µ)u|u⟩ is µ⟨BQ n u|Q n u⟩ and∀u ∈ Q n H , ⟨BQ n u|Q n u⟩ = ∥ ∥P (µ)u∥ H = ∥P QnBQn (µ)u+iµQ n BQ n u∥ H ≤ ∥P QnBQn (µ)u∥ H +|µ|∥Q n BQ n u∥ H . (C.8) Then, we compute, for any u ∈ Q n H ∥Q n BQ n u∥ 2H ≤ ∥BQ n u∥ 2 ⟨BQ n u|Q n u⟩ and using (C.7), we obtain ∥Q n BQ n u∥2

			H ≤	f (µ) µ	∥P (µ)u∥ H + g(µ)∥ √ Bu∥ H .	(C.4)
			∥u∥ H ≤	M (µ) |µ|	∥P QnBQn (µ)u∥ H ,	(C.5)
	where	M (µ) = 3 max	(	f (µ) , f (µ) 2 ∥	√ B∥ 2 L(H) , g(µ) 2	)	.	(C.6)
									√ Bu∥ 2 H ≤	1 |µ|	∥P QnBQn (µ)u∥ H ∥u∥ H .	(C.7)
	We write
				H ≤ ∥ √ L(H) H ≤ ≤ ∥ √ B∥ 2 ∥ √ B∥ 2 L(H) |µ| ∥P QnBQn (µ)u∥ H ∥u∥ H . B∥ 2 L(H) ∥ √ BQ n u∥ 2 H
	Combining this result with (C.4), (C.7) and (C.8), we get
	∥u∥ H ≤	f (µ) |µ|	∥P QnBQn (µ)u∥ H +	f (µ) √ |µ|	∥ √	B∥ L(H) ∥P QnBQn (µ)u∥ 1/2 H ∥u∥	1/2 H
			+	g(µ) √ |µ|	∥P QnBQn (µ)u∥ 1/2 H ∥u∥	1/2

Proof: Let u ∈ Q n D(L). We have ⟨P QnBQn (µ)u|u⟩ = -⟨Lu|u⟩ + µ 2 ⟨u|u⟩ -iµ⟨BQ n u|Q n u⟩ . H .

Proof: The proof of (ii)⇒(i) is detailed in Proposition C.1 below, adding projections on the high-frequencies. The implication stated here is simply the complete case n = 0.

Let us show the converse implication. Take v 1 = 0 in (B.1). We have u 2 = iµu 1 and P B (µ)u 1 = v 2 . If (i) holds, that is ∥(A -iµ) -1 ∥ L(X) ≤ M (|µ|), we must have in particular that

In some cases, to obtain an estimation of ∥(P B (µ)) -1 ∥ L(H) , it is more convenient to prove an observability estimate and to use the following proposition. Notice that this proposition yields a loss due to the term f (µ) 2 in M (µ). This loss may sometimes be avoided but this may require an accurate study, based on particular dynamical properties of the geodesic flow.

Proposition B.3. We set

Assume that there exist two positive functions f and g and µ 0 ≥ 0 such that, for any µ with |µ| ≥ µ 0 and any u ∈ D(L),

Then, for any µ with |µ| ≥ µ 0 and any u ∈ D(L),

where

4)

Proof: In fact, this proposition is simply Proposition C.2 below in the particular case n = 0 that is Q n = Id. We choose to copy this particular case in this section for clarity. □

To finish, let us study the case where L is replaced by L = L + V where V is a bounded non-negative operator, typically a potential or a linearized term. Using the previous propositions, we show that, if M (µ) = o(µ), then the estimates for L are equivalent to the estimates for L. Using Theorems A.3 or A.4, we may obtain a relation between the decays of the semigroups.

Proposition B.4. We use the above notations and set

) also holds for large µ.

We set w = (P QnBQn (µ)) -1 (v 1 ). Since P QnBQn (µ) + L + iµQ n BQ n = µ 2 Id, we have

and so

. We have

where the above estimation are independent of n. The estimate ∥|(P QnBQn (µ)) -1 ∥| L(H,H) is given by Hypothesis (C.1), yielding

Using that M is uniformly positive, M (µ)/µ = o(M (µ) 2 ) and so

Since P QnBQn (µ) defined from D(L -1/2 ) in H is the adjoint of P QnBQn (-µ) defined from H in D(L 1/2 ), we also have

Coming back to (C.