Assia Benabdallah
email: assia.benabdallah@univ-amu.fr

Franck Boyer
email: franck.boyer@math.univ-toulouse.fr

AND Morgan Morancey
email: morgan.morancey@univ-amu.fr

A BLOCK MOMENT METHOD TO HANDLE SPECTRAL CONDENSATION PHENOMENON IN PARABOLIC CONTROL PROBLEMS

Keywords: Control theory, parabolic partial differential equations, minimal null control time, block moment method AMS subject classifications. 93B05, 93C20, 93C25, 30E05, 35K90, 35P10

This article is devoted to the characterization of the minimal null control time for abstract linear control problem. More precisely we aim at giving a precise answer to the following question: what is the minimal time needed to drive the solution of the system starting from any initial condition in a given subspace to zero? Our setting will encompass a wide variety of systems of coupled one dimensional linear parabolic equations with a scalar control.

Following classical ideas we reduce this controllability issue to the resolution of a moment problem on the control and provide a new block resolution technique for this moment problem. The obtained estimates are sharp and hold uniformly for a certain class of operators. This uniformity allows various applications for parameter dependent control problems and permits us to deal naturally with the case of algebraically multiple eigenvalues in the underlying generator.

Our approach sheds light on a new phenomenon: the condensation of eigenvalues (which can cause a non zero minimal null control time in general) can be somehow compensated by the condensation of eigenvectors. We provide various examples (some are abstract systems, others are actual PDE systems) to highlight those new situations we are able to manage by the block resolution of the moment problem we propose.

This paper is concerned with the following abstract linear control system [START_REF] Allonsius | Boundary null-controllability of semi-discrete coupled parabolic systems in some multi-dimensional geometries[END_REF] y (t) + Ay(t) = Bu(t), y(0) = y 0 .

We are more precisely interested in the minimal time issue for null-controllability, which can be roughly expressed as follows: what is the smallest time T 0 ≥ 0 such that, for any T > T 0 , for any initial condition y 0 , there exists a control u such that the associated solution of (1) satisfies y(T ) = 0? Under quite general assumptions, we shall give formulas (that are reasonably explicit and usable) for such a minimal control time. The precise notion of solution as well as the wellposedness result for such system will be detailed below (see Propositions 1.1 and 1.2).

We will consider assumptions on the operator A that will relate (1) to parabolic evolution equations. Thus, due to regularizing properties, one cannot expect to reach any target in the state space and should rather try to reach any trajectory. By linearity, this is equivalent to the aforementioned null-controllability property (see for instance [START_REF] Coron | Control and nonlinearity[END_REF]Secs. 2.3 and 2.5]).

Pioneering works for null-controllability of a scalar one dimensional heat equations are due to H.O. Fattorini and D.L. Russell [START_REF] Fattorini | Exact controllability theorems for linear parabolic equations in one space dimension[END_REF][START_REF] Fattorini | Uniform bounds on biorthogonal functions for real exponentials with an application to the control theory of parabolic equations[END_REF]. For instance, they proved null-controllability of the one dimensional heat equation using a nonhomogeneous boundary condition as a control. For this purpose, they give a direct strategy reducing the null-controllability property to a moment problem that the control should satisfy (see Section 1.4 for a presentation of the moment problem). The moment method they propose consists in solving this problem using a biorthogonal family in L 2 (0, T ) to the family of exponentials associated to the eigenvalues of A * . Let us mention that this idea of reducing a (optimal) control problem to a moment problem is already present in the works [START_REF] Egorov | Some problems in the theory of optimal control[END_REF] by J.V. Egorov and [30] by L.I. Gal'chuk.

Later on, A.V. Fursikov, O.Yu. Imanuvilov [START_REF] Fursikov | Controllability of evolution equations[END_REF] and G. Lebeau, L. Robbiano [START_REF] Lebeau | Contrôle exact de l'équation de la chaleur[END_REF] used Carleman estimates to solve the boundary and internal null-controllability problem of the heat equation in any space dimension. These two papers have generated plenty of null-controllability results for various parabolic problems. The common qualitative behavior of these results is that for scalar parabolic equations nullcontrollability holds in arbitrary time (i.e. T 0 = 0) and without any restriction on the control domain.

Among all of these results let use mention the peculiar work [START_REF] Dolecki | Observability for the one-dimensional heat equation[END_REF] by S. Dolecki.

He considered a one dimensional heat equation, with a scalar control located at one point inside the space interval, and proved that choosing suitably the location of this control point one can achieve any value in [0, +∞] for the minimal null-control time T 0 . Until the years 2000's this work seemed to be considered too peculiar and the possible presence of a positive minimal null-control time (a very natural property in the hyperbolic case) was expected to be generically not possible in a parabolic setting. However, this point of view has been reconsidered recently in various works as for instance: [START_REF] Khodja | Null-controllability of some systems of parabolic type by one control force[END_REF] for abstract control problems, [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF] for the abstract problem (1) with applications to systems of one dimensional coupled parabolic equations, and [START_REF] Beauchard | Null controllability of Grushin-type operators in dimension two[END_REF] for a degenerate parabolic two dimensional equation of Grushin type.

Since then, the minimal null-control time property was investigated on various examples, still in the setting of coupled one dimensional parabolic systems [START_REF] Khodja | New phenomena for the null controllability of parabolic systems: minimal time and geometrical dependence[END_REF][START_REF] Duprez | Controllability of a 2 × 2 parabolic system by one force with space-dependent coupling term of order one[END_REF][START_REF] Ouaili | Minimal time of null controllability of two parabolic equations[END_REF][START_REF] Samb | Boundary null-controllability of two coupled parabolic equations: simultaneous condensation of eigenvalues and eigenfunctions[END_REF] or in the setting of degenerate parabolic scalar equations [START_REF] Beauchard | Heat equation on the Heisenberg group: observability and applications[END_REF][START_REF] Beauchard | Minimal time issues for the observability of Grushin-type equations[END_REF][START_REF] Beauchard | Degenerate parabolic operators of Kolmogorov type with a geometric control condition[END_REF][START_REF] Beauchard | 2D Grushin-type equations: minimal time and null controllable data[END_REF][START_REF] Duprez | Control of the Grushin equation: non-rectangular control region and minimal time[END_REF].

For coupled parabolic systems a geometric control condition may also be needed for approximate controllability to hold [START_REF] Boyer | Approximate controllability conditions for some linear 1D parabolic systems with space-dependent coefficients[END_REF][START_REF] Olive | Boundary approximate controllability of some linear parabolic systems[END_REF], proving once again that hyperbolic-like behavior can be observed in the parabolic setting.

Concerning the study of the abstract control problem [START_REF] Allonsius | Boundary null-controllability of semi-discrete coupled parabolic systems in some multi-dimensional geometries[END_REF], some characterizations are provided in the series of works [START_REF] Gaunard | Problèmes d'Interpolation dans les Espaces de Paley-Wiener et Applications en Théorie du Contrôle[END_REF][START_REF] Jacob | Interpolation by vector-valued analytic functions, with applications to controllability[END_REF][START_REF] Jacob | Weighted interpolation in Paley-Wiener spaces and finite-time controllability[END_REF][START_REF] Jacob | Weighted multiple interpolation and the control of perturbed semigroup systems[END_REF][START_REF] Jacob | Applications of Laplace-Carleson embeddings to admissibility and controllability[END_REF] using the formalism of Carleson measures. However the precise question of an abstract characterization of the minimal null-control time has not been much considered. A formula has been given for the minimal null-control time of system [START_REF] Allonsius | Boundary null-controllability of semi-discrete coupled parabolic systems in some multi-dimensional geometries[END_REF], in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF], in a setting encompassing coupled one dimensional parabolic equations with a scalar control. Its value depends on the condensation index of the eigenvalues of A * (see Section 7.5 for a precise definition)

and the observation of the associated eigenvectors. In this work the authors assume that the eigenvectors form a Riesz basis of the state space. Let us also mention the work [START_REF] Khodja | Quantitative Fattorini-Hautus test and minimal null control time for parabolic problems[END_REF] where the null-control time is studied through a resolvent-like inequality (introduced in [START_REF] Duyckaerts | Resolvent conditions for the control of parabolic equations[END_REF]) that is a quantification of the well-known Fattorini-Hautus test for approximate controllability (see [START_REF] Fattorini | Some remarks on complete controllability[END_REF][START_REF] Olive | Boundary approximate controllability of some linear parabolic systems[END_REF]). It is an abstract characterization that might not be easily computable on actual examples but provides a common setting for all the previous examples exhibiting a positive null-control time. The last two mentioned results also rely on the moment method. Note that, even if the Carleman approach is very powerful, it does not seem to be applicable to all the systems of interest: in many situations (including the ones discussed in Section 5.2) the moment

2
This manuscript is for review purposes only.

method is still the only successful technique up to now.

To highlight the limitations of the existing litterature on such problems and the improvements we propose, let us consider the following control problem

             ∂ t y + -∂ xx 1 0
-∂ xx + c(x) y = 0 0 , (t, x) ∈ (0, T ) × (0, 1), y(t, 0) = 0 u(t) , y(t, 1) = 0 0 , t ∈ (0, T ), y(0, x) = y 0 (x),

x ∈ (0, 1), where c ∈ L 2 (0, 1; R) is a given potential. We insist on the fact that our goal is not to study this particular example but to develop a general characterization. The application to this particular example is detailed in Section 5.2. The study of the minimal null-control time for this system for an arbitrary potential c is not covered by the litterature for several reasons.

• First, depending on c, the underlying operator can have geometrically double eigenvalues. This induces (a finite number of) non-observable modes and thus prevents even approximate controllability. We thus propose to extend the study of the minimal null-control time to a given subspace of initial conditions.

This allows to still analyze the controllability properties in this case.

• Even if the potential c is such that the eigenvalues are geometrically simple it can happen that some of them are algebraically double. In this case, to the best of our knowledge, the only existing results are [START_REF] Fernández-Cara | Boundary controllability of parabolic coupled equations[END_REF][START_REF] Khodja | The Kalman condition for the boundary controllability of coupled parabolic systems. Bounds on biorthogonal families to complex matrix exponentials[END_REF] which ensures null-controllability in arbitrary time if the eigenvalues are well separated (i.e.

satisfy the classical gap condition recalled in [START_REF] Fursikov | Controllability of evolution equations[END_REF]).

• Finally, if the potential c is such that the eigenvalues are geometrically and algebraically simple, to the best of our knowledge, the only existing result can be found in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF]. Under an extra assumption (on the observability of eigenfunctions), it provides null-controllability at any time T satisfying

T > T * = lim sup λ∈σ(A * )
-ln dist λ, σ(A * )\{λ} λ .

However, their arguments to disprove null-controllability at time T < T * cannot be applied in this example when the potential c is such that the family of eigenvectors forms a complete family but is not a Riesz basis. Therefore, the above formula for T * may dramatically overestimate the actual null-control time for the system. We will see in Section 5.2.1 that it may happen that T * = +∞ whereas the system is actually null-controllable at any time T > 0.

We will use quite general assumptions in our analysis answering all these concerns in the case of scalar control (see [START_REF] Boyer | Analysis of non scalar control problems for parabolic systems by the block moment method[END_REF] for an extension to non-scalar control). Doing so, we will prove that the difference between the Riesz basis assumption and the complete family assumption for the eigenvectors is not only technical and a new phenomenon can appear: the condensation of eigenvalues can be compensated by the condensation of eigenvectors.

We continue this introduction by stating more precisely the problem under consideration and the obtained results.

Functional setting.

Let X be a separable Hilbert space, whose inner product and norm are denoted by (•, •) and • respectively. We shall systematically identify X to its dual through the Riesz theorem. Let (A, D(A)) be an unbounded operator in X such that -A generates a C 0 -semigroup in X and (A * , D(A * )) its adjoint in X. Up to a suitable translation, we can assume that 0 is in the resolvent set of A. We denote by X 1 (resp.

X * 1 ) the Hilbert space D(A) (resp. D(A * )) equipped with the norm x 1 := Ax (resp. x 1 * := A * x ). We define X -1 as the completion of X with respect to the norm

y -1 := sup z∈X * 1 (y, z) z 1 * .
Notice that X -1 is isometrical to the topological dual of X * 1 using X as a pivot space (see for instance [49, Proposition 2.10.2]); the corresponding duality bracket will be denoted by •, • -1,1 * .

Let U be an Hilbert space (that we will identify to its dual) and B : U → X -1 be a linear continuous control operator. We denote by B * : X * 1 → U its adjoint in the duality described above.

proposition 1.1. Under the above assumptions, for any T > 0, any y 0 ∈ X -1 , and any u ∈ L 2 (0, T ; U ), there exists a unique y ∈ C 0 ([0, T ]; X -1 ) solution to [START_REF] Allonsius | Boundary null-controllability of semi-discrete coupled parabolic systems in some multi-dimensional geometries[END_REF] in the sense that it satisfies for any t ∈ [0, T ] and any

z t ∈ X * 1 , (2)
y(t), z t -1,1 * -y 0 , e -tA * z t -1,1 * = t 0 u(s), B * e -(t-s)A * z t U ds.
Moreover there exists C T > 0 such that sup t∈[0,T ] y(t) -1 ≤ C T y 0 -1 + u L 2 (0,T ;U ) .

The proof of this result is recalled in Appendix 7.1. Let us mention that this notion of solution is very weak. In most works concerning controllability properties for abstract systems like [START_REF] Allonsius | Boundary null-controllability of semi-discrete coupled parabolic systems in some multi-dimensional geometries[END_REF], an extra admissibility assumption is made on the control operator B to ensure more regularity for the solutions. Note however that this is not mandatory to prove wellposedness of the system in the weak sense above. We will discuss below the regularity properties of the control problem.

Let (X * , . * ) be an Hilbert space such that X * 1 ⊂ X * ⊂ X with dense and continuous embeddings. We assume that X * is stable by the semigroup generated by -A * (see Remark 1.1). We also define X -as the subspace of X -1 defined by

X -:= y ∈ X -1 ; y -:= sup z∈X * 1 y, z -1,1 * z * < +∞ ,
which is also isometrical to the dual of X * with X as a pivot space. The corresponding duality bracket will be denoted by •, • -, . Thus, we end up with the following five functional spaces

X * 1 ⊂ X * ⊂ X ⊂ X -⊂ X -1 .
We say that the control operator B is an admissible control operator for (1) with respect to the space X -if for any T > 0 there exists C T > 0 such that (3)

T 0 B * e -(T -t)A * z 2 U dt ≤ C T z 2
Notice that if (3) holds for some T > 0 it holds for any T > 0. The admissibility condition [START_REF] Khodja | Null-controllability of some systems of parabolic type by one control force[END_REF] implies that, by density, we can give a meaning to the map t → B * e -(T -t)A * z ∈ L 2 (0, T ; U ), for any z ∈ X * .

In this setting, following the lines of [START_REF] Coron | Control and nonlinearity[END_REF]Theorem 2.37] we obtain the following regularity result for the solutions.

proposition 1.2. Assume that (3) holds. Then, for any T > 0, any y 0 ∈ X -, and any u ∈ L 2 (0, T ; U ), there exists a unique y ∈ C 0 ([0, T ]; X -) solution to (1) in the sense that it satisfies for any t ∈ [0, T ] and any z t ∈ X * ,

y(t), z t -, -y 0 , e -tA * z t -, = t 0 u(s), B * e -(t-s)A * z t U ds.
Moreover there exists C T > 0 such that

sup t∈[0,T ] y(t) -≤ C T y 0 -+ u L 2 (0,T ;U ) .
Remark 1.1. Note that a similar regularity result holds if we don't assume that X * is stable by the semigroup generated by -A * except that we need to restrict ourselves to initial data y 0 ∈ X. In that case the solution satisfies for any t ∈ [0, T ] and any z t ∈ X * ,

y(t), z t -, -y 0 , e -tA * z t = t 0 u(s), B * e -(t-s)A * z t U ds, sup t∈[0,T ] y(t) -≤ C T y 0 + u L 2 (0,T ;U ) .
Remark 1.2. The case where X * = X * 1 means that we do not have any additional regularity property for B. Conversely, the case X * = X means that we have the best regularity we can hope for system (1) (this is the usual definition of an admissible control operator as in [START_REF] Coron | Control and nonlinearity[END_REF][START_REF] Tucsnak | Observation and control for operator semigroups[END_REF]).

To give more accurate results, we aim at analyzing the minimal null-control time problem for each specified set of initial data. This is the object of the following definition.

Definition 1.1. Let Y 0 be a closed subspace of X -.

We say that system (1) is null-controllable at time T from Y 0 if for any y 0 ∈ Y 0 there exists u ∈ L 2 (0, T ; U ) such that the associated solution of (1) satisfies y(T ) = 0.

As a specific choice of Y 0 one can think of Y 0 = X -, in which case we recover the classical notion of null-controllability. On the opposite side, if Y 0 is a one dimensional subspace Y 0 = Span{y 0 }, then the notion above amounts to consider only the nullcontrollability of the system for that particular initial condition y 0 .

From now on, we will assume that the space Y 0 is given, and we denote by P Y0 the orthogonal projection onto Y 0 with respect to • -and by P * Y0 ∈ L(X * ) its adjoint in the duality X -, X * . Notice that these definitions yield (4)

P * Y0 z * ≤ z * , ∀z ∈ X * . 5
This manuscript is for review purposes only.

Notations. We give here some notations that will be used throughout this article.

• For any integers a, b, c ∈ N, we shall define the following subsets of N:

a, b := [a, b] ∩ N, a, b =c := a, b \ {c}.
• For any complex number µ ∈ C we define e µ : (0, +∞) → C to be the exponential function

(5) e µ : s → e -µs .

• We shall denote by C γ1,...,γ l > 0 a constant possibly varying from one line to another but depending only on the parameters γ 1 , . . . , γ l .

• For any multi-index α ∈ N n , we denote its length by |α| = n j=1 α j and its maximum by |α| ∞ = max j∈ 1,n α j .

For α, µ ∈ N n , we say that µ ≤ α if and only if µ j ≤ α j for any j ∈ 1, n .

• For any finite subset A ⊂ C, we will make use of the polynomial P A defined by ( 6)

P A (x) := µ∈A (x -µ).
It satisfies in particular, for any λ ∈ A,

P A (λ) = µ∈A µ =λ (λ -µ).
1.3. Presentation of the main results.

1.3.1. Spectral assumptions. In addition to the hypotheses described in the introduction that are necessary for the well-posedness and regularity of our control problem, we shall make now the following structural assumptions.

• First of all, we shall only consider scalar controls in this paper, that is U = R.

• We assume that the spectrum of A * is only made of a countable number of geometrically simple eigenvalues denoted by Λ. We refer to Section 6.3 for a discussion on this assumption.

We shall also assume for simplicity that the eigenvalues are all real (see however the discussion in Section 6.1) and that

(7) Λ ⊂ [1, +∞).
Note that, if [START_REF] Khodja | Quantitative Fattorini-Hautus test and minimal null control time for parabolic problems[END_REF] does not hold, we can replace A by A + γ for γ > 0 large enough and find an associated null-control u. A null-control for the original problem is then given by t → e γt u(t) and we can explicitly bound its cost with respect to the parameters γ and T .

• For any eigenvalue λ ∈ Λ, we denote by α λ ≥ 1 its algebraic multiplicity and we assume that there exists an integer η ≥ 1 such that α λ ≤ η for any λ ∈ Λ.

• The main structural assumptions on the eigenvalues Λ we shall make in this paper are the following:

-Asymptotic behavior:

(8) λ∈Λ 1 λ < +∞.
-Weak gap condition with parameters p ∈ N and ρ > 0:

(9) # Λ ∩ [µ, µ + ρ] ≤ p, ∀µ ∈ [0, +∞).
In the case p = 1, the weak gap condition above simply reduces to

|λ -λ | > ρ, ∀λ, λ ∈ Λ, λ = λ ,
which is the usual gap condition used for instance in [START_REF] Fattorini | Exact controllability theorems for linear parabolic equations in one space dimension[END_REF]. If the spectrum Λ is increasingly indexed as Λ = (λ m ) m≥1 the weak gap condition (9) reads

λ m+p -λ m > ρ, ∀m ≥ 1.
As we will use a different labelling of the spectrum in this paper we shall not use these notations anymore in what follows.

• We denote by (φ 0 λ ) λ∈Λ an associated family of eigenvectors of A * . These eigenvectors are chosen to be normalized in X * .

As we are interested in null-controllability properties of system (1), we will first assume that (10)

B * φ 0 λ = 0, for any λ ∈ Λ.
This is a necessary condition for the approximate controllability of system (1), and is therefore mandatory if we expect null-controllability to hold. In our setting, the assumption ( 10) is also a sufficient condition for approximate controllability (see [START_REF] Fattorini | Some remarks on complete controllability[END_REF][START_REF] Olive | Boundary approximate controllability of some linear parabolic systems[END_REF]).

When the considered set of initial data Y 0 is not the whole space X -, the approximate controllability condition [START_REF] Beauchard | Heat equation on the Heisenberg group: observability and applications[END_REF] can be too strong and we can relax it. We will discuss this point in Section 6.2.

• For each λ ∈ Λ, we denote by (φ l λ ) l∈ 1,α λ -1 a Jordan chain associated with φ 0 λ , that is a family satisfying

A * φ l λ = λφ l λ + φ l-1 λ , ∀l ∈ 1, α λ -1 .
By [START_REF] Beauchard | Heat equation on the Heisenberg group: observability and applications[END_REF], we may uniquely determine such Jordan chain if we impose in addition that the generalized eigenvectors satisfy

(11) B * φ l λ = 0, ∀l ∈ 1, α λ -1 .
This particular choice of the Jordan chain is not mandatory but will simplify the forthcoming computations. In the case were the eigenvalues are algebraically simple (η = 1) we drop the superscipt 0 for the eigenvectors.

• We introduce the notation

Φ := φ l λ , λ ∈ Λ, l ∈ 0, α λ -1 ,
for the family of all the (generalized) eigenvectors of A * . We assume that Φ is complete in X * i.e. for any y ∈ X -, we have [START_REF] Beauchard | Minimal time issues for the observability of Grushin-type equations[END_REF] y, φ -, = 0, ∀φ ∈ Φ =⇒ y = 0.

We emphasize the fact that we will not make any additional assumptions on the family Φ. This is a very important difference with related results in the literature in which, most of the time, it is assumed that Φ forms a Riesz basis of X * . This is discussed in Sections 1.3.4 and 3.

In the forthcoming paper [START_REF] Boyer | Analysis of non scalar control problems for parabolic systems by the block moment method[END_REF], we will study the extension of our analysis to the case of possibly infinite dimensional controls.

1.3.2. Groupings of eigenvalues. To introduce our formula for the minimal null-control time it is convenient to define adapted groupings for the spectrum Λ. We highlight that this notion does not exactly coincide with the condensation groupings introduced by Bernstein [START_REF] Bernstein | Leçons sur les Progrès Récents de la Théorie des Séries de Dirichlet[END_REF], even though it is closely related.

Definition 1.2. Let p ∈ N * and r, ρ > 0. A sequence of sets (G k ) k≥1 ⊂ P(Λ)
is said to be a grouping for Λ with parameters p, r, ρ, and we will write

(G k ) k≥1 ∈ G(Λ, p, r, ρ), if it is a covering of Λ Λ = k≥1 G k ,
with the additional properties that for every k ≥ 1,

g k := #G k ≤ p, sup(G k ) < inf(G k+1 ), (13) dist
(G k , G k+1 ) ≥ r, and
(14) diam G k < ρ.
We prove in Appendix (Proposition 7.1) that such a grouping always exists for any Λ satisyfing the weak gap condition [START_REF] Baiocchi | Ingham-Beurling type theorems with weakened gap conditions[END_REF].

Once we are given such a grouping, we shall always adopt the following labelling of the elements of Λ

G k = {λ k,1 , . . . , λ k,g k } with λ k,1 < • • • < λ k,g k ,
and the (generalized) eigenvectors will be relabelled accordingly

φ l k,j := φ l λ k,j , ∀k ≥ 1, ∀j ∈ 1, g k , ∀l ∈ 0, α k,j -1 ,
where in the same fashion α k,j := α λ k,j . For any k ≥ 1, we gather the multiplicities associated with the elements of G k in the multi-index

α k = (α λ k,1 , . . . , α λ k,g k ) ∈ N g k . Remark 1.3. The condition λ k,1 < • • • < λ k,g k is convenient to treat the abstract problem (1
) but might not be convenient in actual examples. As all the estimates in our analysis will depend on the parameters p and ρ, the eigenvalues inside a same group are mostly interchangeable and thus the increasing labelling is not needed.

Minimal control time definition.

From now on, we assume given a grouping (G k ) k in G(Λ, p, r, ρ). Thanks to the assumption [START_REF] Beauchard | Heat equation on the Heisenberg group: observability and applications[END_REF], we can define the following family of elements in X * (15) ψ l k,j :=

P * Y0 (φ l k,j ) B * φ 0 k,j , ∀k ≥ 1, ∀j ∈ 1, g k , ∀l ∈ 0, α k,j -1 .

8

This manuscript is for review purposes only.

Let (16)

T 0 (Y 0 ) := lim sup k→∞ ln max

µ≤α k ψ[λ (µ1) k,1 , . . . , λ (µg k ) k,g k ] * λ k,1
where the notation ψ[. . . ] stands for the generalized divided differences (see Section

ψ l k,j ; j ∈ 1, g k , l ∈ 0, α k,j -1
whose coefficients can be explicitely computed on actual control problems (see Section 5) and that only depends on the group G k and on the multiplicity multi-index µ.

In the simpler case where the eigenvalues are assumed to be algebraically simple (i.e. η = 1) we can immediately give a more explicit formula for T 0 (Y 0 ). Indeed, in this case one recovers the standard divided differences (whose definition and properties are recalled in Section 7.3.1) and thus

(17) T 0 (Y 0 ) = lim sup k→∞ ln   max m,l∈ 1,g k m≤l ψ[λ k,m , . . . , λ k,l ] *   λ k,1 .
Then, using Corollary 7.1 and ( 14) it comes that the computation of all those divided differences is not needed and the formula reduces to

T 0 (Y 0 ) = lim sup k→∞ ln max l∈ 1,g k ψ[λ k,1 , . . . , λ k,l ] * λ k,1 (18)
= lim sup

k→∞ ln   max l∈ 1,g k l j=1 ψ k,j i∈ 1,l =j (λ k,j -λ k,i ) *   λ k,1 (19)
where the last equality comes from the use of Newton formula (see Proposition 7.3).

Remark 1.4. The definition above corresponds to a given grouping of the spectrum, however the minimal null-control result stated in Theorem 1.1 will show that its value does not depend on this particular choice of a grouping. As a consequence, for specific examples, one can compute the minimal null control time T 0 (Y 0 ) using any convenient such grouping in a class G(Λ, p, r, ρ).

For the sake of simplicity, for any y 0 ∈ X -we denote by T 0 (y 0 ) the quantity T 0 (Span(y 0 )). Of course, we have the following proposition relating T 0 (Y 0 ) and T 0 (y 0 )

for y 0 ∈ Y 0 . proposition 1.3. For any closed subspace Y 0 ⊂ X -, sup y0∈Y0 T 0 (y 0 ) = T 0 (Y 0 ).
This assertion is proved in Section 7.4.

9

This manuscript is for review purposes only.

Remark 1.5. Let us discuss the sign of T 0 (Y 0 ).

• In the case Y 0 = X -(the operator P * Y0 thus reduces to the identity), the minimal time T 0 (Y 0 ) is always non-negative. Indeed, from the case µ =

(1, 0, . . . , 0) in the definition (16) of T 0 we have that

T 0 (X -) ≥ lim sup k→∞ ln φ 0 k,1 * |B * φ 0 k,1 | λ k,1
.

From the admissibility condition (3) applied to z = φ 0 k,1 , we deduce the fol-

lowing upper bound B * φ 0 k,1 ≤ C λ k,1 φ 0 k,1 *
. Thus,

lim sup k→∞ ln φ 0 k,1 * |B * φ 0 k,1 | λ k,1 ≥ 0, which proves that T 0 (X -) ∈ [0, +∞].
• In the general case where Y 0 is a strict closed subspace of X -, it may happen that T 0 (Y 0 ) < 0.

For instance, if we choose y 0 ∈ X 1 to be an eigenvector of A for an eigenvalue λ ∈ Λ, then we have T 0 (y 0 ) = -∞. Indeed, we first observe that

y 0 , φ 0 λ -, = y 0 , φ 0 λ = 0, ∀λ ∈ Λ, λ = λ,
which implies, with Y 0 = Span(y 0 ) that P * Y0 φ 0 λ = 0 for any λ = λ. We deduce that the logarithms in the definition of T 0 (y 0 ) are all equal to -∞ for k large enough.

1.3.4. Null-controllability result. The main result of this paper reads as follows (see also the extension discussed in Section 6.1).

theorem 1.1. Assume that the operators A and B satisfy the assumptions given in Section 1.3.1. Let T > 0 and T 0 (Y 0 ) be defined by [START_REF] Bhandari | Boundary null-controllability of coupled parabolic systems with Robin conditions[END_REF]. Then, i. If T 0 (Y 0 ) < +∞ and T > T 0 (Y 0 ), the system (1) is null-controllable from Y 0 at time T .

ii. If T 0 (Y 0 ) > 0 and T < T 0 (Y 0 ), the system (1) is not null-controllable from Y 0 at time T .

Let us briefly mention that our strategy of proof relies on an adapted block resolution of the associated moment problem (see Theorems 2.1 and 4.1). In the case of spectral condensation this new method of resolution ensures sharper results than the one given by standard biorthogonal families. However, as a by-product, in the case of algebraically simple eigenvalues we recover the known optimal estimates for such biorthogonal families (see Corollary 2.1). In the case of algebraically multiple eigenvalues we provide new estimates for such biorthogonal families (see Corollary 4.1). Before describing with more details this strategy of proof let us make some comments.

• There are settings in which formulas for the minimal null-control time are already known in the literature for instance when the eigenvalues are algebraically simple and:

when the condensation index of Λ (see Appendix 7.5 for a precise definition) is equal to 0 (see [START_REF] Khodja | Quantitative Fattorini-Hautus test and minimal null control time for parabolic problems[END_REF]Remark 1.15]);

or when the family (φ λ ) λ∈Λ of eigenvectors forms a Riesz basis of X *

(see [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF]).

10

This manuscript is for review purposes only.

Obviously, in those settings we recover the known expressions. This is discussed in Section 3.1 and Section 3.2.

• However, we also prove that the Riesz basis assumption considered in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF] is not only technical. More precisely, we show in Proposition 3. • The weak gap condition ( 9) is particularly well adapted to the applications we have in mind, namely coupled one dimensional parabolic equations in which case the spectrum is given by a finite union of sequences satisfying a classical gap condition (see for instance Lemma 2.1).

The restriction to the one dimensional case in those applications comes from the assumption [START_REF] Avdonin | Riesz bases of exponentials and divided differences[END_REF]. Although this assumption can be seen as a restriction due to the use of moment method, as we are considering scalar controls (U = R) it is also a necessary null-controllability condition (see for instance [40, Appendix A]).

• As we precised the space of initial conditions in this study of minimal nullcontrol time, it directly comes that finite linear combination of eigenvectors are null-controllable in arbitrary small time: the existence of positive minimal null-control time is definitely a high-frequency phenomenon as already observed in Remark 1.5.

• In this article we not only prove Theorem 1.1 but we also develop a new strategy to solve moments problem: the block moment method presented in Section 1.4. The resolution of these problems (see Theorems 2.1 and 4.1)

is done with precise estimates. This not only leads to the construction of a control but also to uniform estimates, with respect to Λ in a certain class, on this control (see Corollaries 2.2 and 2.3). Those uniform estimates are important in various contexts when one wants to achieve bounds on the control for parameter-dependent problems (see for instance [START_REF] Allonsius | Boundary null-controllability of semi-discrete coupled parabolic systems in some multi-dimensional geometries[END_REF][START_REF] Allonsius | Spectral analysis of discrete elliptic operators and applications in control theory[END_REF] for an example in numerical analysis of null-controllability problems, or [START_REF] López | Uniform null-controllability for the one-dimensional heat equation with rapidly oscillating periodic density[END_REF] for an example in oscillating coefficient problems). Thus, the block moment method can be of interest to study parameter-dependent problems even in the presence of a strong gap condition, in particular in the case where the strong gap badly behaves with respect to the parameter. Actually, this strategy has already been applied in [START_REF] Bhandari | Boundary null-controllability of coupled parabolic systems with Robin conditions[END_REF] to study the boundary controllability of a coupled parabolic system with Robin boundary conditions, uniformly with respect to the Robin parameters. Moreover, this uniformity property will be crucial in Section 4 to infer the results on multiple eigenvalues from the ones on simple eigenvalues. We end this introduction by describing the global strategy used to prove Theorem 1.1 and giving some further bibliographical comments. Section 2 is dedicated to the proof of Theorem 1.1 in the case of algebraically simple eigenvalues. We provide in Section 3 a comparison of our results with available results of the literature. In

11
This manuscript is for review purposes only.

Section 4 we prove that the uniform estimates obtained in Section 2 allow to prove Theorem 1.1 in the general case of algebraically multiple eigenvalues. To highlight the new cases and phenomenon covered by our analysis we present different examples in Section 5. Then we propose some extensions in Section 6. To ease the reading we gather various technical results in Section 7.

1.4. Strategy of proof.

The proof of the positive controllability result (that is point i. of Theorem 1.1) relies on a block resolution of the moment problem. Let us give more details about this strategy.

Let y 0 ∈ Y 0 and u ∈ L 2 (0, T ; R) given. Using Proposition 1.2, it comes that y(T ) = 0 if and only if the control u satisfies

(20) -y 0 , e -T A * φ -, = T 0 u(t)B * e -(T -t)A * φ dt, ∀φ ∈ X * .
As the family Φ of (generalized) eigenvectors is assumed to form a complete family in X * (see [START_REF] Beauchard | Minimal time issues for the observability of Grushin-type equations[END_REF]) it is in fact sufficient to test [START_REF] Dolecki | Observability for the one-dimensional heat equation[END_REF] against the elements of Φ. Therefore a null control u is characterized by the following countable set of equations

(21) -y 0 , e -T A * φ l λ -, = T 0 u(t)B * e -(T -t)A * φ l λ dt, ∀λ ∈ Λ, ∀l ∈ 0, α λ -1 .
Using the formalism of generalized divided differences, we can give a convenient expression of the action of the semi-group on the generalized eigenvectors as follows [START_REF] Duprez | Control of the Grushin equation: non-rectangular control region and minimal time[END_REF] e -tA * φ l λ = e -λt l p=0

(-t) p p! φ l-p λ = l p=0 e (p) t (λ) p! φ l-p λ = l p=0 e t [λ (p+1) ]φ[λ (l-p+1) ] = (e t φ)[λ (l+1) ],
this last equality coming from Definition 7.3. Then, y(T ) = 0 if and only if for any λ ∈ Λ and any l ∈ 0, α λ -1 ,

T 0 u(T -t)B * (e t φ)[λ (l+1) ] dt = -y 0 , (e T φ)[λ (l+1) ] -,
.

By [START_REF] Beauchard | Null controllability of Grushin-type operators in dimension two[END_REF], and since y 0 ∈ Y 0 , this reduces to find u ∈ L 2 (0, T ; R) such that for any λ ∈ Λ and any l ∈ 0, α λ -1 ,

(B * φ 0 λ ) T 0 u(T -t) (-t) l l! e -λt dt = -y 0 , P * Y0 (e T φ)[λ (l+1) ] -,
, that is, using [START_REF] Beauchard | Heat equation on the Heisenberg group: observability and applications[END_REF] and ( 15),

T 0 u(T -t) (-t) l l! e -λt dt = -y 0 , (e T ψ)[λ (l+1) ] -, , ∀λ ∈ Λ, ∀l ∈ 0, α λ -1 , (23)
To solve this so-called moment problem the classical strategy introduced in [START_REF] Fattorini | Exact controllability theorems for linear parabolic equations in one space dimension[END_REF] consists in designing a biorthogonal family in L 2 (0, T ) to t → t l e -λt ; λ ∈ Λ, l ∈ 0, α λ -1 with associated estimates. Then, thanks to these estimates, a suitable control is defined. Usually in this procedure each biorthogonal element is estimated separately.

Thus, this method is somehow inoperent to analyse the possible condensation of eigenvectors (which is related to possible cancellations in linear combinations of right-hand sides of ( 23)). We will thus propose to solve this moment problem using the grouping introduced in Section 1.3.2, in order to cope with such possible compensations. We then look for a solution u in the form ( 24)

u(t) = - k≥1 q k (T -t)
where each q k will solve the moment problem corresponding to the group G k . More precisely, such a control will formally solve ( 23) if

             T 0 q k (t) (-t) l l ! e -λ k ,j t dt = 0, ∀k = k, ∀j ∈ 1, g k , ∀l ∈ 0, α k ,j -1 , T 0 q k (t) (-t) l l! e -λ k,j t dt = y 0 , (e T ψ)[λ (l+1) k,j ] -, , ∀k ≥ 1, ∀j ∈ 1, g k , ∀l ∈ 0, α k,j -1 . (25)
Then the proof of point i. of Theorem 1.1 reduces to the resolution of such a block moment problem with suitable estimates (see Theorem 4.1). First, we solve in Theorem 2.1 the block moment problem in the case where the eigenvalues are algebraically simple i.e.

(

)          T 0 q k (t)e -λ k ,j t dt = 0, ∀k = k, ∀j ∈ 1, g k , T 0 q k (t)e -λ k,j t dt = e -λ k,j T y 0 , ψ k,j -, , ∀k ≥ 1, ∀j ∈ 1, g k . 26
This construction uses a Laplace transform isomorphism together with a suitable restriction argument (Proposition 2.4). The obtained estimates on q k will allow to prove convergence of the series (24) when T > T 0 (Y 0 ). Those estimates are uniform with respect to Λ in a certain class (see Definition 2.1) which will allow in Section 4

to infer the resolution of (25) in the general case.

Remark 1.6. Contrarily to the classical strategy, notice that the sequence (q k ) k is not a biorthogonal family to

t → t l e -λt ; λ ∈ Λ, l ∈ 0, α λ -1 .
The function q k is only orthogonal to those functions corresponding to groups other than G k . Inside the group G k its definition is adapted to solve the moment problem [START_REF] Duyckaerts | Resolvent conditions for the control of parabolic equations[END_REF]. Through the right-hand side (adapted to each initial condition) we will possibly take into account the unsufficient observation of eigenvectors, the condensation of eigenvalues but also the condensation of eigenvectors. This construction

13
This manuscript is for review purposes only.

can thus be seen as a block moment method. As we consider at the same time the eigenvalues associated to a same group this will lead to sharper estimates than the one coming from the design of a biorthogonal family (i.e. when considering each eigenvalue individually).

However, as already mentioned, our strategy still allows to prove the existence and sharp estimates on biorthogonal families (see Corollary 2.1 and Corollary 4.1). Let us mention that, to the best of our knowledge, the estimate we obtain in Corollary 4.1 for a biorthogonal family in presence of algebraic multiplicity of eigenvalues without the standard gap condition was not known. Even though these biorthogonal families are not always suitable to deal with controllability properties in presence of spectral condensation (this is why we designed this block resolution of the moment problem)

they can be useful for other problems.

Let us mention that, in the context of control problems with a spectrum satisfying the weak-gap condition, divided differences were already used for instance in [START_REF] Avdonin | Riesz bases of exponentials and divided differences[END_REF][START_REF] Baiocchi | Ingham-Beurling type theorems with weakened gap conditions[END_REF].

Among other things, in theses works, the authors give a necessary and sufficient condition for the family of (generalized) divided differences

{t → e t [iλ k,1 ], . . . , t → e t [iλ k,1 , . . . , iλ k,g k ] , k ≥ 1}
to form a Riesz basis of L 2 (0, T ; C). The possible condensation of eigenvalues then appears to deduce properties on the original family of exponentials (t → e iλt ) λ∈Λ .

Their results are then applied to hyperbolic control problems.

The results presented in our work can be seen as the 'parabolic' equivalent (or the 'real-valued' equivalent if one focuses on the exponentials) of these results. Nevertheless, the control problems have really different behaviours as well as the families of exponentials. Indeed, in our setting the considered families never form a Riesz basis.

Thus, neither work can be deduced from the other.

The proof of point ii. of Theorem 1.1, relies on the optimality of the bounds proved in Theorem 2.1 for the resolution of the block moment problem.

As dealing with null-controllability from a proper subspace of initial conditions is not classical let us recall the following lemma that characterizes this controllability through an observability inequality.

Lemma 1.1 (see for instance [1, Lemma 2.1]). Let M > 0.
The following two propositions are equivalent.

1. For any y 0 ∈ Y 0 there exists a u ∈ L 2 (0, T ; U ) such that y(T ) = 0 and u L 2 (0,T ;U ) ≤ M y 0 -.

2.

For any z T ∈ X * , the following partial observability holds:

(27) P * Y0 e -T A * z T 2 * ≤ M 2 T 0 B * e -(T -t)A * z T 2 U dt.
In this case, the best constant M satisfying those properties is called the cost of controllability from Y 0 at time T and is denoted M (Y 0 , T ).

2. The case of simple eigenvalues.

Null-controllability in large time.

The goal of this section is to prove point i. of Theorem 1.1 in the case of algebraically simple eigenvalues. Thus, in all this section we assume that η = 1.

As explained in Section 1.4, we will now focus on the construction of a solution to [START_REF] Fattorini | Exact controllability theorems for linear parabolic equations in one space dimension[END_REF]. Of course as we want to design a control u ∈ L 2 (0, T ; U ) the estimate of q k L 2 (0,T ;R) will play a crucial role to prove that the series (24) makes sense. Actually we will prove sharp estimates that are uniformly valid for Λ in a certain class. These uniform estimates can be used for various applications and will be crucial to deal with algebraic multiplicity of eigenvalues in Section 4. We start by precising the class of Λ we will consider.

Definition 2.1. Let p ∈ N * , ρ > 0 and N : (0, +∞) → R. We say that a countable family Λ belongs to the class L w (p, ρ, N ) if Λ satisfies the weak-gap condition [START_REF] Baiocchi | Ingham-Beurling type theorems with weakened gap conditions[END_REF] with parameters p and ρ and if for any ε > 0 we have

(28) λ∈Λ λ≥N (ε) 1 λ < ε.
This definition is directly inspired by the pioneering work [START_REF] Fattorini | Uniform bounds on biorthogonal functions for real exponentials with an application to the control theory of parabolic equations[END_REF]. More precisely, the class of sequences used in [START_REF] Fattorini | Uniform bounds on biorthogonal functions for real exponentials with an application to the control theory of parabolic equations[END_REF] is similar to L w (1, ρ, N ), but it is however slightly different since in [START_REF] Fernández-Cara | Boundary controllability of parabolic coupled equations[END_REF] the summation condition is given on the value of λ itself whereas in the above reference the condition is on the index of the eigenvalue in Λ (which is supposed to be sorted increasingly). Despite this small difference (whose aim is to simplify some computations) the results we shall take from [START_REF] Fattorini | Uniform bounds on biorthogonal functions for real exponentials with an application to the control theory of parabolic equations[END_REF] that use this definition are also valid with this alternative definition and thus we set L(ρ, N ) := L w (1, ρ, N ).

Remark 2.1 (The usual gap condition). With our definition, a sequence Λ belongs to L(ρ, N ) if it satisfies the classical gap condition

(29) |λ -λ| > ρ, ∀λ, λ ∈ Λ, λ = λ ,
and the asymptotic behavior estimate [START_REF] Fernández-Cara | Boundary controllability of parabolic coupled equations[END_REF].

As we will see in the examples (Section 5), the typical situation where sequences satisfying the weak gap condition appear is when one glues a finite number of sequences, each of them satisfying a standard gap condition as in Remark 2.1. This is formalized in the following lemma.

Lemma 2.1. Let p, p ∈ N * , ρ, ρ > 0 and N , Ñ : (0, +∞) → R given. Then, for any Λ ∈ L w (p, ρ, N ) and Λ ∈ L w (p, ρ, Ñ ), we have

Λ ∪ Λ ∈ L w (p, ρ, N ), with p = p + p, ρ = min(ρ, ρ) and N (ε) = max(N (ε/2), Ñ (ε/2)).
Proof. Let us first prove the weak gap condition. For any µ ≥ 0, we have

[µ, µ + ρ] ∩ (Λ ∪ Λ) = ([µ, µ + ρ] ∩ Λ) ∪ [µ, µ + ρ] ∩ Λ ⊂ ([µ, µ + ρ] ∩ Λ) ∪ [µ, µ + ρ] ∩ Λ ,
and taking the cardinal, we get

#[µ, µ + ρ] ∩ (Λ ∪ Λ) ≤ p + p = p.
For the asymptotic behavior of the sequences, we have

λ∈Λ∪ Λ λ> N (ε) 1 λ ≤ λ∈Λ λ> N (ε) 1 λ + λ∈ Λ λ> N (ε) 1 λ ≤ λ∈Λ λ>N (ε/2) 1 λ + λ∈ Λ λ> Ñ (ε/2) 1 λ ≤ ε 2 + ε 2 = ε.
The claim is proved.

The following straightforward facts will also be useful.

Remark 2.2. Let Λ ∈ L w (p, ρ, N ). • Any Λ ⊂ Λ also satisfies Λ ∈ L w (p, ρ, N ). • For any h ∈ (0, inf Λ), Λ + h ∈ L w p, ρ, N 2 .
Using this class we prove the following theorem.

theorem 2.1. Let T ∈ (0, +∞]. Let p ∈ N * , r, ρ > 0 and N : (0, +∞) → R.
Assume that Λ ∈ L w (p, ρ, N ) and let (G k ) k ∈ G(Λ, p, r, ρ) be an associated grouping.

For any ε > 0, there exists a constant C ε,T,p,r,ρ,N > 0 such that for any k ≥ 1,

for any ω k,1 , . . . , ω k,g k ∈ C, there exists q k ∈ L 2 (0, T ; C) satisfying T 0 q k (t)e -λ k ,j t dt = 0, ∀k = k, ∀j ∈ 1, g k , (30a) T 0 q k (t)e -λ k,j t dt = ω k,j , ∀j ∈ 1, g k , (30b) and (31)
q k L 2 (0,T ;C) ≤ C ε,T,p,r,ρ,N e ελ k,1 max i∈ 1,g k ω[λ k,1 , . . . , λ k,i ] .
Moreover, up to the factor e ελ k,1 , this last estimate is sharp: any solution q k ∈ L 2 (0, T ; C) of (30b) satisfies

(32) q k L 2 (0,T ;C) ≥ C p max i∈ 1,g k ω[λ k,1 , . . . , λ k,i ] ,
for some C p > 0.

The proof of Theorem 2.1 is conducted all along Sections 2.1.1 and 2.1.2.

Before going on with the proof, let us notice that the resolution of the block moment problem [START_REF] Gal'chuk | Optimal control of systems described by parabolic equations[END_REF] for a specific choice of ω k,j allows to prove, as a by-product, the existence and uniform estimates of a biorthogonal family to the exponentials (e λ ) λ∈Λ where e λ is defined by [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF].

Corollary 2.1. Let T ∈ (0, +∞]. Let p ∈ N * , r, ρ > 0 and N : (0, +∞) → R. Assume that Λ ∈ L w (p, ρ, N ) and let (G k ) k ∈ G(Λ, p, r, ρ
) be an associated grouping.

For any k ≥ 1, for any j ∈ 1, g k , there exists q k,j ∈ L 2 (0, T ; R) satisfying

(33) T 0 q k,j (t)e -λ k ,j t dt = δ k,k δ j,j , ∀k, k ≥ 1, ∀j ∈ 1, g k , ∀j ∈ 1, g k ,
where δ denotes the Kronecker symbol. Moreover, for any ε > 0, there exists a constant C ε,T,p,r,ρ,N > 0 such that for any k ≥ 1 and for any j ∈ 1, g k ,

q k,j L 2 (0,T ;R) ≤ C ε,T,p,r,ρ,N e ελ k,1 P G k (λ k,j )
,

where P G k is defined in [START_REF] Khodja | New phenomena for the null controllability of parabolic systems: minimal time and geometrical dependence[END_REF].

Moreover, up to the factor e ελ k,1 , this estimate is optimal since any function q k,j satisfying (33), satisfies the lower bound

q k,j L 2 (0,T ;R) ≥ C p 1 P G k (λ k,j )
for some C p > 0.

Proof. Let k ≥ 1 and j ∈ 1, g k . Let q k,j ∈ L 2 (0, T ; C) be the solution of the block moment problem [START_REF] Gal'chuk | Optimal control of systems described by parabolic equations[END_REF] given by Theorem 2.1 associated with the right-hand side ω k,j = δ j,j for any j ∈ 1, g k . Since those values of ω are real we can change q k,j in its real part without changing its properties. Then, the equalities (33) follow directly.

Moreover we have

q k,j L 2 (0,T ;R) ≤ C ε,T,p,r,ρ,N e ελ k,1 max i∈ 1,g k ω[λ k,1 , . . . , λ k,i ] .
From the Newton formula (see Proposition 7.3) it comes that for any

i ∈ 1, g k , ω[λ k,1 , . . . , λ k,i ] =        0, if i < j, 1 m∈ 1,i =j (λ k,j -λ k,m ) , if i ≥ j.
To conclude the proof of Corollary 2.1 we prove that there exists C p,ρ > 0 such that

for any k ≥ 1, j ∈ 1, g k and any i ∈ j, g k , ( 34
) m∈ 1,i =j |λ k,j -λ k,m | ≥ C p,ρ |P G k (λ k,j )|.
Indeed, we have

m∈ 1,i =j |λ k,j -λ k,m | |P G k (λ k,j )| = m∈ 1,i =j |λ k,j -λ k,m | m∈ 1,g k =j |λ k,j -λ k,m | = 1 m∈ i+1,g k |λ k,j -λ k,m | .
By [START_REF] Beauchard | 2D Grushin-type equations: minimal time and null controllable data[END_REF], we get

|λ k,j -λ k,m | ≤ ρ, ∀m ∈ i + 1, g k . Thus, ( 35
) m∈ 1,i =j |λ k,j -λ k,m | |P G k (λ k,j )| ≥ 1 ρ g k -i
.

As the right-hand side only takes a finite number of values, inequality [START_REF] Jacob | Applications of Laplace-Carleson embeddings to admissibility and controllability[END_REF] proves [START_REF] Jacob | Weighted multiple interpolation and the control of perturbed semigroup systems[END_REF] and ends the proof of Corollary 2.1.

The lower bound directly follows from [START_REF] Jacob | Interpolation by vector-valued analytic functions, with applications to controllability[END_REF] and the inequality max

i∈ 1,g k ω[λ k,1 , . . . , λ k,i ] ≥ ω[λ k,1 , . . . , λ k,g k ] = 1 P G k (λ k,j )
.

Resolution of block moment problems in infinite time.

In this section, we start by proving Theorem 2.1 in the case of simple eigenvalues and with T = +∞. More precisely, we prove the following proposition.

proposition 2.1. Let p ∈ N * , r, ρ > 0 and N : (0, +∞) → R. Assume that Λ ∈ L w (p, ρ, N ) and let (G k ) k ∈ G(Λ, p, r, ρ) be an associated grouping.

For any ε > 0, there exists a constant C ε,p,r,ρ,N > 0 such that for any k ≥ 1, for

any ω k,1 , . . . , ω k,g k ∈ C, there exists q k ∈ L 2 (0, +∞; C) satisfying (36)
         +∞ 0 q k (t)e -λ k ,j t dt = 0, ∀k = k, ∀j ∈ 1, g k , +∞ 0 q k (t)e -λ k,j t dt = ω k,j , ∀j ∈ 1, g k , and
q k L 2 (0,+∞;C) ≤ C ε,p,r,ρ,N e ελ k,1 max i∈ 1,g k ω[λ k,1 , . . . , λ k,i ] .
The proof relies on the construction of an holomorphic function satisfying suitable properties and estimates. The resolution of the block moment problem (36) then comes from the isomorphism induced by the Laplace transform.

Proof. Let us start by recalling classical properties of the Laplace transform (see for instance [47, pp. 19-20] and the references therein). Let H 2 (C + ) the space of

holomorphic functions F on C + = {z ∈ C ; (z) > 0} such that sup σ>0 F (σ + i•) L 2 (R;C) < +∞,
endowed with the norm

F 2 H 2 (C + ) := sup σ>0 F (σ + i•) 2 L 2 (R;C) .
From the properties of H 2 (C + ), it comes that

F 2 H 2 (C + ) = R |F (iτ )| 2 dτ, ∀F ∈ H 2 (C + ).
Then the Laplace transform

L : f ∈ L 2 (0, +∞; C) → F : λ ∈ C + → +∞ 0 e -λt f (t)dt ∈ H 2 (C + )
is an isomorphism.

We shall construct for each k, a function

J k ∈ H 2 (C + ) satisfying J k (λ) = 0, ∀λ ∈ Λ\G k , (37)
J k (λ k,j ) = ω k,j , ∀j ∈ 1, g k , (38)
and such that for any ε > 0, there exists C ε,p,r,ρ,N > 0 such that

R |J k (iτ )| 2 dτ ≤ C ε,p,r,ρ,N e ελ k,1 max i∈ 1,g k ω[λ k,1 , . . . , λ k,i ] , ∀k ≥ 1.
Taking advantage of the isomorphism property of the Laplace transform we will then set q k := L -1 (J k ), to conclude the proof.

Construction of J k .

We define J k as

J k : z ∈ C + → P k (z) (1 + z) p W k (z)
where P k is a polynomial of degree less than p which is precised below and W k is the following Blaschke-type product

W k (z) = p j=1   λ∈Λj \G k λ -z λ + z   , where (39)
Λ j := {λ l,min(j,g l ) , l ≥ 1}.

The sequence Λ j contains the j-th element of each group G l , except if this group contains less than j elements, in which case, we replace it by the largest element of

G l that is λ l,g l .
In particular, we observe that Λ j is a subsequence of Λ.

From ( 8), we deduce that λ∈Λj 1 λ < +∞, so that for any j, the associated infinite product uniformly converges on any compact of C + . As a consequence, W k is well-defined and holomorphic in C + (see for instance [START_REF] Rudin | Real and complex analysis[END_REF]Chapter 15]). It follows that J k is also holomorphic on C + .

We shall need the following property, whose proof is technical and postponed to Section 2.1.4.

proposition 2.2. There exists a constant C ε,p,r,ρ,N > 0 such that for any k ≥ 1, any l ∈ 0, p , and any

θ ∈ Conv(G k ), ( 40
) 1 W k (l) (θ) ≤ C ε,p,r,ρ,N e ελ k,1 .
From the definition of W k it comes that (37) is satisfied. Next, it comes that [START_REF] Lebeau | Contrôle exact de l'équation de la chaleur[END_REF] is equivalent to

P k (λ k,j ) = (1 + λ k,j ) p W k (λ k,j ) ω k,j , ∀j ∈ 1, g k . Let f : s ∈ R → (1 + s) p , and
f k : s ∈ R → f (s) W k (s) .
To satisfy [START_REF] Lebeau | Contrôle exact de l'équation de la chaleur[END_REF], we define P k as the Lagrange interpolating polynomial at points λ k,j

with values (f k ω)[λ k,j ] := f k (λ k,j )ω k,j that is, in Newton form, P k (z) := g k j=1 (f k ω)[λ k,1 , . . . , λ k,j ] j-1 i=1 (z -λ k,i ).
Thus, to conclude it remains to estimate

R |J k (iτ )| 2 dτ . Estimate of J k .
Notice that since the eigenvalues in Λ are real, for any k ≥ 1 and any τ ∈ R, we

have |W k (iτ )| = 1. This implies (41) |J k (iτ )| ≤ g k (|τ | + λ k,g k ) p-1 (1 + τ 2 ) p/2 max j∈ 1,g k (f k ω)[λ k,1 , . . . , λ k,j ]
and thus

J k ∈ H 2 (C + ).
Using Leibniz formula (see Proposition 7.6),

(f k ω)[λ k,1 , . . . , λ k,j ] ≤ j i=1 f k [λ k,i , . . . , λ k,j ] ω[λ k,1 , . . . , λ k,i ] . (42)
Using again Leibniz formula (see Proposition 7.6), ( 43)

f k [λ k,i , . . . , λ k,j ] = j m=i f [λ k,i , . . . , λ k,m ] 1 W k [λ k,m , . . . , λ k,j ].
The two factors in each term of this sum are estimated using Lagrange theorem (see Proposition 7.4):

• First, we have

f [λ k,i , . . . , λ k,m ] = f (m-i) (θ k ) (m -i)! with θ k ∈ [λ k,i , λ k,m ]. It comes that there exists C p > 0 such that (44) f [λ k,1 , . . . , λ k,m ] ≤ C p (1 + λ k,m ) p ≤ C p (1 + λ k,g k ) p .
• Second, we have

1 W k [λ k,m , . . . , λ k,j ] = 1 (j -m + 1)! 1 W k (j-m+1) ( θk ) with θk ∈ [λ k,m , λ k,j ].
By using [START_REF] Miller | On the controllability of anomalous diffusions generated by the fractional Laplacian[END_REF], it follows that

(45) 1 W k [λ k,m , . . . , λ k,j ] ≤ C ε,p,ρ,N e ελ k,g k .
Recall that ( 14) implies λ k,g k -λ k,1 < ρ. Then, using ( 44) and ( 45) into the identity [START_REF] Powell | Approximation theory and methods[END_REF] proves that there exists C ε,p,r,ρ,N > 0 such that for any i, j ∈ 1, g k , i ≤ j, we have

f k [λ k,i , . . . , λ k,j ] ≤ C ε,p,r,ρ,N e ελ k,1 .
Plugging it in [START_REF] Ouaili | Minimal time of null controllability of two parabolic equations[END_REF] we obtain max

j∈ 1,g k |(f k ω)[λ k,1 , . . . , λ k,j ]| ≤ C ε,p,r,ρ,N e ελ k,1 max i∈ 1,g k ω[λ k,1 , . . . , λ k,i ] .
Finally, getting back to estimate [START_REF] Olive | Boundary approximate controllability of some linear parabolic systems[END_REF], and using the isomorphism property of L ends the proof of Proposition 2.1.

20

This manuscript is for review purposes only.

From infinite time horizon to finite time horizon.

In this section we first prove that the estimates on the solution on (0, +∞) of the block moment problem [START_REF] Jensen | Sur une expression simple du reste dans la formule d'interpolation de Newton[END_REF] for simple eigenvalues given in Proposition 2.1 implies the resolution on (0, T ) of the similar block moment problem [START_REF] Gal'chuk | Optimal control of systems described by parabolic equations[END_REF]. More precisely we prove the following.

proposition 2.3. Let p ∈ N * , ρ > 0 and N : (0, +∞) → R. Assume that Λ ∈ L w (p, ρ, N ). For any T > 0, there exists a constant C T,p,ρ,N > 0 such that for any q ∈ L 2 (0, +∞; C) there exists q ∈ L 2 (0, T ; C) satisfying

T 0 q(t)e -λt dt = +∞ 0 q(t)e -λt dt, ∀λ ∈ Λ, and
q L 2 (0,T ;C) ≤ C T,p,ρ,N q L 2 (0,+∞;C) .
For any T ∈ (0, +∞], we set

A(Λ, T ) := Span {e λ ; λ ∈ Λ} L 2 (0,T ;C)
, where e λ is defined in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF]. The proof of Proposition 2.3 mainly relies on the following proposition that gives an estimate on the inverse of the restriction operator.

proposition 2.4. Let p ∈ N * , ρ > 0 and N : (0, +∞) → R. Assume that Λ ∈ L w (p, ρ, N ).
Let T > 0 be fixed. Then, the restriction operator

(46) R Λ,T : q ∈ A(Λ, +∞) → q |(0,T ) ∈ A(Λ, T )
is an isomorphism. Moreover there exists a constant C T,p,ρ,N > 0 such that

(47) R -1 Λ,T ≤ C T,p,ρ,N .
In the case p = 1, this result is due to Fattorini and Russell [27, Theorem 1.3]. Our proof follows closely the strategy developed in this reference and takes advantage of the uniform estimates we established in the previous sections.

Proof. The fact that R Λ,T is an isomorphism is proved in [START_REF] Schwartz | Étude des sommes d'exponentielles réelles, NUMDAM[END_REF] under the sole assumption [START_REF] Avdonin | Riesz bases of exponentials and divided differences[END_REF]. The only thing to prove is thus the bound [START_REF] Schwartz | Étude des sommes d'exponentielles réelles, NUMDAM[END_REF].

The proof is done by contradiction. Assume that the estimate does not hold for given T , p, ρ, and N , then there exists a sequence (Λ m ) m≥1 belonging to the same

class L w (p, ρ, N ), such that (48) R -1 Λ m ,T ----→ m→∞ +∞.
For each m, by Proposition 7.1, we consider a grouping (

G m k ) k ∈ G(Λ m , p, ρ/p, N ),
and from [START_REF] Shackell | Overconvergence of Dirichlet series with complex exponents[END_REF] we know that there exists coefficients a m k,j such that the finite linear combination

P m : t → K m k=1 g m k j=1 a m k,j e -λ m k,j t , satisfies P m L 2 (0,∞;C) = 1, and P m L 2 (0,T ;C) ----→ m→∞ 0.

21

This manuscript is for review purposes only.

Let 0 < ε < T 2 be fixed and let C + 2ε = {z ∈ C ; (z) > 2ε}. We prove that the sequence z → P m (z) is uniformly bounded on any compact of C + 2ε .

Let m ≥ 1 and z ∈ C + 2ε . Then for any k ∈ {1, . . . , K m } the application of Proposition 2.1 to the sequence Λ m yields the existence of q m,z k ∈ L 2 (0, +∞; C) satisfying ( 49)

         +∞ 0 q m,z k (t)e -λ m k ,j t dt = 0, ∀k = k, ∀j ∈ 1, g m k , +∞ 0 q m,z k (t)e -λ m k,j t dt = e -λ m k,j z , ∀j ∈ 1, g m k ,
and

q m,z k L 2 (0,+∞;C) ≤ C ε,p,ρ,N e ελ m k,1 max j∈ 1,g m k e z [λ m k,1 , . . . , λ m k,j ] ,
where e z is defined in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF].

The previous right-hand side is estimated using Lagrange theorem (see Proposition 7.4). As the function e z is complex-valued we apply it on both its real and imaginary parts. This yields

max j∈ 1,g m k e z [λ m k,1 , . . . , λ m k,j ] ≤ C p,ρ |z| p e -(z)λ m k,1 .
Thus,

q m,z k L 2 (0,+∞;C) ≤ C ε,p,ρ,N |z| p e -( (z)-ε)λ m k,1 .
Then, using [START_REF] Tucsnak | Observation and control for operator semigroups[END_REF] it comes that, for m sufficiently large,

P m , q m,z k L 2 (0,∞;C) = K m k =1 g m k j=1 a m k ,j +∞ 0 e -λ m k ,j t q m,z k (t)dt = g m k j=1
a m k,j e -λ m k,j z .

From Cauchy-Schwarz inequality we deduce that

g m k j=1 a m k,j e -λ m k,j z ≤ P m L 2 (0,+∞;C) q m,z k L 2 (0,+∞;C) ≤ C ε,p,ρ,N |z| p e -( (z)-ε)λ m k,1 .
Summing these inequalities we obtain that for any

z ∈ C + 2ε , |P m (z)| ≤ K m k=1 g m k j=1 a m k,j e -λ m k,j z ≤ C ε,p,ρ,N |z| p k≥1 e -( (z)-ε)λ m k,1 ,
From the properties of the groupings (see Definition 1.2), it comes that λ m k,1 ≥ ρ p (k-1).

Thus, for any z

∈ C + 2ε , ( 50
) |P m (z)| ≤ C ε,p,ρ,N |z| p e -ρ p ( (z)-ε) k≥1 e -ε ρ p (k-2) ≤ C ε,p,ρ,N |z| p e -ρ p (z) .
This gives that (P m ) m is a sequence of holomorphic functions uniformly bounded on any compact of C + 2ε . From Montel's theorem it comes that we can extract a subsequence converging uniformly on any compact of C + 2ε to an holomorphic function P .

22

This manuscript is for review purposes only.

Now recall that P m L 2 (0,T ;C) goes to 0 as m goes to infinity. This implies that P (t) = 0 for any t ∈ (2ε, T ). The function P being holomorphic it comes that it vanishes on C + 2ε . Using (50) and the Lebesgue dominated-convergence theorem yields

P m L 2 (0,+∞;C) ----→ m→∞ 0.
This is in contradiction with P m L 2 (0,+∞;C) = 1 and ends the proof of Proposition 2.4.

We now have all the ingredients to prove Proposition 2.3. is a proper subspace of L 2 (0, +∞; C). Let Π Λ the associated orthogonal projection.

Proof
Let q ∈ L 2 (0, +∞, C). Then, by construction, we have

(51) +∞ 0 Π Λ q(t)e -λt dt = +∞ 0 q(t)e -λt dt, ∀λ ∈ Λ.
From Proposition 2.4, the restriction operator R Λ,T defined by ( 46) is an isomorphism.

Thus, setting q := (R -1 Λ,T ) * Π Λ q ends the proof of Proposition 2.3. Indeed, there exists

C T,p,ρ,N > 0 such that q L 2 (0,T ;C) ≤ C T,p,ρ,N q L 2 (0,+∞;C) ,
and, using (51), for every λ ∈ Λ,

T 0 q(t)e -λt dt = (R -1 Λ,T ) * Π Λ q, e λ L 2 (0,T ) = Π Λ q, R -1 Λ,T R Λ,T e λ L 2 (0,+∞) = +∞ 0 q(t)e -λt dt.
We can now conclude the proof of Theorem 2.1 for simple eigenvalues.

Proof (of Theorem 2.1). The resolution of the block moment problem [START_REF] Gal'chuk | Optimal control of systems described by parabolic equations[END_REF] as well as the estimate (31) follow directly from Propositions 2.1 and 2.3.

The only thing left to prove is the lower bound [START_REF] Jacob | Interpolation by vector-valued analytic functions, with applications to controllability[END_REF]. Let q k ∈ L 2 (0, T ; C) be any solution of (30b). Using the linearity of divided differences, equalities (30b) imply that for any i ∈ 1, g k

(52) ω[λ k,1 , . . . , λ k,i ] = T 0 q k (t)e t [λ k,1 , . . . , λ k,i ]dt
where e t is defined by [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF]. From the Lagrange theorem (see Proposition 7.4) and the fact that e (i-1) t is decreasing on (0, +∞), it comes that for any t ∈ (0, T ), we have

|e t [λ k,1 , . . . , λ k,i ]| ≤ t p e -λ k,1 t ≤ t p e -t ,
since we assumed that Λ ⊂ [1, +∞). Thus applying the Cauchy-Schwarz inequality to (52) gives

|ω[λ k,1 , . . . , λ k,i ]| ≤ q k L 2 (0,T ;C) T 0 t 2p e -t dt 1 2 ≤ (2p)! q k L 2 (0,T ;C)
which proves [START_REF] Jacob | Interpolation by vector-valued analytic functions, with applications to controllability[END_REF] and ends the proof of Theorem 2.1.

23

This manuscript is for review purposes only.

Construction of the control.

In this section we gather the previous ingredients to prove the positive controllability result. We also give an upper bound and a lower bound on the cost of controllability from Y 0 at time T as defined by Lemma 1.1.

Proof (of point i. of Theorem 1.1 for simple eigenvalues).

Assume that T 0 (Y 0 ) < +∞ and let us consider an initial data y 0 ∈ Y 0 . Without loss of generality we assume that y 0 -= 1.

Let T ∈ (T 0 (Y 0 ), +∞) and ε > 0 be such that T > T 0 (Y 0 ) + 2ε. From Proposition 1.3, it comes that T > T 0 (y 0 ) + 2ε.

For any k ≥ 1 and j ∈ 1, g k we set ω k,j := e -λ k,j T y 0 , ψ k,j -, .

Let (q k ) k≥1 be the solution of the block moment problem [START_REF] Gal'chuk | Optimal control of systems described by parabolic equations[END_REF] given in Theorem 2.1.

There exists a constant C ε,T,p,r,ρ,N > 0 such that

q k L 2 (0,T ;R) ≤ C ε,T,p,r,ρ,N e ελ k,1 max i∈ 1,g k ω[λ k,1 , . . . , λ k,i ] , ∀k ≥ 1. Let ξ k,j := y 0 , ψ k,j -, . Notice that ω k,j = e T (λ k,j )ξ k,j
, where e T is defined in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF].

From Leibniz formula (see Proposition 7.6), ( 53)

ω[λ k,1 , . . . , λ k,i ] = i m=1 e T [λ k,m , . . . , λ k,i ] ξ[λ k,1 , . . . , λ k,m ].
In this expression, ξ[. . . ] stands for the divided differences associated with the values

(λ k,1 , ξ k,1 ), . . . , (λ k,g k , ξ k,g k ). From Lagrange theorem (see Proposition 7.4) it comes that e T [λ k,m , . . . , λ k,i ] = e (i-m) T (θ k ) (i -m)! = (-T ) i-m (i -m)! e -θ k T with θ k ∈ [λ k,m , λ k,i ]. Using the definition (18) of T 0 (y 0 ), it comes that, (54) ξ[λ k,1 , . . . , λ k,m ] = y 0 , ψ[λ k,1 , . . . , λ k,m ] -, ≤ Ce λ k,1 (T0(y0)+ε) .
Thus, there exists C T,p > 0 such that

ω[λ k,1 , . . . , λ k,i ] ≤ C T,p e λ k,1 (T0(y0)+ε-T ) .
Then, as T > T 0 (y 0 ) + 2ε, the series ( 24) is convergent in L 2 (0, T ; R) and defines a control u that solves the moment problem [START_REF] Fattorini | Exact controllability theorems for linear parabolic equations in one space dimension[END_REF], which implies that the associated solution of (1) satisfies y(T ) = 0.

With the same strategy we can prove a more accurate result. Namely we get the following uniform bound for the cost of controllability.

(55) max l∈ 1,g k ψ[λ k,1 , . . . , λ k,l ] * ≤ C * e λ k,1 T 0 (Y 0 )+T 2 , ∀k ≥ 1.
Then, there exists a constant C T0(Y0),T,p,r,ρ,N > 0 such that for any y 0 ∈ Y 0 , there exists a control u ∈ L 2 (0, T ; R) such that the associated solution y of (1) satisfies y(T ) = 0 and u L 2 (0,T ;R) ≤ C T0(Y0),T,p,r,ρ,N C * y 0 -.

Proof. We follow the same strategy as in the proof of point i. of Theorem 1.1

with ε = T -T0(Y0) 4 but we do not use (54). Instead notice that using (55) we have

ξ[λ k,1 , . . . , λ k,m ] = y 0 , ψ[λ k,1 , . . . , λ k,m ] -, ≤ C * e λ k,1 T 0 (Y 0 )+T 2 y 0 -.
From (53) it comes that

ω[λ k,1 , . . . , λ k,m ] ≤ C T,p,r,ρ,N e -λ k,1 T C * e λ k,1 T 0 (Y 0 )+T 2 y 0 -.
Thus, writing that u ≤ k≥1 q k , we get

u L 2 (0,T ;R) ≤ C T,p,r,ρ,N C * y 0 - k≥1 e λ k,1 T 0 (Y 0 )-T 2 .
From Definition 1.2 it comes that λ k,1 ≥ r(k -1) which ends the proof of Corollary 2.2.

We also provide the following lower bound for the cost of controllability.

Corollary 2.3. Let p ∈ N * , r, ρ > 0 and N : (0, +∞) → R. Assume that Λ ∈ L w (p, ρ, N ) and let (G k ) k ∈ G(Λ, p, r, ρ) be an associated grouping. Let T > 0 and Y 0 ⊂ X -. Assume that η = 1. If system (1) is null-controllable from Y 0 at time T then, M (Y 0 , T ) ≥ sup k≥1 max l∈ 1,g k l! λ k,1 T l (e T ψ)[λ k,1 , . . . , λ k,l ] *
where e T is defined by (5), ψ is defined by [START_REF] Bernstein | Leçons sur les Progrès Récents de la Théorie des Séries de Dirichlet[END_REF] and M (Y 0 , T ) the cost of controllability from Y 0 at time T is defined in Lemma 1.1.

Proof. Let k ≥ 1 and l ∈ 1, g k . If system (1) is null-controllable at time T from Y 0 , we apply Lemma 1.1 with

z T := l j=1 φ k,j B * φ k,j i∈ 1,l =j (λ k,j -λ k,i ) .
By definition of z T we have

e -T A * z T = l j=1 e -λ k,j T φ k,j B * φ k,j i∈ 1,l =j (λ k,j -λ k,i ) ,
and thus

M (Y 0 , T ) 2 ≥ P * Y0 (e -T A * z T ) 2 * T 0 l j=1 e -λ k,j t i∈ 1,l =j (λ k,j -λ k,i ) 2 dt .
From Newton formula (see Proposition 7.3) and Lagrange theorem (see Proposition 7.4), for any t ∈ [0, T ], there exists

ξ t in [λ k,1 , λ k,l ] such that l j=1 e -λ k,j t i∈ 1,l =j (λ k,j -λ k,i ) ≤ t l l! e -ξtt ≤ T l l! e -λ k,1 t .
Using Newton formula (see Proposition 7.3), we have

P * Y0 (e -T A * z T ) = (e T ψ)[λ k,1 , . . . , λ k,l ],
which ends the proof of Corollary 2.3.

Estimates on Blaschke products.

This aim of this section is to prove the technical estimate stated in Proposition 2.2.

This relies on an extension of the following result by Fattorini and Russell.

Lemma 2.2 (see [START_REF] Fattorini | Uniform bounds on biorthogonal functions for real exponentials with an application to the control theory of parabolic equations[END_REF]Theorem 1.1]). Let γ > 0 and J : R + → R. Let L(γ, J ) be the class introduced in Remark 2.1. For any Σ ∈ L(γ, J ) and σ ∈ Σ we define

W Σ σ : z ∈ C + → σ ∈Σ σ =σ σ -z σ + z .
Then, for any ε > 0, there exists C ε,γ,J > 0 such that The generalisation we propose is the following.

W Σ σ (σ) ≥ C ε,γ,J e -
Lemma 2.3. Let γ > 0 and J : R + → R. For any ε > 0, there exists C ε,γ,J > 0 such that, for any sequence Σ ∈ L(γ, J ), for any σ ∈ Σ, we have

W Σ σ (z) ≥ C ε,γ,J e -εσ , ∀z ∈ C + , s.t. |z -σ| ≤ γ 2
Proof. For any σ > 0, since (σ

-(z)) 2 ≤ (σ + (z)) 2 , it comes that σ -z σ + z 2 = (σ -(z)) 2 + (z) 2 (σ + (z)) 2 + (z) 2 ≥ (σ -(z)) 2 (σ + (z)) 2 ,
and thus, (56)

|W Σ σ (z)| ≥ |W Σ σ ( (z))|.
We introduce the family Σ obtained from Σ by replacing σ by (z), that is

Σ := (Σ \ {σ}) ∪ { (z)}.
Since only one value has been modified, Σ also satisfies

σ∈ Σ 1 σ < +∞.
As, (57)

| (z) -σ| ≤ |z -σ| ≤ γ 2 ,
it comes that Σ satisfies the gap condition [START_REF] Fursikov | Controllability of evolution equations[END_REF] with ρ replaced by γ 2 . Notice that

{ (z)} ∈ L 1, ε → 1 ε .
Thus using Remark 2.2 and the arguments of the proof of Lemma 2.1 it comes that Σ ∈ L γ 2 , J with J depending only on J .

Obviously, as the terms σ ∈ Σ that are different from σ have not been modified it comes that

W Σ σ = W Σ (z) .
Applying Lemma 2.2 it comes that for any ε > 0, there is C ε,γ,J > 0 such that

|W Σ (z) ( (z))| ≥ C ε,γ,J e -ε (z) .
Finally, recalling (56) and (57), we obtain

|W Σ σ (z)| = |W Σ (z) (z)| ≥ |W Σ (z) ( (z))| ≥ C ε,γ,J e -ε (z) ≥ C ε,γ,J e -ε γ 2 e -εσ
which ends the proof of Lemma 2.3.

We now turn to the estimates we need for the derivatives of 1

W Σ σ .
proposition 2.5. Let γ > 0 and J : R + → R.

Then, for any l ≥ 0, for any ε > 0, there exists C l,ε,γ,J > 0 such that for any

Σ ∈ L(γ, J ), 1 W Σ σ (l) (σ) ≤ C l,ε,γ,J e εσ , ∀σ ∈ Σ.
Proof. The case l = 0 is nothing but the estimate given in Lemma 2.3.

Let D σ,γ := z ∈ C + ; |z -σ| ≤ γ 2 , C σ,γ := z ∈ C + ; |z -σ| = γ 2 .
As W Σ σ does not vanish in an open neighbourhood of D σ,γ it comes that 1

W Σ σ is
holomorphic on this domain. Thus applying Cauchy formula yields

1 W Σ σ (l) (σ) = l! 2iπ Cσ,γ 1 W Σ σ (z) (z -σ) l+1 dz.
From Lemma 2.3 it comes that for any ε > 0 there exists C ε,γ,J > 0 such that

1 W Σ σ (z) ≤ C ε,γ,J e εσ , ∀z ∈ C σ,γ .
This directly implies

1 W Σ σ (l) (σ) ≤ C ε,γ,J l! γ l e εσ
and ends the proof of Proposition 2.5.

27

This manuscript is for review purposes only.

We shall now move to the proof of Proposition 2.2 which is the main objective of this section.

Proof (of Proposition 2.2). Recall that the function N : R + → R is the one appearing in [START_REF] Fernández-Cara | Boundary controllability of parabolic coupled equations[END_REF] and that the subsequences Λ j are defined in [START_REF] López | Uniform null-controllability for the one-dimensional heat equation with rapidly oscillating periodic density[END_REF].

We recall that the index k is fixed, as well as the value θ ∈ Conv(G k ). We introduce the new sequence Λj obtained from Λ j by replacing the k-th value λ k,min(j,g k ) by θ i.e.

Λj := Λ j \{λ k,min(j,g k ) } ∪ {θ}.

Notice that, using Proposition 7.1, the fact that Λ j is a subsequence of Λ such that each term belong to a different group, and by the assumption on θ, we obtain that Λj satisfies the gap condition [START_REF] Fursikov | Controllability of evolution equations[END_REF] with ρ replaced by γ = ρ p . Notice that {θ} ∈ L 1, ε → 1 ε . Thus using Remark 2.2 and the arguments of the proof of Lemma 2.1 it comes that Λj ∈ L γ, J with J depending only on N .

With these notations and Proposition 7.1 it comes that

1 W k (z) = p j=1 1 W Λj θ (z)
.

Finally, using Leibniz rule (for derivatives), evaluating the result at z = θ and using Proposition 2.5 yield the claim.

Lack of null-controllability in small time.

The goal of this section is to prove the point ii. of Theorem 1.1 in the case of algebraically simple eigenvalues. Thus, in all this section we assume that η = 1. The proof mainly relies on the optimality on the bound obtained in the resolution of the block moments problem in [START_REF] Jacob | Interpolation by vector-valued analytic functions, with applications to controllability[END_REF].

Proof. Let T > 0 and assume that null controllability from Y 0 in time T of (1) holds.

Thus, there exists C T > 0 such that: for any y 0 ∈ X -, there exists u ∈ L 2 (0, T ; R) such that the associated solution of (1) with initial condition P Y0 y 0 satisfies y(T ) = 0 and u L 2 (0,T ;R) ≤ C T y 0 -.

Due to the equivalence between null controllability and the moment problem [START_REF] Duyckaerts | Resolvent conditions for the control of parabolic equations[END_REF] it comes that

T 0 e -λ(T -t) u(t)dt = -y 0 , e -λT ψ λ -, , ∀λ ∈ Λ.
Recall that ψ λ is defined by [START_REF] Bernstein | Leçons sur les Progrès Récents de la Théorie des Séries de Dirichlet[END_REF]. Thus, for any k ∈ N * , the control u(T -•) solves (30b). Using the lower bound [START_REF] Jacob | Interpolation by vector-valued analytic functions, with applications to controllability[END_REF] with ω k,j := y 0 , e -λ k,j T ψ k,j -, , it comes that

ω[λ k,1 , . . . , λ k,l ] ≤ C p u L 2 (0,T ;R) ≤ C p,T y 0 -, ∀k ∈ N * , ∀l ∈ 1, g k .
Due to the definition of ω k,j , this can be rewritten as

y 0 , (e T ψ)[λ k,1 , . . . , λ k,l ] -, ≤ C p,T y 0 -, ∀k ∈ N * , ∀l ∈ 1, g k ,
where e t is defined by [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF]. By the dual characterization of the norms this implies

(58) (e T ψ)[λ k,1 , . . . , λ k,l ] * ≤ C p,T , ∀k ∈ N * , ∀l ∈ 1, g k ,
Notice that ψ k,j = e λ k,j T e -λ k,j T ψ k,j = e -T [λ k,j ](e T ψ)[λ k,j ]. Thus, using Leibniz formula (see Proposition 7.6), we obtain, (59)

ψ[λ k,1 , . . . , λ k,l ] = l j=1 (e T ψ)[λ k,1 , . . . , λ k,j ]e -T [λ k,j , . . . , λ k,l ], ∀k ∈ N * , ∀l ∈ 1, g k .
From Proposition 7.4, it comes that for any j, l ∈ 1, g k , there exists z

∈ [λ k,j , λ k,l ] such that |e -T [λ k,j , . . . , λ k,l ]| = e (l-j) -T (z) (l -j)! ≤ C p e λ k,j T ≤ C p,ρ,T e λ k,1 T .
Finally, plugging this estimate and (58) into (59) we obtain,

ψ[λ k,1 , . . . , λ k,l ] * ≤ C p,ρ,T e λ k,1 T .
Due to the definition of T 0 (Y 0 ), this implies that T ≥ T 0 (Y 0 ) and ends the proof on the lack of null-controllability at time T < T 0 (Y 0 ).

3.

Comparison with some already known results.

In this section, we prove that we actually recover the known formulas for the minimal null-control time when there is no condensation of eigenvalues or when the eigenvectors are assumed to form a Riesz basis of X * . Doing so we will highlight in Proposition 3.2 that the actual minimal null-control time is always smaller than the value predicted by the formula that would be valid under the Riesz basis assumption.

As all these results were proved for albreaically simple eigenvalues we assume in all this section that η = 1.

Notice that the proofs in all this section only rely on the definition of the minimal null-control time [START_REF] Coron | Control and nonlinearity[END_REF] and thus do not depend on Theorem 1.1.

3.1.

When there is no condensation of eigenvalues.

In this section we prove that, if the condensation index of the sequence Λ vanishes

T 0 (X -) = lim sup λ→∞ λ∈Λ -ln |B * φ λ | λ .
This result was already proved in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF] with the additional assumption that the family of eigenvectors Φ = (φ λ ) λ∈Λ forms a Riesz basis of X * or in [START_REF] Khodja | Quantitative Fattorini-Hautus test and minimal null control time for parabolic problems[END_REF]Remark 1.15] in a more general framework encompassing the one studied here.

Proof. Notice that when Y 0 = X -, the operator P * Y0 reduces to the identity.

Thus, considering l = 1 in (19) always lead to

T 0 (X -) ≥ lim sup λ→∞ λ∈Λ -ln |B * φ λ | λ .

29

This manuscript is for review purposes only.

We assume that (60)

T 0 (X -) > lim sup λ→∞ λ∈Λ -ln |B * φ λ | λ ,
and we will prove that c(Λ) > 0.

We shall reason as in the proof of point ii. of Theorem 1.1 (see Section 2.2) but starting with the formula (17) instead of [START_REF] Boyer | Approximate controllability conditions for some linear 1D parabolic systems with space-dependent coefficients[END_REF]. We can find an integer l * ≥ 1, an extraction (κ n ) n≥1 and integers m n such that 1 ≤ m n ≤ m n + l * -1 ≤ g κn and such that if

x n := ψ[λ κn,mn , . . . , λ κn,mn+l * -1 ],

we have

lim n→∞ ln x n * λ κn,1 = T 0 (X -).
Moreover, we can assume that for any l ∈ 1, l * -1 , we have for some ε > 0

T (l) := lim sup n→+∞ ln     max m,r∈ 1,gκ n m≤r r-m<l ψ[λ κn,m , . . . , λ κn,r ] *     λ κn,1 < T 0 (X -) -ε,
since, if it is not the case, we can reduce the value of l * accordingly. Note that, as

P * Y0 φ λ * ≤ 1, by (60), we know that l * > 1.
From the definition of divided differences (see Definition 7.1), it comes that

x n = ψ[λ κn,mn+1 , . . . , λ κn,mn+l * -1 ] -ψ[λ κn,mn , . . . , λ κn,mn+l * -2 ] λ κn,mn+l * -1 -λ κn,mn .

For n sufficiently large, we have

x n * ≥ e λκ n ,1(T0(X-)-ε/2) .

Using the definition of T (l * -1) it comes that, for n large enough, ψ[λ κn,mn+1 , . . . , λ κn,mn+l * -1 ] * + ψ[λ κn,mn , . . . , λ κn,mn+l * -2 ] * ≤ e λκ n,1 ( T (l * -1)+ε/2) .

Thus, since l * ≥ 2, we can combine the last two estimates to obtain

|λ κn,mn+1 -λ κn,mn | ≤ |λ κn,mn+l * -1 -λ κn,mn | ≤e -λκ n ,1 (T0(X-)-ε-T (l * -1))
≤e ρ(T0(X-)-ε-T (l * -1)) e -λκ n ,mn (T0(X-)-ε-T (l * -1)) .

In particular, we have

lim sup n→∞ -ln |P Gκ n (λ κn,m )| λ κn,mn ≥ lim sup n→∞ -ln |λ κn,mn+1 -λ κn,mn | λ κn,mn ≥ T 0 (X -) -ε -T (l * -1) > 0.
Using Proposition 7.12, we conclude that c(Λ) > 0, and the claim is proved.

30

This manuscript is for review purposes only.

When there is a Riesz basis of eigenvectors.

As already mentioned the null-control problem for (1) has been considered in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF] with the additional assumption that the family (φ λ ) λ∈Λ forms a Riesz basis of X * .

Observe that it is equivalent to ask that (φ λ / φ λ ) λ∈Λ is a Riesz basis of X.

With this additional assumption, the minimal null-control time from Y 0 = X - was proved to be equal to (61)

T * := lim sup λ→∞ λ∈Λ ln 1 |B * φ λ | |E Λ (λ)| λ ,
where the interpolating function E Λ is defined in (126).

Remark 3.1. Notice that, since φ λ is normalized in X * , there exists C > 0 such that

1 Cλ ≤ φ λ ≤ C, ∀λ,
so that the value of T * in (61) does not change if one considers the normalization of eigenvectors in X instead of in X * .

In our setting, we prove that the formula above for T * is always an upper bound of the actual minimal null-control time, without assuming the Riesz basis condition.

proposition 3.2. Assume that A and B satisfy the assumptions of Theorem 1.1 with η = 1. Then, T 0 (X -) ≤ T * where T * is defined by (61).

Proof. First step: we begin by proving that the grouping designed in Proposition 7.1 ensures a simpler expression for T * . Let (G k ) k≥1 ∈ G(Λ, p, r, ρ) be a grouping as introduced in Section 1.3.2. For each λ ∈ Λ, we denote by G [λ] the unique group in (G k ) k≥1 that contains λ. Then, we have (62)

T * = lim sup λ→∞ λ∈Λ ln 1 |B * φ λ | |P G [λ] (λ)| λ ,
where, for each group G, the polynomial P G is defined by [START_REF] Khodja | New phenomena for the null controllability of parabolic systems: minimal time and geometrical dependence[END_REF].

• Let G be a group of eigenvalues and λ ∈ G. We prove that, for any finite subset M of Λ \ G, whose cardinal is denoted by n := #M , we have

(63) µ∈M |λ -µ| ≥ r n n 2p !
where r := ρ p . To this end, for any j ∈ 1, p , we define

M j := {µ ∈ M ; ∃k ≥ 1 such that µ = λ k,j }.
Since the groups are covering Λ, we have a disjoint union M = p j=1 M j . It follows that there exists j 0 ∈ 1, p such that #M j0 ≥ n p . From (13) it comes that

|λ -µ| ≥ r, ∀µ ∈ M, and
|µ -µ | ≥ r, ∀j ∈ 1, p , ∀µ, µ ∈ M j , µ = µ .

31

This manuscript is for review purposes only.

Then,

µ∈M |λ -µ| =     p j=1 j =j0 µ∈Mj |λ -µ|       µ∈Mj 0 |λ -µ|   ≥ r #(M \Mj 0 ) r #(Mj 0 ) #M j0 2 !
This proves (63).

• From (63) we apply [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF]Theorem 3.8] to obtain that for any subsequence

(λ n ) n≥1 ⊂ Λ, (64) lim n→∞   ln 1 |E Λ (λn)| λ n - ln 1 |P G [λn ] (λn)| λ n   = 0.
This directly implies (62).

Remark 3.2. Notice that (63) is not the exact assumption required in [5, Theorem 3.8]. For this result the authors assumed

(65) µ∈M |λ -µ| ≥ r n n!,
with the same notation as in the proof above. We claim that with the exact same proof it is sufficient to assume (63). Indeed, in the proof of [5, Theorem 3.8], the only point were assumption (65) is used is the Second step in the middle of page 2097. Then the term n! is estimated asymptotically using Stirling formula to prove that the term Γ k,1 goes to 0 as k goes to ∞. As the rest of the proof is long, technical and remains unchanged when replacing (65) by (63) we do not reproduce it here for the sake of brevity.

Second step: we end the proof of Proposition 3.2. Recall that from [START_REF] Jacob | Weighted multiple interpolation and the control of perturbed semigroup systems[END_REF] we have that there exists C p,ρ > 0 such that for any k ≥ 1, l ∈ 1, g k and any j ∈ 1, l ,

i∈ 1,l =j |λ k,j -λ k,i | ≥ C p,ρ |P G k (λ k,j )|.
As we have considered normalized eigenvectors, and by ( 4), for any k ≥ 1 and any l ∈ 1, g k , we have

l j=1 ψ k,j i∈ 1,l =j (λ k,j -λ k,i ) * ≤ l j=1 φ k,j B * φ k,j i∈ 1,l =j (λ k,j -λ k,i ) * ≤ l max j∈ 1,l φ k,j B * φ k,j i∈ 1,l =j (λ k,j -λ k,i ) * ≤ l max j∈ 1,l 1 |B * φ k,j | i∈ 1,l =j |λ k,j -λ k,i | . 32
This manuscript is for review purposes only.

Using [START_REF] Jacob | Weighted multiple interpolation and the control of perturbed semigroup systems[END_REF] this leads to

l j=1 ψ k,j i∈ 1,l =j (λ k,j -λ k,i ) * ≤ Cl max j∈ 1,l 1 |B * φ k,j ||P G k (λ k,j )| .
Thus, ln max

l∈ 1,g k   l j=1 ψ k,j i∈ 1,l =j (λ k,j -λ k,i ) *   λ k,1 ≤ max j∈ 1,g k ln 1 |B * φ k,j ||P G k (λ k,j )| λ k,j λ k,j λ k,1 + ln(Cl) λ k,1 .
Then, using (62), we obtain

T 0 (X -) ≤ T * .
We now prove that we indeed recover exactly the expression of the minimal time (61) (or (62)) when we assume that the eigenvectors form a Riesz basis.

proposition 3.3. Assume that A and B satisfy the assumptions of Theorem 1.1 with η = 1 and that (φ λ ) λ∈Λ forms a Riesz basis of X * . Then, T 0 (X -) = T * where T * is defined by (62).

Remark 3.3. It will appear clearly in the proof that the Riesz basis assumption is much stronger than what we really need. The only thing that we actually use at the very beginning of the proof, is that the spectral radius of the inverse of the Gram matrix M

k := Gram X * (φ k,1 , . . . , φ k,g k ) satisfies sup k≥1 ρ(M -1 k ) < +∞.
A careful inspection of the proof shows that it is in fact sufficient to assume that

lim sup k→∞ ln ρ(M -1 k ) λ k,1 = 0.
Note in particular that, in practice, estimating such a spectral radius in each group is much simpler than proving that the whole family is a Riesz basis.

Proof. As we assumed that (φ λ ) λ∈Λ is a Riesz basis of X * it comes that there exists C > 0 such that for any k ≥ 1, for any

α k,1 , . . . , α k,g k ∈ R, max j∈ 1,g k |α k,j | ≤   g k j=1 α 2 k,j   1 2 ≤ C g k j=1 α k,j φ k,j * , and thus max j∈ 1,g k |α k,j | ≤ C max l∈ 1,g k l j=1 α k,j φ k,j * . Setting α k,j := 1 B * φ k,j i∈ 1,g k =j (λ k,j -λ k,i ) , 33
This manuscript is for review purposes only.

yield max j∈ 1,g k 1 |B * φ k,j | i∈ 1,g k =j |λ k,j -λ k,i | ≤ C g k j=1 φ k,j B * φ k,j i∈ 1,g k =j (λ k,j -λ k,i ) * ≤ C max 1≤l≤g k    l j=1 φ k,j B * φ k,j i∈ 1,l =j (λ k,j -λ k,i ) *    .
It follows that for any j ∈ 1, g k ,

1 |B * φ k,j ||P G k (λ k,j )| = 1 |B * φ k,j | i∈ 1,g k =j |λ k,j -λ k,i | ≤ C max 1≤l≤g k    l j=1 φ k,j B * φ k,j i∈ 1,l =j (λ k,j -λ k,i ) *    .
Thus, taking the logarithm,

-ln |B * φ k,j ||P G k (λ k,j )| λ k,j ≤ ln max l∈ 1,g k   l j=1 φ k,j B * φ k,j i∈ 1,l =j (λ k,j -λ k,i ) *   λ k,1 + ln C λ k,1 .
Since by definition we have G k = G [λ k,j ] , this ends the proof of Proposition 3.3.

The case of multiple eigenvalues.

In this section we prove Theorem 1.1 in the case where we allow algebraic multiplicity for the eigenvalues i.e. η ≥ 2. As previously, the main issue is the resolution of the block moment problem given in [START_REF] Fattorini | Some remarks on complete controllability[END_REF]. This is detailed in the next subsection.

Resolution of block moment problems.

We prove in this subsection the following theorem which is the generalization of Theorem 2.1. Assume that Λ ∈ L w (p, ρ, N ) and let (G k ) k ∈ G(Λ, p, r, ρ) be an associated grouping.

We also consider an integer η ≥ 1.

For any ε > 0, there exists a constant C ε,T,p,r,ρ,η,N > 0 such that for any k ≥ 1,

for any multi-index α k ∈ N g k with |α k | ∞ ≤ η, any set of values ω α k ∈ C |α k | , there exists q k ∈ L 2 (0, T ; C) satisfying T 0 q k (t) (-t) l l ! e -λ k ,j t dt = 0, ∀k = k, ∀j ∈ 1, g k , ∀l ∈ 0, η , (66a) T 0 q k (t) (-t) l l! e -λ k,j t dt = ω l k,j , ∀j ∈ 1, g k , ∀l ∈ 0, α k,j -1 , (66b)
and the bound (67) q k L 2 (0,T ;C) ≤ C ε,T,p,r,ρ,η,N e ελ k,1 max

µ∈N g k µ≤α k ω[λ (µ1) k,1 , . . . , λ (µ k ) k,g k ] .
Moreover, up to the factor e ελ k,1 , this last estimate is sharp: any solution q k of (66b),

satisfy (68) q k L 2 (0,T ;C) ≥ C p,η max µ∈N g k µ≤α ω[λ (µ1) k,1 , . . . , λ (µ k ) k,g k ] ,
for some C p,η > 0.

In the case p = 1 (usual gap condition), a solution to (66) is given by the biorthogonal family built in [START_REF] Khodja | The Kalman condition for the boundary controllability of coupled parabolic systems. Bounds on biorthogonal families to complex matrix exponentials[END_REF]. Here, we extend this resolution using a weak gap condition (9) and we prove that the obtained estimates are uniform with respect to Λ in a given class Assume that Λ ∈ L w (p, ρ, N ) and let (G k ) k ∈ G(Λ, p, r, ρ) be an associated grouping.

L w (•, •, •).
We consider an integer η ≥ 1 and for any k we suppose given a multi-index

α k ∈ N g k such that |α k | ∞ ≤ η.
Then, for any k ≥ 1, for any j ∈ 1, g k and any l ∈ 0, α k,j -1 there exists

q k,j,l ∈ L 2 (0, T ; C) satisfying T 0 q k,j,l (t) (-t) l l ! e -λ k ,j t dt = δ k,k δ j,j δ l,l ,
for any k, k ≥ 1, any j ∈ 1, g k , j ∈ 1, g k and any l ∈ 0, α k,j -1 , l ∈ 0, α k ,j -1 . Moreover, for any ε > 0, there exists a constant C ε,T,p,r,ρ,η,N > 0 such that for any k ≥ 1, any j ∈ 1, g k and any l ∈ 0, α k,j -1 , we have

q k,j,l L 2 (0,T ;C) ≤ C ε,T,p,r,ρ,η,N e ελ k,1 i∈ 1,g k =j |λ k,j -λ k,i | α k,i 1 min i∈ 1,g k =j |λ k,j -λ k,i | α k,j -l-1 .
The proof of Corollary 4.1 is left to the reader: it follows closely the one of Corollary 2.1 and makes use of the estimate given in Proposition 7.11 instead of the Newton formula for standard divided differences.

Remark 4.1. Contrary to the estimate in Corollary 2.1, the above estimate is not optimal in general, even if we do not consider the exponential factor. Indeed, some cancellations can occur depending on the relative positions and multiplicities of the eigenvalues that are not taken into account in the above general bound. In actual examples, one needs to compute carefully the coefficients of the generalized divided differences introduced in Proposition 7.11 to see whether or not a sharper estimate can be obtained.

Here also, the proof of Theorem 4.1 relies on the resolution of the block moment problem (66) with T = +∞ and then on a restriction argument. For pedagogical resaons (the proof being less technical) let us present first this restriction argument (which is the generalization of Proposition 2.3).

35

This manuscript is for review purposes only. For any T > 0, there exists a constant C T,p,r,ρ,η,N > 0 such that for any q ∈ L 2 (0, +∞; C), there exists q ∈ L 2 (0, T ; C) satisfying

T 0 q(t) (-t) l l! e -λt dt = +∞ 0 q(t) (-t) l l! e -λt dt, ∀λ ∈ Λ, ∀l ∈ 0, η ,
and the estimate q L 2 (0,T ;C) ≤ C T,p,r,ρ,η,N q L 2 (0,+∞;C) .

Proof. For any h > 0, we define

Λ h := η l=0 (Λ + lh).
Using Remark 2.2 and Lemma 2.1 we have that Λ h ∈ L w (pη, ρ, Ñ ) for some Ñ which does not depend on h. We suppose given a fixed q and, for any h > 0, we can apply Proposition 2.3 with the sequence Λ h and obtain the existence of a function

q h ∈ L 2 (0, T ; C) such that (69) T 0 q h (t)e -(λ+lh)t dt = +∞ 0 q(t)e -(λ+lh)t dt, ∀λ ∈ Λ, ∀l ∈ 0, η ,
and satisfying moreover the uniform estimate q h L 2 (0,T ;C) ≤ C T,pη,r,ρ, Ñ q L 2 (0,+∞;C) , ∀h > 0.

We can then find a subsquence (q hn ) n that weakly converges towards some q ∈ L 2 (0, T ; C) such that q L 2 (0,T ;C) ≤ C T,pη,r,ρ, Ñ q L 2 (0,+∞;C) . We will show that q solves the required equations.

Let λ ∈ Λ and l ∈ 0, η -1 be fixed. Combining the equations (69) to make appear divided differences, we have the equality (70) Let us now turn to the resolution of the block moment problem (66) for T = +∞.

The next proposition is the generalization of Proposition 2.1.

proposition 4.2. Let p ∈ N * , r, ρ > 0 and N : (0, +∞) → R. Assume that Λ ∈ L w (p, ρ, N ) and let (G k ) k ∈ G(Λ, p, r, ρ) be an associated grouping. We also consider an integer η ≥ 1.

For any ε > 0, there exists a constant C ε,p,r,ρ,η,N > 0 such that for any k ≥ 1, for

any multi-index α k ∈ N g k with |α k | ∞ ≤ η, and any set of values ω α k ∈ C |α k | , there exists q k ∈ L 2 (0, +∞; C) satisfying          +∞ 0 q k (t) (-t) l l ! e -λ k ,j t dt = 0, ∀k = k, ∀j ∈ 1, g k , ∀l ∈ 0, η , +∞ 0 q k (t) (-t) l l! e -λ k,j t dt = ω l k,j , ∀j ∈ 1, g k , ∀l ∈ 0, α k,j -1 ,
and the bound

(71) q k L 2 (0,+∞;C) ≤ C ε,p,r,ρ,η,N e ελ k,1 max µ∈N g k µ≤α k ω[λ (µ1) k,1 , . . . , λ (µ k ) k,g k ] .
Before getting to the proof let us mention that Propositions 4.1 and 4.2 imply Theorem 4.1. The lower bound (68) is proved in the exact same way as [START_REF] Jacob | Interpolation by vector-valued analytic functions, with applications to controllability[END_REF] and is thus left to the reader.

Proof. As in the previous proof, for h > 0, we define

Λ h := η l=0 (Λ + lh),
that belongs to the class L w (pη, ρ, Ñ ). For any k ≥ 1, we set

G k,h := η l=0 (G k + lh).
For any h < r/(2η), the family (G k,h ) k is a grouping in G(Λ h , pη, r/2, ρ + r/2). Now, we are given a fixed index k. We observe that, there exists a h 0 ∈ (0, r/(2η))

(possibly depending on k) such that, for any h < h 0 , the sets G k , G k + h, . . . , G k + ηh are pairwise disjoint.

Since we need to take into account precisely the multiplicities we are interested in, encoded in the multi-index α k , we introduce the modified k-th group

G k,h = g k j=1 {λ k,j , λ k,j + h, . . . , λ k,j + (α k,j -1)h} ⊂ G k,h ,
and the new family

Λh =     l≥1 l =k G l,h     ∪ G k,h ,
which satisfies Λh ⊂ Λ h and therefore also belongs to the class L w (pη, ρ, Ñ ).

By construction, the family of points in G k,h , that we denote by µ

k,h,1 < • • • < µ k,h,|α k |
is an approximation of the weighted family ((λ k,1 , . . . , λ k,g k ), α k ) in the sense of Definition 7.2. Let F : R → C be a smooth function satisfying the conditions

(72) 1 l! F (l) (λ k,j ) = ω l k,j , ∀j ∈ 1, g k , ∀l ∈ 0, α k,j -1 .
For each h > 0, we apply Proposition 2.1 to the family Λh to find a solution q k,h ∈ L 2 (0, +∞; C) to the following moment problem (73)

         +∞ 0 q k,h (t)e -(λ k ,j +hl)t dt = 0, ∀k = k, ∀j ∈ 1, g k , ∀l ∈ 0, η , +∞ 0 q k,h (t)e -(λ k,j +hl)t dt = F (λ k,j + hl), ∀j ∈ 1, g k , ∀l ∈ 0, α k,j -1 ,
and satisfying the following bound, with a constant uniform with respect to h,

q k,h L 2 (0,+∞;C) ≤ C ε,ηp,r,ρ, Ñ e ελ k,1 max i∈ 1,|α k | F [µ k,h,1 , . . . , µ k,h,i ] .
By Proposition 7.7, we know that the right-hand side in the above estimate converges when h → 0 towards a similar quantity with generalized divided differences instead of the usual divided differences. It follows that we can extract a subsequence (q k,hn ) n that weakly converges in L 2 (0, +∞; C) towards a function q k that satisifies the bound (71).

Finally, by the same argument as in the proof of Proposition 4.1 above, we can combine the equations (73) to make appear divided differences on both side and pass to the weak-strong limit in the integral to finally get

         +∞ 0 q k (t) (-t) l l! e -λ k ,j t dt = 0, ∀k = k, ∀j ∈ 1, g k , ∀l ∈ 0, η , +∞ 0 q k (t) (-t) l l! e -λ k,j t dt = F [λ (l) k,j ], ∀j ∈ 1, g k , ∀l ∈ 0, α k,j -1 ,
which is exactly our claim since, by the computation rule (121) and by (72), we have

F [λ (l) k,j ] = ω l k,j .

Proof of the minimal null-control time property.

In this section we end the proof of Theorem 1.1. The extension of Corollaries 2.2 and 2.3 as well as their proofs to the case η ≥ 2 are straightforward and left to the reader.

Controllability in large time: proof of point i. of Theorem 1.1.

Let T > T 0 (Y 0 ) and y 0 ∈ Y 0 . For any k ≥ 1, let q k ∈ L 2 (0, T ; C) be given by Theorem 4.1 with

ω l k,j := y 0 , (e T ψ)[λ (l+1) k,j ] -,
.

As in Section 2.1.3, since T > T 0 (Y 0 ), the estimates (67) imply that

u := - k≥1 q k (T -•) ∈ L 2 (0, T ; C).
Moreover, as q k solves the block moment problem (66) it comes that u solves the moment problem [START_REF] Duyckaerts | Resolvent conditions for the control of parabolic equations[END_REF] and thus y(T ) = 0.

Lack of null-controllability in small time: proof of point ii. of Theorem 1.1.

The proof follows exactly the lines of Section 2.2 and relies on the lower bound (68)

given for the solution of the block moments problem (66).

Examples.

In The design of these abstract examples is inspired from the work [START_REF] Khodja | Null-controllability of some systems of parabolic type by one control force[END_REF]. Our goal is to illustrate, in particular, the fact that, even if the control operator has no influence on the minimal null-control time, the knowledge of the condensation index of the eigenvalues of the operator A is not sufficient to understand the null-controllability properties of system (1).

Let A be a positive definite self-adjoint operator with compact resolvent in a

Hilbert space H whose eigenvalues (µ k ) k≥1 are assumed to be sorted in increasing order. One can think of A, for instance, as the Laplace operator -∂ xx or any Sturm-Liouville operator with homogeneous Dirichlet boundary conditions.

If we denote by (ϕ k ) k≥1 a corresponding Hilbert basis of eigenvectors, A may be written

A = k≥1 µ k (•, ϕ k ) H ϕ k , D(A) =    x ∈ H ; k≥1 µ 2 k (x, ϕ k ) 2 H < +∞    ,
where (•, •) H denotes the scalar product in H. We assume that (µ k ) k≥1 satisfies ( 8) and ( 9) with p = 1, i.e., satisfies the so-called gap property. Let ρ > 0 be such that

(74) 0 < ρ < inf k≥1 (µ k+1 -µ k )
and f : σ(A) → R a positive function defined on σ(A) the spectrum of A satisfying

(75) 0 < f (µ k ) < ρ, ∀k ≥ 1.
Let f (A) be the operator defined on D(A) by

f (A) := k≥1 f (µ k ) (•, ϕ k ) H ϕ k .

39

This manuscript is for review purposes only.

Let x 0 ∈ H fixed satisfying (76) | (x 0 , ϕ k ) H | ≥ e -√ µ k , ∀k ≥ 1.
Remark 5.1. This vector x 0 will be used to design the control operator B. This assumption will ensure that the terms B * φ λ appearing in the definition [START_REF] Coron | Control and nonlinearity[END_REF] have no influence. This will allow us to really emphasize the role of the condensation of eigenvectors.

5.1.1. Perturbation of a 2 × 2 Jordan block.

Let X = H × H, (77) A = A I 0 A + f (A) , D(A) = D(A) × D(A),

and

(78)

B : u ∈ R → u 0 x 0 ∈ X.
It is easy to see that (-A, D(A)) generates a C 0 -semigroup on X and that B : R → X is bounded. Thus we consider for this example that X * = X = X -and Y 0 = X.

The spectrum of (A * , D(A)) is given by i. For any function f satisfying (75), null-controllability from X holds in any time i.e. T 0 (X) = 0.

Λ = {µ k , µ k + f (µ k ) ; k ≥ 1}.
ii. For any τ ∈ [0, +∞], there exists a function f satisfying (75) such that c(Λ) = τ . This gives a first example in this setting where the minimal time is not related to the condensation index. As it will appear from the proof, see (79), this is due to a condensation of eigenvectors compensating the condensation of eigenvalues.

Proof. The proof of point ii. directly follows from straightforward computations using Proposition 7.12 with the explicit choices f :

s → ρe - √ s , f : s → ρe -cs with c > 0 or f : s → ρe -s 2 .
We now turn to the computation of the minimal null-control time. Using ( 74) and ( 75), it comes that ( 8) and ( 9) are satisfied with p = 2. We define our grouping by setting λ k,1 := µ k and λ k,2 := µ k + f (µ k ). The associated normalized eigenvectors are

φ k,1 := 1 1 + f (µ k ) 2 -f (µ k ) 1 ϕ k , φ k,2 := 0 1 ϕ k ,
which do form a complete family in X. Moreover, for all k ≥ 1,

B * φ k,1 = 1 1 + f (µ k ) 2 (x 0 , ϕ k ) H , and B * φ k,2 = (x 0 , ϕ k ) H , 40
This manuscript is for review purposes only.

so that, with (15), we have

ψ k,1 = 1 (x 0 , ϕ k ) H -f (µ k ) 1 ϕ k , ψ k,2 = 1 (x 0 , ϕ k ) H 0 1 ϕ k .
From Definition [START_REF] Boyer | Approximate controllability conditions for some linear 1D parabolic systems with space-dependent coefficients[END_REF], we have

T 0 (X) = lim sup k→∞ ln max ψ k,1 , ψ k,2 -ψ k,1 f (µ k ) µ k .
One has

(79) ψ k,2 -ψ k,1 f (µ k ) = 1 | (x 0 , ϕ k ) H | 1 0 ϕ k = 1 | (x 0 , ϕ k ) H | .
Using ( 76) and ( 75) we easily deduce that

T 0 (X) = lim sup k→∞ ln 1 |(x0,ϕ k ) H | µ k = 0. Remark 5.2. Notice that, φ k,2 -φ k,1 2 = 2 1 + f (µ k ) 2 -1 + f (µ k ) 2 1 + f (µ k ) 2 -→ k→∞ 0,
thus the eigenvectors of A * do not form a Riesz basis of X. If this family were a Riesz basis, then we would deduce from [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF] that the minimal null-control time would be equal to the condensation index c(Λ). In this case, the operator A * has spectrum σ(A * ) = {µ k ; k ≥ 1} with algebraically double eigenvalues satisfying the gap property and an associated Hilbert basis of (generalized) eigenvectors given by

φ 0 k = 0 1 ϕ k , and
φ 1 k = 1 0 ϕ k .
Notice that from (79) one has

ψ k,2 -ψ k,1 f (µ k ) = 1 (x 0 , ϕ k ) H 1 0 ϕ k = φ 1 k B * φ 0 k .
Thus, the analysis of (77)-( 78), is unchanged if ones sets f = 0.

Algebraically multiple eigenvalues.

Let X = H × H × H. Let β > 0 and g : σ(A) → R be such that

(80) g(µ k ) = ρe -βµ k , with ρ satisfying (74). Let (81) A =   A I 0 0 A 0 0 0 A + g(A)   , D(A) = D(A) × D(A) × D(A),

and

(82)

B : u ∈ R → u   0 x 0 x 0   .
Again B is a bounded control operator and we also set for this example X * = X = X - and Y 0 = X.

The spectrum of (A * , D(A)) is given by Remark 5.4. In this case, the family of (generalized) eigenvectors do form a Hilbert basis in X. However due to the presence of algebraically multiple eigenvalues one cannot compute the value of the minimal null-control time using [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF]. Its value is still related to the condensation index of Λ but also depends on the multiplicity of each eigenvalue in the system.

Λ = {µ k , µ k + g(µ k ) ; k ≥ 1}.
Proof. From (80), we see that the eigenvalues are geometrically simple. Then, it comes that (8) and ( 9) are satisfied with p = 2. We define our grouping by setting

λ k,1 := µ k and λ k,2 := µ k + g(µ k ).
In this setting, the eigenvalue λ k,1 is algebraically double and λ k,2 is algebraically simple. The associated eigenvectors and generalized eigenvectors of A * are

φ 0 k,1 :=   0 1 0   ϕ k , φ 1 k,1 :=   1 0 0   ϕ k , φ 0 k,2 :=   0 0 1   ϕ k ,
which obviously form a complete family in X. Moreover, for all k ≥ 1,

B * φ 0 k,1 = B * φ 0 k,2 = (x 0 , ϕ k ) H , leading to ψ 0 k,1 = 1 (x 0 , ϕ k ) H   0 1 0   ϕ k , ψ 0 k,2 = 1 (x 0 , ϕ k ) H   0 0 1   ϕ k , and
ψ 1 k,1 = 1 (x 0 , ϕ k ) H   1 0 0   ϕ k .
To compute the minimal time T 0 (X), let us estimate the different terms appearing in [START_REF] Boyer | Approximate controllability conditions for some linear 1D parabolic systems with space-dependent coefficients[END_REF]. We have

ψ[λ k,1 ] = ψ 0 k,1 and ψ[λ k,1 , λ k,1 ] = ψ 1 k,1 implying ψ[λ k,1 ] = ψ[λ k,1 , λ k,1 ] = 1 | (x 0 , ϕ k ) H | .
Using Proposition 7.10, it only remains to compute and estimate the generalized

divided difference ψ[λ k,1 , λ k,1 , λ k,2
]. This comes from ( 121) and (122) as follows

ψ[λ k,1 , λ k,2 ] = 1 (x 0 , ϕ k ) H 1 λ k,2 -λ k,1   0 -1 1   ϕ k = 1 (x 0 , ϕ k ) H 1 g(µ k )   0 -1 1   ϕ k , and
ψ[λ k,1 , λ k,1 , λ k,2 ] = ψ[λ k,1 , λ k,2 ] -ψ[λ k,1 , λ k,1 ] λ k,2 -λ k,1 = 1 (x 0 , ϕ k ) H 1 g(µ k ) 2   -g(µ k ) -1 1   ϕ k .
Thus, using (80) we obtain

ψ[λ k,1 , λ k,1 , λ k,2 ] = 1 (x 0 , ϕ k ) H 1 g(µ k ) 2 g(µ k ) 2 + 2 = 1 (x 0 , ϕ k ) H ρ -2 e 2βµ k ρ 2 e -2βµ k + 2.
Then, for k large enough, we have

max 1, ρ -2 e 2βµ k ρ 2 e -2βµ k + 2 = ρ -2 e 2βµ k ρ 2 e -2βµ k + 2,
and, using (76), this leads to (83).

Competition between different perturbations.

Let

X = H × H × H. Let α, β > 0 with α = β and f, g : σ(A) → R be such that f (µ k ) = ρe -αµ k , g(µ k ) = ρe -βµ k , with ρ satisfying (74). Let (84) A =   A I 0 0 A + f (A) 0 0 0 A + g(A)   , D(A) = D(A) × D(A) × D(A), and (85)
B : u ∈ R → u   0 x 0 x 0   .
Again B is a bounded control operator and we still set for this example X * = X = X - and Y 0 = X.

proposition 5.3. Let us consider the control system (1) with A and B given by (84)-(85). Then,

T 0 (X) = β + min{α, β}.
Proof. The spectrum of (A * , D(A)) is given by

Λ = {µ k , µ k + f (µ k ), µ k + g(µ k ) ; k ≥ 1}.
By construction, these eigenvalues are geometrically simple. Then, it comes that [START_REF] Avdonin | Riesz bases of exponentials and divided differences[END_REF] and ( 9) are satisfied with p = 3. We define our grouping by setting

λ k,1 := µ k , λ k,2 := µ k + f (µ k ), and λ k,3 = µ k + g(µ k ).
Notice that the eigenvalues are not necessarily increasingly sorted inside the k th group depending on the relative positions of α and β but, due to the invariance of divided differences with respect to permutations, this does not change our analysis.

These eigenvalues are algebraically and geometrically simple and the associated eigenvectors are

φ k,1 := 1 1 + f (µ k ) 2   -f (µ k ) 1 0   ϕ k , φ k,2 :=   0 1 0   ϕ k , φ k,3 :=   0 0 1   ϕ k ,
which do form a complete family in X. Moreover, for all k ≥ 1,

B * φ k,1 = 1 1 + f (µ k ) 2 (x 0 , ϕ k ) H , and B * φ k,2 = B * φ k,3 = (x 0 , ϕ k ) H , leading to ψ k,1 = 1 (x 0 , ϕ k ) H   -f (µ k ) 1 0   ϕ k , ψ k,2 = 1 (x 0 , ϕ k ) H   0 1 0   ϕ k , and
ψ k,3 = 1 (x 0 , ϕ k ) H   0 0 1   ϕ k .
To compute the minimal time T 0 (X), let us determine the different terms appearing in [START_REF] Boyer | Approximate controllability conditions for some linear 1D parabolic systems with space-dependent coefficients[END_REF]. We have

ψ[λ k,1 ] = ψ k,1 , ψ[λ k,1 , λ k,2 ] = 1 (x 0 , ϕ k ) H 1 λ k,2 -λ k,1   f (µ k ) 0 0   ϕ k = 1 (x 0 , ϕ k ) H   1 0 0   ϕ k , ψ[λ k,2 , λ k,3 ] = 1 (x 0 , ϕ k ) H 1 λ k,3 -λ k,2   0 -1 1   ϕ k = 1 (x 0 , ϕ k ) H 1 g(µ k ) -f (µ k )   0 -1 1   ϕ k , 44
This manuscript is for review purposes only.

and finally

ψ[λ k,1 , λ k,2 , λ k,3 ] = ψ[λ k,1 , λ k,2 ] -ψ[λ k,2 , λ k,3 ] λ k,1 -λ k,3 = 1 (x 0 , ϕ k ) H 1 g(µ k ) (g(µ k ) -f (µ k ))   g(µ k ) -f (µ k ) -1 1   ϕ k .
Since lim k→+∞ g(µ k ) = 0, we immediately see that, for k large enough, we have

max ψ[λ k,1 ] , ψ[λ k,1 , λ k,2 ] , ψ[λ k,1 , λ k,2 , λ k,3 ] = ψ[λ k,1 , λ k,2 , λ k,3 ] ,
so that using (76) and ( 18) we get

(86) T 0 (X) = lim sup k→∞ ln ψ[λ k,1 , λ k,2 , λ k,3 ] µ k = lim sup k→∞ -ln |g(µ k )(g(µ k ) -f (µ k ))| µ k .
The analysis is now split into two cases:

1. Assume first that (87) β < α.

We deduce from (86) that

T 0 (X) = lim sup k→∞ -ln e -2βµ k 1 -e -(α-β)µ k µ k = 2β. 2. Assume now that (88) β > α.
We deduce from (86) that

T 0 (X) = lim sup k→∞ -ln e -(α+β)µ k 1 -e -(β-α)µ k µ k = α + β.
Remark 5.5. As previously the family of eigenvectors does not form a Riesz basis since for instance we have φ k,1 -φ k,2 -→ k→∞ 0. Thus, one cannot apply the results of [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF] that would give that the minimal null-control time is (see Appendix 7.5)

c(Λ) = lim sup k→∞ -ln (f (µ k )g(µ k )) µ k = α + β.
Yet, in the case (88) we still have T 0 (X) = c(Λ). However, in the case (87) we have

0 < T 0 (X) = 2β < c(Λ).
Notice that, in this case, setting f = 0 one recovers the system studied in subsection 5.1.2 for which the minimal time is exactly 2β.

Condensation in partial differential equations.

We provide in this section actual PDE examples covered by our analysis. First of all, let us emphasize that our setting naturally covers a wide range of coupled one dimensional parabolic equations. Indeed if there exists p ∈ N * such that the spectrum of A is given by the union of p families

Λ j = λ j k ; k ≥ 1
such that each family satisfies ( 9) and ( 8), then the structural assumptions on Λ are automatically satisfied (see Lemma 2.1).

A system with two different potentials.

Let us consider the following boundary control system (89)

             ∂ t y + -∂ xx + c 1 (x) 1 0 -∂ xx + c 2 (x) y = 0 0 , (t, x) ∈ (0, T ) × (0, 1), y(t, 0) = 0 u(t)
, y(t, 1) = 0 0 , t ∈ (0, T ),

y(0, x) = y 0 (x),
where c 1 , c 2 ∈ L 2 (0, 1; R). Without loss of generality we assume that c 1 and c 2 are nonnegative. The operator A appearing in this system is defined in

X = (L 2 (0, 1; R)) 2 with domain X * 1 = D(A) = H 2 (0, 1; R) ∩ H 1 0 (0, 1; R) 2 .
The control operator B is defined in a weak sense as in [START_REF] Tucsnak | Observation and control for operator semigroups[END_REF]. The expression of its adjoint is easier to rule out and is given by

B * : f g ∈ X * 1 → 0 B * g = 0 -g (0)
.

Here we denoted by B * the (scalar) normal derivative operator at x = 0 defined on H 2 (0, 1; R). Standard parabolic regularity properties show that, if we define

X * = (H 1 0 (0, 1; R)) 2 , then the operator B is admissible with respect to X -= (H -1 (0, 1; R)) 2
, in the sense of (3). Therefore, for any u ∈ L 2 (0, T ; R), (89) is well-

posed in C 0 ([0, T ], (H -1 (0, 1; R)) 2 ).
For any non-negative potential c ∈ L 2 (0, 1; R), we denote by A c the definite posi-

tive self-adjoint operator in L 2 (0, 1; R) with domain H 2 (0, 1; R) ∩ H 1 0 (0, 1; R) defined by A c y = -∂ xx y + c(x)y. Its spectrum is denoted by Λ c ⊂ (π 2 , +∞) and satisfies (90) inf λ,µ∈Λ c λ =µ | √ λ - √ µ| > 0.
We choose associated eigenfunctions denoted by ϕ c λ that are normalized in L 2 (0, 1; R)

and that satisfy (see for instance [START_REF] Kirsch | An introduction to the mathematical theory of inverse problems[END_REF]Theorem 4.11])

(91) ϕ c λ (x) = √ 2 sin( √ λx) + O 1 √ λ , uniformly in x, (92)
∂ x ϕ c λ (x) = √ 2 √ λ cos( √ λx) + O (1) , uniformly in x.
In particular, there exist C, C > 0 such that,

(93) C√ λ ≤ |B * ϕ ci λ | = |∂ x ϕ ci λ (0)| ≤ C√ λ, ∀λ ∈ Λ ci , ∀i = 1, 2.
The analysis will be based on the careful inspection of spectral properties of the adjoint operator

A * = A c1 0 1 A c2 .
It is easily seen that the spectrum of A * is given by Λ = Λ c1 ∪ Λ c2 . We will often use the following straightforward property

(94) (A ci ϕ cj λ , ϕ cj λ ) ≤ Cλ, ∀λ ∈ Λ cj , ∀i, j ∈ {1, 2},
where C depends only on c 1 and c 2 .

Our controllability result concerning system (89) is the following. theorem 5.1. For any non-negative potentials c 1 , c 2 , there exists a closed subspace Y 0 of (H -1 (0, 1; R)) 2 of finite codimension such that:

• For any y 0 ∈ Y 0 , system (89) is not approximately controllable.

• For any y 0 ∈ Y 0 , system (89) is null-controllable at any time T > 0.

Remark 5.6. The set Y 0 can be equal to the whole space (H -1 (0, 1; R)) 2 , for instance if c 1 and c 2 are close enough.

Before proving this theorem, we would like to emphasize the fact that for a system like (89), the condensation index of its spectrum can be arbitrary. Therefore, Theorem 5.1 gives another example of a system which is null-controllable at any time T > 0 (for well-prepared initial data) despite the fact that the condensation index of the spectrum is non zero. Proof. This follows from inverse spectral theory. Indeed, it is proven in [44,

Chapter 3] for instance, that for any α ∈ R and any sequence (ν k ) k≥1 ∈ l 2 , one can find a potential c ∈ L 2 (0, 1; R) such that the spectrum of A c is given by (k

2 π 2 + α + ν k ) k .
It is then clear that we can choose c 1 and c 2 such that the spectrums of A c1 and A c2 are asymptotically as close as we want and then generate an arbitrary condensation index for the spectrum of A * . Note that such potentials are not necessarily nonnegative, but this is actually not really needed in our analysis (we simply need that the spectrum of A c is made of positive eigenvalues).

In the context of parabolic control problems, this was already noticed and used in [START_REF] Ouaili | Minimal time of null controllability of two parabolic equations[END_REF].

We can now move to the proof of the main result of this section.

Proof (of Theorem 5.1). The first part of the proof consists in a precise description of the spectral properties of A * .

• For any λ ∈ Λ c2 , we have a first eigenfunction of A * given by (95)

φ 0 λ := 0 ϕ c2 λ .
Moreover, by (93), we have

B * φ 0 λ = B * ϕ c2 λ = 0,
so that all those eigenfunctions are observable.

-If λ ∈ Λ c1 , this eigenfunction is algebraically and geometrically simple.

-However, if λ ∈ Λ c2 ∩ Λ c1 , this eigenvalue is (algebrically or geometrically) double. As detailed in Section 6.3, we can deal with geometric multiplicity of eigenvalues with an adequate choice of the space of initial conditions Y 0 . This choice will be precised in (101).

Let us define (96)

β λ := (ϕ c1 λ , ϕ c2 λ ).
By (91), we see that there exists λ 0 such that

(97) 1 2 ≤ β λ ≤ 1, ∀λ ∈ Λ c2 ∩ Λ c1 , s.t. λ > λ 0 . * If β λ = 0 then there exists a solution of (A c2 -λ)ϑ λ = -ϕ c1 λ ,
that we can choose to satisfy B * ϑ λ = 0 in such a way that

φ0 λ := ϕ c1 λ ϑ λ ,
is another independent eigenfunction of A * associated with λ that satisfy B * φ0 λ = 0. Note that, by (97), we know that β λ can vanish only for a finite number of values of λ. * Assume now that β λ = 0. In that case, λ is geometrically simple but there exists a generalized eigenfunction φ 1 λ associated with φ 0 λ of the following form

φ 1 λ := 1 β λ ϕ c1 λ χ λ ,
where χ λ is the unique solution of

(A c2 -λ)χ λ = β λ ϕ c2 λ -ϕ c1 λ ,
that satisfy B * χ λ = 0.

We can express χ λ in the basis ϕ c2

• as follows

χ λ = a λ ϕ c2 λ - µ∈Λ c 2 µ =λ (ϕ c1 λ , ϕ c2 µ ) λ -µ ϕ c2 µ , with a λ = 1 B * ϕ c2 λ µ∈Λ c 2 µ =λ (ϕ c1 λ , ϕ c2 µ ) λ -µ B * ϕ c2 µ .
• Consider now λ ∈ Λ c1 \Λ c2 . We obtain another family of eigenfunctions given by (98)

φ 0 λ := ϕ c1 λ ξ λ ,
where ξ λ satisfies (A c2 -λ)ξ λ = -ϕ c1 λ . This last equation has a unique solution since λ ∈ Λ c2 and it can be expressed as follows (99)

ξ λ = µ∈Λ c 2 (ϕ c1 λ , ϕ c2 µ ) λ -µ ϕ c2 µ .
We now state the following lemma, whose proof is postponed at the end of this subsection.

Lemma 5.1. There exists C 1 , C 2 > 0 depending only on c 1 , c 2 such that

(100) |B * ξ λ | 2 ≥ C 1 λ -C 2 , ∀λ ∈ Λ c1 \ Λ c2 .
This lemma shows in particular that B * ξ λ can only vanish for a finite number of values of λ.

It is straightforward to prove that the family of (generalized) eigenfunctions we just computed is complete in X. We can now define Y 0 to be the set of initial data

y 0 ∈ X -such that (101)    y 0 , φ0 λ -, = 0, ∀λ ∈ Λ c1 ∩ Λ c2 , s.t. β λ = 0, see (96), y 0 , φ 0 λ -, = 0, ∀λ ∈ Λ c1 \ Λ c2 , s.t. B * ξ λ = 0, see (99)
.

By construction, this set is closed and of finite codimension, moreover it is clear that initial data not belonging to this set are not approximately controllable. Note that this definition actually excludes the influence of the possible presence of a geometrically double eigenvalue in the system.

We will now endow the space X * with the following norm f g

2 * := A c1 f, f H -1 ,H 1 0 + A c2 g, g H -1 ,H 1 0 ,
which is equivalent to the usual H 1 -norm and more comfortable for the following computations. Note that, if f, g ∈ H 2 (0, 1; R), this quantity is simply equal to

(A c1 f, f ) + (A c2 g, g).
From (90), we can find a ρ > 0 such that

| √ λ - √ µ| > ρ, ∀λ = µ ∈ Λ ci , ∀i = 1, 2.
This implies, in particular, that (102) |λ -µ| > ρ(

√ λ + √ µ) > ρ, ∀λ = µ ∈ Λ ci , ∀i = 1, 2.
Without loss of generality, we can assume that ρ < C 2C1 where C and C 1 are respectively defined in (93) and (100).

It follows that Λ satisfies the summability condition [START_REF] Avdonin | Riesz bases of exponentials and divided differences[END_REF], as well as the weak gap condition [START_REF] Baiocchi | Ingham-Beurling type theorems with weakened gap conditions[END_REF] with p = 2. We can thus consider a grouping (G k ) k ∈ G(Λ, 2, r, ρ) for a suitable r > 0. We will now use the formula [START_REF] Boyer | Analysis of non scalar control problems for parabolic systems by the block moment method[END_REF] we obtained for T 0 (Y 0 ) to prove that the system is null-controllable from Y 0 at any time T > 0. For that we will consider one of the groups G (we drop the index k which is not important here) and give estimates of the corresponding divided differences.

• Case 1: G = {λ} is of cardinal 1.

-If λ ∈ Λ c2 we need to estimate the quantity

ψ[λ] 2 * := P * Y0 φ 0 λ B * φ 0 λ 2 * ≤ φ 0 λ B * φ 0 λ 2 * except if λ ∈ Λ c1 ∩ Λ c2
and β λ = 0. The computations above, and (93), show that

ψ[λ] 2 * ≤ 1 |B * ϕ c2 λ | 2 (A c2 ϕ c2 λ , ϕ c2 λ ) = λ |B * ϕ c2 λ | 2 ≤ 1 C2 .
-If λ ∈ Λ c1 \ Λ c2 , recall that φ 0 λ is given by ( 98) and that we need to estimate the same quantity ψ[λ] * , in the case where B * ξ λ = 0. Since λ is the only element in the group G, we know that |λ -µ| ≥ r for any other eigenvalue µ. With this remark, we can deduce that

ψ[λ] 2 * ≤ 1 |B * ξ λ | 2 (A c1 ϕ c1 λ , ϕ c1 λ ) + (A c2 ξ λ , ξ λ ) = 1 |B * ξ λ | 2   λ + µ∈Λ c 2 µ(ϕ c1 λ , ϕ c2 µ ) 2 (λ -µ) 2   ≤ 1 |B * ξ λ | 2   λ + 1 r 2 µ∈Λ c 2 (ϕ c1 λ , ϕ c2 µ )(ϕ c1 λ , A c2 ϕ c2 µ )   .
Using Parseval's identity, (94) and then (100), we finally obtain

ψ[λ] 2 * ≤ 1 |B * ξ λ | 2 λ + 1 r 2 (ϕ c1 λ , A c2 ϕ c1 λ ) ≤ λ |B * ξ λ | 2 1 + C r 2 ≤ C.
-Finally, if λ ∈ Λ c1 ∩ Λ c2 , then we need to estimate the contribution of the generalized eigenvector ψ[λ, λ]

2 * := P * Y0 φ 1 λ /(B * φ 0 λ ) 2 * . A com-
putation similar to the one above, for λ > λ 0 , leads to

ψ[λ, λ] 2 * ≤ 1 |B * φ 0 λ | 2 β 2 λ ((A c1 ϕ c1 λ , ϕ c1 λ ) + (A c2 χ λ , χ λ )) = 1 |B * ϕ c2 λ | 2 β 2 λ     λ(1 + a 2 λ ) + µ∈Λ c 2 µ =λ (ϕ c1 λ , ϕ c2 µ ) 2 (λ -µ) 2 µ     ≤ 1 |B * ϕ c2 λ | 2 β 2 λ λ(1 + a 2 λ ) + 1 r 2 (A c2 ϕ c1 λ , ϕ c1 λ ) ≤ λ |B * ϕ c2 λ | 2 β 2 λ 1 + a 2 λ + C r 2 ≤ C(1 + a 2 λ ).
Here, we have used (97) to bound from below the term β λ . It remains to bound a λ . We proceed as follows, by using (92), (102), and (94)

|a λ | ≤ 1 |B * ϕ c2 λ | µ∈Λ c 2 µ =λ (ϕ c1 λ , ϕ c2 µ ) λ -µ |B * ϕ c2 µ | ≤ C √ λ µ∈Λ c 2 µ =λ |(ϕ c1 λ , ϕ c2 µ )| √ µ √ λ + √ µ ≤ 1 √ λ     µ∈Λ c 2 µ =λ |(ϕ c1 λ , ϕ c2 µ )| 2 µ     1 2   µ∈Λ c 2 1 µ   1 2 ≤ C √ λ (A c2 ϕ c1 λ , ϕ c1 λ ) 1 2 ≤ C.
This concludes the proof of the uniform bound of ψ[λ, λ] * .

• Case 2: G = {λ 1 , λ 2 } is of cardinal 2. Since the diameter of G is smaller than ρ, we can choose the numbering such that

λ 1 ∈ Λ c1 \ Λ c2 and λ 2 ∈ Λ c2 \ Λ c1 .
In particular, we have B * φ 0 λ1 = 0, B * φ 0 λ2 = 0 and there is no generalized eigenvector associated to this group G. Therefore, the only new quantity we need to estimate is the contribution of the following divided difference

ψ[λ 1 , λ 2 ] 2 * ≤ 1 |λ 1 -λ 2 | 2 φ 0 λ1 B * φ 0 λ1 - φ 0 λ2 B * φ 0 λ2 2 * .
Using formulas (95) and (98), we find

ψ[λ 1 , λ 2 ] 2 * ≤ 1 |λ 1 -λ 2 | 2 1 B * ξ λ1 ϕ c1 λ1 ξ λ1 - 1 B * ϕ c2 λ2 0 ϕ c2 λ2 2 * .
Since λ 1 and λ 2 can be arbitrarily close it is not clear that this estimate does not blow up. In particular, if we use the triangle inequality, we will not be able to take benefit of compensations that occur in the divided difference.

We will thus make appear from (99) the principal part of ξ λ1 as follows

ξ λ1 = β λ1,λ2 λ 1 -λ 2 ϕ c2 λ2 + ζ λ1 ,
with β λ1,λ2 := (ϕ c1 λ1 , ϕ c2 λ2 ) and

ζ λ1 := λ 1 -λ 2 β λ1,λ2 µ∈Λ c 2 µ =λ2 (ϕ c1 λ1 , ϕ c2 µ ) λ 1 -µ ϕ c2 µ .
Thus,

ξ λ1 B * ξ λ1 - ϕ c2 λ2 B * ϕ c2 λ2 = ϕ c2 λ2 + ζ λ1 B * (ϕ c2 λ2 + ζ λ1 ) - ϕ c2 λ2 B * ϕ c2 λ2 = 1 B * (ϕ c2 λ2 + ζ λ1 ) - 1 B * ϕ c2 λ2 ϕ c2 λ2 + ζ λ1 B * (ϕ c2 λ2 + ζ λ1 ) (103)
Since we are only interested in the asymptotic behavior when λ 1 and λ 2 are large, we see that we can assume from (91) that the following properties hold (104)

|β λ1,λ2 | ≥ 1/2, λ 2 ≥ λ 1 /2 ≥ 1.
Using that |λ 1 -µ| ≥ r for all µ ∈ Λ c2 , with µ = λ 2 , we can find with ( 104) and (94) the following bound

(105) (A c2 ζ λ1 , ζ λ1 ) ≤ |λ 1 -λ 2 | 2 (A c2 ϕ c1 λ1 , ϕ c1 λ1 ) r 2 |β λ1,λ2 | 2 ≤ C * |λ 1 -λ 2 | 2 λ 1 .
Moreover, we have

|B * ζ λ1 | ≤ |λ 1 -λ 2 | |β λ1,λ2 | µ∈Λ c 2 µ =λ2 |(ϕ c1 λ1 , ϕ c2 µ )| |λ 1 -µ| |B * ϕ c2 µ |, ≤ C C|λ 1 -λ 2 | µ∈Λ c 2 µ =λ2 |(ϕ c1 λ1 , ϕ c2 µ )| |λ 1 -µ| √ µ ≤ C C|λ 1 -λ 2 |     µ∈Λ c 2 µ =λ2 (ϕ c1 λ1 , ϕ c2 µ ) 2 µ     1/2     µ∈Λ c 2 µ =λ2 1 |λ 1 -µ| 2     1/2 .
We use Parseval's identity and (94) to bound the second factor by C √ λ 1 . Moreover, by using (102), we have for any µ ∈ Λ c2 , µ = λ 2 ,

|λ 1 -µ| ≥ |λ 2 -µ| -|λ 1 -λ 2 | ≥ ρ( √ µ + λ 2 ) -ρ ≥ ρ √ µ,
so that the value of the series in the last factor is uniformly bounded. Hence, we have proved

(106) |B * ζ λ1 | ≤ C 1 |λ 1 -λ 2 | λ 1 .
From this last estimate, (100) and (93), we deduce that

|B * (ϕ c2 λ2 + ζ λ1 )| ≥ |B * ϕ c2 λ2 | -C 1 |λ 1 -λ 2 | λ 1 ≥ C λ 1 -C 1 |λ 1 -λ 2 | λ 1 .
Recall that ρ < C/(2C 1 ). Since λ 1 and λ 2 belong to the same group G, we have |λ 1 -λ 2 | ≤ ρ and thus, we obtain the estimate

(107) |B * (ϕ c2 λ2 + ζ λ1 )| ≥ C 2 λ 1 .
Coming back to the definition of ψ[λ 1 , λ 2 ] and using (103) and the triangle 52

This manuscript is for review purposes only.

inequality, we write

ψ[λ 1 , λ 2 ] 2 * = 1 |λ 1 -λ 2 | 2 1 B * (ϕ c2 λ2 + ζ λ1 ) λ1-λ2 β λ 1 ,λ 2 ϕ c1 λ1 ϕ c2 λ2 + ζ λ1 - 1 B * ϕ c2 λ2 0 ϕ c2 λ2 2 * ≤ ϕ c1 λ1 0 2 * β 2 λ1,λ2 |B * (ϕ c2 λ2 + ζ λ1 )| 2 + 2 0 ζ λ1 2 * |B * (ϕ c2 λ2 + ζ λ1 )| 2 |λ 1 -λ 2 | 2 + 2 0 ϕ c2 λ2 2 * |λ 1 -λ 2 | 2 1 B * (ϕ c2 λ2 + ζ λ1 ) - 1 B * ϕ c2 λ2 2 =:S 1 + S 2 + S 3 .
We now analyze each of the three terms.

-Using (94), ( 104) and (107), we can obtain

S 1 ≤ 16 C2 .
-Using ( 105) and (107), we get

S 2 ≤ 8C * C2 .
-Finally, with (94), we write

S 3 = λ 2 |λ 1 -λ 2 | 2 (B * ζ λ1 ) 2 (B * (ϕ c2 λ2 + ζ λ1 )) 2 (B * ϕ c2 λ2 ) 2 ,
so that, with (93), ( 107) and (106), we get

S 3 ≤ C 2 1 C4 .
All in all, we have obtained a uniform bound for ψ[λ 1 , λ 2 ] * , which is exactly the compensation phenomenon we were expecting for this particular system.

As a conclusion, we finally proved that, whatever the group G is, all the divided differences ψ[λ], ψ[λ, λ] or ψ[λ 1 , λ 2 ] remain bounded uniformly. It follows from [START_REF] Boyer | Analysis of non scalar control problems for parabolic systems by the block moment method[END_REF] that T 0 (Y 0 ) ≤ 0, so that our main Theorem 1.1 show that (89) is null-controllable at any time T > 0 for any initial data y 0 ∈ Y 0 .

It remains to prove the lemma.

Proof (of Lemma 5.1). By definition, the function ξ λ satisfies (108)

-∂ xx ξ λ + c 2 (x)ξ λ = λξ λ -ϕ c1 λ , in (0, 1).
Using [2, Lemmas 2.2 and 2.3], and the fact that ϕ c1 λ is normalized in L 2 (0, 1; R), we have

|ξ λ (x)| 2 + 1 λ |∂ x ξ λ (x)| 2 ≤ C |ξ λ (y)| 2 + 1 λ |∂ x ξ λ (y)| 2 + C λ , ∀x, y ∈ [0, 1].
We take y = 0 in this inequality and we integrate with respect to x to obtain

ξ λ 2 ≤ C λ |∂ x ξ λ (0)| 2 + C λ = C λ |B * ξ λ | 2 + C λ .
It remains to bound from below the L 2 norm of ξ λ . To this end, we multiply (108) by ϕ c1 λ and integrate over (0, 1). After integration by parts, and using the equation satisfied by ϕ c1 λ , we get

-1 = 1 0 (c 2 -c 1 )ξ λ ϕ c1 λ .
The Cauchy-Schwarz inequality gives

1 ≤ c 1 -c 2 ξ λ ϕ c1 λ L ∞ ,
and since by (91), we have a uniform L ∞ bound on ϕ c1 λ , the proof is complete.

5.2.2.

A system with different diffusions and a non constant coupling term. Let us briefly describe another example of a coupled parabolic system with boundary control which has motivated our study. This example is analyzed in details in [START_REF] Samb | Boundary null-controllability of two coupled parabolic equations: simultaneous condensation of eigenvalues and eigenfunctions[END_REF]. We consider the following system (109)

             ∂ t y + -∂ xx q(x) 0 -ν∂ xx y = 0 0 , (t, x) ∈ (0, T ) × (0, 1), y(t, 0) = 0 u(t) , y(t, 1) = 0 0 , t ∈ (0, T ), y(0, x) = y 0 (x),
where q ∈ L ∞ (0, 1) and ν > 0.

The spectrum of

A * = -∂ xx 0 q(x) -ν∂ xx is Λ = {k 2 , ν k 2 , k ≥ 1}.
• System (109) in the case where q(x) = 1 and ν = 1 was studied in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF] where the influence of the condensation of eigenvalues in the system was first pointed out. It was proved that the minimal null-control time was exactly the condensation index of Λ, provided that √ ν ∈ Q.

• System (109) with a non constant q but with the same diffusions, that is ν = 1, was studied in [START_REF] Khodja | New phenomena for the null controllability of parabolic systems: minimal time and geometrical dependence[END_REF]. The picture is different since in that case, there is no condensation of eigenvalues but there may however exist a minimal nullcontrol time (depending on the coupling term q) due to very weak observation properties of the eigenfunctions.

• In the general case, assuming that √ ν ∈ Q, the eigenvalues are algebraically and geometrically simple and it is proved in [START_REF] Samb | Boundary null-controllability of two coupled parabolic equations: simultaneous condensation of eigenvalues and eigenfunctions[END_REF] that the associated family of eigenfunctions is complete in X * = (H -1 (0, 1)) 2 , and that, moreover, there exist functions q and values of ν, √ ν ∈ Q, such that this family (properly normalized) is not a Riesz basis of X * . Therefore the abstract results in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF][START_REF] Khodja | New phenomena for the null controllability of parabolic systems: minimal time and geometrical dependence[END_REF] do not apply.

Inspirated by the block moment method presented in the present paper, a suitable value of T q,ν 0 is defined in [START_REF] Samb | Boundary null-controllability of two coupled parabolic equations: simultaneous condensation of eigenvalues and eigenfunctions[END_REF] (taking into account both effects of condensation of eigenvalues and weak observation of eigenfunctions) such that T q,ν 0 is the minimal null-control time of (109).

6. Extensions.

6.1. Dealing with complex eigenvalues.

In the previous sections, we decided to state our results in the framework of real eigenvalues to simplify the presentation. However, most of them still hold for complex eigenvalues satisfying assumptions largely inspired from [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF]. More precisely, for a function N : R + → R, we will consider the class L w (δ, p, ρ, N ) of the families Λ ⊂ C satisfying

• Parabolicity condition: λ ≥ δ|λ|, ∀λ ∈ Λ.

• Asymptotic behavior: for any ε > 0, we have

λ∈Λ |λ|>N (ε) 1 |λ| ≤ ε.
• Weak gap condition with parameters ρ > 0 and p ∈ N * :

#Λ ∩ ([µ, µ + ρ] + iR) ≤ p, ∀µ > 0. In that case, a grouping (G k ) k should satisfy Λ = k≥1 G k , #G k ≤ p, diam(G k ) < ρ, inf( G k+1 ) -sup( G k ) > r.
The corresponding formula the minimal time T 0 (Y 0 ) will be now given by

T 0 (Y 0 ) := lim sup k→∞ ln   max µ∈N g k µ≤α k ψ[λ (µ1)
k,1 , . . . , λ

(µg k ) k,g k ] *   λ k,1 ,
and our results (namely Theorems 1.1, 2.1 and 4.1) still hold in that case.

Most of the proofs are very similar by taking care of the following points:

• The divided differences associated with pairwise distinct points x 0 , . . . , x n in the complex plane do not satisfy the Lagrange theorem but instead the following slightly weaker result, due to Jensen [START_REF] Jensen | Sur une expression simple du reste dans la formule d'interpolation de Newton[END_REF].

. , x n }) such that |f [x 0 , . . . , x n ]| ≤ f (n) (z) n! .
-For any z ∈ U , we have

f [x 0 , . . . , x n ] - f (n) (z) n! ≤ C f,n diam(U ).
• The Blaschke product W k should be replaced by

W k (z) = p j=1 λ∈Λj λ -z λ + z .
• Finally, in the restriction argument of Section 2.1.2, the holomorphy domain C + 2ε should be replaced by a sector {z ∈ C, z > 2ε, | z| ≤ δ 2 |λ|}.

6.2.

Weakening the assumptions on the control operator.

In this article, we not only study the classical null-controllability property (i.e.

Y 0 = X -), we also provide a more accurate description depending on the space of initial conditions Y 0 one wants to drive to 0. In this setting, the assumption (10) can be too strong.

It is easily seen that a necessary approximate null-controllability condition in that case is the following: for any λ ∈ Λ and any l ∈ 0, α λ -1 we have

(110) B * φ j λ = 0, ∀j ∈ 0, l ⇒ P * Y0 φ j λ = 0, ∀j ∈ 0, l ,
where, in this formula, (φ j λ ) j is a Jordan chain associated with the eigenvalue λ. Note that such a Jordan chain is not unique but (110) does not depend on the particular chain we choose. Note also that the assumption ( 12) can be verified using any Jordan chain.

From now on, we assume that (110) holds. For any λ ∈ Λ, two cases have to be considered:

• Case 1: We have B * φ j λ = 0, for all j ∈ 0, α λ -1 .

From (110), it follows that for any y 0 ∈ Y 0 , any T > 0, all the moment equation ( 21) corresponding to this eigenvalue are automatically satisfied. It follows that we can simply remove this eigenvalue from the family Λ when studying the control problem at time T from Y 0 .

• Case 2:

(111) there exists j * ∈ 0, α λ -1 s.t. B * φ j λ = 0, ∀j < j * , and B * φ j * λ = 0.

In that case, for j > j * we set

β j := - B * φ j λ B * φ j * λ ,
and then by induction, we define

(112) φj λ =        φ j λ , for j ≤ j * , φ j λ + j-1 k=j * β j+j * -k φk λ , for j > j * .
This construction ensures that ( φj λ ) j and (φ j λ ) j span the same space, that (113) B * φj λ = 0, if and only if j = j * ,

56

This manuscript is for review purposes only.

and finally satisfy the equations

A * φj λ = λ φj λ + φj-1 λ + γ j φ j * -1 λ ,
for some γ j ∈ R whose precise value is unimportant in the sequel.

A straightforward computation shows that the semi-group generated by -A * satisfy e -tA * φj λ ∈ (e t φ)[λ (j+1) ] + V j * , where V j * := Span(φ 0 λ , . . . , φ j * -1 λ

). We shall prove that the term in V j * does not contribute to the moment problem. Indeed, from (110) and (111), we have V j * ⊂ Ker B * ∩ Ker P * Y0 . Thus:

-Concerning the control term, we have ) ], and by (113), it simply remains

B * e -tA * φj λ = B * (e t φ)[λ (j+1
B * e -tA * φj λ = 0, if j < j * , B * φj * λ e t [λ (j-j * ) ], if j ≥ j * .
-Concerning the contribution of the source term, we have

P * Y0 e -T A * φj λ = P * Y0 (e T φ)[λ (j+1) ] = P * Y0 j+1 k=1 e T [λ (j+2-k) ] φ[λ (k) ] = P * Y0 j+1 k=j * +1 e T [λ (j+2-k) ] φ[λ (k) ],
with the convention that the sum is 0 as soon as j < j * .

We may now adapt the definition of our null-control time by setting (114)

ψ l λ = P * Y0 φj * +l λ B * φj * λ , ∀l ∈ 0, α λ -1 -j * ,
so that the moment problem associated with this eigenvalue becomes

T 0 u(T -t) (-t) l l! e -λt dt = -y 0 , (e T ψ)[λ (l+1) ] -, , ∀l ∈ 0, α λ -1 -j * .
This is formally exactly the same as [START_REF] Duyckaerts | Resolvent conditions for the control of parabolic equations[END_REF] except that the multiplicity of the eigenvalue have been changed into α λ -j * and the associated values of ψ • λ have been constructed as explained above by ( 112) and (114).

As a conclusion, to obtain the definition of the minimal null-control time from Y 0 assuming that (110) holds, we simply need to ignore the eigenvalues corresponding to case 1, and to modify the multiplicity and the Jordan chain as explained above for the eigenvalues that are in case 2. Then, we define formally T 0 (Y 0 ) by the same formula as [START_REF] Bhandari | Boundary null-controllability of coupled parabolic systems with Robin conditions[END_REF] and we obtain exactly the same result as Theorem 1.1.

Moreover, it clearly appears from the proof that (110) is actually a necessary and sufficient condition to solve the moment problem associated to any finite number of eigenvalues. Thus (110) is a necessary and sufficient condition for the approximate null-controllability from Y 0 .

57

This manuscript is for review purposes only.

Dealing with geometrical multiplicities.

In the considered setting of scalar controls, if one wants to be able to control every initial data, that is to take Y 0 = X -, then it is absolutely necessary to assume, as we did in Section 1.3.1, that all the eigenvalues are geometrically simple. Indeed, as soon as dim ker(A * -λ) ≥ 2, there necessarily exists an eigenfunction that belongs to ker B * and the system is thus not even approximately controllable.

However, if we are willing to restrict the space of initial data we consider, then our result can apply when some eigenvalues are not geometrically simple. More precisely, we will show that our results can be adapted under the condition (115) ker(A * -λ) ∩ ker B * ⊂ ker P * Y0 , ∀λ ∈ Λ, replacing the geometrical simplicity condition and the observation condition [START_REF] Beauchard | Heat equation on the Heisenberg group: observability and applications[END_REF].

We will also assume, for simplicity, that geometrically multiple eigenvalues are algebraically simple. Let us give the main arguments.

• First, it is clear that (115) is a necessary condition for the null-controllability from initial data in Y 0 . Indeed, if this condition does not hold, there exists a φ ∈ ker(A * -λ) ∩ ker B * and a y 0 ∈ Y 0 such that y 0 , φ -, = 0. For this particular y 0 , it cannot exist a u such that y(T ) = 0 because it would contradict the equality [START_REF] Dolecki | Observability for the one-dimensional heat equation[END_REF].

• Assume now that (115) holds. For each λ ∈ Λ, there are two cases:

-If ker(A * -λ) ⊂ ker B * , then by (115), we also have ker(A * -λ) ⊂ ker P * Y0 . In that case, the controllability condition (20) automatically holds for any φ ∈ ker(A * -λ) and any control u since both terms in the equality are equal to 0.

Hence, everything happens as if λ were not an eigenvalue of A * , and we can essentially not consider it in the moment problem to be solved.

-If ker(A * -λ) ⊂ ker B * , we fix any φ 0 λ in ker(A * -λ) such that B * φ 0 λ = 0. As B * is a linear form, we observe that ker(A * -λ) = Span φ 0 λ + (ker(A * -λ) ∩ ker B * ) .

It follows from (115) that (20) holds for any φ ∈ ker(A * -λ) if and only if [START_REF] Dolecki | Observability for the one-dimensional heat equation[END_REF] holds for φ = φ 0 λ .

Therefore, everything happens as if λ were geometrically simple with φ 0 λ as a unique (up to a scalar) eigenvector. Our analysis can then be pushed forward without change.

Appendices.

We gather in this final section some definitions or intermediate results that we used in this paper.

Wellposedness.

This section is dedicated to the proof of Proposition 1.1.

First of all, let us notice that the problem (2) admits at most one solution y ∈ C 0 ([0, T ]; X -1 ) and that the continuous dependancy directly follows from (2). Thus, it remains to prove the existence of a function y ∈ C 0 ([0, T ]; X -1 ) satisfying [START_REF] Allonsius | Spectral analysis of discrete elliptic operators and applications in control theory[END_REF].

From [49, Proposition 2.10.3] it comes that A can be uniquely extended to an operator à ∈ L(X, X -1 ). Moreover it comes from [49, Proposition 2.10.4] that -à generates a C 0 -semigroup in X -1 satisfying e -t à = Ãe -tA Ã-1 , ∀t ≥ 0.

Thus, for any T > 0, any y 0 ∈ X -1 and any u ∈ L 2 (0, T ; U ) the problem y (t) + Ãy(t) = Bu(t), y(0) = y 0 admits a unique mild solution y ∈ C 0 ([0, T ], X -1 ) given by ( 116) y(t) = e -t Ãy 0 + t 0 e -(t-s) ÃBu(s)ds.

We prove now that this function satisfies (2). To do so, we simply prove that the semigroup e -t à is the adjoint of e -tA * in the duality between X * 1 and X -1 .

Let x ∈ X and z ∈ X 1 such that x = Az. As à is an extension of A it also comes that x = Ãz. Then, as e -tA (X 1 ) ⊂ X 1 it comes that e -t Ãx = Ãe -tA Ã-1 Ãz = Ãe -tA z = Ae -tA z = e -tA Az = e -tA x.

Then, for any x ∈ X and any z ∈ X * Thus, the density of X in X -1 implies (117) e -t Ãy, z -1,1 * = y, e -tA * z -1,1 * , ∀y ∈ X -1 , ∀z ∈ X * 1 .

Finally, the duality pairing of (116) with any z t ∈ X * 1 with the computation rule (117) directly gives (2). By [START_REF] Baiocchi | Ingham-Beurling type theorems with weakened gap conditions[END_REF], we know that one of this sets is empty since if it not the case, there is at least p + 1 elements in Λ ∩ [µ 1 , µ 1 + ρ] because µ 1 ∈ Λ and pr = ρ. Let j ∈ 1, p such that Λ ∩ (µ 1 + (j -1)r, µ 1 + jr] = ∅. We define

G 1 := Λ ∩ [µ 1 , µ 1 + (j -1)r], 59
This manuscript is for review purposes only.

whose cardinal is, by [START_REF] Baiocchi | Ingham-Beurling type theorems with weakened gap conditions[END_REF], less or equal than p and diameter is less than ρ. Moreover, by construction, we have inf(Λ \ G 1 ) -sup G 1 > r.

This allows to build G 2 by the same construction applied on Λ \ G 1 while ensuring the required properties, and following this process we construct the sequence (G k ) k .

7.3. About divided differences.

In this section we give all the properties concerning divided differences that are used all along this article. This notion is a key technical tool in our analysis as it drastically eases the computations and the formulation of the results. The definition and results given in Section 7.3.1 are classical in the field of interpolation (see for instance [START_REF] Powell | Approximation theory and methods[END_REF]Chap. 5]). To deal with algebraic multiplicity we use a generalization of divided differences where the 'interpolation points' are not necessarily distincts.

Let us mention that there exists generalizations in this direction (see for instance [START_REF] Powell | Approximation theory and methods[END_REF]Chap. 5]) in the context of Hermite interpolation. However as we are not directly dealing with interpolation, we propose such a generalization adapted to our purposes. This is detailed in Section 7.3.2.

Definitions and basic properties.

Let V be a real vector space, n ∈ N and x 1 , . . . , x n ∈ R. Assume that x 1 , . . . , x n are pairwise distinct (see Section 7.3.2 for a generalization). Let f 1 , . . . , f n ∈ V be given. In all what follows, if f : R → V is a given function it will be implicitely assumed It is clear that (z h ) h is contained in a compact set so that, up to a subsequence, we may find a limit z of (z h ) h that belongs to Conv({x 1 , . . . , x n }) and satisfies the required property.

that f i = f [x i ] = f (x i ).
proposition 7.10. Let (µ p ) p∈ 0,|α| ⊂ N n be a family of multi-indices satisfying (123). For any multi-index µ such that µ ≤ α we have

f [x (µ) ]] ≤ N 2 N -1 (1 + R) N -1 max p∈ 1,|α| f [x (µ p ) ] .
Proof. We proceed as in the proof of Proposition 7.8 by passing to the limit in the similar result for standard divided differences (Corollary 7.1).

For generalized divided differences, there is no simple equivalent to the Newton formula (Proposition 7.3). However, we can state the following result. Proof. Since the divided differences are clearly linear with respect to the data f α , the existence of the coefficients θ µ j,l is straightforward. Let us prove the claimed estimates. From now on we assume that l is fixed. Moreover, for any j ∈ 1, n we introduce the notation d j := min i∈ 1,n =j |x i -x j | and we define δ j ∈ N n to be the Kronecker multi-index, that is δ j i = 0 for i = j and δ j j = 1.

• When µ j < l + 1, it is clear from the recurrence formulas (121) and (122) that the value of f [x (µ) ] does not dependent on the value f l j , and therefore θ µ j,l = 0. .

-Assume first that N = l + 1 and let µ such that |µ| = N . If µ j < l + 1 we have already seen that θ µ j,l = 0 which obviously implies (124). If µ j = l + 1, since |µ| = l + 1, we necessarily have µ i = 0 for any i = j, so that (121) gives f [x (µ) ] = f l j , which implies that θ µ j,l = 1, that is exactly (124) with C |µ| = 1 in that case.

65

This manuscript is for review purposes only.

-Assume now that, for some N ≥ l + 1, (124) holds and let µ such that |µ| = N + 1.

If µ i = 0 for any i = j, then we have

f [x (µ) ] = f µj -1 j
, which implies that θ µ j,l = 0 since l = µ j -1 and (124) is obvious.

If there is a i 0 = j such that µ i0 ≥ 1 then we use (122) to get

f [x (µ) ] = f [x (µ-δ j ) ] -f [x (µ-δ i 0 ) ] x i0 -x j ,
which implies the formula θ µ j,l = θ µ-δ j j,l -θ µ-δ i 0 j,l T 0 (y 0 ) = T 0 (Y 0 ).

66

This manuscript is for review purposes only.

To ease the reading let us do the computations in the simpler case η = 1; the extension to the case η ≥ 2 being straightforward. Let us introduce

x k,l := • Since y 0 is normalized, we have y 0 , x k,l -, ≤ x k,l * , for any k and l, it immediately comes that T 0 (y 0 ) ≤ T 0 (Y 0 ) and thus sup y0∈Y0 T 0 (y 0 ) ≤ T 0 (Y 0 ).

• Conversely, let T be such that sup y0∈Y0 y0 -=1

T 0 (y 0 ) < T.

Setting xk,l := e -λ k,1 T x k,l , it comes that for any y 0 ∈ Y 0 , we have sup k≥1 l∈ 1,g k y 0 , xk,l -, < +∞, and this property is in fact true for any y 0 ∈ X -since P * Y0 xk,l = xk,l , so that we have y 0 , xk,l -, = y 0 , P * Y0 xk,l -, = P Y0 y 0 , xk,l -, ,

67

This manuscript is for review purposes only. This ends the proof of Proposition 1.3.

7.5. On the condensation index.

In this appendix we give some useful properties concerning the condensation index of a sequence. Let Σ be a family of positive real numbers. We start by recalling the definition of c(Σ). This definition (and also its extension to complex sequences) is given in [START_REF] Shackell | Overconvergence of Dirichlet series with complex exponents[END_REF]. Notice that due to the assumption (125), both functions

z ∈ C → σ∈Σ 1 - z σ , z ∈ C → σ∈Σ 1 + z σ ,
are entire. Since Σ ⊂ (0, +∞), the function E Σ has simple roots corresponding exactly to Σ. Thus, the condensation index of such sequences is well defined. The fact that it belongs to [0, +∞] is proved in [START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF].

In the case where the considered sequence satisfies the weak gap condition (9) the computation of the condensation index can be simplified: the grouping introduced in Proposition 7.1 is an optimal condensation grouping in the following sense. Recall that G [σ] is the element of (G k ) k≥1 containing σ.

68

This manuscript is for review purposes only.

Proof. The proof follows directly from (64).

Using this result, we compute easily the condensation index of the particular sequence used in Section 5.1.

proposition 7.13. Let (µ k ) k≥1 be a real increasing sequence such that Thus, as 2β < α + β, we obtain c(Θ) = α + β.

1. Introduction. 1 . 1 .

 11 Problem under study and state of the art.

1. 3 . 5 .

 35 Structure of the article.

Corollary 2 . 2 .

 22 Let p ∈ N * , r, ρ > 0 and N : (0, +∞) → R. Assume that Λ ∈ L w (p, ρ, N ) and let (G k ) k ∈ G(Λ, p, r, ρ) be an associated grouping. Assume that η = 1 and T 0 (Y 0 ) < +∞. For any T > T 0 (Y 0 ), let C * > 0 be such that

( 1 with η = 1 .

 11 the definition of c(Λ) is recalled in Appendix 7.5) then the expression (19) coincides with the known expression relating the minimal time for null-controllability to the observation of the eigenvectors φ λ through the operator B * . proposition 3.1. Assume that A and B satisfy the assumptions of Theorem 1.If c(Λ) = 0, then, we have

theorem 4 . 1 .

 41 Let T ∈ (0, +∞]. Let p ∈ N * , r, ρ > 0 and N : (0, +∞) → R.

Corollary 4 . 1 .

 41 Let T ∈ (0, +∞]. Let p ∈ N * , r, ρ > 0 and N : (0, +∞) → R.

proposition 4 . 1 .

 41 Let p ∈ N * , r, ρ > 0 and N : (0, +∞) → R. Assume that Λ ∈ L w (p, ρ, N ) and let (G k ) k ∈ G(Λ, p, r, ρ) be an associated grouping. We also consider an integer η ≥ 1.

5 . 1 .

 51 this section we study various examples. In Section 5.1, we design 'abstract examples' to highlight the phenomenon described in Section 1.3.4: the condensation of eigenvectors can compensate the condensation of eigenvalues. More precisely we design an example which is null-controllable in arbitrary time but with an arbitrary condensation of the eigenvalues. We also give examples to illustrate the new settings covered by our analysis when the eigenvalues are algebraically multiple in the absence of a gap condition. The interest of these abstract examples is to highlight the different phenomena as the computations are straightforward. Finally, we provide in Section 5.2, actual examples of one dimensional coupled parabolic control systems that have motivated the present study. The precise analysis of null-controllability for those systems was not possible using existing results in the literature. Abstract examples: a possible compensation of condensation of eigenvalues.

proposition 5 . 1 .

 51 Let us consider the control system (1) with A and B given by (77)-(78).

Remark 5 . 3 .

 53 Let us consider in the same setting the evolution problem[START_REF] Allonsius | Boundary null-controllability of semi-discrete coupled parabolic systems in some multi-dimensional geometries[END_REF]

proposition 5 . 2 .

 52 Let us consider the control system (1) with A and B given by (81)-(82). Then, (83) T 0 (X) = 2β = 2 c(Λ).

proposition 5 . 4 .

 54 For any τ ∈ [0, +∞] there exist c 1 , c 2 ∈ L 2 (0, 1; R) such that the condensation index of the spectrum Λ of the operator A * satisfies c(Λ) = τ .

proposition 6 . 1 .

 61 Let U ⊂ C be a convex open set and x 0 , . . . , x n ∈ U be pairwise distinct. For any holomorphic function f : U → C, there exists a z ∈ Conv({x 0 , . .

1 e

 1 -t Ãx, z -1,1 * = e -tA x, z -1,1 * = e -tA x, z = x, e -tA * z = x, e -tA * z -1,1 * .

7. 2 .proposition 7 . 1 .

 271 Existence of a grouping for sequences satisfying the weak gap condition. For any Λ satisyfing (9), there exists at least one grouping in G Λ, p, ρ p , ρ .Proof. Let r = ρ/p. We set µ 1 = inf Λ. and we consider the p disjoint sets Λ ∩ (µ 1 , µ 1 + r], Λ ∩ (µ 1 + r, µ 1 + 2r], . . . , Λ ∩ (µ 1 + (p -1)r, µ 1 + pr].

Definition 7 . 1 .

 71 The divided differences are defined byf [x i ] := f i , ∀i ∈ 1, n ,and then recursively for any k ∈ 2, n , for any pairwise distinct i 1 , . . . , i k ∈ 1, n , byf [x i1 , . . . , x i k ] := f [x i1 , . . . , x i k-1 ] -f [x i2 , . . . , x i k ] x i1 -x i k .

proposition 7 . 2 .proposition 7 . 3 .

 7273 The divided differences are symmetric with respect to their arguments: for any k ∈ 1, n , for any pairwise distinct i 1 , . . . , i k ∈ 1, n and anyσ ∈ S({i 1 , . . . , i k }), f [x σ(i1) , . . . , x σ(i k ) ] = f [x i1 , . . . , x i k ].The following property states another (equivalent) definition of divided differences known as Newton formula. For any k ∈ 1, n , for any pairwise distinct i 1 , . . ., i k ∈ 1, n f [x i1 , . . . , x i k ] = k j=1 f [x ij ] l∈ 1,k =j (x ij -x i l ).difference, we can apply Lagrange theorem (see Proposition 7.4) to get the existence of a point z h ∈ Conv({y h 1 , . . . , y h N }) such that f [y h 1 , . . . , y h N ] = 1 (N -1)! f (N -1) (z h ).

proposition 7 . 11 .0

 711 For any multi-index µ ≤ α, there exists coefficients (θ µ j,l ) j,l depending only on x and µ, such thatf [x (µ) if µ j < l + 1, C |µ| i∈ 1,n =j |x i -x j | µi 1 (min i∈ 1,n =j |x i -x j |) µj -l-1 , if µ j ≥ l + 1.

•i∈ 1 ,n =j |x j -x i | µi 1 d

 11 Let us show by induction on N = |µ| that, for all j ∈ 1, n , with µ j ≥ l + 1µj -l-1 j

x 1 i∈ 1 ,n =j |x j -x i | µi 1 d 1 d 1 . 7 . 4 .

 1111174 i0 -x j , and thus|θ µ j,l | ≤ |θ µ-δ j j,l | |x i0 -x j | + |θ µ-δ i 0 j,l | |x i0 -x j | . Since |µ -δ j | = |µ -δ i0 | = N ,we can apply the induction hypothesis to bound the two terms in the right-hand side as follows|θ µ-δ j j,l | |x i0 -x j | ≤ C N -1 |θ µ-δ j i -x j | ≤ C N -1 |θ µ-δ i 0 j,l | |x j -x i0 | ≤ C N -1 i∈ 1,n =j |x j -x i | µi-δ i 0 i µj -l-1 j 1 |x j -x i0 | .Summing those two inequalities gives (124) with C N = 2C N -The supremum of T 0 (y 0 ).We prove here Propositon 1.3, that issup y0∈Y0 T 0 (y 0 ) = T 0 (Y 0 ).Proof. Since by definition T 0 (y 0 ) only depends on Span(y 0 ), it is actually equivalent to prove sup y0∈Y0 y0 -=1

P

  ,j -λ k,i ) , ∀k ≥ 1, ∀l ∈ 1, g k . with ψ k,j := P * Y 0 φ k,jB * φ k,j as defined in[START_REF] Bernstein | Leçons sur les Progrès Récents de la Théorie des Séries de Dirichlet[END_REF].Notice that, since y 0 -= 1, for any z ∈ X * , Span(y0) y, z -,= sup y∈Xy -=1 (y, y 0 ) -y 0 , z -, = y 0 , z -, .Thus, with those notations, we haveT 0 (y 0 ) = lim sup k→+∞ ln max l∈ 1,g k y 0 , x k,l -,

and P Y0 y 0 ∈

 0 Y 0 . Applying the Banach-Steinhaus theorem, this implies that sup k≥1 l∈ 1,g k xk,l * < +∞. Thus there exists C > 0 such that x k,l * ≤ Ce λ k,1 T for any k ≥ 1 and any 1 ≤ l ≤ g k . Finally this yields, T 0 (Y 0 ) = lim sup k→∞ ln max l∈ 1,g k x k,l * λ k,1 ≤ T.

Definition 7 . 4 . 2 .

 742 Assume that Σ satisfies The condensation index c(Σ) ∈ [0, +∞] is defined by c(Σ) = lim sup σ∈Σ σ→∞ -ln |E Σ (σ)| σ .

proposition 7 . 12 .

 712 Assume that Σ satisfies the assumptions of Definition 7.4 as well as the weak gap condition (9). Denote by (G k ) k≥1 a grouping satisfying the conditions of Definition 1.2. Then, c(Σ) = lim sup σ∈Σ σ→∞ -ln |P G [σ] (σ)| σ .

  Let α > β > 0 and Θ = {µ k , µ k + e -αµ k , µ k + e -βµ k ; k ∈ N * }. Then, c(Θ) = α + β.Proof. One can directly verify that the grouping defined byG k := µ k , µ k + e -αµ k , µ k + e -βµ ksatisfies the requirements given in Proposition 7.1. Then, direct computations lead to|P G k (µ k )| = |µ k -(µ k + e -αµ k )| |µ k -(µ k + e -βµ k )| = e -(α+β)µ k , |P G k (µ k + e -αµ k )| = |µ k + e -αµ k -µ k | |µ k + e -αµ k -(µ k + e -βµ k )| = e -(α+β)µ k 1 -e -(α-β)µ k ,and|P G k (µ k + e -βµ k )| = |µ k + e -βµ k -µ k | |µ k + e -βµ k -(µ k + e -αµ k )| = e -2βµ k 1 -e -(α-β)µ k .

  (of Proposition 2.3). This proof follows closely the one of [4,Section 4] and[START_REF] Khodja | Minimal time for the null controllability of parabolic systems: The effect of the condensation index of complex sequences[END_REF] Lemma 4.2]. From [4, Corollary 4.3], as Λ satisfies[START_REF] Avdonin | Riesz bases of exponentials and divided differences[END_REF], it comes that A(Λ, +∞)

This manuscript is for review purposes only.

The next result about divided differences is crucial to obtain the different estimates we need. It is known as Lagrange theorem.

proposition 7.4. Assume that V = R and that f ∈ C n-1 (Conv{x 1 , . . . , x n }).

For any k ∈ 1, n , for any pairwise distinct i 1 , . . . , i k ∈ 1, n , there exists a z ∈ Conv{x i1 , . . . , x i k } such that f [x i1 , . . . , x i k ] = f (k-1) (z) (k -1)! .

The divided differences naturally appear in polynomial interpolation problems as recalled in the following classical result.

proposition 7.5. The polynomial function P : R → V defined by

is the unique polynomial of degree less than n -1 such that

We recall a simple way to compute divided differences of a product which is known as the Leibniz rule.

proposition 7.6. Let g : R → R and (gf )

Finally, we deduce from the results above the following useful corollary.

Corollary 7.1. Assume that V is equipped with a norm • V . For any k ∈ 1, n and any pairwise distinct i 1 , . . . , i k ∈ 1, n , we have

where R = diam({x 1 , . . . , x n }).

Proof. Let P be the Lagrange interpolation polynomial defined in (118) and let i 1 , . . . , i k be fixed.

By the Hahn-Banach theorem, there exists φ ∈ V , such that φ V = 1 and

Additionally, by (119) and by linearity of φ, we know that

Applying Proposition 7.4 to x → φ, P (x) V ,V ∈ R we find that for some z ∈ Conv{x 1 , . . . , x n }, we have

.

This manuscript is for review purposes only.

Combining those identities, we arrive at

Let us compute the derivatives of P . Let C be the circle of center z and radius R in the complex plane. The Cauchy formula leads to

Then, the triangle inequality implies that for any w ∈ C,

which finally gives the result.

Generalization of divided differences.

Assume that V is a normed vector space.

Let x = (x 1 , . . . , x n ) ∈ R n be pairwise distinct real numbers and let α ∈ N n a multi-index such that α > 0. To such a multi-index we associate elements of V that we gather in a f α ∈ V |α| and that are indexed as follows

Definition 7.2. We set N = |α|. We say that a family of points (y h p ) p∈ 1,N , depending on a small parameter h > 0, is an approximation of the weighted family

• For each h > 0, the points y h 1 , . . . , y h N are pairwise distinct.

• There exist disjoint subsets P j ⊂ 1, N such that for any j ∈ 1, n , #P j = α j , and y h p ---→ h→0

x j , ∀p ∈ P j .

proposition 7.7. With the notation above, let F : R → V be any smooth function satisfying

For any approximation of the weighted family (x, α), the (usual) divided difference

weakly converges when h → 0 towards an element in V that depends only on x, α and f α . In particular it does not depend on the particular choice of F nor or the approximation families (y h p ) p .

This limit is called the generalized divided difference associated with the points x, the multi-index α and the values f α and is denoted by

or, in a more compact way, f [x (α) ].

Moreover, we extend this definition if some of the α j are 0, simply by not considering the corresponding points.

62

This manuscript is for review purposes only.

Remark 7.1. If the function F is chosen to take its values in a finite dimension space then the above convergence is actually strong. It is always possible to make this assumption, for instance by chosing F that takes its values in the subspace of V spanned by the elements f α .

Proof (of Proposition 7.7). The proof is done by recurrence on N .

• If N = 1, then we necessarily have n = 1 and α 1 = 1. The result is just a consequence of the continuity of F and we simply have f [x 1 ] = f 0 1 .

• Assume that the result holds for a given value of N and let us prove it for the value

-First case: If there is only one point x 1 . It means that n = 1 and α 1 = N + 1.

In this case, for any h > 0, and any ψ ∈ V , we use the Lagrange theorem to get the existence of a z ψ,h ∈ Conv({y h 1 , . . . , y h N +1 }) such that

Since, by assumption, all the points y h p converge to the same point x 1 , we have z ψ,h → x 1 and thus

-Second case: We assume that n > 1. By assumption there exists two distinct indices j 1 , j 2 ∈ 1, n and two distinct indices p 1 , p 2 ∈ 1, N + 1 such that y h p1 → x j1 and y h p2 → x j2 . By symmetry of the usual divided differences, we can always assume that p 1 = N and p 2 = N + 1. It follows that we can write

The recurrence assumption shows that the two terms in the numerator have weak limits that only depends on the points x, the multiplicities α and on the values f α , whereas the denominator y h N +1 -y h N converges to x j2 -x j1 which is not zero. The result follows.

The above construction also shows, as a by-product, the following rules to compute the generalized divided differences: for any µ ∈ N n such that µ ≤ α

, if µ j = 0 for all j = j, and for all j 1 = j 2 and µ j1 > 0,

Let us now give some useful properties that are the extension of the classical properties recalled in Section 7.3.1.

Definition 7.3. Let α ∈ N n be a multi-index, g α ∈ R |α| a set of real values associated with α and f α ∈ V |α| a set of elements of V associated with α.

63

This manuscript is for review purposes only.

We define (gf ) α ∈ V |α| to be the product set of values as follows:

(gf ) l j := l k=0 g k j f l-k j , ∀j ∈ 1, n , ∀l ∈ 0, α j -1 .

we have the Leibniz formula

Proof. By assumption, for each p ∈ 1, |α| , the multi-index µ p is obtained from µ p-1 by incrementing exactly one of its element. We denote by i p ∈ 1, n this index, and we define y h p := x ip + ph. It is easily seen that, for h > 0 small enough, those points are pairwise distinct.

Let F : R → V be a function satisfying (120) and G : R → R be a function satisfying (120) but with the values g α instead of f α . The usual Leibniz formula as well as the Definition 7.3 shows that the product GF exactly satisfies 1 l! (GF ) (l) (x j ) = (gf ) l j , ∀j ∈ 1, n , ∀l ∈ 0, α j -1 .

We can thus apply the Leibniz formula from Proposition 7. and then pass to the limit as h → 0 to obtain the claim.

proposition 7.9 (Lagrange theorem). Let x, α as before. We set N = |α|. With any f : R → R of class C N -1 , we associate the set of values f α ∈ R N by f l j := 1 l! f (l) (x j ), ∀j ∈ 1, n , ∀l ∈ 0, α j -1 .

Then, there exists a z ∈ Conv({x 1 , . . . , x n }) such that the generalized divided difference built on these data satisfies f [x (α) ] = 1 (N -1)! f (N -1) (z).

Proof. Let y h 1 , . . . , y h N be an approximation of the weighted family of points (x, α) as in Definition 7.2. By definition, the generalized divided difference f [x (α) ] is the limit as h goes to 0, of the usual divided difference f [y h 1 , . . . , y h N ]. For this last divided