Camille Laurent
email: laurent@ann.jussieu.fr

Matthieu Léautaud
email: matthieu.leautaud@math.u-psud.fr

On uniform observability of gradient flows in the vanishing viscosity limit

Keywords: Transport equation, gradient flow, vanishing viscosity limit, parabolic equation, minimal control time, semiclassical Schrödinger operator. 2010 Mathematics Subject Classification: 93B07, 93B05, 35B25, 35F05, 35K05, 93C73

We consider a transport equation by a gradient vector field with a small viscous perturbation -ε∆g. We study uniform observability (resp. controllability) properties in the (singular) vanishing viscosity limit ε → 0 + , that is, the possibility of having a uniformly bounded observation constant (resp. control cost). We prove with a series of examples that in general, the minimal time for uniform observability may be much larger than the minimal time needed for the observability of the limit equation ε = 0. We also prove that the two minimal times coincide for positive solutions. The proofs rely on a semiclassical reformulation of the problem together with (a) Agmon estimates concerning the decay of eigenfunctions in the classically forbidden region [HS84] (b) fine estimates of the kernel of the semiclassical heat equation [LY86].

Introduction and main results

Introduction

Given a smooth connected compact manifold M without boundary (the case of a bounded domain of R n is also discussed in Section 1.3 below), a smooth real valued vector field X on M and a real valued potential q(x), we consider the question of observability/detectability for the autonomous transport equation (∂ t -Xq)u = 0, in R × M, u| t=0 = u 0 , on M, (1.1) from an observation (open) set ω ⊂ M through the time interval (0, T ). More precisely, the question is whether there exists a constant C 0 = C 0 (T, ω) > 0 such that

C 2 0 T 0 ω |u(t, x)| 2 ds(x)dt ≥ u(T ) 2 L 2 (M) ,
for all u 0 ∈ L 2 (M) and u solution of (1.1). (1.2)

Here, ds(x) denotes any positive density measure1 on M, and the L 2 norm is defined accordingly. The observability question (1.2) is naturally solved by introducing an appropriate Geometric Control Condition (recall ∂M = ∅): we say that (M, X, ω, T ) satisfies (GCC) if for all x ∈ M, there is t ∈ (0, T ) such that φ -t (x) ∈ ω, where (φ t ) t∈R denotes the flow of X (see Section 2.3 for precise statements and proofs). We also say that (M, X, ω) satisfies (GCC) if (M, X, ω, T ) does for some T > 0; and if so, we denote by T GCC (M, X, ω) the infimum of times for which (M, X, ω, T ) satisfies (GCC).

On the other hand, endowing M with a Riemannian metric g, one may want to investigate the observability question for the viscously damped transport equation:

(∂ t -X -q -ε∆ g )u = 0, in R + * × M, u| t=0 = u 0 , on M, (1.3)
and its link with the minimal observation time T GCC (M, X, ω) associated to the limit problem (1.1) is left open. In particular, the formulation of the results in [START_REF] Guerrero | Singular optimal control for a transport-diffusion equation[END_REF] (see e.g. Theorem 2 and the discussion thereafter in that reference) suggests the possible existence of a universal constant K ≥ 1 such that T unif (ω) ≤ K T GCC (M, X, ω).

(1.5)

The present article investigates this question in a very particular case, namely assuming the vector field X is a gradient vector field, i.e. X = ∇ g f for a function f ∈ W 2,∞ (M; R) (note that the gradient is taken with respect to the Riemannian metric g). Hence, Equation (1.3) becomes

(∂ t -∇ g f • ∇ g -q -ε∆ g )u = 0, in R + * × M, u| t=0 = u 0 , on M, (1.6)
Here, given two vector fields Y 1 and Y 2 , we have denoted

Y 1 • Y 2 = g(Y 1 , Y 2 ) or (Y 1 • Y 2 )(x) = Y 1 (x) • Y 2 (x) = g x (Y 1 (x), Y 2 (x)
) for all x ∈ M. We denote similarly |Y | g = √ Y • Y the associated Riemannian norm of a vector (or a vector field). Note that the vector field ∇ g f is canonically identified with the derivation ∇ g f • ∇ g .

In this context, the first consequence of our main results can be (loosely) stated as follows.

Theorem 1.2. 1. There are geometries (M, g) such that for all Λ > 0, one can find f ∈ C ∞ (M) and ω open such that (M, ∇ g f, ω) satisfies (GCC) and T unif (ω) ≥ Λ T GCC (M, ∇ g f, ω).

2. There are (M, f, X, ω) such that for all Λ > 0, one can find a metric g Λ on M such that

• X = ∇ gΛ f,
• (M, X, ω) satisfies (GCC),

• T unif (ω) ≥ Λ T GCC (M, X, ω).

In particular, Theorem 1.2 states that there is no K such that (1.5) holds for all (M, X, ω). The second item in Theorem 1.2 stresses the importance of the viscosity one chooses. Namely, with the same vector field X, changing the metric g, that is the viscous perturbation, may change the minimal uniform observability time. We also obtain related results for domains of R n (see Section 1.3.1).

Our second main result in this setting concerns the uniform observability of positive solutions to (1.6). Recall that nonnegative data u 0 ≥ 0 give rise to positive solutions to (1.6). We define C + 0 (T, ε) the observability constant for positive solutions, that is for which (1.2) holds for all u 0 ≥ 0, and accordingly set T + unif (ω) = inf {T > 0 for which there exist K 0 , ε 0 > 0 such that (1.4) holds for all u 0 ≥ 0, with C + 0 (T, ε) ≤ K 0 for all ε ∈ (0, ε 0 ) . (1.7)

Theorem 1.3 (Positive solutions). For all f ∈ C 3 (M; R), and ω ⊂ M such that (M, ∇ g f, ω) satisfies (GCC), we have T + unif (ω) = T GCC (M, ∇ g f, ω). As usual, these uniform observability/non-observability results can be reformulated in terms of uniform controllability/non-controllability statements for an adjoint controlled equation, see Section 2.1.1 below.

Background and motivation

Uniform controllability problems for singular perturbations of partial differential equations already appeared in the reference book of Lions [START_REF] Lions | Contrôlabilité exacte, perturbations et stabilisation de systèmes distribués[END_REF]Chapter 3]. In the context of transport/heat equation in vanishing viscosity limit, this study was initiated by Coron and Guerrero on the 1D problem with constant speed in [START_REF] Coron | Singular optimal control: A linear 1-D parabolic-hyperbolic example[END_REF], where the authors make a conjecture on the minimal time needed to achieve uniform controllability. Then, the estimates on this minimal time have been improved successively in [START_REF] Glass | A complex-analytic approach to the problem of uniform controllability of a transport equation in the vanishing viscosity limit[END_REF][START_REF] Lissy | A link between the cost of fast controls for the 1-D heat equation and the uniform controllability of a 1-D transport-diffusion equation[END_REF][START_REF] Lissy | An application of a conjecture due to Ervedoza and Zuazua concerning the observability of the heat equation in small time to a conjecture due to Coron and Guerrero concerning the uniform controllability of a convection-diffusion equation in the vanishing viscosity limit[END_REF][START_REF] Lissy | Explicit lower bounds for the cost of fast controls for some 1-D parabolic or dispersive equations, and a new lower bound concerning the uniform controllability of the 1-D transport-diffusion equation[END_REF]. We also refer to the articles [START_REF] Münch | Numerical estimations of the cost of boundary controls for the equation ytεyxx + M yx = 0 with respect to ε[END_REF][START_REF] Amirat | On the controllability of an advection-diffusion equation with respect to the diffusion parameter: asymptotic analysis and numerical simulations[END_REF][START_REF] Amirat | Asymptotic analysis of an advection-diffusion equation and application to boundary controllability[END_REF] proposing numerical experiments to find the optimal minimal time. Such uniform control properties in singular limits are also addressed for vanishing dispersion in [START_REF] Glass | Some exact controllability results for the linear KdV equation and uniform controllability in the zero-dispersion limit[END_REF] and for vanishing dispersion and viscosity in [START_REF] Glass | Uniform controllability of a transport equation in zero diffusiondispersion limit[END_REF].

Whereas the one dimensional problem with a constant vector field has received a lot of attention in the past fifteen years, there are very few results in higher dimension or for non-constant vector fields. Besides [START_REF] Guerrero | Singular optimal control for a transport-diffusion equation[END_REF] we are only aware of the results of [START_REF] Asier | Uniform controllability of a Stokes problem with a transport term in the zero-diffusion limit[END_REF] and [START_REF] Asier | Cost of null controllability for parabolic equations with vanishing diffusivity and a transport term[END_REF] for the flat Laplace operator and the vector field ∂ x1 , with several boundary conditions.

Note that controllability problems for nonlinear conservation laws in vanishing viscosity have also been studied in [START_REF] Glass | On the uniform controllability of the Burgers equation[END_REF], [START_REF] Léautaud | Uniform controllability of scalar conservation laws in the vanishing viscosity limit[END_REF], and [START_REF] Marbach | Small time global null controllability for a viscous Burgers' equation despite the presence of a boundary layer[END_REF].

There are several motivations for studying the vanishing viscosity limit. A first motivation comes from the theory of conservation laws, for which the vanishing viscosity criterium is a selection principle for the physical (called entropy) solution, see [START_REF] Stanislav | First order quasilinear equations with several independent variables[END_REF] or [START_REF] Dafermos | Hyperbolic conservation laws in continuum physics[END_REF]Chapter 6]. It is therefore very natural, when considering control problems for conservation laws, to study the cost of the viscosity, that is, to determine if known controllability properties for the hyperbolic equation are still valid for the model with small viscosity, and how the size of the control evolves as the viscosity approaches 0. So far the only known results in this directions seem to be [START_REF] Glass | On the uniform controllability of the Burgers equation[END_REF] and [START_REF] Léautaud | Uniform controllability of scalar conservation laws in the vanishing viscosity limit[END_REF].

Another important motivation for studying singular limits in control problems is the seek of controllability properties for the perturbated system itself. This is well-illustrated by the papers [START_REF] Coron | On the controllability of the 2-D incompressible Navier-Stokes equations with the Navier slip boundary conditions[END_REF][START_REF] Coron | Global exact controllability of the 2-D Navier-Stokes equations on a manifold without boundary[END_REF][START_REF] Chapouly | On the global null controllability of a Navier-Stokes system with Navier slip boundary conditions[END_REF][START_REF] Coron | Small-time global exact controllability of the Navier-Stokes equation with Navier slip-with-friction boundary conditions[END_REF], where the authors investigate the Navier-Stokes system with Navier slip or slip-withfriction boundary conditions. They use a global controllability result for the inviscid equation (in this case, the Euler equation) to deduce global approximate controllability of the Navier-Stokes system.

On the other hand, the study of gradient fields naturally arises as the simplest dynamical situation among all vector fields. The importance of gradient vector fields with a vanishing viscosity coefficient also appears in theoretical physics and differential topology, through the Witten-Helffer-Sjöstrand theory [START_REF] Witten | Supersymmetry and Morse theory[END_REF][START_REF] Helffer | Puits multiples en mécanique semi-classique. IV. Étude du complexe de Witten[END_REF]. See e.g. the monographs [START_REF] Helffer | Semi-classical analysis for the Schrödinger operator and applications[END_REF][START_REF] Hans | Schrödinger operators with application to quantum mechanics and global geometry[END_REF]. In that theory, the operator -∇ g f•∇ g -ε∆ g (and its analogues acting on forms) is conjugated to a particular semiclassical Schrödinger operator, namely

P ε = -ε 2 ∆ g + |∇ g f| 2 g 4 + ε 2 ∆ g f, (1.8)
sometimes called the Witten Laplacian. Topological properties of the couple (M, f) (e.g. the Morse inequalities, linking the number of critical points of the Morse function f with the Betti numbers of M) are deduced from spectral properties of the Witten Laplacian. We also refer to [START_REF] Viet | Pollicott-Ruelle spectrum and Witten Laplacians[END_REF] for the understanding of other links between the spectral theories of the Witten Laplacian and the vector field ∇ g f • ∇ g (in appropriate spaces linked to the dynamics of the gradient flow), in the semiclassical limit ε → 0 + .

Viscous perturbations of gradient dynamics also arise naturally in molecular dynamics. Indeed, in R n , the operator -∇f • ∇ -ε∆ is the infinitesimal generator of the so-called overdamped Langevin process

dX t = ∇f(X t )dt + √ 2εdB t ,
where X t ∈ R n and (B t ) t≥0 is a standard Brownian motion of dimension n. This stochastic process is a classical model in statistical physics. It is used in particular for the simulation of molecular dynamics at low temperature (proportional to ε), see [START_REF] Chandresekhar | Stochastic problems in physics and astronomy[END_REF][START_REF] Schuss | The exit problem: a new approach to diffusion across potential barriers[END_REF]. The possible convergence to equilibrium, as well as the so-called metastability phenomenon are closely related to the low-lying eigenvalues (and associated eigenfunctions) of -∇f•∇-ε∆, or equivalently of the Witten Laplacian (1.8). For a very precise asymptotic study of the exponentially small eigenvalues of this operator, we refer e.g. to [START_REF] Helffer | Quantitative analysis of metastability in reversible diffusion processes via a Witten complex approach[END_REF][START_REF] Michel | About small eigenvalues of the Witten Laplacian[END_REF] in the case of a compact manifold and [HN06, LP10, DLLN19] in the case with boundary (see also the references therein). We finally remark that the above-mentioned works concerning the Witten Laplacian mostly study the behavior of the bottom of the spectrum P ε (thus linked to critical points of f). In the present work, we rely on a similar conjugation.

Main results

As already seen in the end of Section 1.1, the results of this article go in two different directions.

In a first part (Section 3), we prove some general lower bounds on the time T unif (ω) for a general class of domains and vector fields. This implies in particular that the quite natural idea that T unif (ω) is T GCC (M, X, ω) or even KT GCC (M, X, ω) for a universal constant K is false in general. This might be interpreted by the fact that in the vanishing viscosity limit, some strong oscillations can be responsible for concentration phenomena. The latter are not only described by the flow of X, but other parameters where an Agmon distance plays a crucial role. We also study (in Section 4) the particular case of surfaces of revolutions where we obtain refined lower bounds. This analysis also shows that the global geometry of the Riemannian manifold we consider has an effect on the vanishing viscosity limit. In particular, with a fixed vector field, we show that the choice of the Laplacian ∆ g can change drastically the time T unif (ω) of uniform controllability. This shows definitely that the flow of the vector field is not the unique parameter defining T unif (ω).

In a second part (Section 5), our results go exactly in the opposite direction, but for a specific class of solutions, namely positive solutions. As announced in Theorem 1.3, we prove that T unif (ω) + = T GCC (M, ∇ g f, ω) for positive solutions. This shows that the dynamics of positive solutions are actually well represented by the sole flow of ∇ g f.

In both cases, using the change of unknown v = e f 2ε u, see Section 2.2, the problem is reduced (modulo lower order terms, and in weighted spaces) to observability of solutions of a semiclassical heat equation

ε∂ t v -ε 2 ∆ g v + V v = 0, where V = |∇g f| 2 g 4 .
Note that most of the results we obtain are of interest for this particular question as well.

A general lower bound

The first family of results in this paper concern the general setting of a compact connected Riemannian manifold (M, g), with or without boundary ∂M, and the associated internal/boundary observability question. Namely, we consider the parabolic-transport problem with small viscosity ε > 0 and Dirichlet boundary conditions:

     (∂ t -X -q -ε∆ g )u = 0, in (0, T ) × Int(M), u = 0 on (0, T ) × ∂M, u| t=0 = u 0 , in M.
(1.9)

Moreover, we assume that the vector field X is a gradient vector field for the metric g, that is: there is

f ∈ W 2,∞ (M) (at least) such that X = ∇ g f • ∇ g .
For the Dirichlet problem (1.9) as well as for the case ∂M = ∅ discussed in Section 1.1, one may discuss the behavior of the observability constant (and in particular its possible uniform boundedness in the limit ε → 0 + ) in the internal observability inequality (1.4). Also, a boundary observability problem for (1.9) can be formulated as follows (see [START_REF] Guerrero | Singular optimal control for a transport-diffusion equation[END_REF] and Section 2.1.2 below). Given θ ∈ C ∞ (∂M), there exist a constant C 0 (T, ε) > 0 such that

C 0 (T, ε) 2 T 0 θε∂ ν u| ∂M (t) 2 H 1/2 (∂M) dt ≥ u(T ) 2 H 1 0 (M) ,
for all u 0 ∈ H 1 0 (M ) and u solution of (1.9). (1.10) Here ∂ ν denotes a unit normal (for the metric g) vector field to ∂M. Then, the uniform observability question is whether C 0 (T, ε) remains uniformly bounded in the limit ε → 0 + , and the associated minimal uniform observation time is defined again by T unif (θ) = inf {T > 0 for which there exist K 0 , ε 0 > 0 such that (1.10) holds for all u 0 ∈ H 1 0 (M), with C 0 (T, ε) ≤ K 0 for all ε ∈ (0, ε 0 ) . Before going further, let us first give the definition of an analogue of the condition (GCC) in case of a manifold with boundary ∂M = ∅ (called here Flushing Condition), as used in the Guerrero-Lebeau result [START_REF] Guerrero | Singular optimal control for a transport-diffusion equation[END_REF]. For this, we need to extend (M, g) in a slightly bigger Riemannian manifold ( M, g), i.e. such that M ⊂ Int( M) and g any Riemannian metric on M such that g| M = g. In the case of a bounded domain of R n , one may choose M = R n . We also extend f ∈ W 2,∞ (M) as a compactly supported function f ∈ W 2,∞ ( M) such that f| M = f. Definition 1.4. For x ∈ M, denote by γ x the maximal solution to

γx (t) = -∇ gf (γ x (t)), γ x (0) = x.
Note that this solution is defined globally in time since f is compactly supported.

We say that (M,

∇ g f, T ) (resp. (M, ∇ g f, ω, T )) satisfies the Flushing Condition (FC) if for all x ∈ M there is t ∈ (0, T ) such that γ x (t) / ∈ M (resp. γ x (t) / ∈ M or γ x (t) ∈ ω). We also say that (M, ∇ g f) (resp. (M, ∇ g f, ω)) satisfies (FC) if there is T > 0 such that (M, ∇ g f, T ) (resp. (M, ∇ g f, ω, T )) does. We then define accordingly the time T F C (M, ∇ g f) (resp. T F C (M, ω, ∇ g f))
as the infimum of T > 0 for which this property holds.

Remark finally that these definitions do not depend on the extensions ( M, g) and f.

Guerrero-Lebeau [START_REF] Guerrero | Singular optimal control for a transport-diffusion equation[END_REF] prove an analogue of Theorem 1.1 in the present setting (and for general vector fields), namely: if (M, ∇ g f) (resp. (M, ∇ g f, ω)) satisfies the Flushing Condition, given θ any nontrivial observation function, then there is T u > 0 and K 0 > 0 such that for all T ≥ T u , (1.10) holds (resp. (1.4) holds for all solutions to (1.9)), with C 0 (T, ε) ≤ K 0 for all ε ≤ 1.

Two important geometric quantities in our results are the potential associated to the function f, defined by

V (x) := |∇ g f(x)| 2 g 4 , (1.11)
and the Agmon distance (see e.g. [Hel88, Chapter 3]) to the minimum of this potential, namely, with

E 0 = min M V , d A (x, y) = inf 1 0 (V (γ(t)) -E 0 ) + | γ(t)| g dt, γ ∈ U 1 (x, y) , U 1 (x, y) = γ ∈ W 1,∞ ([0, 1]; M), γ(0) = x, γ(1) = y , d A (x) = d A (x, V -1 (E 0 )) = inf y∈V -1 (E0) d A (x, y).
(1.12)

Here V -1 (E 0 ) is the classically allowed region at the potential minimum,

(V (x) -E 0 ) + = max (V (x) -E 0 , 0) ,
and d A (x)
is the Agmon distance of x to the set V -1 (E 0 ) for the (pseudo-)metric (V -E 0 ) + g. Remark that the index (•) + is not needed at the bottom energy E 0 = min M V ; however, we keep it here since the definition (1.12) will also be useful for a general energy level.

Our main result in this general setting formulates as follows.

Theorem 1.5. We assume that

f ∈ W 2,∞ (M) (hence V ∈ W 1,∞ (M)) and q ∈ L ∞ (M). We let E 0 = min M V , set W E0 (x) = f(x) 2 + d A (x),
and fix ω ⊂ M (resp. θ ∈ C ∞ (∂M)), and in the case of boundary estimates, we further assume f, q ∈ C ∞ (M). For any δ > 0, there is ε 0 > 0 such that for all ε ∈ (0, ε 0 ) the observability inequality (1.4)

(resp. (1.10)) with constant C 0 (T, ε) implies C 0 (T, ε) ≥ exp 1 ε min ω W E0 -max V -1 (E0) f 2 -E 0 T -δ , resp. C 0 (T, ε) ≥ exp 1 ε min supp θ W E0 -max V -1 (E0) f 2 -E 0 T -δ .
In particular, we have

E 0 T unif (ω) ≥ min ω W E0 -max V -1 (E0) f 2 , (1.13) (resp. E 0 T unif (θ) ≥ min supp θ W E0 -max V -1 (E0) f 2 ).
(1.14)

Note that the quantity in the right hand-side of (1.13)-(1.14) as well as E 0 are invariant under the change f → f+C for C constant on M. This is consistent with the fact that the equations remain unchanged by such a modification of f. Note also that if

E 0 = 0 and V -1 (E 0 ) ∩ ω = ∅ (resp. V -1 (E 0 ) ∩ supp(θ) = ∅),
a precised version of this result (see Theorem 3.1 below) actually shows that C 0 (T, ε) ≥ e c/ε for one c > 0 and all time T > 0 (in particular, uniform observability never holds). This is consistent with (and a particular case of) the Guerrero-Lebeau [START_REF] Guerrero | Singular optimal control for a transport-diffusion equation[END_REF] result (first part of Theorem 1.1 above) for in this case, (M, ∇ g f, ω) does not satisfy (GCC). Indeed, a point x 0 ∈ V -1 (E 0 ) satisfies ∇ g f(x 0 ) = 0 and is thus a stationary point of the gradient dynamics.

We refer to Remark 3.5 concerning the additional smoothness assumption for the boundary estimate.

Theorem 1.6. Assume M = Ω where Ω ⊂ R n is any smooth bounded connected open set endowed with g = Eucl the Euclidean metric and q ∈ L ∞ (M). There exist ω ⊂ Ω and constants c ω , δ > 0 such that for any λ > 0, there is a function f λ ∈ C ∞ (Ω) such that:

• (Ω, ∇f λ , ω) satisfies (FC) and T F C (Ω, ∇f λ , ω) ≤ diam(Ω); • T unif (ω) ≥ c ω λ, • λ 2 δ ≤ ∇f λ 2 L ∞ (Ω) ≤ λ 2 diam(Ω) 2 + n.
In particular, for all Λ > 0, there is f ∈ C ∞ (Ω) such that

T unif (ω) ≥ ΛT F C (Ω, ∇f, ω).
The result of Theorem 1.6 is already of interest in dimension one. In this case Ω = (-L, L), the vector field we consider is f ′ (x)∂ x with f ′ ≥ 1 on [-L, L] and f ′ (0) = 1 and the observation set ω is a neighborhood of the boundary point L (note that this would correspond to the case M < 0 in the Coron-Guerrero problem with the notation of [START_REF] Coron | Singular optimal control: A linear 1-D parabolic-hyperbolic example[END_REF]). Note that the function f λ in this result satisfies max f λmin f λ ≃ λ. As a consequence, one cannot even hope to have existence of a constant K > 0 depending only on min x∈[0,L] f ′ (x) (a uniform flushing time) such that T unif (ω) ≤ K T F C (Ω, ∇f λ , ω). However, at this point, it does not seem hopeless that such a constant K depends only on ∇f λ L ∞ (Ω) , at least for a fixed metric.

Remark 1.7. In the case ∂M = ∅, Theorem 1.5 does not seem to suffice to construct functions f, ω having

T unif (ω)
TGCC (M,∇g f,ω) arbitrarily large. In a domain of R n , Theorem 1.5 is however enough to provide counterexamples.

Another drawback of Theorem 1.5 is that it does not produce any useful lower bound in case ω is a whole neighborhood of ∂M (or in the boundary observation case from the whole boundary ∂M).

We remedy these issues in the next section on surfaces of revolution.

Lower bounds on surfaces of revolution

In Theorems 1.5 and 1.6 above, the lower bound of the minimal uniform observability time is essentially due to the contribution of the potential

V (x) = |∇g f(x)| 2 g 4
(and the difference between its maximal and minimal values on M). In this section, we consider a family of geometric settings, namely surfaces of revolution, for which the contribution of the geometry of (M, g) plays an important role. This leads in particular to explicit versions of Theorem 1.2.

The precise description of the geometry of the surfaces we consider is given in Section 4.1 and we only describe here features required to state the result. We may consider either:

1. M = S ⊂ R 3 a smooth compact surface diffeomorphic to the sphere S 2 ; 2. M = S ⊂ R 3 a smooth compact surface diffeomorphic to the disk D;

3. M = S ⊂ R 3 a smooth compact surface diffeomorphic to the cylinder [0, 1] × S 1 ; 4. M = S a smooth compact surface diffeomorphic to the torus T 2 = S 1 × S 1 .

We assume moreover that it has revolution invariance around an axis. In particular, the axis may intersect S in two points (sphere), one point (disk) or no point (cylinder or torus). Except near these points, S has a global coordinate chart (s, θ) ∈ (0, L) × S 1 for some L > 0. In the first three cases, the surface is endowed with the metric g inherited from the Euclidean metric on R 3 which writes

g = ds 2 + R(s) 2 dθ 2 ,
(1.15) on account to the rotation invariance. Here the function R : (0, L) → (0, ∞) describes the shape of S (distance to the revolution axis). In the torus case, we simply endow T 2 with the metric (1.15). We further assume that the function f and q are themselves rotationally invariant, that is f = f(s) and q = q(s) in these coordinates. In this setting (and as opposed to results presented in above Section 1.3.1), the relevant quantities for our analysis are the following.

We define for any c > 0 (that can be chosen) the (θ-invariant) effective potential

V c (s) := c 2 R(s) 2 + |f ′ (s)| 2 4 . (1.16)
Note that, as opposed to the potential appearing in (1.11), this potential V c is different from

|∇g f| 2 g 4 = |f ′ (s)| 2 4
. Moreover, it depends explicitly on the geometry (namely, on R). We shall make the simplifying assumption that V -1 c (min V c ) = {s min } consists in a single point s min .

(1.17)

Note that in case S has a boundary, one may have s min at the boundary (see Section 4 for more precise statements). As in the previous section, we introduce the associated Agmon distance, which simply writes in the first three cases:

d c A (s) = s smin V c (y) -V c (s min )dy .
(1.18)

In the fourth case S = S 1 L × S 1 , an analog of (1.12) still makes sense on S 1 L when V c is defined on S 1 L . We may choose a representation of S 1 L = R/LZ in which s min = 0 + LZ, a definition of the Agmon distance then reads

d c A (s) = min smin∈LZ s smin V c (y) -V c (s min )dy .
(1.19)

We also set

W c (s) = d c A (s) + f(s) 2 , for s ∈ (0, L).
(1.20)

Then, our main result in this geometric context can be (loosely) stated as follows.

Theorem 1.8. Let c > 0, assume that V c satisfies (1.17) and that, in the coordinates (s, θ) ∈ (0, L) × S 1 , we have ω = (0, δ) × S 1 ∪ (Lδ, L) × S 1 (this can be rewritten in an intrinsic way on S). Assuming the observability inequality (1.4) (resp. (1.10)) with constant C 0 (T, ε, ω) (resp. C 0 (T, ε, {L} × S 1 )), there is a sequence ε k → 0 + such that for any δ > 0, there is k 0 (δ) > 0 such that

C 0 (T, ε k , ω) ≥ e 1 ε k W c ω -W c m -Vc(smin)T -δ for all k ≥ k 0 (δ), W c ω = min ω W c , resp. C 0 (T, ε k , {L} × S 1 ) ≥ e 1 ε k W c (L)-W c m -Vc(smin)T -δ for all k ≥ k 0 (δ),
where W c m = inf (0,L) W c . In particular, we have

V c (s min )T unif (ω) ≥ W c ω -W c m , (1.21) V c (s min )T unif ({L} × S 1 ) ≥ W c (L) -W c m .
(1.22) Theorem 1.8 differs from Theorem 1.5 in several respects. First notice that the potential appearing in Theorem 1.5 is |f ′ (s)| 2 4 , that is V 0 (s) with the definition of V c in (1.16). In particular, it does not depend on R: neither does its minimal value, nor the associated Agmon distance and function W 0 . Therefore, in this very particular geometric context, the results of Theorem 1.5 do not depend on the geometry of R, and hence only formulate as a one dimensional result in the s variable. As such, they do not care about the "transverse dynamics" in the θ-variable. Theorem 1.8 overcomes this lack and shows that both have to be taken into account.

Another difference with the estimates of Theorem 1.5 is thatmax V -1 (E0)

f 2 = -max V -1 (E0)
W is here replaced byinf (0,L) W c . This improvement is due to the "one dimensional" underlying framework (in which localization properties of eigenfunctions are better understood).

Again, we remark that the initial problem is invariant by changing f by a constant f + C 0 . In Estimate (1.21), both the potential V c and the quantity W c ω -W c m are as well invariant by this change of the function f.

We now state three particular examples of application of Theorem 1.8. The latter imply Theorem 1.2.

Corollary 1.9. Assume S is a surface of revolution in R 3 diffeomorphic to S 2 (resp. D) and such that R -1 (max R) is a single point (R has a unique maximum). Denote by N, S (resp. N only) the north/south poles of S, which are the only two (resp. the unique) invariant points under the revolution symmetry. Then, for any δ > 0, there exists f δ ∈ C ∞ (S) invariant by rotation such that with ω δ = B g (N, δ) ∪ B g (S, δ) (resp. ω δ = B g (N, δ)) we have 1. (S, ∇ g f δ , ω δ ) satisfies (GCC) (resp. satisfies (FC)) and T GCC (S,

∇ g f δ , ω δ ) = L-2δ ≤ L = dist g (N, S) (resp. T F C (S, ∇ g f δ , ω δ ) = L -δ ≤ L = dist g (N, ∂S));
2. For all c > 0, there is δ 0 , C > 0 such that for all δ ∈ (0, δ 0 )

T unif (ω δ ) ≥ c 2 R(s min ) 2 + 1 4 -1 c log( 1 δ ) -C .
This result proves the first item in Theorem 1.2. In particular, the limit δ → 0 + prevents from the existence of a universal constant K ≥ 1 such that T unif (ω) ≤ K T GCC (M, X, ω). Note that in this construction, the vector fields ∇ g f δ • ∇ g are rotationally invariant and independent of the metric g. Moreover, for δ < δ ′ , the two functions f δ and f δ ′ coincide on M \ ω δ ′ .

In our second result, the geometry is close to that of Corollary 1.9. However, we consider fixed ω and f (and even a fixed vector field), but let the metrics g vary. We denote S 1 L = R/LZ and S 1 = S 1 2π . Corollary 1.10. Assume S = S 1 L × S 1 with coordinates (s, θ) and let f ∈ C ∞ (S 1 L ). Let ω = I ω × S 1 with I ω a nonempty interval such that I ω = S 1 L . Then, there is a constant C > 0 such that for any δ ∈ (0, 1), there exists a function R δ ∈ C ∞ (S 1 L ; R + * ) such that 1. the vector field ∇ g δ f•∇ g δ = f ′ (s)∂ s (defined by (1.15) associated to R δ ) does not depend on the metric g δ , the triple (S, ∇ g δ f, ω) satisfies (GCC) if and only if f ′ = 0 on S 1 L \ I ω , and, in this case, we have

T GCC (S, ∇ g δ f, I ω × S 1 ) = S 1 L \Iω ds f ′ (s) < +∞,
2. for the transport equation (1.6) with viscosity given by the metric g δ under the form (1.15) associated to R δ , we have

T unif (I ω × S 1 ) ≥ Cδ -1/2 , 3. δ 1/2 (1 -Cδ) ≤ min S 1 R δ ≤ δ 1/2 for all δ ∈ (0, 1).
This result implies the second item in Theorem 1.2. Another application is given by the following result, which is an analogue of Corollary 1.10 for the boundary observability problem in the cylindrical geometry.

Corollary 1.11. Assume S = [0, L] × S 1 (whence ∂S = ({0} ∪ {L}) × S 1 ), and let f ∈ C ∞ ([0, L]). Then, for any γ > 2, δ 0 > 0 there is a constant C > 0 such that for all δ ∈ (0, δ 0 ], there exists a function R δ : [0, L] → R + * such that 1. the vector field X = ∇ g δ f • ∇ g δ = f ′ (s)∂ s does not depend on the metric g δ (defined by (1.15) with R δ ); 2. (S, ∇ g δ f) satisfies (FC) (in the sense of Definition 1.4) if and only if f ′ = 0 on [0, L] and in this case,

T F C (M, ∇ g δ f) = L 0 ds f ′ (s)
< +∞;

3. for the transport equation with viscosity and Dirichlet boundary conditions (1.6), and with metric g δ under the form (1.15) associated to R δ , we have

T unif (∂S) ≥ max [0,L] |f ′ | 2 4 -1 1 δ γ/2-1 1 γ/2 -1 -C . 4. (s + δ) γ 2 (1 -C(s + δ)) ≤ R δ (s) ≤ (s + δ) γ 2 for s ∈ [0, L/4] and (L -s + δ) γ 2 (1 -C(L -s + δ)) ≤ R δ (s) ≤ (L -s + δ) γ 2 for s ∈ [3L/4, L];
Note the link between the asymptotic singularity of the metric R δ (s) ∼ (s + δ) γ 2 (i.e. S becomes close to a "conical geometry" for small δ) and the blowup of the minimal time T unif 1 δ γ/2-1 . Note also that the limit case γ = 2, all calculations can be done as well and lead to R δ (s) ∼ s + δ together with T uniflog δ. This corresponds to the case where the geometry of the cylinder degenerates towards that of the disk, leading to the same blowup estimate as in Corollary 1.9.

Observability for positive solutions

As already mentioned, our last result concerns the uniform observability question for positive solutions of (1.6), and is restricted to the case ∂M = ∅. We also assume f ∈ C 3 . Note that if u 0 ∈ L 2 (M; R + ), then the associated solution u to (1.6) satisfies u(t, x) ≥ 0 for a.e. (t, x) ∈ R + × M (see e.g. Theorem 9 in Section 7.1 p369 together with Problem 7 in Section 7.5 in [START_REF] Evans | Partial differential equations[END_REF], or Chapter III, Theorem 7.1, p181 in [START_REF] Ol | Linear and quasilinear equations of parabolic type[END_REF]).

We consider the observability inequality for nonnegative solutions:

C + 0 (T, ε) 2 T 0 ω |u(t, x)| 2 ds(x)dt ≥ u(T ) 2 L 2 (M) ,
for all u 0 ∈ L 2 (M; R + ) and u solution of (1.6). (1.23) and the associated minimal time T + unif (ω) of uniform observability for positive solution, already defined in (1.7). The main result we obtain in this context is the above Theorem 1.3, stating that T + unif (ω) = T GCC (M, ∇ g f, ω). As a byproduct of our analysis, we also obtain a lower estimate on the blow up of the control cost when the Geometric Control Condition is not satisfied. It involves the definition of a quantity that roughly speaking, measures how two points are far from being the image of a trajectory at time t, namely

d ∇g f (x, y, t) = 1 4 inf t 0 | γ(s) -∇ g f(γ(s))| 2 g ds, γ ∈ U t (x, y) , (1.24) U t (x, y) = γ ∈ W 1,∞ ([0, t]; M), γ(0) = x, γ(t) = y .
Note that we have d ∇g f (x, y, t) ≈ d(x, φ -t (y)) for bounded t, where d denotes the Riemannian distance and φ t the flow of the vector field ∇ g f (see Lemma A.2 where d ∇g f (x, y, t) is interpreted as a control cost from x to φ -t (y) with time varying metric).

Proposition 1.12. Assume that (M, ∇ g f, ω) does not satisfy (GCC). Then, we have

d ([0,T ],ω) := sup y∈M inf x∈ω,t∈[0,T ] d ∇g f (x, y, t) > 0.
Moreover, for any δ > 0, there is ε 0 > 0 such that for all ε ∈ (0, ε 0 ), we have

C 0 (T, ε) ≥ C + 0 (T, ε) ≥ e d ([0,T ],ω) -δ ε .
(1.25)

Note that this exponential blowup is a refinement of the Guerrero-Lebeau [START_REF] Guerrero | Singular optimal control for a transport-diffusion equation[END_REF] result (first part of Theorem 1.1 above). However, we provide here with a precise geometric rate (namely d ([0,T ],ω) ) quantifying this blowup phenomenon.

The proofs of Theorem 1.3 and Proposition 1.12 rely on estimates on the kernel of the associated equation. Note that kernel estimates have already been used in [START_REF] Miller | Geometric bounds on the growth rate of null-controllability cost for the heat equation in small time[END_REF] to prove lower bounds for the cost of controllability of the usual heat equation in the short time asymptotics, and in [START_REF] Laurent | Observability of the heat equation, geometric constants in control theory, and a conjecture of Luc Miller[END_REF] to prove observability of positive solutions to the heat equation with optimal constants.

Further remarks

In this section, we collect several remarks and comments related to our results.

1. The principal interest of working with gradient vector fields X = ∇ g f is that the associated operator -∇ g f -ε∆ g can be conjugated to a selfadjoint Schödinger operator (1.8). And the limit ε → 0 + then corresponds to the semiclassical limit, which has been the object of many studies (see e.g. [Wit82, Sim83, HS84, HS85, Hel88, CFKS87, All98, HKN04, HN06]. This conjugation does no longer hold in case X is not a gradient vector field. One could also consider that giving counterexamples with gradient flows is "stronger" than general counter-examples. We do not know wether an analogue of Theorem 1.3 for positive solutions remains true for general vector fields X. This seems to be an interesting open problem.

2. In the context of surfaces of revolution, as presented in Section 1.3.2, we only provide with lower bounds of T unif . It would of course be interesting to obtain related upper bounds on this uniform control time. This would require to provide a very precise description of several spectral quantities (spectral gaps, localization of eigenfunctions at all energy levels...) for the semiclassical Schrödinger operator P ε in (1.8), and seems to be a difficult question. See the companion paper [START_REF] Laurent | On uniform controllability of 1D transport equations in the vanishing viscosity limit[END_REF] for an upper bound in a related one dimensional situation.

3. The one dimensional one well problem is considered in [START_REF] Laurent | On uniform controllability of 1D transport equations in the vanishing viscosity limit[END_REF]. In this very particular situation, we are able to provide with

• improved lower bounds on the minimal time when compared to Theorem 1.5;

• an upper bound on the minimal time.

This requires the knowledge of precise information on the spectral gap and the localization of eigenfunctions at all energy levels E ≥ E 0 (whereas Theorem 1.5 is only concerned with the bottom energy level E 0 ). See also the discussion at the beginning of Section 3.

4. Notice that if one is not interested in null-controllability (i.e. driving the solution exactly to zero at time T ), but rather in approximate controllability with a reasonable cost (and a precision depending on the viscosity ε), one might be satisfied by the following statement.

Proposition 1.13. Suppose ∂M = ∅ and (M, X, ω, T ) satisfies (GCC) (resp. ∂M = ∅ and (M, X, ω, T ) satisfies (FC)). Then, there exist C, C 0 > 0 such that for all

y 0 ∈ L 2 (M), ε ∈ (0, 1] there is h = h ε ∈ L 2 ((0, T ) × ω) with h ε L 2 ((0,T )×ω) ≤ C y 0 L 2 (M) ,
such that the associated solution to (2.1) satisfies

y(T ) L 2 (M) ≤ Ce -C 0 ε y 0 L 2 (M) .
That is to say, one can drive the solutions e -C 0 ε close to zero with a uniformly bounded cost. This result follows from Proposition 5.7 below (a particular case of [GL07, Proposition 3]) together with [START_REF] Laurent | The cost function for the approximate control of waves[END_REF]Appendix]. This can be particularly useful for numerical purposes, since e -C 0 ε = 0 numerically for ε small enough.

In the situation of Theorems 1.5 or 1.8, this means that for intermediate times 5. Note that in the context of revolution surfaces of Section 1.3.2, we prove a complementary result compared to [LL21b, Theorem 1.9/Corollary 1.10]. We prove in Proposition 4.9 that in all cases of Section 1.3.2, for any rotationally invariant set ω, we have (with the notation of [START_REF] Laurent | Observability of the heat equation, geometric constants in control theory, and a conjecture of Luc Miller[END_REF], the constant K eig (ω) being the smallest constant K in the inequality

T ∈ T GCC (M, ∇ g f, ω), T unif (ω) (resp. for T ∈ T F C (M, ∇ g f, ω), T unif (ω) if ∂M = ∅),
ψ λ L 2 (M) ≤ Ce K √ λ ψ λ L 2 (ω) where -∆ g ψ λ = λψ λ ) K eig (ω) = R(s min )d A (ω), with d A (ω) = inf x∈ω d A (x).
In [START_REF] Laurent | Observability of the heat equation, geometric constants in control theory, and a conjecture of Luc Miller[END_REF], we only proved K eig (ω) ≥ R(s min )d A (ω) (and only in case S is diffeomorphic to a sphere). This result is close to that of Allibert [START_REF] Allibert | Contrôle analytique de l'équation des ondes et de l'équation de Schrödinger sur des surfaces de révolution[END_REF], which already proves this in case S is diffeomorphic to a cylinder and the function R has a single local maximum which is non-degenerate.

Acknowledgements. The authors are partially supported by the Agence Nationale de la Recherche under grants SRGI ANR-15-CE40-0018 (for C.L.), SALVE ANR-19-CE40-0004 (for M.L.) and ISDEEC ANR-16-CE40-0013 (for C.L. and M.L.). Part of this research was done when the second author was in École Polytechnique, Centre de Mathématiques Laurent Schwartz UMR7640, 91128 Palaiseau cedex France.

We would like to thank Olivier Glass (who suggested Proposition 1.13), Franck Sueur and Dorian Le Peutrec for interesting discussions related to this work. We are also grateful to Jon Asier Bárcena-Petisco for his comments on a first version of the article. Finally we would like to thank both referees for their careful reading of the manuscript and their constructive remarks that helped to improve the presentation of the paper.

2 Preliminaries: duality, conjugation of gradient flows and (GCC)

Uniform controllability problems and dual formulation

In this section, we reformulate the above uniform observability questions in terms of uniform controllability statements. This follows essentially the classical duality approach, see [START_REF] Dolecki | A general theory of observation and control[END_REF] or [Cor07, Chapter 2.3].

Duality between internal control and observation problems

In this section, we present the controllability problems associated to the above observability questions, and we briefly describe the duality between the control and the observation problems. We introduce the internal control problem

     (∂ t + X + div g (X) -q -ε∆ g )y = 1 ω h, in (0, T ) × M, y = 0, on (0, T ) × ∂M, y| t=0 = y 0 , on M.
(2.1)

Notice that, as opposed to (1.3), the operator appearing in these control problems is

X * = -X -div g (X)
, where the adjoint is taken in the space L 2 (M, d Vol g ).

That the appropriate dual observation problem is (1.3) is a consequence of the following lemma.

Lemma 2.1 (Duality equation). For all solutions u ∈ C 0 ([0, T ]; L 2 (M)) of (1.3) on [0, T ] with Dirichlet boundary conditions and all y ∈ C 0 ([0, T ]; L 2 (M)) solution to (2.1) with h ∈ L 2 ((0, T ) × M), we have (u(T ), y 0 ) L 2 (M) -(u 0 , y(T )) L 2 (M) + T 0 1 ω u(t), h(T -t) L 2 (M) dt = 0. (2.2)
Notice that one passes from the observed evolution to the controlled evolution by changing (X, q) into (-X, qdiv g (X)). The interest of adding a potential term q(x) in these equations is that the free equation (1.9) and the controlled equation (2.1) then have the same form (i.e. the adjoint of a vector field is not a vector field but the adjoint of a first order differential operator is a first order differential operator). Definition 2.2 (Controllability and cost). Given (ω, ε, T ), we say that (2.1) is null-controllable from (ω, T ) if for any y 0 ∈ L 2 (M), there is h ∈ L 2 ((0, T ) × M) such that the associated solution to (2.1) satisfies y(T ) = 0. If (2.1) is null-controllable from (ω, T ), we define for y 0 ∈ L 2 (M) the set U (y 0 ) = ∅ of all such controls h ∈ L 2 ((0, T ) × M), and the cost function

C 0 (T, ε) := sup y0∈L 2 (M), y0 L 2 (M) ≤1 inf h∈U(y0) h L 2 ((0,T )×M) .
As a corollary of Lemma 2.1, together with classical arguments (see e.g. [START_REF] Dolecki | A general theory of observation and control[END_REF] or [Cor07, Chapter 2.3]) we deduce the following statement.

Corollary 2.3 (Observability constant = control cost). Given (ω, ε, T ), Equation (2.1) is null-controllable from (ω, T ) if and only if the observability inequality (1.4) holds. Moreover, we then have

C 0 (T, ε) = C 0 (T, ε).
As a consequence, all lower bounds on C 0 (T, ε) formulated in Theorems 1.5 and 1.8 translate into lower bounds on C 0 (T, ε). The time T unif (ω) is equal to the minimal time of uniform controllability, and all lower bounds on the time T unif (ω) obtained in Theorems 1.5 and 1.8 and their corollaries apply.

The uniform observability result for positive solutions of the heat equation in Theorem 1.3 also has a controllability counterpart. This fact was indeed proved by Le Balc'H [LB20, Theorem 4.1] for the classical heat equation. In the present context, the uniform observability estimate for positive solutions, associated to Theorem 1.3, implies the following controllability result.

Corollary 2.4. Let M be a compact Riemannian manifold with ∂M = ∅, X = ∇ g f where f ∈ C 3 (M), and ω ⊂ M an open subset. Assume that (M, ∇ g f, ω) satisfies (GCC), and T > T GCC (M, ∇ g f, ω). Then, there exist C, ε 0 > 0 so that for any

y 0 ∈ L 2 (M) and 0 < ε ≤ ε 0 , there exists a control h ∈ L 2 ([0, T ], L 2 (ω)) with h L 2 ([0,T ],L 2 (ω)) ≤ C y 0 L 2 (M)
such that the solution of (2.1) satisfies y(T ) ≥ 0.

We refer to Section 5.4 for a proof.

Duality between boundary control and observation problems

We now briefly discuss the boundary case and we refer to [START_REF] Guerrero | Singular optimal control for a transport-diffusion equation[END_REF] for the details. The boundary control problem under interest is

     (∂ t + X + div g (X) -q -ε∆ g )y = 0, in (0, T ) × Int(M), y = θh, on (0, T ) × ∂M, y| t=0 = y 0 , on M, (2.3) where θ ∈ C ∞ (∂M; R) is meant to be a smooth version of 1 Γ , Γ ⊂ ∂M. Solutions of (2.
3) are defined in the sense of transposition, and a well-posedness statement can be written as follows.

Lemma 2.5 (Guerrero-Lebeau [GL07] pp 1814-1815). Assume X is a L ∞ vector field on M with div g (X) ∈ L ∞ (M), q ∈ L ∞ (M)
, and let T > 0. Then, there exists C > 0 such that for all y 0 ∈ H -1 (M), all h ∈ L 2 (0, T ; H -1/2 (∂M)) and all ε > 0, there is a unique solution y to (2.3) in the sense of transposition, which satisfies

y ∈ L 2 ((0, T ) × M) ∩ C 0 ([0, T ]; H -1 (M)) ∩ H 1 (0, T ; H -2 (M)) with y L 2 ((0,T )×M) + y L ∞ (0,T ;H -1 (M)) + ∂ t y L 2 (0,T ;H -2 (M)) ≤ C ε y 0 H -1 (M) + ε 1/2 h L 2 (0,T ;H -1/2 (∂M)) .
Such solutions in particular solve the first equation of (2.3) in the sense of distributions.

Definition 2.6 (Controllability and cost). Given (θ, ε, T ), we say that (2.3) is null-controllable from (θ, T ) if for any y 0 ∈ H -1 (M), there is h ∈ L 2 (0, T ; H -1/2 (∂M)) such that the associated solution to (2.3) satisfies y(T ) = 0. If (2.3) is null-controllable from (θ, T ), we define for y 0 ∈ H -1 (M) the set U (y 0 ) = ∅ of all such controls h ∈ L 2 (0, T ; H -1/2 (∂M)), and the cost function

C 0 (T, ε) := sup y0∈H -1 (M), y0 H -1 (M) ≤1 inf h∈U(y0) h L 2 (0,T ;H -1/2 (∂M)) .
We now describe the link with the boundary observation problem (1.9). We start with the duality identity.

Lemma 2.7 (Duality equation). For all solutions

u ∈ C 0 ([0, T ]; H 1 0 (M)) of (1.9) on [0, T ] and all y ∈ C 0 ([0, T ]; H -1 (M)) solution to (2.1) with h ∈ L 2 (0, T ; H -1/2 (M)), we have u(T ), y 0 H 1 0 ,H -1 -u 0 , y(T ) H 1 0 ,H -1 - T 0 θε∂ ν u| ∂M (t), h(T -t) H 1/2 (∂M),H -1/2 (∂M) dt = 0.
The proof is omitted here and only consists in an integration by parts for smooth solutions, and then a density argument. As in the internal case, classical duality arguments (see [START_REF] Dolecki | A general theory of observation and control[END_REF] or [Cor07, Chapter 2.3]) yield the following statement.

Corollary 2.8 (Observability constant = control cost). Given (θ, ε, T ), Equation (2.3) is null-controllable from (θ, T ) if and only if the observability inequality (1.10) holds. Moreover, we then have

C 0 (T, ε) = C 0 (T, ε).
Again, all lower bounds on C 0 (T, ε) formulated in Theorems 1.5 and 1.8 translate into lower bounds on C 0 (T, ε). The time T unif (θ) is equal to the minimal time of uniform controllability, and all lower bounds on the time T unif (θ) obtained in Theorems 1.5 and 1.8 and their corollaries apply.

The vanishing viscosity limit for gradient flows. Conjugation and reformulation

We focus in this article on the very particular case (1.6) where X is a gradient vector field (with respect to the same metric g defining the viscous perturbation ε∆ g ) of a weight function f : M → R, that is X = ∇ g f • ∇ g . In this case, it is classical (see e.g. [START_REF] Witten | Supersymmetry and Morse theory[END_REF][START_REF] Helffer | Puits multiples en mécanique semi-classique. IV. Étude du complexe de Witten[END_REF]) that the operator -ε∆ g -∇ g f • ∇ gq can be conjugated to a "semiclassical selfadjoint operator". Here, X • Y is the inner product of the two vector fields X and Y given by the metric g in T M.

The first basic computation is the following:

e -f 2ε ∆ g e f 2ε = ∆ g + 1 ε ∇ g f • ∇ g + |∇ g f| 2 g 4ε 2 + ∆ g f 2ε .
We denote by

1 ε 2 P ε := -∆ g + |∇ g f| 2 g 4ε 2 + ∆ g f 2ε - q ε , that is P ε := -ε 2 ∆ g + |∇ g f| 2 g 4 + εq f , (2.4)
where q f = ∆gf 2q. The above computation implies that

e -f 2ε 1 ε 2 P ε e f 2ε = -∆ g - 1 ε ∇ g f • ∇ g - q ε . (2.5)
The interest of this conjugation is that the operator

P ε is selfadjoint in L 2 (M, d Vol g ) endowed with domain D(P ε ) = H 2 (M) ∩ H 1 0 (M). Henceforth, the operator -∆ g -1 ε ∇ g f • ∇ g -q ε is also selfadjoint in L 2 (M, e f 2ε d Vol g ).
Let us now reformulate the uniform observability problem (1.4) in terms of the heat equation involving the operator P ε defined in (2.4).

Note that the constant coefficient one dimensional problem introduced in [CG05] enters the "gradient flow" setting with

M = (0, L) ⊂ R, g = 1, ∆ g = ∂ 2 x , q = 0, f = M x for M ∈ R, and thus ∇ g f • ∇ g = M ∂ x .
In that context, this form together with its formulation (2.4) have already been used in [CG05, Gla10, Lis14, Lis15].

Lemma 2.9. Given T 0 , C 0 , ε > 0 and a function u, the following statements are equivalent.

The function u solves

(∂ t -∇ g f • ∇ g -q -ε∆ g )u = 0, in (0, T 0 ) × Int(M), u = 0 on (0, T 0 ) × ∂M, (2.6) resp. u(T 0 ) 2 L 2 (M) ≤ C 2 0 T0 0 u 2 L 2 (ω) dt, (2.7) resp. u(T 0 ) 2 H 1 0 (M) ≤ C 2 0 T0 0 θε∂ ν u| ∂M 2 H 1/2 (∂M) dt. (2.8)
2. The function v(t, x) = e f(x)/2ε u(t, x) solves

ε∂ t v + P ε v = 0, in (0, T 0 ) × Int(M), v = 0 on (0, T 0 ) × ∂M, (2.9) resp. e -f 2ε v(T 0 ) 2 L 2 (M) ≤ C 2 0 T0 0 e -f 2ε v 2 L 2 (ω) dt, (2.10) resp. e -f 2ε v(T 0 ) 2 H 1 0 (M) ≤ C 2 0 T0 0 θe -f 2ε ε∂ ν v| ∂M 2 H 1/2 (∂M)
dt.

(2.11)

3. The function w(t, x) = v(t/ε, x) = e f(x)/2ε u(t/ε, x) solves    ∂ t w + 1 ε 2 P ε w = 0, in (0, T 0 ) × Int(M), w = 0 on (0, T 0 ) × ∂M,
(2.12)

resp. e -f 2ε w(εT 0 ) 2 L 2 (M) ≤ C 2 0 εT0 0 e -f 2ε w 2 L 2 (ω) dt, resp. e -f 2ε w(εT 0 ) 2 L 2 (M) ≤ C 2 0 εT0 0 θe -f 2ε ε∂ ν w 2 H 1/2 (∂M)
dt.

(2.13)

Proof. Start e.g. with u satisfying (2.6). Using (2.5) and the definition of P ε , Equation (2.6) rewrites equivalently as

∂ t u + e -f 2ε 1 ε P ε e f 2ε u = 0, t ∈ [0, T 0 ].
The function v = e f 2ε u then satisfies (2.9) (and conversely). Setting w(t, x) = v(t/ε, x) then satisfies (2.12), and conversely.

The proof that (2.8)⇔(2.11)⇔(2.13) uses additionally that, on account to the Dirichlet boundary condition, we have ∂ ν (e

f 2ε u)| ∂M = e f 2ε ∂ ν u| ∂M .

(GCC) and controllability of the limit equation ε = 0

In this section, we characterize the observability inequality (1.2) for solutions of (1.1) in terms of the Geometric Control Condition (GCC). In this section, M is always assumed without boundary.

We denote by (φ t ) t∈R the flow of the vector field X, namely

φt (x) = X(φ t (x)), φ 0 (x) = x ∈ M. (2.14)
This flow is globally defined on account to the compactness of M. We consider the following definition of the geometric control condition in the manifold M for the vector field X and the set ω, which we denote (GCC).

Definition 2.10. Let M be a compact manifold without boundary, X a Lipschitz vector field on M, ω ⊂ M, χ ∈ C 0 (M), I ⊂ R and interval and T > 0. We say that

• (M, X, ω, I) satisfies (GCC) if for all x ∈ M, there is t ∈ I such that φ -t (x) ∈ ω;

• (M, X, ω, T ) satisfies (GCC) if (M, X, ω, (0, T )) satisfies (GCC); • (M, X, ω) satisfies (GCC) if there is T > 0 such that (ω, T ) satisfies (GCC); • (M, X, χ, T ) satisfies (GCC) if (M, X, {χ = 0}, T ) satisfies (GCC); • (M, X, χ) satisfies (GCC) if (M, X, {χ = 0}) satisfies (GCC).
In this section, the manifold M is fixed. To lighten notation, we omit the dependence on M in (M, X, ω, I) and we simply write (X, ω, I) instead of (M, X, ω, I) (with a similar notation for the other definitions).

Note in particular that an open set ω satisfying (GCC) must contain all singular points of the vector field X (i.e. all points x ∈ M such that X(x) = 0). We now provide with different reformulations of this property.

Lemma 2.11. Let M be a compact manifold and X a Lipschitz vector field on M. Given ω ⊂ M and T > 0, the following properties are equivalent:

1. (X, ω, T ) satisfies (GCC); 2. t∈(0,T ) φ t (ω) ⊃ M; 3. t∈(0,T ) φ -t (ω) ⊃ M; 4. (-X, ω, T ) satisfies (GCC).
Proof. The definition of (X , ω, T ) satisfying (GCC) is equivalent to: for all x ∈ M, there is t ∈ (0, T ) such that x ∈ φ t (ω). Equivalence between the Items 1 and 2 follows. Item 3 is equivalent to t∈(0,T ) φ T -t (ω) ⊃ φ T (M) = M after having applied φ T , which itself is equivalent to Item 2. Equivalence between Item 4 and Item 3 finally follows from the fact that the flow of -X is (φ -t ) t∈R .

Proposition 2.12. Assume M is a compact manifold, X is a Lipschitz vector field on M, ds is a positive density on M, and q ∈ L ∞ (M). Given ω ⊂ M, χ ∈ C 0 (M) and T > 0, the following statements hold true:

1. If (X, ω, T ) satisfies (GCC), then the observability inequality (1.2) for solutions of (1.1) is true; 2. The observability inequality (1.2) for solutions of (1.1) implies that (X, ω, [0, T ]) satisfies (GCC);

The observability inequality

C 2 0 T 0 M |χ(x)u(t, x)| 2 ds(x)dt ≥ u(T ) 2 L 2 (M) ,
for all u 0 ∈ L 2 (M) and u solution of (1.1). (2.15)

holds true if and only if (X, χ, T ) satisfies (GCC);

4. In all the above observability statements, u(T ) 2 L 2 (M) can be equivalently replaced by u(0) 2 L 2 (M) . The proof below is inspired by that in [START_REF] Dehman | Analysis of the HUM control operator and exact controllability for semilinear waves in uniform time[END_REF][START_REF] Laurent | Uniform observability estimates for linear waves[END_REF] for the wave equation. It is constructive and would also yield a characterisation of the HUM control operator (see e.g. [START_REF] Han-Kwan | Geometric analysis of the linear Boltzmann equation I. Trend to equilibrium[END_REF] or [Léa18, Section 1.2] for more on controllability/stabilization properties for transport equations).

Proof. First notice that for u 0 ∈ L 2 (M), the unique solution to (1.1) is explicitly given by

u(t, x) = e t 0 q•φt-τ (x)dτ u 0 • φ t (x) ∈ C 0 (R; L 2 (M)).
A first direct consequence is the existence of a constant C T,q > 1 such that

C -1 T,q u(0) 2 L 2 (M) ≤ u(T ) 2 L 2 (M) ≤ C T,q u(0) 2 L 2 (M)
, for all solutions to (1.1), which proves Item 4. Next, we write the observation term in (1.2) (the same holds for (2.15) if we replace 1 ω by χ) as

T 0 ω |u(t, x)| 2 ds(x)dt = T 0 M 1 ω (x)|u(t, x)| 2 ds(x)dt = T 0 M 1 ω (x)u 2 0 • φ t (x)e 2 t 0 q•φt-τ (x)dτ ds(x)dt.
Using the change of variable y = φ t (x) (see e.g. [Lee13, Proposition 16.42 p432]), we obtain

T 0 ω |u(t, x)| 2 ds(x)dt = T 0 M 1 ω (φ -t (y))u 2 0 (y)e 2 t 0 q•φ-τ (y)dτ (φ * -t ds)(y)dt
(note that div ds (X) is defined by d dt (φ * t ds)| t=0 = div ds (X)ds, so that this expression simplifies slightly in case div ds (X) = 0). Using that the density is positive on the compact M, we get the existence of C T > 1 such that

C -1 T ds(y) ≤ e 2 t 0 q•φ-τ (y)dτ (φ * -t ds)(y) ≤ C T ds(y) uniformly for (t, y) ∈ [0, T ] × M.
As a consequence, we obtain

C -1 T T 0 M 1 ω (φ -t (y))u 2 0 (y)ds(y)dt ≤ T 0 M 1 ω (x)|u(t, x)| 2 ds(x)dt ≤ C T T 0 M 1 ω (φ -t (y))u 2 0 (y)ds(y)dt.
Hence setting

g ω,T (y) := T 0 1 ω (φ -t (y))dt ∈ L ∞ (M) resp. g χ,T (y) := T 0 χ 2 (φ -t (y))dt ∈ C 0 (M),
we deduce

C -1 T M g ω,T (y)u 2 0 (y)ds(y) ≤ T 0 M 1 ω (x)|u(t, x)| 2 ds(x)dt ≤ C T M g ω,T (y)u 2 0 (y)ds(y). (2.16)
Recalling that ω is an open set and M compact, together with Definition 2.10, we deduce that if (X, ω, T ) satisfies (GCC), then we have the existence of c > 0 such that g ω,T (y) ≥ c for a.e. y ∈ M. The lower bound in (2.16) then implies the observability inequality (1.2), and Item 1 follows. Concerning Item 2, if (X, ω, [0, T ]) does not satisfy (GCC), then there is a point

x 0 ∈ M such that φ -t (x 0 )∩ω = ∅ for all t ∈ [0, T ]. The set ω×[0, T ] being compact, there is a neighborhood U of x 0 such that φ -t (U ) ∩ ω = ∅ for all t ∈ [0, T ]. Setting u 0 = 1 U ,
we have on the one hand that u 0 L 2 (M) > 0. On the other hand, we have 1 ω •φ -t (y) = 0 for all y ∈ U and t ∈ (0, T ). This implies that g ω,T u 0 = g ω,T (y)1 U = 0 a.e. and, according to the upper bound in (2.16), that T 0 ω |u(t, x)| 2 ds(x)dt = 0. This contradicts (2.16) and concludes the proof of Item 2.

Finally, the proof of Item 3 is split in two parts. That (X, χ, T ) satisfies (GCC) implies the observability inequality (2.15) follows as in the proof of Item 1. Now assume that (X, χ, T ) does not satisfy (GCC). Then there is a point x 0 ∈ M such that φ -t (x 0 ) ∩ {χ = 0} = ∅ for all t ∈ (0, T ). Hence, we have g χ,T (x 0 ) = T 0 χ 2 (φ -t (x 0 ))dt = 0. We now choose a sequence of continuous real-valued initial data (u n 0 ) n∈N such that u n 0 L 2 (M) = 1 and (u n 0 ) 2 (x)ds(x) ⇀ δ x0 in the sense of measures on M. The fact that g χ,T is continuous on M together with the upper bound in (2.16) implies that, denoting by u n the solution of (1.1) with initial datum u n 0 , we have

T 0 M |χ(x)u n (t, x)| 2 ds(x)dt ≤C T M g χ,T (y)(u n 0 ) 2 (y)ds(y) → C T δ x0 , g χ,T Meas,C 0 = 0,
which contradicts the observability inequality (2.15), and concludes the proof of Item 3.

General lower bounds without geometric assumption

In this section, we consider a general manifold (with or without boundary) M, and prove the lower bound for the minimal time of uniform controllability provided in Theorem 1.5. We also give a proof of Theorem 1.6 as a corollary. To do this, we use the semiclassical reformulation (2.9)-(2.10) of the problem in Lemma 2.9, as well as exponential decay properties of eigenfunctions of the operator P ε . We rely on the Helffer-Sjöstrand theory as developed in [START_REF] Helffer | Multiple wells in the semiclassical limit[END_REF][START_REF] Helffer | Semi-classical analysis for the Schrödinger operator and applications[END_REF][START_REF] Dimassi | Spectral asymptotics in the semi-classical limit[END_REF]. All results presented in this section apply as well for the semiclassical heat equation. The result of Theorem 1.5 is stated for a potential minimum. However, we shall prove a seemingly more general result, at any energy level in V (M). We shall then explain why this latter result is not more general, and how it can be improved in dimension one. We recall the definition of V in terms of f in (1.11), and define the classically allowed region at any energy level E:

K E = {x ∈ M, V (x) ≤ E}.
We then define the Agmon distance (see e.g. [Hel88, Chapter 3]) to the set K E at the energy level E:

d A,E (x, y) = inf 1 0 (V (γ(t)) -E) + | γ(t)| g dt, γ ∈ U 1 (x, y) , U 1 (x, y) = γ ∈ W 1,∞ ([0, 1]; M), γ(0) = x, γ(1) = y , d A,E (x) = inf y∈KE d A,E (x, y). (3.1)
That is to say, d A,E (x) is the distance of x to the set K E for the (pseudo-)metric (V -E) + g. Here again (V (x) -E) + = max (V (x) -E, 0). We will use, as in (1.12) the notation d A = d A,E0 where E 0 = min M V for the Agmon distance at the bottom energy. Note that d A,E vanishes identically on K E (and only on this set). Finally, an important function in the estimates below is given by

W E (x) = d A,E (x) + f(x) 2 . (3.2)
We shall prove in this section the following result.

Theorem 3.1. Assume the observability estimate (2.7) (resp. the boundary observability inequality (2.8)) for all solutions to (2.6)

with constant C 0 = C 0 (T 0 , ε). Then, for all E ∈ V (M) = [min M V, max M V
] and all δ > 0, there is ε 0 > 0 such that we have for all ε ∈ (0, ε 0 )

C 0 (T 0 , ε) ≥ exp 1 ε min ω W E -max KE W E -δ -ET 0 ,
in the internal observation case, and

C 0 (T 0 , ε) ≥ exp 1 ε min Γ W E -max KE W E -δ -ET 0 ,
in the boundary observation case.

In particular, we have for all E ∈ V (M), for each respective case,

T unif (ω) ≥ 1 E min ω W E -max KE W E , T unif (Γ) ≥ 1 E min Γ W E -max KE W E .
Theorem 1.5 is then the particular case

E = E 0 = min M V in Theorem 3.1. Unfortunately, the function E → 1 E (min ω W E -max KE W E ) is a decreasing function of E. Indeed, • the sets K E are increasing in E, hence the function E → max KE W E = max KE f 2 increases; • E → d A,E (x) is decreasing in E, hence the function E → min ω W E decreases.
Therefore, the estimate of Theorem 3.1 simply reduces to that for E = E 0 , that is Theorem 1.5 in the introduction. This comes from the fact that the estimate involving the term max KE W E is very rough (see Section 3.2 below for a more precise discussion). This can be improved in the one dimensional one well case, see [START_REF] Laurent | On uniform controllability of 1D transport equations in the vanishing viscosity limit[END_REF]. Moreover, it is interesting to notice that the proof of Theorem 3.1 is not more involved than the direct proof of Theorem 1.5, and we shall re-use part of it in case of revolution surfaces in Section 4.

Eigenfunctions of semiclassical Schrödinger operators

In this section, we collect classical results concerning eigenfunctions of semiclassical Schrödinger operators, and some of their decay properties. Recall that

P ε = -ε 2 ∆ g + |∇g f| 2 g 4 + εq f = -ε 2 ∆ g + V + εq f is defined in (2.4
). We first need to prove existence of eigenfunctions near any energy level.

Lemma 3.2 (Existence of eigenfunctions). Assume V ∈ W 1,∞ (M) and q f ∈ L ∞ (M) are both real valued. For all E ∈ V (M) = [min M V, max M V ] and all ε ∈ (0, 1], there is E ε = E + O ε 2/3 and ψ ε ∈ H 2 (M) ∩ H 1 0 (M) such that P ε ψ ε = E ε ψ ε . Note that the O ε 2/3
precision is relatively poor, and can be improved in different situations (e.g. if there is a critical point of V at energy E). These refinements are however not needed here.

Proof. The proof consists in constructing a (very rough) quasimode. Assume first that E is reached by an interior point, i.e. there is x 0 ∈ Int(M) such that V (x 0 ) = E. We then only work in a local chart near x 0 , centered at x 0 (hence we work in R n in a neighborhood of 0).

We take a cutoff function χ ∈ C ∞ c (R) such that χ = 1 in a neighborhood of 0. We set u ε (x) = ε -n 3 χ ε -2 3 |x| , so that u ε is smooth and moreover supported in the chart for all ε < ε 0 with ε 0 sufficiently small, and thus u ε ∈ H 2 (M) ∩ H 1 0 (M). Notice also that

u ε 2 L 2 (M) = |u ε (x)| 2 |g|(x)dx = χ 2 (|y|) |g|(ε 2/3 y)dy = c 0 + O ε 2/3 ,
with c 0 > 0. We now estimate (P ε -E)u ε . For this, we first have εq f u ε L 2 ≤ Cε. Second, we always have the rough estimate

V (x) -E = V (x) -V (0) = O (|x|) so that we have (V -E)u ε 2 L 2 ≤ C |x| 2 χ 2 ε -2 3 |x| dx ε 2n 3 ≤ Cε 4/3 . (3.3)
Third, we have

ε 2 ∆ g u ε 2 L 2 = ε 2 |g| i,j ∂ i g ij |g|∂ j u ε 2 L 2 ≤ Cε 4 i,j ∂ i g ij |g| x j |x| ε -2 3 χ ′ ε -2 3 |x| + g ij |g| x j |x| x i |x| ε -4/3 χ ′′ ε -2 3 |x| 2 dx ε 2n 3 ≤ Cε 4 3 . (3.4)
Combining the above four estimates yields the existence of D, ε 0 > 0 such that for all ε < ε 0 , we have,

(P ε -E)u ε L 2 ≤ Cε 2 3 ≤ Dε 2 3 u ε L 2 . Hence, if E / ∈ Sp(P ε ), this implies (P ε -E) -1 L 2 →L 2 ≥ (Dε 2
3 ) -1 . Finally, the operator P ε being selfadjoint, we have, for

z ∈ C \ Sp(P ε ), (P ε -z) -1 L 2 →L 2 = 1 d(z,Sp(Pε)) , so that, if E / ∈ Sp(P ε ), 1 d(E, Sp(P ε )) ≥ (Dε 2 3 ) -1 .
In any case, this implies d(E, Sp(P ε )) ≤ Dε 2 3 , and using that the spectrum of P ε is purely pointwise, this proves the sought result.

Assume now that E is not reached by an interior point, i.e. E / ∈ V (Int(M)). This means in particular that there is x 0 ∈ ∂M such that V (x 0 ) = E. Then, we again work in a local chart near x 0 , centered at x 0 . In this chart, M is given by R n-1 × R -and x 0 by 0. We denote (x ′ , x n ) ∈ R n-1 × R -local coordinates. We then take χ as above and further define χ ∈ C ∞ c (R), non-identically vanishing, such that supp( χ) ⊂ (-1, 0). We define

u ε (x) = ε -n 3 χ(ε -2 3 x n )χ(ε -2 3 |x ′ |).
One can check that all above properties of u ε are still satisfied, and in particular (3.3)-(3.4). In addition, we have supp u ε ⊂ R n-1 × R - * , and thus

u ε ∈ H 2 (M) ∩ H 1 0 (M).
The remainder of the proof then follows the same as in the first case.

Remark 3.3. Note that near a noncritical value of V , or near the boundary of M, the appropriate local model is -ε 2 ∂ 2 x ± x. Considering concentrating quasimodes of the form χ x ε α leads to

ε 2 ∂ 2 x χ x ε α = ε 2 ε 2α χ ′′ x ε α , and xχ x ε α = ε α x ε α χ x ε α .
Henceforth, the right scaling is given by 2 -2α = α, that is α = 2/3. The quasimode we construct are then O ε 2/3 . If one wants to obtain a better remainder, one could replace χ by an Airy function, as one should replace χ by a Hermite function in the case of bottom of potential (in which case the precision of the quasimode is improved). Also, the remainder ε 2/3 is actually the worst possible case.

The next result states the decay estimates for eigenfunctions in the classically forbidden region, and is a consequence of so-called Agmon estimates (see [START_REF] Helffer | Multiple wells in the semiclassical limit[END_REF][START_REF] Helffer | Semi-classical analysis for the Schrödinger operator and applications[END_REF][START_REF] Dimassi | Spectral asymptotics in the semi-classical limit[END_REF]). Here, it is a particular case of [Hel88, Propositions 3.3.1 and 3.3.4]. Note that with respect to [Hel88, Proposition 3.3.1], our operator P ε contains an additional term, namely multiplication by εq f . However, this contribution is of lower order and can be absorbed in the proof of the Agmon estimates.

Theorem 3.4 (Decay of eigenfunctions in the classically forbidden region).

Assume V ∈ W 1,∞ (M) and q f ∈ L ∞ (M). Let E ∈ V (M) = [min M V, max M V ]
and assume ψ ε ∈ H 2 (M) ∩ H 1 0 (M) and E ε satisfy

P ε ψ ε = E ε ψ ε , ψ ε L 2 (M) = 1, max(E ε -E, 0) = o(1) as ε → 0 + . (3.5)
Then for all δ > 0, there exist C = C(δ), ε 0 = ε 0 (δ) > 0 such that, for all ε < ε 0 , we have

e 1 ε dA,E ψ ε L 2 (M) ≤ Ce δ ε . (3.6)
Assuming further that V, q f ∈ C ∞ (M), we have ψ ε ∈ C ∞ (M) and for all δ > 0 and all smooth vector field Y on M, there exist C = C(δ), ε 0 = ε 0 (δ) > 0 such that, for all 0 < ε < ε 0 , we have

|ψ ε (x)| + |(Y ψ ε )(x)| ≤ Ce -1 ε (dA,E(x)-δ) ,
for all x ∈ M.

(3.7)

Remark 3.5. Note that the smoothness assumption V, q f ∈ C ∞ (M) (as well as the smoothness assumptions in Theorem 1.5) is essentially only used in [HS84, Hel88, DS99] to perform elliptic regularity estimates so that the pointwise estimate (3.7) makes sense. A finer (much less demanding) regularity assumption can be formulated.

As a direct corollary, we have that most of the norm of ψ ε is near

K E , see [Hel88, Corollary 3.3.2].
Corollary 3.6 (Most of the norm is in the classically allowed region).

Let E ∈ V (M) = [min M V, max M V ]
and assume ψ ε , E ε satisfy (3.5). For any open set U containing K E , there is δ, ε 0 > 0 such that for all 0 < ε < ε 0 , we have

ψ ε 2 L 2 (U) = 1 + O e -δ ε .
3.2 Rough localization of eigenfunctions, and a proof of Theorem 3.1

From the decay estimates in the classically forbidden region (Theorem 3.4) and the rough localization of the L 2 mass of eigenfunctions (Corollary 3.6), together with the existence of eigenfunctions at any energy level (Lemma 3.2), we may now deduce a proof of Theorem 3.1. Recall that W E is defined in (3.2). We first prove the following proposition, from which Theorem 3.1 will follow.

Proposition 3.7. Let E ∈ V (M) = [min M V, max M V ] and assume ψ ε ∈ H 2 (M) ∩ H 1 0 (M)
and E ε satisfy (3.5). Then for all δ > 0, there exists ε 0 = ε 0 (δ) > 0 such that, for any open set ω ⊂ M and for all ε < ε 0 , we have

e -f 2ε ψ ε L 2 (M) ≥ e -1 ε (maxK E WE +δ) , (3.8) e -f 2ε ψ ε L 2 (M) ≤ e -1 ε (minM WE -δ) , (3.9) e -f 2ε ψ ε L 2 (ω) ≤ e -1 ε (minω WE -δ) .
(3.10)

Assuming also that V, q f ∈ C ∞ (M) and Γ ⊂ ∂M, we have

e -f 2ε ∂ ν ψ ε L 2 (Γ) ≤ e -1 ε (min Γ WE -δ) . (3.11)
Note that Estimate (3.8) is very rough, due to our lack of knowledge on the localization of ψ ε in the classically allowed region K E . In the one dimensional one well case, this bound can be refined, see the companion paper [START_REF] Laurent | On uniform controllability of 1D transport equations in the vanishing viscosity limit[END_REF].

We first prove Proposition 3.7, and then deduce a proof of Theorem 3.1.

Proof of Proposition 3.7. First, setting

U δ = {x ∈ M, f(x) < max KE f + δ} ⊃ K E , we have e -f 2ε ψ ε L 2 (M) ≥ e -f 2ε ψ ε L 2 (U δ ) ≥ e -1 ε (maxK E f 2 + δ 2 ) ψ ε L 2 (U δ ) ≥ 1 2 e -1 ε (maxK E f 2 + δ 2 ) ,
after having used Corollary 3.6 in the last inequality for ε < ε 0 (δ). This proves (3.8), recalling that d A,E = 0 on K E . Second, notice that (3.9) is a consequence of (3.10). Third, to prove (3.10), we use (3.6) as follows

e -f 2ε ψ ε L 2 (ω) ≤ e -f 2ε e -1 ε dA,E e 1 ε dA,E ψ ε L 2 (ω) ≤ e -1 ε (minω( f 2 +dA,E )) e 1 ε dA,E ψ ε L 2 (ω) ≤ e -1 ε (minω( f 2 +dA,E)) Ce δ ε .
Finally, to prove (3.11), we proceed similarly using (3.7) (instead of the sole (3.6)) with Y = ∂ ν . We have

e -f 2ε ∂ ν ψ ε L 2 (Γ) ≤ e -f 2ε Ce -1 ε (dA,E-δ) L 2 (Γ) ≤ Ce -1 ε (min Γ ( f 2 +dA,E)-δ) ,
which concludes the proof of the proposition.

Proof of Theorem 3.1 from Proposition 3.7 and Lemma 3.2. We use the reformulation in Lemma 2.9 and consider the observability estimate (2.10) for solutions to the evolution equation (2.9). More precisely, we select E ∈ V (M), and we let v ε be the solution to (2.9) associated to the initial condition v ε (0) = ψ ε , where ψ ε is given by Lemma 3.2. That is to say, v ε (t, x) = e -Eε ε t ψ ε (x). We estimate both sides of (2.10).

Firstly, using (3.8), we have

e -f 2ε v ε (T 0 ) L 2 (M) = e -Eε ε T0 e -f 2ε ψ ε L 2 (M) ≥ 1 2 e -Eε ε T0 e -1 ε (maxK E WE + δ 2 ) , for ε < ε 0 (δ). Recalling that E ε = E + O ε 2/3 , this implies the existence of ε 0 (δ) > 0 such that for ε < ε 0 (δ), e -f 2ε v ε (T 0 ) L 2 (M) ≥ 1 2 e -1 ε (ET0+maxK E WE +δ) ,
Secondly, concerning the internal case, using (3.10), we have

T0 0 e -f 2ε v ε 2 L 2 (ω) dt = T0 0 e -2 Eε ε t e -f 2ε ψ ε 2 L 2 (ω) dt = ε 2E ε (1 -e -2 Eε ε T0 ) e -f 2ε ψ ε 2 L 2 (ω) ≤ Ce -2 ε (min ω WE -δ) .
As a consequence of these two estimates, applying (2.10) implies

e -2 ε (ET0+maxK E WE +δ) ≤ C 0 (T 0 , ε) 2 Ce -2 ε (minω WE -δ) ,
which concludes the proof.

In the case of observation from part of the boundary, we simply replace the above use of (3.10) by that of (3.11).

An explicit counter-example for a domain of R n

The purpose of this section is to prove Theorem 1.6. Here M = Ω where Ω ⊂ R n is an open set, endowed with the Euclidean metric. For δ > 0, we may assume, up to an appropriate translation of the domain Ω, that there is η > 0 such that B(0, η) ⊂ Ω.

(3.12)

ω ⊂ (Ω ∩ {x 1 > 0, x 2 > 0, • • • , x n > 0}) , and
B(0, η) ∩ ω = ∅. (3.13)
We let f λ (x 1 , • • • , x n ) be defined as follows:

f λ (t) := t 0 λ 2 s 2 + 1ds and f λ (x 1 , • • • , x n ) := n i=1 f λ (x i ). (3.14)
With this definition, the associated gradient vector field and potential are given respectively by

∇f λ (x 1 , • • • , x n ) = n i=1 f ′ λ (x i )e i = n i=1 λ 2 x 2 i + 1e i , V λ (x) = |∇f λ | 2 4 = λ 2 |x| 2 + n 4 ,
where (e 1 , • • • , e n ) denotes the canonical basis of R n . The proof of Theorem 1.6 now directly follows from the following Lemmata 3.8 and 3.9, when taking λ large enough.

Lemma 3.8. In the above setting, recalling (3.12), (3.13) and that f λ is defined in (3.14), we have T unif (ω) ≥ λη 2 n for all λ > 0. Lemma 3.9. In the above setting, recalling (3.12), (3.13) and that f λ is defined in (3.14), (Ω, ∇f λ , ω) and (Ω, ∇f λ ) both satisfy (FC) and we have for all λ ≥ 0

T F C (Ω, ∇f λ , ω) ≤ T F C (Ω, ∇f λ ) ≤ min v∈R n ,|v|=1,v•ei≥0∀i max x∈Ω x • v -min x∈Ω x • v ≤ diam(Ω). (3.15)
In particular, both are bounded by a constant uniformly in λ ≥ 0.

We recall that (FC) and T F C are defined in this context in Definition 1.4. Notice that the quantity max x∈Ω x • vmin x∈Ω x • v represents the minimal Euclidean distance between two parallel hyperplanes (normal to v) such that Ω is contained between the two hyperplanes.

Proof of Lemma 3.8. The minimum of V λ is reached at x min = 0 and V λ (x min ) = n 4 . The Agmon distance (3.1) at the bottom energy V λ (x min ) can be explicitly computed for points x ∈ B(0, η) ⊂ Ω \ ω. Indeed, for x ∈ Ω we have V λ (x) -V λ (x min ) = λ 2 |x| and thus for x ∈ B(0, η),

d A (x) = λ 2 inf 1 0 |γ(t)|| γ(t)|dt, γ ∈ W 1,∞ ([0, 1]; M), γ(0) = 0, γ(1) = x = λ|x| 2 2 1 0 tdt = λ|x| 2 4 ,
where, by symmetry arguments, we have noticed that the straight line γ(t) := xt reaches the infinimum for x ∈ B(0, η). Recalling that V λ (0) = n 4 and f(0) = 0, application of (1.13) in Theorem 1.5 implies

n 4 T unif (ω) = V λ (0)T unif (ω) ≥ min ω f 2 + d A . (3.16)
By a connectedness argument, for any x ∈ ω, we have from (3.13) that d A (x) ≥ min y∈∂B(0,η) d A (y) ≥ λη 2 4 . Moreover, the condition (3.13) together with the definition of f λ in (3.14) imply that f λ (x) ≥ 0 for x ∈ ω. We thus have

min ω f λ 2 + d A ≥ min ω (d A ) ≥ λη 2 4 .
Combined with (3.16), this concludes the proof of the lemma.

Proof of Lemma 3.9. Notice first that given v ∈ R n with |v| = 1 such that e i • v ≥ 0 for all i ∈ {1, • • • , n}, we have for all x ∈ Ω,

∇f λ (x) • v = n i=1 λ 2 x 2 i + 1e i • v ≥ n i=1 e i • v ≥ n i=1 (e i • v) 2 = |v| = 1.
(3.17)

Second, following Definition 1.4, we extend f λ as a smooth compactly supported function fλ in R n Given x ∈ Ω, we denote by γ x the maximal (global) solution to γx (t) = -∇ fλ (γ x (t)) with γ x (0

) = x, defined in R n for t ∈ R. Given v ∈ R n with |v| = 1 such that e i •v ≥ 0, Estimate (3.17) thus implies that γx (t)•v = -∇ fλ (γ x (t))• v ≤ -1. Integrating this between 0 and t, we obtain γ x (t) • v -x • v ≤ -t. Assuming that T > max y∈Ω y • v -min y∈Ω y • v thus implies γ x (T ) • v -x • v ≤ -T < -max y∈Ω y • v -min y∈Ω y • v , that is x • v -γ x (T ) • v > max y∈Ω y • v -min y∈Ω y • v.
Since x ∈ Ω, this implies γ x (T ) / ∈ Ω. This holds true for any x ∈ Ω. Recalling the definition of (FC) in Definition 1.4, we have obtained that this condition is satisfied by both (Ω, ∇f λ , T ) and (Ω, ∇f λ , ω, T ). Moreover, given the definition of T F C as an infimum, we have also obtained that T F C (Ω, ∇f λ , ω) ≤ T F C (Ω, ∇f λ ) ≤ max y∈Ω y • vmin y∈Ω y • v. Since this holds for all v, we have proved (3.15), which concludes the proof of the lemma.

Surfaces of revolution

General setting

In this section we introduce the geometric setting for the results presented in Section 1.3.2. We are concerned with a revolution surface S ⊂ R 3 being either 1. Case 1: diffeomorphic to a sphere S 2 (in which case ∂S = ∅); 2. Case 2: diffeomorphic to a disk (in which case ∂S is a circle embedded in R 3 ).

3. Case 3: diffeomorphic to a cylinder [0, 1] × S 1 ⊂ R 3 (in which case ∂S consists in two disjoint circles embedded in R 3 , and belonging to two parallel hyperplanes);

We follow [START_REF] Laurent | Observability of the heat equation, geometric constants in control theory, and a conjecture of Luc Miller[END_REF] and [Bes78, Chapter 4B p95] for the precise geometric description of such manifolds. At some places, we also consider the case of the torus T 2 , endowed with a metric invariant in one direction.This setting does not strictly speaking enter the framework of the present section, but is much simpler to describe (and we thus do not mention it in the present section).

Definition and differentiable structure. We assume that (S, g) is an embedded 2D submanifold of R 3 (endowed with the induced Euclidean structure), having S 1 = (R/2πZ) ∼ SO(2) as an effective isometry group. The action of S 1 on S, denoted by θ → R θ (such that R θ S = S) has:

1. exactly two fixed points denoted by N, S ∈ S called North and South poles in Case 1; we write ∂ N S = {N } and ∂ S S = {S};

2. exactly one fixed point denoted by N ∈ S called North pole in Case 2; in this case, we write ∂ N S = {N }, and ∂ S S = ∂S has a single connected component (called "south boundary") which is also invariant by R θ ;

3. no fixed point in Case 3; in this case ∂S has two connected components denoted ∂ N S and ∂ S S (called "north and south boundaries") which are both invariant by R θ .

We denote by P the set of poles, that is P = {N, S} in Case 1, P = {N } in Case 2 and P = ∅ in Case 3 and set

U = S \ P. (4.1)
We now describe a nice parametrization of (S, g), and, in particular, useful coordinates on the set U . We set L = dist g (∂ N S, ∂ S S) and denote by

I L = (0, L) in Case 1, I L = (0, L] in Case 2 , I L = [0, L] in Case 3. (4.2)
We let γ 0 be a geodesic curve of S joining N (resp. ∂ N S in case 3) to S (resp. ∂ S S in Cases 2 and 3). Note in particular that length(γ 0 ) = L. For any θ ∈ S 1 , the isometry R θ transforms the geodesic γ 0 into R θ (γ 0 ), which is another geodesic joining N (resp. ∂ N S) to S (resp. ∂ S S). For every m ∈ U (defined in (4.1)), there exists a unique θ ∈ S 1 such that m belongs to R θ (γ 0 ). The geodesic R θ (γ 0 ) can be parametrized by arclength

ρ : [0, L] → R θ (γ 0 ), ρ(0) ∈ ∂ N S, ρ(L) ∈ ∂ S S, s = dist g (ρ(s), ∂ N S) = L -dist g (ρ(s), ∂ S S),
and there exists a unique s ∈ I L such that ρ(s) = m. We use (s, θ) as a parametrization of U ⊂ S:

ζ : U = S \ P → I L × S 1 m → ζ(m) = (s, θ).
In Case 3, P = ∅ and thus the whole surface S = U is diffeomorphic to the cylinder I L × S 1 via ζ. In Cases 3 and 1, we further need to describe coordinate charts around the poles. In cases 1 and 2, we define another exponential chart (U N , ζ N ) centered at the pole N by

U N = {N } ∪ ζ -1 0, L 2 × S 1 = B g N, L 2 ⊂ S, ζ N : U N → B R 2 0, L 2 , ζ N (N ) = 0.
with the transition map

ζ N • ζ -1 : ζ U ∩ U N = 0, L 2 × S 1 → ζ N U ∩ U N = B R 2 0, L 2 \ {0} (s, θ) → s cos(θ), s sin(θ) .
In Case 1, we add similarly a last exponential chart (U S , ζ S ) centered at the pole S by

U S = {S} ∪ ζ -1 L 2 , L × S 1 = B g S, L 2 ⊂ S, ζ S : U S → B R 2 0, L 2 , ζ S (S) = 0,
with the transition map

ζ S • ζ -1 : ζ U ∩ U S = L 2 , L × S 1 → ζ S U ∩ U S = B R 2 0, L 2 \ {0} (s, θ)
→ (Ls) cos(θ), (Ls) sin(θ) .

We shall need the following notation. For a subset J ⊂ I L , we denote by

C J = ζ -1 (J × S 1 ) = m = ζ -1 (s, θ) ∈ U ; s ∈ J ⊂ U ⊂ S
the R θ invariant set which projects downto J. We will also extend this definition to sets J ⊂ [0, L] by adding the point N if 0 ∈ J (in Cases 1 and 2) and the point S if L ∈ J (in Case 1).

Riemannian structure and operators involved. On the cylinder I L × S 1 , the metric g is given by

(ζ -1 ) * g = ds 2 + R(s) 2 dθ 2 , (4.3)
for some smooth function R : I L → R + * (the function R can be interpreted as the Euclidean distance in R 3 of the point parametrized on S by s to the symmetry axis, see e.g. [LL21b, Section 3]). Since g is a smooth metric on S, [Bes78, Proposition 4.6] gives that R extends to a C ∞ function [0, L] → R + satisfying moreover R(0) = 0, R ′ (0) = 1, R (2p) (0) = 0 for any p ∈ N, in Cases 1 and 2, (4.4)

R(L) = 0, R ′ (L) = -1, R (2p) (L) = 0 for any p ∈ N, in addition in Case 1.
For other parametrizations of surfaces of revolution, or direct computations on the sphere S 2 and the disk D, we refer to [LL21b, Section 3].

Example 4.1. In particular, we consider here the following three examples:

• the unit sphere of R 3 is given by case 1 L = π, s ∈ (0, π), R(s) = sin(s) (and the maximum of R is reached at s 0 = π 2 ) • the unit disk of R 2 is given by case 2 with L = 1 and R(s) = s;

• flat cylinder of length L > 0 and radius R 0 > 0 is given by case 3 with R(s) = R 0 .

In these coordinates, the Riemannian volume form is hence R(s)dsdθ, the Riemannian gradient of a function is

∇ g u = ∂ s u ∂ ∂s + 1 R(s) 2 ∂ θ u ∂ ∂θ , with g(∇ g u, ∇ g u) = |∂ s u| 2 + 1 R(s) 2 |∂ θ u| 2 , (4.5)
and the Laplace-Beltrami operator is given by

∆ s,θ = 1 R(s) ∂ s R(s)∂ s + 1 R(s) 2 ∂ 2 θ .
We define by L 2 (S) := L 2 (S, d Vol g ) the space of square integrable functions, which is also invariant by the action of (R θ ) θ∈S 1 . We will sometime also use the same definition for

L 2 (C J ) := L 2 (C J , d Vol g ) for J ⊂ [0, L].
Another important operator is the infinitesimal generator X θ of the group (R θ ) θ∈S 1 , defined, for u ∈ C ∞ (S), by

X θ u = lim ϑ→0 ϑ -1 (u • R ϑ -u).
(4.6)

In the chart (U, ζ), the action of R θ is given by

(ζ -1 ) * R θ (u, θ ′ ) = (u, θ ′ + θ), so that (ζ -1 ) * X θ = ∂ θ . It is proved in [LL21b, Section 3.2
] that X θ is a smooth vector field on S. Note also that X θ (N ) = X θ (S) = 0 and that its norm is given by g(X θ , X θ )(s, θ) = R(s) (in the coordinates of U ). Now, remark that (R θ ) θ∈S 1 acts as a (periodic) one-parameter unitary group on L 2 (S) by f → f • R θ . The Stone Theorem (see e.g. [RS80, Theorem VIII-8 p266]) hence implies that its infinitesimal generator is iA, where A is a selfadjoint operator on L 2 (S) with domain D(A) ⊂ L 2 (S). Since iAf = X θ f for f ∈ C ∞ (S) (which is dense in D(A)) according to (4.6), we have that A is the selfadjoint extension of X θ i . From now on, we slightly abuse the notation and still denote X θ i for its selfadjoint extension A.

Gradient vector field and conjugated operator. We finally introduce a function f : S → R, at least of class W 2,∞ to define the gradient flow. Throughout this section, we assume that X θ f = 0, i.e. f is invariant by rotation and the same property holds for q. In the coordinates of U , we shall thus simply write f = f(s). These regularity assumptions can be written in the coordinate of U by s → f(s) ∈ W 2,∞ (0, L), with f ′ (0) = 0 (in Cases 1 and 2), and f ′ (L) = 0 (in Case 1). (4.7)

We may now define as in (2.4) the conjugated operator P ε as

P ε = -ε 2 ∆ g + |∇ g f| 2 g 4 + ε ∆ g f 2 -εq = -ε 2 1 R(s) ∂ s R(s)∂ s + 1 R(s) 2 ∂ 2 θ + |f ′ (s)| 2 4 + ε 2 1 R(s) ∂ s (R(s)f ′ (s)) -εq, (4.8)
where the second writing, in the coordinates of U , uses the invariance of f. Note that the last term in this expression acts as a multiplication operator by a function in L ∞ (S) with size ε. We shall often consider it as a lower order term, and keep the shorter notation ∆ g f in place of 1 R(s) ∂ s (R(s)f ′ (s)). Since both g and f are invariant by the action of R θ , we have

[X θ , P ε ] = 0.
(4.9)

Moreover, P ε is selfadjoint in L 2 (S, d Vol g ) with domain H 2 (S) ∩ H 1 0 (S) (= H 2 (S)
in Case 1), and has compact resolvent. Therefore, the operators P ε and X θ share a common basis of eigenfunctions (see e.g. [LL21b, Section 3.2] for a proof). If λ ∈ R is an eigenvalue of P ε , then (in the coordinates of U ) the associated eigenfunction can be written as

e ikθ v(s) with k ∈ Z, v ∈ H 2 loc (I L ) ∩ L 2 ((0, L), R(s)ds) satisfying - ε 2 R(s) ∂ s (R(s)∂ s v) + ε 2 k 2 R(s) 2 v + |f ′ (s)| 2 4 + εq f v = λv, (4.10)
together with v(L) = 0 in Case 2 and v(0) = v(L) = 0 in Case 3.

Restoring the dependence of the eigenelements in the parameter ε, we call the normalized eigenfunctions of P ε : ϕ ε k,n = e ikθ v ε k,n (s) with eigenvalues λ ε k,n , where n ∈ N. In particular, for all ε > 0, we can write

L 2 (S) = ⊕ ⊥ (k,n)∈Z×N span(ϕ ε k,n ). We further denote L 2 k = ker(X θ -ik) = ϕ ∈ L 2 (S); ϕ |U = e ikθ f (s), f ∈ L 2 ((0, L), R(s)ds) , and
H 2 k = D(P ε ) ∩ L 2 k = H 2 (S) ∩ H 1 0 (S) ∩ L 2 k .
The commutation property (4.9) implies that for all ε > 0, P ε H 2 k ⊂ L 2 k , so we can define the operator

P (k) ε = P ε|L 2 k , with domain H 2 k , (4.11)
which is selfadjoint. This can be seen for instance directly on the simultaneous diagonalization which implies an isometry

L 2 (S) ≈ ℓ 2 (Z × N) where L 2 k ≈ {(k, n) |n ∈ N } as a closed subspace of ℓ 2 (Z × N).
The fact that P ε has compact resolvent implies that this is also the case for P

P (k) ε w = - ε 2 R(s) ∂ s (R(s)∂ s w) + ε 2 k 2 R(s) 2 + |f ′ (s)| 2 4 + εq f w. (4.12)

The conditions (GCC) and (FC) on surfaces of revolution

In this section, we characterize the conditions (GCC) (see Definition 2.10 if ∂M = ∅) and (FC) (see Definition 1.4 if ∂M = ∅) in the above very particular geometry, and further assuming that the observation region ω is rotationally invariant as well. In case ∂M = ∅, we consider the two cases of internal and boundary observation, and describe the associated minimal times T F C .

Proposition 4.2. Let δ > 0 and recall that f is assumed to be θ-invariant.

1. In Case 1, we consider the set ω = B g (N, δ) ∪ B g (S, δ); then (S, ∇ g f, ω) satisfies (GCC) if and only

if f ′ (s) = 0 for all s ∈ [δ, L -δ] and T GCC (S, ∇ g f, ω) = L-δ δ ds f ′ (s) .
2. In Case 2, we consider the set ω = B g (N, δ); then (S, ∇ g f, ω) satisfies (FC) if and only if f ′ (s) = 0

for all s ∈ [δ, L] and T F C (S, ∇ g f, ω) = L δ ds f ′ (s) .
3. In Case 3, (S, ∇ g f) satisfies (FC) if and only if f ′ (s) = 0 for all s ∈ [0, L] and

T F C (S, ∇ g f) = L 0 ds f ′ (s) .
Note that in Case 2, the situations f ′ > 0 and f ′ < 0 play two different roles (see the proof below). Indeed, in case f ′ > 0, all trajectories of -f ′ enter the controlled region ω, whereas, in case f ′ < 0, all trajectories of -f ′ flow out of the domain S through ∂ S (without passing into ω). However, the definition of (FC) in Definition 1.4 does not make a distinction between these two situations.

Proof. We only prove the second item; the other two items are proved similarly. According to Definition 1.4, (4.5) and the θ-invariance of f, it suffices to check under which conditions the solutions to ṡ(t) = -f ′ (s(t)) all enter ω = B g (N, δ) (resp. all exit S, that is satisfy s(T ) > L). If there is s 0 ∈ [δ, L] such that f ′ (s 0 ) = 0, then the associated solution satisfies s(t) = s 0 ∈ [δ, L] for all t ∈ R, and (S, ∇ g f, ω) does not satisfy (FC).

If f ′ > 0 on [δ, L], then s(t) is decreasing, and for any σ 0 , σ 1 ∈ R, one has

s(σ1) s(σ0) ds f ′ (s) = σ 0 -σ 1 . The longest trajectory that does not enter ω is such that s(0) = L and s(T ) = δ so that T = L δ ds f ′ (s) . This proves T F C (S, ∇ g f, ω) = L δ ds f ′ (s) in this case. Finally, if f ′ < 0 on [δ, L],
then s(t) is increasing, and for any σ 0 , σ 1 ∈ R, one has

s(σ1) s(σ0) ds f ′ (s) = σ 0 -σ 1 .
The longest trajectory that does not enter ω is such that s(0) = δ and s(T ) = L so that -T = L δ ds f ′ (s) . This proves T F C (S, ∇ g f, ω) = -L δ ds f ′ (s) in this case, and hence the proposition.

Existence of eigenfunctions

One may consider different asymptotic regimes in the parameters ε → 0 + and k → +∞. Note that the case k bounded would correspond to the one-dimensional situation treated in the companion paper [START_REF] Laurent | On uniform controllability of 1D transport equations in the vanishing viscosity limit[END_REF]. Here, we shall consider the limit k → +∞ and make the following choice of the parameter ε:

ε = ε k = ck -1 (4.13)
considered as a semiclassical parameter, where c > 0 is a fixed parameter (i.e. which does not depend on k) that will be chosen but fixed. All constants that appear below might depend on c. The analysis of the asymptotic of the constant involved as c → 0 (low level of rotation) or c → ∞ (high level of rotation) would be interesting but would require much more work. In view of (4.12), the choice (4.13) naturally leads to consider

s → V c (s) := c 2 R(s) 2 + |f ′ (s)| 2 4 , (4.14)
as the effective potential of the operator P (k) ε k in the semiclassical limit ε = ε k = ck -1 → 0 + . In particular, the operator P (k) ε is now a semiclassical operator with small parameter ε k and (4.12) can be rewritten

P (k) ε k w = - ε 2 k R(s) ∂ s (R(s)∂ s w) + (V c (s) + ε k q f ) w. (4.15)
In the present section, we recall the existence of eigenfunctions (Analogue of Lemma 3.2 above) associated to any value of the effective potential V c . More precisely, in the chosen regime (4.13), for any s 0 ∈ I L (recall that I L is defined in (4.2)), we construct a sequence ψ k such that P ε k ψ k = (V c (s 0 ) + r(k))ψ k , with r(k) → 0 as k → +∞. As in Section 3.1, the precision r(k) might depend on whether V ′ c (s 0 ) = 0 or V ′ c (s 0 ) = 0 but we will only state the worst estimate, which is sufficient for our needs. We shall later on prove localization properties of the ψ k 's assuming further that V c (s 0 ) = min IL V c (which is a global assumption). We recall the choice (4.13) and the definition (4.14).

Lemma 4.3 (Existence of eigenfunctions). For all s 0 ∈ I L , there is k 0 > 0 such that for all k ∈ N, k ≥ k 0 , there exists

ψ k ∈ H 2 (S) ∩ L 2 k in case 1 (resp. ψ k ∈ H 2 (S) ∩ H 1 0 (S) ∩ L 2 k in cases 2-3), and µ k ∈ R such that µ k = V c (s 0 ) + O 1 k 2/3 = c 2 R(s 0 ) 2 + |f ′ (s 0 )| 2 4 + O 1 k 2/3 , P ε k ψ k = µ k ψ k , ψ k L 2 (S) = 1, ψ k (s, θ) = e ikθ ϕ k (s),
with P ε defined in (4.8) and ε k in (4.13).

The proof is very similar to Lemma 3.2. Indeed, the study of the operator P ε k in L 2 k reduces to the 1D problem for the operator P (k) ε defined in (4.11) and (4.15). The proof consists first in constructing quasimodes exactly as in the proof of Lemma 3.2. Deducing existence of an exact eigenfunction from a quasimode requires the use of the right selfadjoint extension P (k) ε k . This issue is however treated in detail in [START_REF] Laurent | Observability of the heat equation, geometric constants in control theory, and a conjecture of Luc Miller[END_REF]Lemma 3.6]. Note that s 0 ∈ I L implies that it cannot be a pole so that R(s 0 ) -1 is finite.

Geometric assumptions and the Agmon distance

The next step is to study the behavior of the eigenfunction ψ k constructed in Lemma 4.3. This will require some global assumptions on the effective potential V c . Recall that V c is defined in (4.14), is continuous on I L , and tends to +∞ near to the poles. Indeed, in Cases 1 and 2 we have for instance

V c (s) ∼ s→0 + c 2 R(s) 2 ∼ s→0 + c 2 s 2 → s→0 + +∞,
as a consequence of (4.4) and (4.7) (and similarly when s → L -in Case 1). As a consequence, V c admits a minimum on the interval I L , which we denote by

V min = min s∈IL V c (s) ∈ R * + .
In the following, we make the following assumption (a precised version of (1.17)) on the set where V c reaches its global minimum.

Assumption 4.4. The set V -1 c (V min ) = {s min } ⊂ I L consists in a single point. Note that this is assumption is generic. Here it is not strictly needed to prove the main results, but simplifies the presentation and statements slightly. We again introduce the relevant Agmon distance at the minimal energy level V min , defined in the coordinates of U by the eikonal equation

(d c A ) ′ (s) 2 -(V c (s) -V c (s min )) = 0, d c A (s min ) = 0, sgn((d c A ) ′ (s)) = sgn(s -s min ),
or, more explicitly, for s ∈ I L , by (1.18). In view of the W 2,∞ regularity of f on S and the definition of V c in (4.14), the function d c A is of class C 2 away from s min , 0 and L and is locally Lipschitz on I L . Note that this includes Lipschitz regularity up to the boundary s = L in Case 2 and to both boundaries s = 0, L in Case 3. We also consider d c A as a θ-invariant function on the surface S.

d c A (s) = -c log(L -s) + O (1) , as s → L -, in Case 1. (4.17)
Proof. We only consider the asymptotics as s → 0 + , that is, prove (4.16). The proof of (4.17) follows the same. Remark that according to (4.4), we have 1 R(y) → +∞ as y → 0 + with

R(s) = s + O s 3 , when s → 0 + .
As a consequence, since f ∈ W 2,∞ (I L ), we have

V c (s) = c 2 s 2 + O (1) , when s → 0 + .
With (1.18), we obtain d c A (s) =

Upper bounds for eigenfunctions: Agmon estimates

As far as upper bounds on ψ k are concerned, we have the following Agmon type estimate.

Proposition 4.6. Under Assumption 4.4, assume that

µ k = V c (s min ) + r(k) with r(k) → 0 and ψ k ∈ H 2 (S) ∩ L 2 k (S) solves P ε k ψ k = µ k ψ k , on S, ψ k | ∂S = 0, ψ k L 2 (S) = 1,
with P ε defined in (4.8) and ε k in (4.13). Then for all 0 < δ < 1, there exist C = C(δ), k 0 = k 0 (δ) > 0 such that, for all k ∈ N, k ≥ k 0 , the following integral is well defined with the estimate

S e 2(1-δ) d c A (m) ε k |ψ k | 2 (m)d Vol g (m) ≤ Ce 2 δ ε k .
Also, if ∂S = ∅ (that is, in Cases 2 and 3), for all 0 < δ < 1, there exist

C = C(δ), k 0 = k 0 (δ) > 0 such that, for all k ∈ N, k ≥ k 0 , ∂ s ψ k (L, •) 2 H 1 (S 1 ) ≤ Ce -2(1-δ) d c A (L) ε k , ∂ s ψ k (0, •) 2 H 1 (S 1 ) ≤ Ce -2(1-δ) d c A (0) ε k , (4.18)
where the last estimate (at s = 0) holds true in Case 3 only.

Note that given the asymptotic expansion of d c

A in Lemma 4.5, this estimate implies that ψ k vanishes strongly near the poles of S. The proof is made with classical Agmon type identity with some care with respect to the degeneracy at the poles. It is very similar to the one performed in [LL21b, Theorem 3.9] and we omit it. Note that, as opposed to [LL21b, Theorem 3.9], we do not assume here that the minimum be non-degenerate, and only deduce an estimate with loss (δ > 0), which is sufficient for our needs. In the non-degenerate case, one can take δ = 0 in this estimate and replace the right hand-side by a polynomial bound of the type ε -M k , see [LL21b, Theorem 3.9]. The proof of the boundary estimate also requires a bootstrap argument to estimate higher H s norms and the use of a trace estimate, see [Hel88, Propositions 3.3.1 and 3.3.4].

We obtain the following two direct Corollaries.

Corollary 4.7. Under the assumptions of Proposition 4.6, for all 0 < δ < 1, there exist C = C(δ), k 0 = k 0 (δ) > 0 such that, for all k ∈ N, k ≥ k 0 , and for all rotationally invariant set ω, we have

ω |ψ k | 2 (m)d Vol g (m) ≤ Ce -2 ε k ((1-δ)d c A (ω)-δ) , with d c A (ω) = inf s∈ω d c A (s).
Proof. This is a direct consequence of the following estimate:

e 2(1-δ) d c A (ω) ε k ω |ψ k | 2 (m)d Vol g (m) ≤ ω e 2(1-δ) d c A (m) ε k |ψ k | 2 (m)d Vol g (m) ≤ S e 2(1-δ) d c A (m) ε k |ψ k | 2 (m)d Vol g (m) ≤ Ce 2 δ ε k ,
where we have used Proposition 4.6 in the last inequality.

Corollary 4.8 (most of the norm is close to the minimum). Under the assumptions of Proposition 4.6, for any ρ > 0, there exists k 0 ∈ N so that

ψ k 2 L 2 (C (s min -ρ,s min +ρ) ) ≥ 1/2, for all k ∈ N, k ≥ k 0 .
Proof. Applying Corollary 4.7 with ω = S \ C (smin-ρ,smin+ρ) ensures that for any δ > 0, there is k 0 ∈ N such that for all k ≥ k 0 , we have

ω |ψ k | 2 (m)d Vol g (m) ≤ C(δ)e -2 ε k ((1-δ)d c A (ω)-δ) , with d c A (ω) = inf s∈ω d c A (s).
From Assumption 4.4, we have d c A (s) > 0 for all s ∈ I L \ {s min }. That ω is closed and does not contain s min implies that d c A (ω) > 0. Then we fix δ > 0 small enough so that

(1 -δ)d c A (ω) -δ > 0. There is k 0 ∈ N such that we have ψ k 2 L 2 (ω) → 0 for k ≥ k 0 ,
which implies the result.

Lower bounds for eigenfunctions: Allibert estimates

In Corollary 4.7, we proved that the family of eigenfunctions ψ k decays on ω at least like e

- d c A (ω) ε k
. The purpose of this section is to prove the converse, i.e. that the ψ k 's decay at most like e

- d c A (ω) ε k
up to δ loss. This comes from the particular one-dimensional underlying context. We follow in this section the method of Allibert [START_REF] Allibert | Contrôle analytique de l'équation des ondes et de l'équation de Schrödinger sur des surfaces de révolution[END_REF]. More precisely, we prove the following estimates. Proposition 4.9. Under Assumption 4.4, assume further that

µ k = V c (s min ) + r(k) with r(k) → 0, (4.19)
and

ψ k ∈ H 2 (S) ∩ L 2 k (S) solves P ε k ψ k = µ k ψ k , on S, ψ k | ∂S = 0, ψ k L 2 (S) = 1, ψ k (s, θ) = e ikθ ϕ k (s),
with P ε defined in (4.8) and ε k in (4.13).

Then, for any η, δ > 0, there exist k 0 , C > 0 so that

ψ k L 2 (C B(s,η) ) ≥ Ce -1 ε k (d c A (s)+δ) ,
for all k ≥ k 0 and s ∈ I L s.t. dist(s, P) > η.

Note that η is a safety distance to the set of poles P defined in Section 4.1. The proof of Proposition 4.9 relies on two lemmata, which we give in the next section.

Two preliminary lemmata

In this section, we assume that the assumptions of Proposition 4.9 are satisfied. In particular, the eigenfunctions under consideration are of the form ψ k (s, θ) = e ikθ ϕ k (s). We define the following "semiclassical energy densities" of the eigenfunctions ψ k , for s ∈ I L by

E k (s) := ε 2 k |∂ s ψ k | 2 (s) + (V c (s) -µ k + 1)|ψ k | 2 (s) (4.20) = ε 2 k |ϕ ′ k | 2 (s) + (V c (s) -µ k + 1)|ϕ k | 2 (s), E + k (s) := ε 2 k |∂ s ψ k | 2 (s) + (V c (s) -µ k )|ψ k | 2 (s) (4.21) = ε 2 k |ϕ ′ k | 2 (s) + (V c (s) -µ k )|ϕ k | 2 (s).
Note that according to elliptic regularity, we have ϕ k ∈ H 2 loc (I L ) and, due to Sobolev embeddings,

ϕ ′ k , E k , E + k ∈ C 0 (I L )
and in particular E k , E + k are defined everywhere on I L . For s, t ∈ I L , we define I s,t to be the interval between the real numbers s and t, that is, either [s, t] or [t, s]. We also set E α = {s ∈ I L ; |s -p| ≥ α, for all p ∈ P, and |ss min | ≥ α}.

(4.22)

Recall that P is the set of poles defined in Section 4.1 (and is aimed at covering all Cases 1-3 at the same time) and s min is the point at which V c reaches its minimum.

Lemma 4.10. Assume the assumptions of Proposition 4.9 and recall that E + k is defined in (4.21). Then, for any α, δ > 0, there exists k 0 > 0 so that for all s, t ∈ I L such that I s,t ⊂ E α , we have

E + k (t) ≤ e 2 ε k (|d c A (s)-d c A (t)|+δ) E + k (s), for all k ≥ k 0 .
Lemma 4.10 provides with a Grönwall type estimate on the energy E + k , with a precise description of the constant, under the condition that we remain finitely away from s min . It is an analogue of [All98, Lemma 12] in our setting (with an additional uniform dependence).

Note that

|d c A (s) -d c A (t)| = d c A (s, t)
is the Agmon distance between s and t at the lowest energy level. Proof of Lemma 4.10. On the interval z ∈ I s,t ⊂ E α , we differentiate E + k , yielding

(E + k ) ′ (z) = 2ε 2 k Re(ϕ k ′ ϕ ′′ k ) + V ′ c (z)ϕ 2 k + 2(V c (z) -µ k ) Re(ϕ k ϕ k ′ ).
We recall the choice of ε k in (4.13) from the definition of P ε k in (4.8), and the definition of

P (k)
ε k in (4.10)-(4.15) that we have

µ k ϕ k = P ε k ϕ k = -ε 2 k ϕ ′′ k -ε 2 k R ′ R ϕ ′ k + V c ϕ k + ε k q f ϕ k .
Replacing ε 2 k ϕ ′′ k in the above identity yields

(E + k ) ′ (z) = 2 (V c (z) -µ k + ε k q f )) Re(ϕ k ϕ k ′ ) -2ε 2 k R ′ R |ϕ ′ k | 2 + V ′ c (z)|ϕ k | 2 + 2(V c (z) -µ k ) Re(ϕ k ϕ k ′ ) = 4(V c (z) -µ k ) + 2ε k q f Re(ϕ k ϕ k ′ ) -2ε 2 k R ′ R |ϕ ′ k | 2 + V ′ c (z)|ϕ k | 2 . (4.23)
First, using the continuity of V c on I L and the compactness of E α in I L , we see that V c reaches its minimum on E α . This, together with the fact that

s min / ∈ E α , implies that C -2 α ≤ V c (s) -V c (s min ) ≤ C 2 α
uniformly for s ∈ E α . Recalling (4.19), this yields the existence of k 0 (α) such that for k ≥ k 0 (α), we have V cµ k ≥ 1 2Cα > 0 on E α . We now estimate each of the terms in (4.23). We first have

|4(V c -µ k ) Re(ϕ k ϕ k ′ )| ≤ 4ε -1 k V c -µ k (ε k |ϕ ′ k |) V c -µ k |ϕ k | ≤ 2ε -1 k V c -µ k ε 2 k |ϕ ′ k | 2 + (V c -µ k )|ϕ k | 2 = 2ε -1 k V c -µ k E + k .
Moreover, according to (4.19), there exists a constant C α > 0 such that we have

√ V c -µ k ≤ V c -V c (s min )+ C α |r(k)| uniformly for s ∈ E α .
Together with the above inequality, this implies

|4(V c -µ k ) Re(ϕ k ϕ k ′ )| ≤ ε -1 k (2 V c -V c (s min ) + C α |r(k)|)E + k , on E α .
Second, we have

V ′ c |ϕ k | 2 = |V ′ c | V c -µ k (V c -µ k )|ϕ k | 2 ≤ 2C α V ′ c L ∞ (Eα) E + k , on E α .
Third, we have

|ε k q f Re(ϕ k ϕ k ′ )| ≤ C α E + k , on E α . Finally, since R ′ R is bounded on E α , we have ε 2 k R ′ R |ϕ ′ k | 2 ≤ C α ε 2 k |ϕ ′ k | 2 ≤ C α E + k on E α .
Combining the last four estimates in (4.23) yields for another constant C α > 0 and for all k ≥ k 0 (α)

(E + k ) ′ (z) ≤ 2ε -1 k V c (z) -V c (s min ) + C α |r(k)| + C α ε k E + k (z), for all z ∈ E α .
Applying the Grönwall Lemma on the interval I s,t contained in E α yields

E + k (s) ≤ e 2 ε k s t √ Vc(z)-Vc(smin)dz +Cα|r(k)|+Cαε k E + k (t),
which is the sought result.

The next Lemma is aimed at giving a rough Grönwall type estimate for the energy E k , without precise constants. The interest of this less precise result is that it remains true close to the minimum s min . This allows to compensate the fact that Lemma 4.10 is not uniform when s is close to s min . Similarly as E α in (4.22), we define F α = {s ∈ I L ; |s -p| ≥ α, for all p ∈ P}.

Lemma 4.11. Assume the assumptions of Proposition 4.9 and recall that E k is defined in (4.20). For any α > 0, there exist C α , D α > 0 so that for all s, t ∈ F α , we have

E k (s) ≤ C α e 2 ε k Dα|s-t| E k (t), for all k ≥ k 0 .
Lemma 4.11 is an analogue of [All98, Lemma 11] in our setting. Recall that P is the set of poles defined in Section 4.1.

Proof. The proof is quite close to that of Lemma 4.10. We only use the fact that there exists C α so that

C -1 α ≤ V c -µ k +1 ≤ C α on F α if k ≥ k 0 (α). This gives a constant D α > 0 such that D -1 α ε 2 k |ϕ ′ k | 2 + |ϕ k | 2 ≤ E k ≤ D α ε 2 k |ϕ ′ k | 2 + |ϕ k | 2 on F α .
The same computation as in (4.23) gives

E ′ k = (E + k ) ′ + 2 Re(ϕ k ϕ k ′ ) = (4(V c -µ k ) + r k + 1) Re(ϕ k ϕ k ′ ) -2ε 2 k R ′ R |ϕ ′ k | 2 + V ′ c (z)|ϕ k | 2 ,
with sup Fα r k → 0 as k → +∞. As a consequence, for k ≥ k 0 (α), we have constants

C ′ α , C ′′ α such that for all z ∈ F α |(E k ) ′ (z)| ≤ ε -1 k C ′ α ε 2 k |ϕ ′ k | 2 + |ϕ k | 2 ≤ ε -1 k C ′′ α E k (z),
which allows to conclude as in the proof of Lemma 4.10 above by Grönwall estimates.

4.6.2 Proof of Proposition 4.9 from Lemmata 4.10 and 4.11

The sketch of the proof of Proposition 4.9 is as follows:

1. the total mass of ϕ k is dominated by its mass near the minimum s min via Corollary 4.8;

2. the mass near s min is dominated by the energy at s minρ via Lemma 4.11 (with a small loss if ρ is small); Application of Lemma 4.11 gives for any u ∈ (s minρ, s min + ρ) (so that dist(u, P) > η)

E k (s min -ρ) ≥ C -1 e -4 ε k Dρ E k (u) ≥ C -1 e -δ ε k E k (u).
Integrating in u ∈ (s minρ, s min + ρ) gives

E k (s min -ρ) ≥ C 2ρ e -δ ε k smin+ρ smin-ρ E k (u) du ≥ C 1 e -δ ε k , (4.24)
where we have used E k (u) ≥ 1 2 |ϕ k | 2 (u) (for k large enough) and Corollary 4.8. Taking y ∈ [s, s + η/2], we still have y ≤ s minρ from the definition of ρ. Choosing α = ρ/2 ≤ η, we can check that [y, s minρ] ⊂ E α , so that Lemma 4.10 applies on E α and gives

E + k (y) ≥ e -2 ε k (|d c A (y)-d c A (smin-ρ)|+δ) E + k (s min -ρ) ≥ e -2 ε k (d c A (y)+δ) E + k (s min -ρ). (4.25)
where we have noticed

|d c A (y) -d c A (s min -ρ)| = d c A (y) -d c A (s min -ρ) ≤ d c A (y) -d c A (s min ) = d c A (y). Since s min is a strict minimum, there are constants k 0 (ρ, δ) = k 0 (δ, η) > 0 and C(ρ, η) = C(δ, η) > 0 such that for k ≥ k 0 (δ), V c (s min -ρ) -µ k = V c (s min -ρ) -V c (s min ) + r(k) ≥ C(δ, η) -1 .
This implies

E k (s min -δ) = E + k + |ϕ k | 2 ≤ E + k + C(δ, η)(V c -µ k )|ϕ k | 2 ≤ 1 + C(δ, η) E + k (s min -δ)
, where all functions are taken at the point (s minδ). Combining this estimate together with (4.25) and (4.24) yields Now, turning to the global manifold and recalling that ψ k (s, θ) = e ikθ ϕ k (s), (4.5) and (4.13), we have

E k (y) ≥ E + k (y) ≥ Ce -2 ε k (d c A (y)+δ) E + k (s min -ρ) ≥ Ce -2 ε k (d c A (y)+δ) E k (s min -ρ) ≥ Ce -2 ε k (d c A (y)+2δ) ≥ Ce
C [s,s+η/2] E k (y)R(y)dy = ε 2 k C [s,s+η/2] |∇ g ψ k | 2 d Vol g + C [s,s+η/2] |∇ g f| 2 4 -µ k + 1 |ψ k | 2 d Vol g ≤ C ψ k 2 H 1 (C [s,s+η/2] ) . (4.27)
Finally, an interpolation estimates together with P ε k ψ k = µ k ψ k and the definition of P ε k in (4.8) gives

ψ k 2 H 1 (C [s,s+η/2] ) ≤ C η ψ k 2 L 2 (C [s,s+η/2] ) + ψ k L 2 (C [s-η,s+η] ) ∆ g ψ k L 2 (C [s-η,s+η] ) ≤ C η ε -2 k ψ k 2 L 2 (C [s-η,s+η] ) .
(4.28)

Now combining (4.26)-(4.27)-(4.28) gives

ψ k 2 L 2 (C [s-η,s+η] ) ≥ Cε 2 k e -2 ε k (d c A (s)+2δ) , k ≥ k 0 (δ, η).
Finally noticing that ε 2 k ≥ e -δ ε k for k ≥ k 0 (δ) ends the proof of Proposition 4.9 up to replacing δ by δ/3. 4.7 Minimal time for uniform controllability in the limit ε → 0 +

The main purpose of this section is the proof Theorem 1.8 and its corollaries.

Proof of Theorem 1.8

Recall that we consider the following situation:

• S is a surface of revolution as described in Section 4.1, c > 0 is a fixed constant. Moreover, Assumption (1.17) is fulfilled.

• For k ∈ N, ε = ε k = ck -1 and ψ k is the set of solutions of P ε k ψ k = µ k ψ k defined in Lemma 4.3 associated to s 0 = s min is the minimum of V c .

• the function

v k (t, x) = e - µ k ε k t ψ k (x) (4.29)
is the solution to (2.9), namely

(ε k ∂ t + P ε k )v k = 0, v k | ∂S = 0, and v k | t=0 = ψ k .
Here, ψ k denotes the eigenfunction constructed in Section 4.3 above (in particular ψ k L 2 (S) = 1) and studied in Sections 4.5-4.6.

We now want to test Inequality (2.10) on v k , and thus estimate both sides of this inequality. This is achieved in Lemmata 4.12 and 4.13. Theorem 1.8 is then a direct consequence of these two lemmata. We recall that V c is defined in (1.16) and W c in (1.20).

Lemma 4.12. For any δ > 0, there exist C, k 0 > 0 such that for all k ≥ k 0 and all 0 ≤ T 0 ≤ δ -1 , we have

e -f 2ε v k (T 0 ) L 2 (S) ≥ Ce - Vc(s min )T 0 +W c m +δ ε k , W c m = min IL W c , with v k defined in (4.29) and W c (s) = d c A (s) + f(s) 2 .
Proof. Note first that the function W c = d c A + f 2 is continuous on I L and converges to +∞ close to the poles p ∈ P, according to the asymptotics of d c A in Lemma 4.5. Hence, it reaches its minimum in (at least) one point that we denote s 1 ∈ I L , that is W c (s 1 ) = W c m = min IL W c . We take 0 < η < dist(s 1 , P) small enough so that |f(s)f(s 1 )| ≤ δ for |ss 1 | ≤ η. We have for k ≥ k 0 large enough

e -f 2ε v k (T 0 ) L 2 (S) = e - µ k ε k T0 e -f 2ε ψ k L 2 (S) ≥ e - µ k ε k T0 e -f 2ε ψ k L 2 (C (s 1 -η,s 1 +η) ) ≥ e - µ k ε k T0 e - f(s 1 )+δ 2ε k ψ k L 2 (C (s 1 -η,s 1 +η) ) ≥ Ce - µ k ε k T0 e - f(s 1 )+δ 2ε k e -1 ε k (d c A (s1)+δ) ,
where we have used Proposition 4.9 for the last estimate. Since µ k → V c (s min ), we have µ k T 0 ≤ V c (s min )T 0 + δ for k large enough, which, together with the above estimate, concludes the proof (up to changing 3δ into δ).

Lemma 4.13. For any ω ⊂ S and δ > 0, there exist C, k 0 > 0 such that for all k ≥ k 0 , and all 0 ≤ T 0 ≤ δ -1 , we have

T0 0 e -f 2ε v k (t, •) 2 L 2 (ω) dt ≤ Ce -2 W c ω -δ ε k , W c ω = min ω W c , T0 0 e -f 2ε ∂ s v k (t)| s=L 2 H 1 (S 1 ) dt ≤ Ce -2 W c (L)-δ ε k
, in Cases 2 and 3,

T0 0 e -f 2ε ∂ s v k (t)| s=0 2 H 1 (S 1 ) dt ≤ Ce -2 W c (L)-δ ε k
, in Case 3, with v k defined in (4.29).

Proof. Since we need an upper bound, we can assume without loss of generality that ω is invariant by rotation. Also, W c ω is finite except in the trivial case ω ⊂ P. Let δ > 0. We first estimate the contribution close to the poles in case P = ∅. There, the function ψ k (hence v k ) is supposed to be very small since d c A is large. More precisely, using the asymptotics of d c A close to P given by Lemma 4.5, there exists η > 0 so that

d c A (s) 2 ≥ W c ω + f L ∞ 2
+ 1; for all s ∈ N η = {s ∈ I L such that dist(s, P) < η}.

(4.30)

We start with the estimate

T0 0 e -f 2ε v k 2 L 2 (ω∩CN η ) dt ≤ e f L ∞ ε ψ k 2 L 2 (CN η ) T0 0 e -2 µ k ε k t dt. (4.31)
We have

T0 0 e -2 µ k ε k t dt = ε k 2µ k T0 2µ k ε k 0 e -s ds ≤ C ε k 2µ k ≤ 1 for k large.
Corollary 4.7 applied to the set C Nη and the constant δ = 1/2, together with (4.30) implies

ψ k L 2 (CN η ) ≤ Ce -1 ε k (W c ω + f L ∞ 2
) .

With (4.31), this gives

T0 0 e -f 2ε v k 2 L 2 (ω∩CN η ) dt ≤ Ce f L ∞ ε e -2 ε k (W c ω + f L ∞ 2 ) = Ce -2 ε k W c ω ,
which is the expected bound for this part.

Let us now treat the contribution of the norm away from the poles (which is the whole I L in Case 3). Since f is uniformly continuous on [0, L] and d c A is uniformly continuous on the compact set N c η := I L \ N η , there exists η > 0 so that

s, s ′ ∈ N c η and |s -s ′ | ≤ η =⇒ |f(s) -f(s ′ )| ≤ δ and |d c A (s) -d c A (s ′ )| ≤ δ. (4.32)
We now select a finite sequence s i ∈ N c η , i = 1, • • • , N so that {|ss i | ≤ η} is a finite covering of N c η . This property gives the estimate

T0 0 e -f 2ε v k 2 L 2 (ω∩C N c η ) dt ≤ T0 0 e -2 µ k ε k t dt i∈J e -f 2ε ψ k 2 L 2 (C (s i -η,s i +η) )
where

J = {i = 1, • • • , N ; (s i -η, s i + η) ∩ ω = ∅}. Using |f (s) -f (s i )| ≤ δ and |d c A (s) -d c A (s i )| ≤ δ for s ∈ (s i -η, s i + η) and then Corollary 4.7 with some 0 < δ ′ ≤ δ/(1 + d c A (s i )) so that (1 -δ ′ )d c A (s i ) -δ ′ ≥ d c A (s i ) -δ, we obtain e -f 2ε ψ k L 2 (C (s i -η,s i +η) ) ≤ e -f (s i )+δ 2ε ψ k L 2 (C (s i -η,s i +η) ) ≤ Ce -f (s i )+δ 2ε e -1 ε k ((1-δ ′ )d c A (si)-δ ′ ) e δ ε ≤ Ce - W c (s i )-4δ ε k
.

We finally obtain

T0 0 e -f 2ε v k 2 L 2 (ω∩C N c η ) dt ≤ Ce 8δ ε k max i∈J e -2 W c (s i ) ε k
.

We remark from (4.32) that min i∈J W c (s i ) ≥ W c ω -2δ. This finishes the proof since δ is arbitrary and k can be chosen large enough.

Finally, the proof of the boundary estimates simply consists in replacing the use of Corollary 4.7 by that of the Inequality (4.18).

We may now conclude the proof of Theorem 1.8 from Lemmata 4.12 and 4.13 Proof of Theorem 1.8. Using Lemma 2.9, if uniform observability holds for T 0 , then, we have the inequality

e -f 2ε v(T 0 ) 2 L 2 (S) ≤ C 2 0 T0 0 e -f 2ε v 2 L 2 (ω)
dt.

for any solution v of (2.9). In particular, this inequality holds true for the sequence v k described above. So, combining Lemmata 4.12 and 4.13, we obtain

e - Vc (s min )T 0 +W c m +δ ε k ≤ C 0 e - W c ω -δ ε k , k ≥ k 0 (δ), for all k ≥ k 0 .
This implies V c (s min )T 0 ≥ W c ω -W c m -2δ when letting k → +∞, which gives the expected result since δ is arbitrary. The proof of the boundary observability estimate (1.22) follows the same. We also notice that all the previous results also apply in the case with the alternative Definition 1.19 of the Agmon distance which is also Lipschitz with the same properties that we used.

Proof of Corollaries of Theorem 1.8

Corollaries 1.9, 1.10 and 1.11, stated in the introduction, are significant examples of application of Theorem 1.8. In this section, we prove these three results.

Proof of Corollary 1.9. We first consider the case S diffeomorphic to S 2 , i.e. Case 1. The case S diffeomorphic to D, i.e. Case 2, is discussed at the end of the proof. We define

f δ (s) = s 0 χ δ (t)dt where χ δ ∈ C ∞ c ((0, L); [0, 1]), χ δ (s) = 1 in a neighborhood of [δ, L -δ]
. Such a function f δ is constant near 0 and L, and hence can be extended by continuity as a C ∞ function on S (see e.g. [Bes78, Proposition 4.6]). We notice that f ′ δ (s) = 1 for s ∈ [δ, Lδ], so the statements about T GCC are direct consequences of Proposition 4.2.

Notice now that we have

V c (s) := c 2 R(s) 2 + |χ δ (s)| 2 4 . Let us call s min ∈ (0, L) the unique point such that R(s min ) = max R, that is 1 R(smin) 2 = min 1 R 2 .
Claim: For all c > 0, there is δ 0 > 0 such that for all δ ∈ (0, δ 0 ), we have V c (s min ) = min V c and V -1 c (V c (s min )) = {s min }.

To prove the claim, we let δ 0 > 0 be such that

s min ∈ [δ 0 , L -δ 0 ], and
c 2 R 2 (s) > c 2 R(s min ) 2 + 1 4 for s / ∈ [δ 0 , L -δ 0 ]
(note that δ 0 thus depends on c). This is possible since R(s) → 0 as s → 0 + and s → L -. Hence, recalling the definitions of χ δ and V c , for δ < δ 0 we have

V c (s min ) = c 2 R(s min ) 2 + 1 4 , V c (s) > c 2 R(s min ) 2 + 1 4 for all s ∈ [δ 0 , L -δ 0 ] \ {s min }, together with V c (s) ≥ c 2 R 2 (s) > c 2 R(s min ) 2 + 1 4 for s / ∈ [δ 0 , L -δ 0 ].
As a consequence, V c reaches its minimum at s min only, which proves the claim. Now Assumption 4.4 is satisfied and s min does not depend on δ. We compute the Agmon distance (1.18)

d c A (s) = s smin V c (y) -V c (s min )dy = s smin c 2 R(y) 2 + |χ δ (y)| 2 4 - c 2 R(s min ) 2 - 1 4
dy .

37

Notice that 0 ≥ |χ δ (y)| 2

4

-1 4 ≥ -1 4 uniformly with respect to δ so that the asymptotic expansion of Lemma 4.5 is valid uniformly in δ: there are C, γ > 0 such that for all δ ∈ (0, δ 0 ),

|d c A (s) + c log(s)| ≤ C for s ∈ (0, γ], |d c A (s) + c log(L -s)| ≤ C, for s ∈ [L -γ, L). Now we recall the definition of W c = d c A + f δ 2 , notice that 0 ≤ f ′ δ ≤ 1 so that 0 ≤ f δ (s) ≤ L for s ∈ I L .
As a consequence, using that d c A ≥ 0 and d c A (s min ) = 0, we have for δ ∈ (0, γ],

W c m = min IL W c = min IL d c A + f δ 2 ≤ L 2 + min IL d c A ≤ L 2 W c ω = min ω W c ≥ min ω d c A ≥ -c log(δ) -C for ω = B g (N, δ) ∪ B g (S, δ).
The bound (1.21) of Theorem 1.8 then yields

c 2 R(s min ) 2 + 1 4 T unif (ω) ≥ -c log(δ) -C - L 2 ,
and hence concludes the proof in the case S diffeomorphic to S 2 , i.e. Case 1.

In the case S diffeomorphic to D, i.e. Case 2, we instead define

f δ (s) = s 0 χ δ (t)dt where χ δ ∈ C ∞ c ((0, L]; [0, 1]), χ δ (s) = 1 in a neighborhood of [δ, L].
Then, the remainder of the proof is the same except that the minimum can be achieved at s = L, and all sets of the form [δ 0 , Lδ 0 ] have to be replaced by [δ 0 , L] (i.e., only a neighborhood of zero is avoided, and not a neighborhood of L).

Proof of Corollary 1.10. For the sake of simplicity, we may identify S 1 L = [-L/2, L/2] and I ω = (-α, α) for α ∈ (0, L/2). We first choose χ δ ∈ C ∞ (S 1 L ) even in this identification, and

s min = ±L/2 / ∈ (-α, α) such that χ δ = 1 δ on - 1 2 α + L 2 , 1 2 α + L 2 , χ δ (s min ) = 1, χ ′ δ (s) = 0 if and only if s ∈ - 1 2 α + L 2 , 1 2 α + L 2 or s = s min .
Note that χ δ = 1 δ in a neighborhood of (-α, α), and that for δ < 1, χ δ reaches at s min = ±L/2 a unique global minimum and in particular χ δ ≥ 1 on S 1 L . We then set

V δ (s) := χ δ + M, M := max S 1 L |f ′ | 2 4 , so that V δ (s) > |f ′ (s)| 2 4 on S 1 L , and
R δ (s) := V δ (s) - |f ′ (s)| 2 4 -1 2 , R δ ∈ C ∞ (S 1 L ; R + * ),
where we have used that f ∈ C ∞ (S 1 L ). Notice that with these definitions, we have (where

V 1 denotes V c with c = 1) V 1 (s) = 1 R δ (s) 2 + |f ′ (s)| 2 4 = V δ (s).
Notice then that V δ admits a unique global minimum at s min = ±L/2. On the one hand, we have

(χ δ (s) + M ) -1 2 ≤ R δ (s) ≤ χ -1 2 δ (
s). On the other hand, since χ δ is even and with the appropriate definition (1.19), we have for s ∈ [0, L/2],

d c A (s) = - s L/2 V δ (y) -V δ (L/2)dy = - s L/2 χ δ (y) -1dy.
As a consequence, we have (recall ω = (-α, α) × S 1 )

W c m = min [0,L] W c = min [0,L] d c A + f 2 ≤ d c A + f 2 (L/2) = f(L/2) 2 , W c ω ≥ min ω d c A + min ω f 2 = d c A (α) + min ω f 2 ,
where

d c A (α) = L/2 α χ δ (y) -1dy ≥ 1 2 (α+ L 2 ) α χ δ (y) -1dy = 1 2 (α+ L 2 ) α 1 δ -1dy ≥ 1 2 ( L 2 -α)(δ -1/2 -1)
for δ < 1. Applying Theorem 1.8 concludes the proof of the corollary.

We next prove Corollary 1.11, stated in the introduction. The proof is close to that of Corollary 1.10.

Proof of Corollary 1.11. We first take χ ∈ C ∞ (R + ; [0, 1]) such that χ = 1 on [0, L/4] and supp χ ⊂ [0, L/2).

Next, with

M := max [0,L] |f ′ | 2 4 , we set, V δ (s) := χ(s) (s + δ) γ + (1 -χ(s)) s - L 2 2 + M, for s ∈ [0, L 2 ], V δ (s) := V δ (L -s) for s ∈ [ L 2 , L].
This function is symmetric about L/2, smooth on [0, L], and satisfies V δ (s)

> |f ′ (s)| 2 4 on [0, L]. Hence, defining R δ as R δ (s) := V δ (s) - |f ′ (s)| 2 4 -1 2 , and using that f ∈ C ∞ ([0, L]), we deduce that R δ ∈ C ∞ ([0, L]; R + * ).
Notice that with these definitions, we have (where

V 1 denotes V c with c = 1) V 1 (s) = 1 R δ (s) 2 + |f ′ (s)| 2 4 = V δ (s).
Notice then that V δ admits a unique global minimum at s min = L/2.

We have on the one hand that for all s ∈ [0, L/4] and δ ∈ (0, δ 0 ],

(s + δ)

γ 2 1 1 + M (s + δ) γ = V δ (s) -1 2 ≤ R δ (s) ≤ V δ (s) -M -1 2 = (s + δ) γ 2 ,
which proves the Item 4 (on account to the symmetry V δ (s) := V δ (Ls)).

On the other hand, recalling that s min = L/2, we have for s

∈ [0, L/2] V 1 (s min ) = V δ (L/2) = M, d c A (s) = s L/2 V δ (y) -V δ (s min )dy = s L/2 χ(s) (s + δ) γ + (1 -χ(s)) s - L 2 2 dy .
As a consequence, we have

W c m = min [0,L] W c = min [0,L] d c A + f 2 ≤ d c A + f 2 (L/2) = f(L/2) 2 , W c {0} = W c (0) ≥ d c A (0) + f(0) 2 ,
where

d c A (0) = L/2 0 χ(s) (s + δ) γ + (1 -χ(s)) s - L 2 2 dy ≥ L/4 0 1 (s + δ) γ dy = δ 1-γ/2 γ/2 -1 - (δ + L/4) 1-γ/2 γ/2 -1 = δ 1-γ/2 γ/2 -1 + O (1) , for γ > 2.
By symmetry of V δ about L/2, we also have d c A (L) = d c A (0). Applying Theorem 1.8 concludes the proof of the corollary.

Uniform time of observability for positive solutions

The proofs of Theorem 1.3 and Proposition 1.12 rely on fine estimates on the semiclassical heat kernel, which we borrow from [START_REF] Li | On the parabolic kernel of the Schrödinger operator[END_REF]. The latter are first presented in Section 5.1. Then, in Section 5.2, we deduce L 1 observability statements and finally conclude the proofs of Theorem 1.3 and Proposition 1.12 in Section 5.3. Throughout this section, we assume ∂M = ∅.

Estimates on the semiclassical heat kernel

The main tool we use to estimate the heat kernel in the semiclassical limit is the following theorem taken from Li-Yau [START_REF] Li | On the parabolic kernel of the Schrödinger operator[END_REF] (see also [START_REF] Simon | Instantons, double wells and large deviations[END_REF], for a similar result on R n ).

Theorem 5.1 (Theorem 6.1 of [START_REF] Li | On the parabolic kernel of the Schrödinger operator[END_REF]). Let M be a compact manifold without boundary. Suppose V ε = V + εq with V, q ∈ C 2 (M). For any ε > 0, we consider H ε , the fundamental solution of

∂ t w -∆ g w + 1 ε 2 V ε (x)w = 0, on (0, +∞) × M.
Then, we have

lim ε→0 ε log H ε (x, y, εt) = -ρ(x, y, t) (5.1) where ρ(x, y, t) = inf t 0 1 4 | γ(s)| 2 g + V γ(s) ds, γ ∈ W 1,∞ ([0, t]; M), γ(0) = x, γ(t) = y .
(5.2)

Moreover, the limit in (5.1) is uniform on any compact set of M 2 × (0, +∞).

We recall that H ε (x, y, t) is defined to be the unique solution to

   ∂ t -∆ g + 1 ε 2 V ε (x) H ε (x, y, t) = 0, for (t, x) ∈ R + * × M, H ε (x, y, t)| t=0 = δ x=y , for x ∈ M, (5.3)
where y ∈ M is fixed, and the differential operator -∆ g + 1 ε 2 V ε (x) acts in the x-variable. We also recall that H ε (x, y, t) is (well-defined and) continuous in M 2 × (0, +∞) as soon as V ε ∈ L ∞ (M), see e.g. [START_REF] Simon | Schrödinger semigroups[END_REF]Theorem B.7.1 (a ′′′ )], so that pointwise estimates like (5.1) make sense.

The statement of Theorem 5.1 is not strictly speaking a consequence of [LY86, Theorem 6.1] for the following two reasons: with d ∇g f (x, y, t) := ρ(x, y, t) + f(x)f(y) 2 , (5.5)

where ρ(x, y, t) is defined by (5.2) with V (x) = |∇g f| 2 g 4 . Moreover, the limit is uniform on any compact set of M 2 × (0, +∞).

In particular, for any δ > 0 and any compact subset I ⋐ (0, +∞), there exists ε 0 > 0 such that e - for any (x, y, t) ∈ M 2 × I and 0 < ε ≤ ε 0 .

Note that the definition of d ∇g f in (5.5) is not the same as that given in (1.24) in the introduction. Equivalence between these two definitions is proved in Lemma A.2.

Note that although the kernel H ε (x, y, t) is symmetric (with respect to the Riemannian volume measure d Vol g ) since -∆ g + 1 ε 2 V ε is, this is no longer the case for the kernel K ε (x, y, t) (since the operator in (1.6) is not symmetric in L 2 (M, d Vol g )). Similarly, ρ(x, y, t) is symmetric whereas d ∇g f (x, y, t) is not.

Proof. Setting u y (t, x) = K ε (x, y, t) and w y (t, x) = e f(x)/2ε u y (t/ε, x), we have from Equation (5.4) and Lemma 2.9 that w y (t, x) solves (2.12). Moreover, we have w y (0, x) = e f(x)/2ε u y (0, x) = e f(x)/2ε δ x=y = e f(y)/2ε δ x=y .

This implies that w y (t, x) = e f(y)/2ε H ε (x, y, t) where H ε (x, y, t) is defined in (5.3) with V ε = |∇g f| 2 g 4 + ε ∆gf 2q . Finally, we have proved that H ε (x, y, εt) = e f(x)/2ε e -f(y)/2ε K ε (x, y, t),

L 1 observability estimates for positive solutions

We first prove intermediate observability statements in L 1 . The following elementary abstract lemma shows that concerning positive solutions, observability in a (possibly weighted) L 1 norm is equivalent to the "observability of the Kernel".

Proposition 5.3. Let M be a compact Borel space (on which we denote by dx a distinguished measure) and T > 0. Take K = K(x, y, t) ∈ C 0 (M × M × (0, T ]) be a nonnegative kernel. Assume further that (t, y) → K(•, y, t) L 1 (M) = M K(x, y, t)dx is uniformly bounded for (t, y) ∈ (0, T ] × M. Define for t ∈ (0, T ] the operator Note first that under the assumption of the theorem, both terms in (5.7) are well-defined. Indeed, the Tonelli theorem (all functions/measures are nonnegative) implies From (5.6), the fact that K ε (x, y, t) ≥ 0 for (x, y, t) ∈ M 2 × (0, ∞), and s > δ, we deduce O T (y) ≥ (5.10) Also, for any s > δ, using that d ∇g f ≥ 0 (see Proposition A.4), we have

I s (y) ≤ Vol g (M)e δ ε .
When combined with (5.10), we obtain O T (y) ≥ Ce -3 δ ε I s (y). By Proposition 5.3, this gives (5.8) which concludes the proof of the proposition (up to changing 4δ into δ).

5.3

From L 1 to L 2 observability estimates for positive solutions

In this section, we conclude the proofs of Proposition 1.12 and Theorem 1.3. We first prove the negative result of Proposition 1.12 (uniform observability of positive solutions does not hold for T < T GCC (M, ∇ g f, ω), with an exponential lower bound of the cost).

Proof of Proposition 1.12. Let us check the first part of the proposition (geometric statement). Since (M, ∇ g f, ω, T ) does not satisfy (GCC), there is y 0 ∈ M, so that for all t ∈ [0, T ], φ -t (y 0 ) / ∈ ω. In particular, for any (t, x) ∈ [0, T ] × ω, we have φ t (x) = y 0 , which implies d ∇g f (x, y 0 , t) > 0 by Proposition A.4. By compactness of [0, T ] × ω, inf x∈ω,t∈[0,T ] d ∇g f (x, y 0 , t) > 0. Therefore, d ([0,T ],ω) > 0 as expected.

For the second part, for any δ > 0, select y 0 ∈ M so that inf By uniform continuity of d ∇g f (x, y 0 , t) defined on [0, 2T ]×ω, we can also find η > 0 so that inf (t,x)∈[0,T +η]×ω d ∇g f (x, y 0 , t) ≥ d ([0,T ],ω)δ/4 and d ∇g f (x, y 0 , η) ≤ δ/4 for x ∈ B(y 0 , η). We take as initial datum u 0 (x) = K ε (x, y 0 , η), yielding u(t, x) = K ε (x, y 0 , t + η) (see the definition of K ε in (5.4)) as the associated solution of (1.6). We have u 0 ∈ L 2 (M) together with a lower bound coming from (5.6) with a sufficiently small δ (replaced by δ/4) u 0 L 2 (M) ≥ u 0 L 2 (B(y0,η)) ≥ C(η)e -δ 2ε .

Concerning the observation term, we deduce from the upper bound in (5.6) that for ε small enough, To conclude the proof of Theorem 1.3, we need the following dissipation result taken from Guerrero-Lebeau [START_REF] Guerrero | Singular optimal control for a transport-diffusion equation[END_REF]. In that reference, it is written on an open subset Ω ⊂ R n with the flat metric; however, it can be checked that the result also applies to the case of a Riemannian manifold (M, g) without boundary and with an additional potential q. Lemma 5.9. Let V be a closed convex set of L 2 (M) with 0 ∈ V and V ⊂ V so that V ⊂ Vv for any v ∈ V. Assume moreover that there exists C V > 0 so that

C V (T, ε) 2 T 0 ω |u(t)| 2 ds(x)dt ≥ u(T ) 2 L 2 (M) ,
for all u 0 ∈ V and u solution of (1.3). (5.16)

Then, for any y 0 ∈ L 2 (M) and 0 < ε ≤ ε 0 , there exists a control h ∈ L 2 ([0, T ], L 2 (M)) with

h L 2 ([0,T ],L 2 (M)) ≤ C V (T, ε) y 0 L 2 (M)
so that the solution of (2.1) satisfies (y(T ), u 0 ) L 2 (M) ≥ 0 for any u 0 ∈ V.

For the proof of Corollary 2.4, we apply this lemma to the sets V = V = L 2 (M; R + ). Notice that Lemma 5.9 also contains one implication (namely Observability =⇒ Controllability) in Corollary 2.3 when applied to V = V = L 2 (M; R).

Proof. For any α > 0, we consider the functional J α defined for any u 0 ∈ V by

J α (u 0 ) = 1 2 T 0 ω |u(t, x)| 2 dtdx + α 2 u 0 2 L 2 (M) + (u(T ), y 0 ) L 2 (M) ,
where u is the solution of (1.3). The functional J α is continuous, convex and coercive. Therefore, J α admits a minimum u 0,α ∈ V (se e.g. [ET74, Chapter II, Proposition 1.2]). The minimality condition gives (see e.g. [ET74, Chapter II, Proposition 2.1]) for any p 0 ∈ L 2 (M) that can be written p 0 = v 0u 0,α , with v 0 ∈ V, (5.17)

where we have denoted p (resp. u α ) the solution of (1.3) with p(0) = p 0 (resp. u(0) = u 0,α ). Now, let y α be the solution of (2.1) with control function h α (t, x) = u α (Tt, x) and initial datum y α (0) = y 0 . The duality equation (2.2) gives for any p 0 ∈ L 2 (M) initial datum for p solution of (1.3) Combined with (5.17), this implies (p 0 , y α (T )) L 2 (M) + α (p 0 , u 0,α ) L 2 (M) ≥ 0, (5.18) for every p 0 ∈ Vu 0,α . This also holds for any p 0 ∈ V since V ⊂ Vu 0,α by assumption.

To obtain an estimate of the control, we apply (5.17) to p 0 = 0u 0,α . After an application of the Cauchy-Schwarz inequality, we have The observability inequality (5.16) applies to u 0,α ∈ V, so that

1

  Introduction and main results 1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1.2 Background and motivation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1.3 Main results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1.4 Further remarks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 Preliminaries: duality, conjugation of gradient flows and (GCC) 2.1 Uniform controllability problems and dual formulation . . . . . . . . . . . . . . . . . . . . . 2.2 The vanishing viscosity limit for gradient flows. Conjugation and reformulation . . . . . . . 2.3 (GCC) and controllability of the limit equation ε = 0 . . . . . . . . . . . . . . . . . . . . . . 3 General lower bounds without geometric assumption 3.1 Eigenfunctions of semiclassical Schrödinger operators . . . . . . . . . . . . . . . . . . . . . . 3.2 Rough localization of eigenfunctions, and a proof of Theorem 3.1 . . . . . . . . . . . . . . . 3.3 An explicit counter-example for a domain of R n . . . . . . . . . . . . . . . . . . . . . . . . .

  controlling the solution e -C 0 ε close to zero costs ≈ 1, whereas controlling the solution exactly to zero costs ≈ e C ε .

  a slight abuse of notation, we shall still denote by P (k) ε the one dimensional operator (ζ -1 ) * P (k) ε ζ * defined on I L , namely

  Lemma 4.5 (Properties of d c A ). Under Assumption 4.4, we have d c A ∈ C 2 (I L \ {s min }) together with d c A (s) = -c log(s) + O (1) , as s → 0 + , in Cases 1 and 2, (4.16)

  O (y))dy = -c log(s) + O (1) as s → 0 + .

3

  . the energy at s minρ is dominated by the energy near s via Lemma 4.10 (with a geometric constant e energy near s is dominated by the L 2 norm of ϕ k (or ψ k ) near s via elliptic regularity.Proof of Proposition 4.9. Without loss of generality, we can assume s < s minη. Indeed, the case s > s min + η is treated similarly, and the case s ∈ [s minη, s min + η] is a direct consequence of Corollary 4.8 applied for ρ = η/2. Now, since s min ∈ I L , notice that we have necessarily d := dist(s min , P) > 0, and we may also assume that δ < d/4. Lemma 4.11 can be applied with α = d/4 and produces some constants C = C α , D = D α . Let us now choose ρ = ρ δ = min(δ/(4D), d/4, η/2).

A

  (s)+2δ) , where C is a new constant depending only on δ and η. Note that in the last inequality, we have used that d c A (s) ≥ d c A (y). Integrating for y ∈ [s, s + η/2], we get [s,s+η/2] E k (y)dy ≥ Ce -

d

  ∇g f (x,y,t)+δ ε ≤ K ε (x, y, t) ≤ e - d ∇g f (x,y,t)-δ ε (5.6)

ε

  log K ε (x, y, t) = f(y)f(x) 2 + lim ε→0 ε log H ε (x, y, εt) = f(y)f(x) 2 ρ(x, y, t),after having applied (5.1).

S

  (t) : Meas(M) → C 0 (M), [S(t)µ] (x) = M K(x, y, t)dµ(y) = µ, K(x, •, t) ,where Meas(M) denotes the space of Radon measures on M. Let w 1 ∈ L ∞ ([0, T ] × ω) and w 2 ∈ L ∞ (M) be two nonnegative weight functions.Then, for all T, s, C 0 > 0, the following statements are equivalent 1. Observability of positive measures:w 2 [S(s)µ] L 1 (M) ≤ C 0 w 1 [S(•)µ] L 1 ((0,T ]×ω) (5.7)for any µ ∈ Meas + (M) nonnegative Radon measure.2. Observability of positive L 1 functions:w 2 [S(s)u 0 ] L 1 (M) ≤ C 0 w 1 [S(•)u 0 ] L 1 ([0,T ]×ω)for any u 0 ∈ L 1 (M, dx) with nonnegative value.3. Observability of Dirac distributions:w 2 [S(s)δ y ] L 1 (M) ≤ C 0 w 1 [S(•)δ y ] L 1 ([0,T ]×ω)for any y ∈ M. 4. Observability of the Kernel: I s (y) ≤ C 0 O T (y), for all y ∈ M, where O T (y) = T 0 ω w 1 (t, y)K(x, y, t)dxdt, and I s (y) = M w 2 (x)K(x, y, s)dx.

  w 1 S(•)µ L 1 ([0,T ]×ω) = T 0 ωw 1 (t, x)[S(t)µ](x)dxdt = t, x)K(x, y, t)dx dµ(y)dt≤ T w 1 L ∞ ([0,T ]×M) µ T V sup t∈[0,T ] sup y∈M ω K(x, y, t)dx < +∞,by assumption. Here µ T V denotes the total variation of the measure µ. Note also that u(t) = [S(t)(µ)] is a continuous nonnegative function for any nonnegative measure µ and t > 0. Remark that O T (y) is essentially the observation of solutions starting from δ y while I s (y) is the weighted norm of this solution at time s.Proof. Remark first thatO T (y) = T 0 ω w 1 (t, y)K(x, y, t)dxdt = T 0 ω w 1 (t, y) δ y , K(x, •, t) dxdt = w 1 [S(•)δ y ] L 1 ([0,T ]×ω) I s (y) = M w 2 (x)K(x, y, s)dx = M w 2 (x) δ y , K(x, •, s) dx = w 2 [S(s)δ y ] L 1 (M) ,so that Item 3 ⇔ Item 4. Moreover, applying the Fubini Theorem, we getw 1 [S(•)µ] L 1 ([0,T ]×ω) = T 0 x∈ω y∈M w 1 (t, y)K(x, y, t)dµ(y)dxdt = µ, O T w 2 [S(s)µ] L 1 (M) =y∈M x∈M w 1 (t, x)K(x, y, t)dµ(y)dxdt = µ, I s . Therefore, Item 1 is equivalent to µ, C 0 O T -I s ≥ 0 for any µ ∈ Meas + (M) and Item 2 to M u 0 (y)(C 0 O T -I s )(y)dy ≥ 0 for any u 0 ∈ L 1 + (M, dx). That Item 4 ⇔ Item 1 ⇔ Item 2 follows from the general fact that if f ∈ C 0 (M), one has f ≥ 0 on M ⇔ µ, f ≥ 0 for all µ ∈ Meas + (M) ⇔ M u 0 (y)f (y)dy ≥ 0 for all u 0 ∈ L 1 + (M, dx).

d

  ∇g f (x,y,t)+δ ε and I s (y) ≤ x∈M e - d ∇g f (x,y,s)-δ ε ,where δ ′ > 0 is chosen sufficiently small so that M, ∇ g f, ω, (δ ′ , T ) still satisfies (GCC), which is possible since T > T GCC (M, ∇ g f, ω). Using Lemma 5.6 (where U y and C δ are defined), we now haveO T (y) ≥ (x,t)∈Uy e - d ∇g f (x,y,t)+δ ε dtdx ≥ (x,t)∈Uy e -2δε dtdx ≥ C δ e -2δ ε .

(

  t,x)∈[0,T ]×ω d ∇g f (x, y 0 , t) ≤ d ([0,T ],ω) + δ/8.

K≤.

  ε (x, y 0 , t + η) 2 dxdt ≤ Vol g (ω)T sup x∈ω,t∈[η,T +η] e -2 d ∇g f (x,y 0 ,t)-δ/4 ε Ce -2 d ([0,T ],ω) -δ/2 εApplying the observability inequality (1.23) to u thus implies thatC + 0 (T, ε) 2 e -2 d ([0,T ],ω) -δ/4 ε ≥ C,uniformly in ε ∈ (0, ε 0 ], which concludes the proof of the proposition.

  x)u α (t, x)dtdx + α (p 0 , u 0,α ) L 2 (M) + (p(T ), y 0 ) L 2 (M) ≥ 0.

  , x)u α (t, x)dtdx = (p 0 , y α (T )) L 2 (M) -(p(T ), y 0 ) L 2 (M) .

  |u α (t, x)| 2 dtdx + α u 0,α 2 L 2 (M) ≤ u α (T ) L 2 (M) y 0 L 2 (M) .

  |u α (t, x)| 2 dtdx + α u 0,α 2 L 2 (M) ≤ C V (T, ε) T 0 ω |u α (t, x)| 2 dtdx 1/2 y 0 L 2 (M) .

See e.g. [Lee13, Chapter 16 p427]: given a local chart (Uκ, κ) of M, we have Uκ u ds = κ(Uκ) u • κ -1 (y)ϕ κ (y)dy for an appropriate smooth positive function ϕ κ , and for any u ∈ C 0 c (Uκ).

• the potential V ε is assumed independent of ε in [START_REF] Li | On the parabolic kernel of the Schrödinger operator[END_REF];

• the uniformity of the limit on any compact subset is not explicitly written in [START_REF] Li | On the parabolic kernel of the Schrödinger operator[END_REF].

However, let us explain why the proof of [LY86, Theorem 6.1] actually contains these two points. The limit (5.1) is proved in two steps, a lower bound and an upper bound. The lower bound lim ε→0 ε log H ε (x, y, εt) ≥ -ρ(x, y, t) is obtained as a consequence of the explicit estimate ([LY86, Theorem 2.1]) ε log H ε (x, y, εt 1 ) ≤ ε log H ε (x, y, εt) + ε 2 A 1/ε (tt 1 ) + ρ α,R (x, y, tt 1 ), for all 0 < t 1 < t, α > 1 and where the constant A λ (λ = ε -1 in our context) only depends on ∆ g V ε L ∞ (M) , ∇ g V ε L ∞ (M) and α. In particular, the limit lim ε→0 εA 1/ε = 0 holds uniformly. The proof finally proceeds by taking the limit R → +∞, and then α → 1 + (in which case ρ Vε α,R (x, y, tt 1 )) → ρ Vε (x, y, t-t 1 ), uniformly on compact sets; here ρ Vε denotes the function ρ defined by (5.2) with V replaced by V ε ). Taking finally the limit ε → 0 + , and then t 1 → 0 + , and noticing that lim t1→0 + lim ε→0 + ε log H ε (x, y, εt 1 ) ≥ 0 concludes the proof of the lower bound. This last argument relies only on a uniform upper bound on V ε L ∞ (M) together with a comparison argument, namely Inequality [START_REF] Li | On the parabolic kernel of the Schrödinger operator[END_REF](6.3)]. It can be checked that the convergence is uniform on any compact set since it involves the asymptotics on the diagonal of kernels of heat equations on large balls and without potential, which are known to be uniform on compact sets, see for instance [START_REF] Srinivasa | Diffusion processes in a small time interval[END_REF]Theorem 4.6]. It only remains to notice that the limit lim t1→0 + lim ε→0 + ρ Vε (x, y, tt 1 ) = ρ(x, y, t) is uniform on compact sets.

The upper bound lim ε→0 ε log H ε (x, y, εt) ≤ -ρ(x, y, t)

follows from [LY86, Theorem 3.3]. As for the lower bound, this result also furnishes an explicit and uniform bound involving another constant A λ enjoying the same type of convergence properties as for the lower bound.

Note that we have chosen to use the estimate of the semiclassical limit of the Kernel of [START_REF] Li | On the parabolic kernel of the Schrödinger operator[END_REF] but it is likely that we could have obtained the observability inequality starting directly from Harnack type inequalities like [LY86, Theorem 2.1], as we did in [START_REF] Laurent | Observability of the heat equation, geometric constants in control theory, and a conjecture of Luc Miller[END_REF].

Here, we are mostly interested in the case

In this situation, we can reformulate the result in terms of the transport equation with vanishing viscosity (1.6).

We next define K ε (x, y, t), the fundamental solution of (1.6) on M 2 × (0, +∞) by the unique solution to

where y ∈ M is fixed, and the differential operator -∇ g f • ∇ gq -ε∆ g acts in the x-variable. Recall that for u 0 ∈ L 1 (M), the function

is the unique solution of (1.6) on (0, +∞) × M issued from u| t=0 = u 0 (This uses the choice of volume form d Vol g in the embedding L 1 (M) ֒→ D ′ (M)).

Corollary 5.2. Let M be a compact manifold without boundary. Suppose X = ∇ g f where f is a C 3 function defined on M. For any ε > 0, we consider K ε , the fundamental solution of (1.6) on (0, +∞) × M. Then, we have

We now give the L 1 observability estimate for positive solutions to (1.6).

Proposition 5.4. Assume that T > T GCC (M, ∇ g f, ω). Then, for any δ, s > 0, there exists ε 0 > 0 so that we have

for all ε ∈ (0, ε 0 ), (5.8)

for any u 0 ∈ L 1 (M) with non-negative values and u solution of (1.6).

Remark that we only use the case s = T below. It is however remarkable that the stronger result for s > 0 small holds as well. This is linked to the L 1 setting here. Note also that in L 1 , we have a "converse inequality", which we state for the sake of the comparison. Proposition 5.4 is proved afterwards.

Lemma 5.5. Assume ∂M = ∅. For all T > 0, there is C T > 0 such that for all u 0 ∈ L 1 (M; R + ) and u the associated solution of (1.6), we have

for all t ∈ [0, T ] and ε > 0.

(5.9)

In particular, this implies that one cannot hope to replace the loss e δ ε by a gain e -δ ε in (5.8).

Proof of Lemma 5.5. Assume first that u 0 ∈ W 2,1 (M) with u 0 ≥ 0 a.e. on M. Then, notice that u(t, x) ≥ 0 for a.e. (t, x) ∈ (0, T ) × M. Integrating (1.6) on M, we obtain after an integration by parts (using that ∂M = ∅)

). The Grönwall inequality yields (5.9). The conclusion for a general u 0 ∈ L 1 follows from a density argument.

We now turn to the proof of the L 1 observability estimate of Proposition 5.4, which will use the following lemma.

Lemma 5.6. Assume M, ∇ g f, ω, (δ ′ , T ) satisfies (GCC) (see Definition 2.10) for some δ ′ > 0.Then, for all δ > 0 there is C δ > 0 such that for all y ∈ M there is an open set U y ⊂ ω × (δ ′ , T ) with |U y | ≥ C δ and for all (x, t) ∈ U y , we have d ∇g f (x, y, t) ≤ δ.

Proof of Lemma 5.6. The assumption (GCC) implies that for any y 0 ∈ M, there is t y0 ∈ (δ ′ , T ) and x y0 = φ -ty 0 (y 0 ) ∈ ω where (φ t ) t∈R is the flow of ∇ g f. The trajectory γ(s) = φ s (x y0 ) satisfies γ(s) = ∇ g f(s) with γ(0) = x y0 and γ(t y0 ) = y 0 so that Proposition A.4 implies d ∇g f (x y0 , y 0 , t y0 ) = 0. In particular, we obtain that for any y 0 ∈ M, there is t y0 ∈ (δ ′ , T ) and x y0 ∈ ω such that d ∇g f (x y0 , y 0 , t y0 ) = 0.

By uniform continuity of

By compactness, we can cover M by M = i∈I B g (y i , ν yi ) where I is finite. Then, for any y ∈ M, there is i ∈ I such that y ∈ B g (y i , ν yi ), and the set

satisfies the sought properties.

As a consequence of this lemma together with Corollary 5.2 and Proposition 5.3, we may now deduce a proof of the L 1 observability estimate of Proposition 5.4.

Proof of Proposition 5.4. Without any loss of generality, we can assume 0 < δ < s. According to Proposition 5.3, it is enough to study the "observability of the Kernel". Using Corollary 5.2, for any δ > 0, there exists ε 0 so that (5.6) holds for all (x, y)

Proposition 5.7 (Proposition 3 of [START_REF] Guerrero | Singular optimal control for a transport-diffusion equation[END_REF], m = 1). Assume that ∂M = ∅ and that (M, X, ω, T ) satisfies (GCC). Then, there exist C, C 0 > 0 such that

(5.11) for all ε ∈ (0, 1] and any (not necessarily positive) solution u to (1.3) (and a fortiori for all solutions u to (2.6)). The same statement holds true if ∂M = ∅ and (M, X, ω, T ) satisfies (FC).

We shall also need the following lemma in the proof of Theorem 1.3.

Lemma 5.8. For any δ ′ > δ > 0, there exists ε 0 > 0 so that

(5.12)

for any solution u to (2.6) and 0 < ε ≤ ε 0 .

Proof. Since the manifold is compact, we only need to prove the L ∞ bound which follows from the bound

on the Kernel. This estimate follows from Corollary 5.2 (e.g. (5.6) together with the fact that d ∇g f ≥ 0, see Proposition A.4).

Proof of Theorem 1.3. Inequality (1.25) of Proposition 1.12 directly yields T + unif (ω) ≥ T GCC (M, ∇ g f, ω). Note that it was mostly proved in Guerrero-Lebeau [GL07, Theorem 1] since one can check that the counterexample they build is a nonnegative solution. Now, we prove

For any δ > 0 (we will later need 2δ < C 0 where C 0 is the constant in (5.11)), and for T ≥ T GCC (M, ∇ g f, ω) + 2δ, we prove the observability inequality for positive solution

(5.13)

The combination of (5.11) on the time interval (2δ, T ) together with (5.12) on the time interval (δ, 2δ) implies

(5.14) Now, applying (5.8) with λ > 0 such that λT = δ, we obtain

(5.15)

Combining (5.14), (5.15) together with the Hölder inequality, we deduce

Choosing δ ∈ (0, C 0 /2) implies (5.13) uniformly for ε ∈ (0, ε 0 (δ)], and hence concludes the proof of the theorem.

From observability of positive solutions to a controllability statement

This section is devoted to the proof of the controllability result of Corollary 2.4 from the observability of positive solutions. It relies on the following lemma. The result and its proof follow [LB20, Theorem 4.1].

We obtain successively

(5.20)

We obtain that 1 ω u α is bounded in L 2 ([0, T ] × M) uniformly in α > 0. Take a sequence α n → 0 so that

The associated solutions y αn with control

) and, again up to a subsequence, converges weakly- * to a solution y of (2.1) with control h(t, x) = 1 ω u(Tt, x) and initial datum y 0 . Moreover, up to a subsequence, we can impose y αn (T ) ⇀ y(T ) in L 2 (M). Passing to the limit in (5.18) using (5.20), we finally obtain

for any p 0 ∈ V. We finally get the expected estimate on h(t, x) = u(Tt, x) passing to the limit in (5.19).

We may now conclude the proof of Corollary 2.4 from Lemma 5.9.

Proof of Corollary 2.4. We apply Lemma 5.9, with V = V = L 2 (M; R + ). Note that the Lemma applies because for any

This gives a control h with the expected uniform bound and so that (y(T ), u 0 ) L 2 (M) ≥ 0 for any u 0 ∈ L 2 (M; R + ). This implies y(T ) ≥ 0.

A About the distances

In this section, (M, g) is a compact Riemannian manifold without boundary.

A.1 A general lemma

We start with a general lemma.

Lemma A.1. Let V ∈ W 1,∞ (M) with nonnegative value. Then, for all x, y ∈ M 2 , we have

where

This lemma is particularly useful for V = (V -E) + in which case the (pseudo-) distance defined is the Agmon distance at energy level E.

Proof. We denote by d 1 , d 2 , d 3 respectively the three (pseudo-)distances defined in the statement of the lemma. Then, we notice that the last two quantities are invariant by reparametrization, so that d 2 = d 3 after a change of variable in the integral. Then, the inequality ab ≤ 1 2 (a 2 + b 2 ) directly yields d 2 ≤ d 1 . Let us now prove the converse inequality, namely d 3 ≥ d 1 . For ε > 0 there exist δ > 0 and a path γ : [0, 1] → M such that Using an approximation argument, we can further assume that γ is smooth up to replacing ε by 2ε. We now define the following reparametrization γ(s) = γ(ϕ(s)) where ϕ solves φ(s) = √

which gives d 1 ≤ d 3 + 2ε, and concludes the proof of the lemma.

A.2 Equivalence between the two definitions of d ∇gf

In this section, we prove equivalence between the two definitions of d ∇g f , respectively given in (1.24) in the introduction and in (5.5). We further give an equivalent quantity in terms of the Riemannian distance and the flow φ t of ∇ g f. The function f is assumed to be C 3 throughout the section.

We recall that ρ is defined in (5.2) with V = |∇g f| 2 g 4 , that is to say

and, for all t > 0 fixed, Lipschitz continuous as a function of (x, y) ∈ M 2 . These quantities are related to the Agmon distance but in finite time, see Section A.3 below. Note that the quantity ρ is symmetric, ρ(x, y, t) = ρ(y, x, t), and remains unchanged under the change of f by -f. This is not the case for d ∇g f .

Lemma A.2. The function d ∇g f defined as

is continuous on M 2 × (0, +∞) and, for all t > 0 fixed, Lipschitz continuous as a function of (x, y) ∈ M 2 . Moreover, we have

and g s is the time varying metric defined by

In particular, for any T ≥ 0, there exists some constant C T > 0 so that

where d denotes the Riemannian distance (associated to g).

Proof. The continuity property directly follows from that of ρ proved in [LY86, Appendix]. To prove (A.3), we remark that for any path γ so that γ(0

In particular, from the definition of d ∇g f in (A.2) and ρ in (A.1), we can rewrite d ∇g f (x, y, t) as

The statement (A.4) is obtained thanks to the change of path γ(s) = γ(ts). Now, we compute d ∇g f (x, φ t (y), t) according to formula (A.4). To this aim, let γ ∈ W 1,∞ ([0, t]; M) so that γ(0) = x, γ(t) = φ t (y). Let γ(s) = φ -s (γ(s)) so that γ(s) = ∇ g f(γ(s)) + Dφ s ( ˙ γ(s)). In particular,

. It gives (A.5) since any path γ ∈ W 1,∞ ([0, t]; M) so that γ(0) = x, γ(t) = y can be written γ(s) = φ s ( γ(s)) with γ(0) = x, γ(t) = φ -t (y), and conversely.

A.3 Further links between the different distances

In this section, we relate the above quantities ρ(x, y, t), d ∇g f (x, y, t) with the Agmon distance to the bottom energy (see 1.12 for V =

and E 0 = min M V = 0), that is to say

The results of this section are not explicitly used in the proofs of the main part of the paper; however we believe these links are interesting and enlightening. Indeed, they relate the quantity d A (x, y) appearing in all general bounds of Section 3 together with the quantities ρ(x, y, t), d ∇g f (x, y, t) appearing in results of Section 5 concerning positive solutions.

Lemma A.3. For all (x, y) ∈ M 2 , we have Proof. Equality in (A.8) is a consequence of Lemma A.1 applied to V = |∇ g f| 2 g . Then, (A.9) is a direct consequence of the expression of W and d ∇g f in terms of d A and ρ in (A.7), (A.2), together with (A.8).

Finally, if ∇ g f(y) = 0, then the function t → ρ(x, y, t) is non-increasing. Indeed, taking t 1 ≤ t 2 , from a path γ 1 : [0, t 1 ] → M such that γ 1 (0) = x and γ 1 (t 1 ) = y, we can construct the path γ

]; M), and thus the set of admissible paths on [0, t 2 ] is larger than the set of admissible paths on

g ds = 0 we deduce that ρ(x, y, t 2 ) ≤ ρ(x, y, t 1 ). This proves that the inf is actually a lim t→+∞ in this case.

Note that related properties are proved in the Appendix of [START_REF] Helffer | Puits multiples en mécanique semi-classique. IV. Étude du complexe de Witten[END_REF]. For instance, [HS85, Lemma A2.2] with our notations can be loosely stated as follows: If W (x, y) = 0, then every minimizing geodesic of d A is a generalized integral curve of ∇ g f. Finally, we state a last result that explains that d ∇g f (x, y, t) measures how far x is the final state of a path of the vector field at time t and starting at y. Part of this result is contained in the last statement of Lemma A.2; we here give a different proof, which, we believe, is interesting in itself.

Proposition A.4. With d ∇g f (x, y, t) defined in (5.5), we have for all (x, y, t) ∈ M 2 × (0, +∞), 1. d ∇g f (x, y, t) ≥ 0; 2. d ∇g f (x, y, t) = 0 if and only if there exists a trajectory of γ(s) = ∇ g f(γ(s)) with γ(0) = x, γ(t) = y, that is if and only if y = φ t (x).

In particular, (M, ∇ g f, ω, T ) satisfies (GCC) if and only if for any x ∈ M, there exist y ∈ ω and t ∈ (0, T ) so that d ∇g f (x, y, t) = 0.

Recall that the flow (φ t ) t∈R is defined in (2.14) and the Geometric Control Condition (GCC) is defined in Definition 2.10.

Proof. Statement 1 follows from the definition of d ∇g f in (A.4). Let us now consider Statement 2. Assume first that there exists a trajectory of γ(s) = ∇ g f(γ(s)) with γ(0) = x and γ(t) = y. Then, by definition of the infimum in (A.4), this yields d ∇g f (x, y, t) ≤ 0 and hence d ∇g f (x, y, t) = 0. Conversely, assume d ∇g f (x, y, t) = 0. Take a minimizing sequence in (A.4), that is to say that γ n ∈ W 1,∞ ([0, t]; M) such that γ n (0) = x, γ n (t) = y and γn -∇ g f(γ n ) = R n (bounded continuous with values in the tangent bundle of M) with Since ∇ g f is bounded on M, the sequence t 0 | γn (s)| 2 g ds is then uniformly bounded in R. As a consequence, the sequence of paths (γ n ) n∈N is equicontinuous. From Ascoli's theorem, we may extract a subsequence (which we do not relabel) (γ n ) n∈N which converges strongly for the topology C 0 ([0, t]; M) to a limit γ ∈ C 0 ([0, t]; M). The latter thus has the same and endpoints γ(0) = x and γ(t) = y. It is solution of γ = ∇ g f(γ) in the distributional sense according to (A.10). Bootstrapping in the differential equation implies γ ∈ W 1,∞ ([0, 1]; M) and γ is a strong solution to γ = ∇ g f(γ). This concludes the proof of Statement 2.

According to Lemma 2.11 Item 4, that (M, X, ω, T ) satisfies (GCC) is equivalent to the fact that for any x ∈ M, there exist t ∈ (0, T ) such that y := φ t (x) ∈ ω. As a consequence of Item 2, this is equivalent to having d ∇g f (x, y, t) = 0.

Note that an analogue statement for the function W in (A.9) is proved in [HS85, Lemma A2.2], and could also be deduced from Proposition A.4.