Antoine Mouzard

Weyl law for the Anderson Hamiltonian on a two-dimensional manifold

Keywords: Anderson Hamiltonian, Paracontrolled calculus, White noise, Schrödinger operator

Weyl law for the Anderson Hamiltonian on a two-dimensional manifold

Antoine Mouzard

Introduction

The study of singular stochastic Partial Differential Equations (PDEs) has rapidly grown over the last decade. Following the theory of Lyons' rough paths and Gubinelli's controlled paths developed for singular stochastic ordinary differential equations, new tools have appeared to describe solutions of such PDEs that share the same philosophy. One first constructs a random space of functions/distributions from the noise through a renormalisation step; this is purely probabilistic. One then solves the PDE with classical methods on this random space; this is purely deterministic. The litterature is also growing and two different approaches have emerged. The first is based on a local description of distributions which satisfies a precise algebraic structure in order to reassemble into global objects; this is the theory of regularity structures as devised by Hairer in [START_REF] Hairer | A theory of regularity structures[END_REF]. The second approach works directly with global objects and uses tools from harmonic analysis to study products; this is the paracontrolled calculus designed by Gubinelli, Imkeller and Perkowski in [START_REF] Gubinelli | Paracontrolled distributions and singular PDEs[END_REF]. In both cases, the equation dictates via a fixed point a space of solutions built from the rough source term of the PDE. There exists interesting relations between the local and the global points of views, see for example the works [START_REF] Bailleul | Paracontrolled calculus and regularity structures (1)[END_REF][START_REF] Bailleul | Paracontrolled calculus and regularity structures[END_REF][START_REF] Martin | A Littlewood-Paley description of modelled distributions[END_REF]. As far as the renormalisation step is concerned, one has to give a meaning to a number of ill-defined functionnals of the noise; this is how singular products are dealed with. If the list of such terms is given by the equation, their construction can be performed independantly of the resolution of the PDE.

To work on manifolds, one has to adapt these methods. For the local approach, Dahlqvist, Diehl and Driver have adapted regularity structures using local charts to study the parabolic Anderson model on Riemann surfaces, see [START_REF] Dahlqvist | The parabolic Anderson model on Riemann surfaces[END_REF]. For the global approach, Bailleul, Bernicot and Frey in [START_REF] Bailleul | Heat semigroup and singular PDEs[END_REF][START_REF] Bailleul | Spacetime paraproducts for paracontrolled calculus, 3d-PAM and multiplicative Burgers equations[END_REF] used harmonic analysis tools built from the heat semigroup instead of Fourier analysis and developed paracontrolled calculus on manifolds. As in the initial work [START_REF] Gubinelli | Paracontrolled distributions and singular PDEs[END_REF] of Gubinelli, Imkeller and Perkwoski, this was only a first order calculus and it constrained the roughness one could deal with. Bailleul and Bernicot then generalised it to a high order paracontrolled calculus in [START_REF] Bailleul | High order paracontrolled calculus[END_REF] and extended the range of regularity one can consider, as far as the analytical step of the problem is concerned, again working on manifolds.

The Anderson Hamiltonian is given by

H := ∆ + ξ
where ξ is a space white noise. It is for example involved in the study of evolution equations such as the heat equation with random multiplicative noise ∂ t u = ∆u + uξ called the Parabolic Anderson model. It first appeared in [START_REF] Anderson | Absence of Diffusion in Certain Random Lattices[END_REF] as a description of a physical phenomena involving quantum-mechanical motion with an effect of mass concentration called Anderson localization. Since then, a large class of Schrödinger operators with random potential have beed investigated, for example as a description of random dynamics in random environment. For particular models, see the book [START_REF] König | The parabolic Anderson model : Random walk in random potential[END_REF] of König for a complete survey in a discrete space setting. One of the most important question on this type of operator lies in the line of the first work by Anderson, that is whether it exhibits a localisation phenomenon. This corresponds to localisation in space of eingenfunctions associated to a part of the spectrum. While this was mainly answered in the discrete setting, very accurate results were obtained by Dumaz and Labbé in [START_REF] Dumaz | Localization of the continuous Anderson Hamiltonian in 1-d[END_REF] for the continuous Anderson Hamiltonian in one dimension, see also their recent works [START_REF] Dumaz | The delocalized phase of the Anderson Hamiltonian in 1-d[END_REF][START_REF] Dumaz | Localization crossover for the continuous Anderson Hamiltonian in 1-d[END_REF]. The litterature about localisation of such operators is now quite vast both in mathematics and physics. Formally, the Anderson Hamiltonian with white noise potential is the continuum limit of the discrete Anderson Hamiltonian with independent and identically distrubuted random potential with finite variance. See for example the work [START_REF] Martin | Paracontrolled distributions on Bravais lattices and weak universality of the 2d parabolic Anderson model[END_REF] by Martin and Perkowski. In dimension 1, the noise is regular enough for the multiplication to make sense and the operator has been constructed by Fukushima and Nakao in [START_REF] Fukushima | On Spectra of the Schrödinger Operator with a White Gaussian Noise Potential[END_REF] without renormalisation using Dirichlet space methods. In two dimensions using paracontrolled calculus, Allez and Chouk were the first to construct the operator on the torus, see [START_REF] Allez | The continuous Anderson Hamiltonian in dimension two[END_REF]. They introduced the space of strongly paracontrolled distributions to get an operator from L 2 to itself with a renormalisation procedure and proved self-adjointness with pure point spectrum. They gave bounds on its eigenvalues and a tail estimate for the largest one. They also studied the large volume limit and gave a bound on the rate of divergence. Then Labbé constructed the operator in dimension ≤ 3 in [START_REF] Labbé | The continuous Anderson Hamiltonian in d ≤ 3[END_REF] with different boundary coundition using regularity structures. It relies on a reconstruction theorem in Besov spaces from his work [START_REF] Hairer | The reconstruction theorem in Besov spaces[END_REF] with Hairer. He also proved self-adjointess with pure point spectrum and gave tail estimate for all the eigenvalues as well as bounds for the large volume limit. Chouk and van Zuijlen also studied the large volume limit in two dimensions, see [START_REF] Chouk | Asymptotics of the eigenvalues of the Anderson Hamiltonian with white noise potential in two dimensions[END_REF]. Finally Gubinelli, Ugurcan and Zachhuber constructed in [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF] the operator in dimension 2 and 3 on the torus using a different approach. With a finer description of the paracontrolled structure, they showed density of the domain in L 2 before studying the operator. They also proved self-adjointness with pure point spectrum considering the bilinear form associated to H and considered evolution PDEs associated to the Anderson Hamiltonian such as the Schrödinger equation or the wave equation. Zachhuber used this approach in [START_REF] Zachhuber | Strichartz estimates and low-regularity solutions to multiplicative stochastic NLS[END_REF] to prove Strichartz estimate in two dimensions, the problem for d = 3 being the use of a Hopf-Cole type transformation to construct the domain.

To the best of our knowledge, the present work is the first to deal with the construction of the Anderson Hamiltonian on a manifold. In particular, the paracontrolled approach with the heat semigroup deals naturaly with Sobolev spaces on a manifold while we are not aware of any adaptation of the work [START_REF] Hairer | The reconstruction theorem in Besov spaces[END_REF] of Hairer and Labbé in a manifold setting. We are able to recover geometric information on the manifold from the spectral properties of the Anderson Hamiltonian as one can do from the Laplacian. For example, we recover the volume of M via a Weyl law from the estimates on the spectrum. This raises many interesting associated questions. As far as PDEs are concerned, it also appears in a number of singular SPDEs of interest. For example, the Schrödinger equation has been studied on manifolds by Burq, Gérard and Tzvetkov in [START_REF] Burq | Strichartz inequalities and the nonlinear Schrödinger equation on compact manifolds[END_REF] where they prove Strichartz inequalities. Similar questions for the stochastic version of the equation are natural to ask on manifold and has been investigated in the flat case by Zachhuber in [START_REF] Zachhuber | Strichartz estimates and low-regularity solutions to multiplicative stochastic NLS[END_REF]. This question was investigated in [START_REF] Mouzard | Strichartz inequalities with white noise potential on compact surfaces[END_REF] by Mouzard and Zachhuber on a two-dimensional manifold with or without boundary, both for the Schrödinger and the wave equations. In particular, the result for the wave equation strongly relies on the Weyl-type law obtained in the present work, even in the flat case of the torus.

Main results

In this work, we construct the Anderson Hamiltonian on a two-dimensional manifold using the high order paracontrolled calculus. We adapt the space-time construction [START_REF] Bailleul | Spacetime paraproducts for paracontrolled calculus, 3d-PAM and multiplicative Burgers equations[END_REF][START_REF] Bailleul | High order paracontrolled calculus[END_REF] of Bailleul, Bernicot and Frey to the spatial setting and work with Sobolev spaces; in particular this work is self-contained and can serve as an introduction to the work [START_REF] Bailleul | High order paracontrolled calculus[END_REF] on the high order paracontrolled calculus. The simpler spatial setting forms a gentle introduction to grasp the space-time paracontrolled calculus, the only technical difficulty being the use of Sobolev spaces in addition to the Hölder spaces. We emphasize that these tools are of interest on their own in the study of singular elliptic PDEs on manifolds and somewhat flexible to use. As application, it yields existence and uniqueness to the nonlinear Schrödinger equation with multiplicative white noise on a two-dimensional manifold using a Brezis-Gallouët type inequality. In particular, this work removes the need of the "strongly paracontrolled distributions" introduced by Allez and Chouk and used by Gubinelli, Ugurcan and Zachhuber with a second order expansion rather than an ad-hoc modification of the first order expansion.

Given a regularisation ξ ε of the noise, we consider the renormalised Hamiltonians

H ε := -∆ + ξ ε -c ε with c ε := E[Π(X ε , ξ ε )] and X ε := -L -1 ξ ε .
While in the flat case of the torus c ε is a renormalisation constant, it is only a renormalisation function in the general setting of a two-dimensional manifold. The renormalisation of the noise yields a enhanced noise Ξ ∈ X α , this is done in Section 2.1. The Anderson Hamiltonian is constructed as the limit of this operator in the resolvent norm sense and we are able to get upper and lower bounds for its eigenfunctions from the eigenfunctions of the Laplacian. This is the content of Theorems 2.11, 2.12 and 2.14.

Theorem. There exists a self-adjoint operator H on L 2 (M ) with pure point spectrum λ n (Ξ) n≥1 such that

H = lim ε→0 -∆ + ξ ε -c ε
in the resolvent sense. Moreover, there exists constants m - δ (Ξ), m + δ (Ξ) such that

λ n -m - δ (Ξ) ≤ λ n (Ξ) ≤ (1 + δ)λ n + m + δ (Ξ)
for any δ ∈ (0, 1) and n ≥ 1.

This implies the following two corollaries, that is the upper bouds on the repartition functions of the eigenfunctions and the almost sure Weyl-type law. This corresponds to corollaries 2.13 and 2.16.

Corollary.

For any n ≥ 1 and λ ∈ R, we have

1 -me -h(λ-2λn) 1 12 ≤ P(λ n (Ξ) ≤ λ) ≤ me -h(λn-λ) 1 5
where m = E e h Ξ X α and h > 0 a constant.

Corollary. We have almost surely

lim λ→∞ λ -1 |{λ n (Ξ) ≤ λ} = Vol(M ) 4π .

Organisation of the paper

In the first section, we introduce the approximation theory based on the heat semigroup and use it to build the paracontrolled calculus. The second section is devoted to the construction and study of the Anderson Hamiltonian H on a manifold in two dimensions. We show self-adjointness with pure point spectrum and provide lower and upper bounds for the eigenvalues. We finally study the cubic Schrödinger equation in Section 2.4. Appendix A contains all the technical details of the approximation theory and Appendix B gives the proof of different continuity estimates for the paracontrolled calculus.

The main ingredients for this work are the following. We adapt the work of Bailleul and Bernicot [START_REF] Bailleul | Spacetime paraproducts for paracontrolled calculus, 3d-PAM and multiplicative Burgers equations[END_REF] and give the tools of the spatial high order paracontrolled calculus. These are of interest in themselves to solve other elliptic problems, on manifold or not, and are very flexible. As Gubinelli, Ugurcan and Zachhuber in [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF], our method relies on the almost duality property between the resonant term and the paraproduct. Finally, we introduce a truncated paraproduct P s to describe product on adapted scales with its companion P s that describes associated mild formulation.

Acknowledgements :

The author would like to thank Massimiliano Gubinelli for the invitation to the Haussdorf Institute of Mathematics (HIM) where he enjoyed a lot of fruitful mathematical discussions with him and others. In particular, the author is grateful to Immanuel Zachhuber for interesting discussions that motivated the present work. The author would also like to thank Nikolay Tzvetkov and Ismaël Bailleul for readings of the preliminary drafts of this version.

-Heat semigroup and paracontrolled calculus

On a manifold M , the heat semigroup P = (P t ) t>0 := (e tL ) t>0 associated to a nice enough second order differential operator L can be used to regularize distributions in D (M ) while the tools of Fourier analysis of the flat case can not be used. One can then consider the Calderón decomposition as an analog of the Paley-Littlewood decomposition with a continuous scaling parameter and Q t := -t∂ t P t acting like a localizer on "frequency" of order t -1 2 . After giving the geometric framework, we introduce the standard families of operators we shall use to define the Besov spaces on M . We then construct the paraproducts P and P with a number of tools of the high order paracontrolled calculus needed to study elliptic singular PDEs.

-Geometric framework

Let (M, d, µ) be a complete volume doubling measured Riemannian manifold. We assume M compact so spatial weight are not needed; everything in this section should work in the unbounded setting of [START_REF] Bailleul | Spacetime paraproducts for paracontrolled calculus, 3d-PAM and multiplicative Burgers equations[END_REF]. All the kernels we consider are with respect to this measure µ. Let (V i ) 1≤i≤v be a family of smooth vector fields identified with first order differential operators on M . Consider the associated second order operator L given by

L = - v i=1 V 2 i .
We assume that L is elliptic. In particular, it implies that the vector fields (V i ) 1≤i≤v span smoothly at every point of M the tangent space and the existence of smooth functions (γ i ) 1≤i≤v such that for any f ∈ C 1 (M, R) and x ∈ M , we have

∇f (x) = v i=1 γ i (x)V i (f )(x)V i (x).
It also implies that L is sectorial in L 2 with kernel the constant functions, it has a bounded H ∞calculus on L 2 and -L generates a holomorphic semigroup (e -tL ) t>0 on L 2 , see [START_REF] Duong | Functional calculi of second-order elliptic partial differential operators with bounded measurable coefficients[END_REF]. Given any collection I = (i 1 , . . . , i n ) ∈ {1, . . . , v} n , we denote by V I := V in . . . V i1 the differential operator of order |I| := n. Under the smoothness and ellipticity conditions, the semigroup has regularity estimate at any order, that is (t |I| 2 V I )e -tL and e -tL (t |I| 2 V I ) have kernels K t (x, y) for any t > 0 and x, y ∈ M that satify the Gaussian estimates

K t (x, y) µ B(x, √ t) -1 e -c d(x,y) 2 t and for x ∈ M K t (x, y) -K t (x , y) d(x, x ) √ t µ B(x, √ t) -1 e -c d(x,y) 2 t
for d(x, x ) ≤ √ t and a constant c > 0. The range of application contains the case of a bounded domain with its Laplacian associated with periodic or Dirichlet boundary conditions if the boundary is sufficiently regular, see again [START_REF] Duong | Functional calculi of second-order elliptic partial differential operators with bounded measurable coefficients[END_REF]. Note that the Laplacian can indeed be written in the Hörmander form, see Strook's book [START_REF] Stroock | An Introduction to the Analysis of Paths on a Riemanian Manifold[END_REF] for example. The operator L : H 2 ⊂ L 2 → L 2 is not invertible since its kernel contains constant function however it is invertible up to a smooth error term. Indeed, setting

L -1 := 1 0 e -tL dt,
we have L • L -1 = Id up to the regularizing operator e -L .

-Approximation theory

All computations below make sense for a choice of large enough integers b and that are fixed in any application, we also assume b even. Given x, y ∈ M and t ∈ (0, 1], we define the Gaussian kernel

G t (x, y) := 1 µ B(x, √ t) 1 + c d(x, y) 2 t
with c > 0 a constant. We do not emphasize the dependance on the postive constant c and abuse notation by writing the same letter G t for two functions corresponding to two different values of the constant. We have for any s, t ∈ (0, 1]

M G t (x, y)G s (y, z)dy G t+s (x, z).
A choice of constant large enough ensure that

sup t∈(0,1] sup x∈M M G t (x, y)dy < ∞.
This implies that any linear operator with a kernel pointwisely bounded by G t is bounded in L p (M ) for every p ∈ [1, ∞]. The family (G t ) t∈(0,1] is our reference kernel for Gaussian operator; this is the letter 'G' in the following definition.

Definition. We define G as the set of families (P t ) t∈(0,1] of linear operator on M with kernels pointwisely bounded by

|K Pt (x, y)| G t (x, y)
given any x, y ∈ M .

We consider two such families of operators (Q

(Q t ) t∈(0,1]
is a kind of localizer on 'frequency' of order t -1 2 as one can see with the parabolic scaling of the Gaussian kernel. In the flat framework of the torus, this can be explicitly written using Fourier theory. These tools also enjoy cancellation properties as Fourier projectors however it is not as precise since the operators involved here are not locally supported. For example, the following simple computation show that the composition

Q (b) t • Q (b) s ts (t + s) 2 b Q (2b) t+s
is small for s t or t s but not equal to 0. The importance of the parameter b appears here as a 'degree' of cancellation. One can also see that in the fact that for any polynomial function p of degree less than 2b, we have P t p = 0 for any t ∈ (0, 1]. We now define the standard family of Gaussian operators with cancellation that we shall use in this work.

Definition. Let a ∈ 0, 2b . We define the standard collection of operators with cancellation of order a as the set StGC a of families (t

|I| 2 V I )(tL) j 2 P (c) t t∈(0,1]
with I, j such that a = |I| + j and c ∈ 1, b . These operators are uniformly bounded in L p (M ) for every p ∈ [1, ∞] as functions of the parameter t ∈ (0, 1]. In particular, a standard family of operator Q ∈ StGC a can be seen as a bounded map t → Q t from (0, 1] to the space of bounded linear operator on L p (M ). We also set

StGC [0,2b] := 0≤a≤2b StGC a .
Since the first order differential operators V i do not a priori commute with each other, they do not commute with L and we introduce the notation

V I φ(L) • := φ(L)V I
for any function φ such that φ(L) is defined in order to state the following cancellation property. This is not related to any notion of duality in general. In particular, L is not supposed self-adjoint here.

Proposition . . Given a, a ∈ 0, 2b , let Q 1 ∈ StGC a and Q 2 ∈ StGC a . Then for any s, t ∈ (0, 1], the composition

Q 1 s • Q 2•
t has a kernel pointwisely bounded by Proof : Let t ∈ (0, 1]. We have

K Q 1 s •Q 2• t (x, y)   s t a 2 1 s<t + t s a 2 1 s≥t   G t+s (x, y) ts (t + s)
Q 1 t = t a 2 V I L j 2 P (c) t and Q 2 t = t a 2 V I L j 2 P (c ) t
with c, c ∈ 1, b , a = |I| + j and a = |I | + j . For any t, s ∈ (0, 1], the composition is given by

Q 1 s • Q 2• t = s a 2 t a 2 V I L j+j 2 P (c) s P (c ) t V I = s a 2 t a 2 (t + s) a+a 2 (t + s) a+a 2 V I L j+j 2 P (c) s P (c ) t V I
and this yields

K Q 1 s •Q 2• t (x, y) s a 2 t a 2 (t + s) a+a 2 G t+s (x, y) s t a 2 1 s<t + t s a 2
1 s≥t G t+s (x, y).

The last estimate follows from a direct computation.

Operators with cancellation but not in this standard form also appear in the description of solutions to PDEs. This is the role of the set GC a of the following definition.

Definition. Let a ∈ 0, 2b . We define the subset GC a ⊂ G as families (Q t ) t∈(0,1] of operators with the following cancellation property. For any s, t ∈ (0, 1] and standard family S ∈ StGC a with a ∈ a, 2b , the operator Q s • S • t has a kernel pointwisely bounded by

K Qs•S • t (x, y) ts (t + s) 2 a 2 G t+s (x, y).
The set StGC can be used to define Besov spaces on a manifold. For any f ∈ L p (M ) with p ∈ [1, ∞[ or f ∈ C(M ), we have the following reproducing Calderón formula

f = lim t→0 P (b) t f = 1 0 Q (b) t f dt t + P (b) 1 f.
We interpret it as an analog to the Paley-Littlewood decomposition of f on a manifold but with a continuous parameter. Indeed, the measure dt t gives unit mass to the dyadic intervals [2 -(i+1) , 2 -i ] with the operator Q (b) t as a kind of multiplier roughly localized at frequencies of size t -1 2 . This motivates the following definition.

Definition. Given any p, q ∈ [1, ∞] and α ∈ (-2b, 2b), we define the Besov space B α p,q (M ) as the set of distribution f ∈ D (M ) such that

f B α p,q := e -L f L p (M ) + sup Q∈StGC k |α|<k≤2b t -α 2 Q t f L p x L q (t -1 dt) < ∞.
Remark : As far as regularity is concerned, a limitation appears with this definition of B α p,q since we can only work with regularity exponent α ∈ (-2b, 2b). This is only technical and b can be taken as large as needed.

The Hölder spaces C α := B α ∞,∞ and Sobolev spaces H α := B α 2,2 are of particular interest with

f C α := e -L f L ∞ + sup Q∈StGC k |α|<k≤2b sup t∈(0,1] t -α 2 Q t f L ∞ x and f H α := e -L f L 2 + sup Q∈StGC k |α|<k≤2b 1 0 t -α Q t f 2 L 2 x dt t 1 2
. This is indeed a generalisation of the classical Hölder spaces as stated in the following proposition. We shall denote C α the classical spaces of Hölder functions with the norm

f C α := f L ∞ + sup x =y |f (x) -f (y)| d(x, y) α
for 0 < α < 1. Note that for any integer regularity exponent, C α = C α since C 1 is the space of Lipschitz functions. The proof of the following proposition is left to the reader, it works exactly as Proposition 5 in [START_REF] Bailleul | Spacetime paraproducts for paracontrolled calculus, 3d-PAM and multiplicative Burgers equations[END_REF].

Proposition. For any α ∈ (0, 1), we have C α = C α and the norms • C α and • C α are equivalent.

We have an analog result for Sobolev spaces however one has to be careful in the case of a manifold with boundary. The semigroup is obtained with Dirichlet conditions hence the associated Sobolev spaces are the analog of the classical H α 0 spaces. We keep the notation H α but the reader should keep that in mind.

Given a distribution f ∈ C α and Q ∈ StGC k , we have by definition a bound for Q t f ∞ only for |α| < k. If f is a distribution and not a function, the quantity diverges and we still have the estimate for all k; this will be important to keep an accurate track of the regularity. The same holds for negative Sobolev spaces.

Proposition . . Let -2b < α < 0 and P ∈ StGC k with k ∈ 0, b . For f ∈ C α , we have sup t∈(0,1] t -α 2 P t f L ∞ 1 k -α f C α .
For f ∈ H α , we have

t -α 2 P t f L 2 x L 2 (t -1 dt) 1 k -α f H α . Proof : Since P ∈ StGC k with k ∈ 0, 2b , there exist I = (i 1 , . . . , i n ), j ∈ N and c ∈ 1, b such that k = |I| + j and P t = (t |I| 2 V I )(tL) j 2 P (c) t . If |α| < k, the result holds by definition of C α . If |α| ≥ k, we have P t f = (t |I| 2 V I )(tL) j 2 1 t Q (c) s f ds s + P (c) 1 f = 1 t t s k 2 (s |I| 2 V I )(sL) j+c 2 P (1) s f ds s + R t f = 1 t t s k 2 Q s f ds s + R t f with Q s := (s |I| 2 V I )(sL) j+c 2 P (1) s ∈ StGC k+c and R t := (t |I| 2 V I )(tL) j 2 P (c) 1 . The term R t f is bounded because of the smoothing operator P (c) 1 . Since c ≥ 1, Q belongs at least to StGC k+1 hence if |α| < k + 1 we have t -α 2 P t f L ∞ ≤ t -α 2 1 t t s k 2 Q s f L ∞ ds s ≤ f C α 1 t t s k-α 2 ds s ≤ f C α 2 k -α
and this yields the result using that α < 0 ≤ k hence k -α > 0. If |α| ≥ k + 1, using the same integral representation for Q and an induction completes the proof of the L ∞ -estimate. For the L 2 -estimate, we interpolate between L 1 and L ∞ as in Appendix A to get

t -α 2 P t f L 2 L 2 (t -1 dt) ≤ t -α 2 1 t t s k 2 Q s f L 2 ds s L 2 (t -1 dt) ≤ 2 k -α f H α .
One can see that the bound diverges as α goes to 0 if the operator does not encode any cancellation, that is k = 0. In the case α = 0, we have

P t f L ∞ f L ∞ hence the L ∞ -bound holds. However the L 2 -bound is not satisfied since P t f L 2 f L 2 only implies 1 0 P t f 2 L 2 dt t ≤ f 2 L 2 1 0 dt t = ∞.
This will explain an important difference for paraproducts on negative Hölder and Sobolev spaces as one can see with Propositions 1.3 and 1.4.

-Intertwined paraproducts

We use the standard family of Gaussian operators to study the product of distributions as one can do using the Paley-Littlewood decomposition in the flat case; this lead to the definition of the paraproduct P and the resonant term Π that describe products. Then we introduce the paraproduct P intertwined with P to describe solutions of elliptic PDEs.

-Paraproduct and resonant term

One can define the product of a distributions f ∈ D (M ) with a smooth function g ∈ D(M ). If however the distribution f belongs to a Hölder space C α with α < 0, one might hope to do better.

It is indeed the case as we can see with the next theorem which is nothing more than Young's integration condition.

Theorem. The multiplication (f, g) → f g extends in a unique bilinear operator from

C α × C β to C α∧β if and only if α + β > 0.
We are however interested in the case α + β < 0 when dealing with singular stochastic PDEs, as we are interested to stochastic ODEs where Young's condition is not verified. Following [START_REF] Gubinelli | Paracontrolled distributions and singular PDEs[END_REF], Bailleul, Bernicot and Frey in [START_REF] Bailleul | Heat semigroup and singular PDEs[END_REF][START_REF] Bailleul | Spacetime paraproducts for paracontrolled calculus, 3d-PAM and multiplicative Burgers equations[END_REF][START_REF] Bailleul | High order paracontrolled calculus[END_REF] have defined two bilinear operators P f g and Π(f, g) such that we have the formal decomposition of the product of two distributions as

f g = P f g + Π(f, g) + P g f
where the paraproducts P f g and P g f are well-defined for any distibutions f, g ∈ D (M ). Of course, this means that Π(f, g) does have a meaning for f ∈ C α and g ∈ C β if and only if α + β > 0; this is the resonant term. We want this decomposition to keep an accurate track of the regularity of each terms. More precisely, P f g and Π(f, g) should belong to C α+β if α < 0 while P g f to the less regular space C α as it is the case for the torus. We construct in this work such paraproduct and resonant term for space distributions on our manifold M , we mainly follow [START_REF] Bailleul | Spacetime paraproducts for paracontrolled calculus, 3d-PAM and multiplicative Burgers equations[END_REF] in the simpler spatial setting.

Let f, g ∈ D (M ). Formally, we have

f g = lim t→0 P (b) t P (b) t f • P (b) t g = 1 0 Q (b) t P (b) t f • P (b) t g + P (b) t Q (b) t f • P (b) t g + P (b) t P (b) t f • Q (b) t g dt t + P (b) 1 P (b) 1 f • P (b) 1 g .
The last term being smooth, it does not bother us. Remark that the choice of the constant "1" is arbitrary and it might be useful to change it, as one can see with the construction of the Anderson Hamiltonian. The family P (b) does not encode any cancellation while Q (b) encodes cancellation of order 2b so each terms in the integral have one operator with a lot of cancellations and two with none. Since we do not have nice estimates for these terms, we want to transfer some of the cancellation from Q (b) to the P (b) in each term. To do so, we use the Leibnitz rule

V i (f g) = V i (f )g + f V i (g).
For example, we have

1 0 P (b) t (tV 2 i )Q (b-1) t f • P (b) t g dt t = 1 0 P (b) t ( √ tV i ) ( √ tV i )Q (b-1) t f • P (b) t g dt t - 1 0 P (b) t ( √ tV i )Q (b-1) t f • ( √ tV i )P (b) t g
dt t so if we denote by (c 1 , c 2 , c 3 ) the cancellation of the three operators in the integral, we have

(0, 2b, 0) = (1, 2b -1, 0) + (0, 2b -1, 1).
This shows that we will not be able to have cancellation for all three operators at the same time but at least two. This is where the notation Q • comes into play and multiple uses of this trick allows to decompose the product as

f g = a∈A b Q∈StGC a b Q 1 0 Q 1• t Q 2 t f • Q 3 t g dt t where Q = (Q 1 , Q 2 , Q 3 ), StGC a = StGC a1 × StGC a2 × StGC a3 , A b = (a 1 , a 2 , a 3 ) ∈ N 3 ; a 1 + a 2 + a 3 = 2b and a 1 , a 2 or a 3 = b and b Q ∈ R is a real coefficient associated to Q.
In particular, only one of the a i in a ∈ A b can be less than b 2 and this gives us three terms P f g, P g f and Π(f, g) such that

f g = P f g + Π(f, g) + P g f + P (b) 1 P (b) 1 f • P (b) 1 g .

Definition.

Given two distributions f, g ∈ D (M ), we define the paraproduct and the resonant term as

P f g := a∈A b ;a2< b 2 Q∈StGC a b Q 1 0 Q 1• t Q 2 t f • Q 3 t g dt t .
and

Π(f, g) := a∈A b ;a2,a3≥ b 2 Q∈StGC a b Q 1 0 Q 1• t Q 2 t f • Q 3 t g dt t .
In particular, P f g is a linear combination of

1 0 Q 1• t P t f • Q 2 t g dt t
and Π(f, g) of

1 0 P • t Q 1 t f • Q 2 t g dt t with Q 1 , Q 2 ∈ StGC b 2 and P ∈ StGC [0,b]
. We insist that in the following P will denote an operator with possibly no cancellations while Q will denote an operator with cancellations of order at least b 2 .

These operators enjoy the same continuity estimates as their Fourier counterparts from which one can recover Young's condition as stated in the following Proposition.

Proposition . . Let α, β ∈ (-2b, 2b) be regularity exponents.

If α ≥ 0, then (f, g) → P f g is continuous from C α × C β to C β . If α < 0, then (f, g) → P f g is continuous from C α × C β to C α+β . If α + β > 0, then (f, g) → Π(f, g) is continuous from C α × C β to C α+β .
Proof : Let us first consider the case α < 0 and let Q ∈ StGC r with r > |α + β|. Recall that P f g is a linear combination of terms of the form

1 0 Q 1• t P t f • Q 2 t g dt t with Q 1 , Q 2 ∈ StGC b 2 and P ∈ StGC [0,b] . Since α < 0, 1.2 gives 1 0 Q s Q 1• t P t f • Q 2 t g dt t 1 0 ts (t + s) 2 r 2 f C α g C β t α+β 2 dt t s α+β 2 f C α g C β for any s ∈ (0, 1) hence P f g ∈ C α+β . For α ≥ 0, we consider Q ∈ StGC r with r > |β|. In this case, we have |P t f | ≤ f C α for all t ∈ (0, 1) so 1 0 Q s Q 1• t P t f • Q 2 t g dt t s β 2 f C α g C β hence P f g ∈ C β .
For the resonant term, let Q ∈ StGC r with r > |α + β|. We have

1 0 Q s P • t Q 1 t f • Q 2 t g dt t f C α g C β s 0 t α+β 2 dt t + 1 s s t r 2 t α+β 2 dt t s α+β 2 f C α f C β using that α + β > 0 hence Π(f, g) ∈ C α+β .
We also have estimates for the Sobolev spaces whose proofs are given in Proposition B.1 from Appendix B.

Proposition . . Let α, β ∈ (-2b, 2b) be regularity exponents.

If α > 0, then (f, g) → P f g is continuous from C α × H β to H β and from H α × C β to H β . If α < 0, then (f, g) → P f g is continuous from C α × H β to H α+β and from H α × C β to H α+β . If α + β > 0, then (f, g) → Π(f, g) is continuous from H α × C β to H α+β .
In particular, this implies that (f, g) → P f g is continuous from L 2 × C β to H β-δ for all δ > 0. For Sobolev spaces, there is a small loss of regularity and one does not recover the space H β while this does not happen for Hölder spaces. This comes from the remark following Proposition 1.2.

As in the works [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF][START_REF] Zachhuber | Strichartz estimates and low-regularity solutions to multiplicative stochastic NLS[END_REF] of Gubinelli, Ugurcan and Zachhuber, one last property of P and Π in terms of Sobolev spaces is that P is almost the adjoint of Π when L is self-adjoint in the sense that the difference is more regular. A careful track of the previous computation show that for all a ∈ {(0, b, b), (b, 0, b), (b, b, 0)} and

Q ∈ StGC a , we have b Q = 0 except for Q = (P (b) t , Q (b/2) t , Q (b/2) t ), (Q (b/2) t , P (b) t , Q (b/2) t ) or (Q (b/2) t , Q (b/2) t , P (b) t ) where b Q = 1. Define the correc- tor for almost duality as A(a, b, c) := a, Π(b, c) -P a b, c . Proposition . . Assume L self-adjoint. Let α, β, γ ∈ (-2b, 2b) such that β+γ < 1 and α+β+γ ≥ 0. If α < 1, then (a, b, c) → A(a, b, c) extends in a unique trilinear operator from H α × C β × H γ to R. Proof : A(a, b, c) is a linear combination of 1 0 a, P 1• t Q 1 t b • Q 2 t c -Q 3• t P 2 t a • Q 4 t b , c dt t with P 1 , P 2 ∈ StGC [0,b] and Q 1 , Q 2 , Q 3 , Q 4 ∈ StGC b 2
. We first consider P 1 , P 2 ∈ StGC 0 . By construction of the paraproduct and the resonant term, we have P 1 = P 2 = P (b) =: P and

Q 1 = Q 2 = Q 3 = Q 4 = Q (b/2) =: Q hence we consider 1 0 a, P t Q t b • Q t c -Q t P t a • Q t b , c dt t .
Since L is self-adjoint, P t and Q t are too and we have

1 0 a, P t Q t b • Q t c dt t = 1 0 P t a, Q t b • Q t c dt t = 1 0 P t a • Q t b, Q t c dt t = 1 0 Q t P t a • Q t b , c dt t
hence the difference is equal to 0. Let us now consider the terms with P 1 , P 2 ∈ StGC [1,b] and bound each of them independently. Since α + β + γ ≥ 0, we have

1 0 a, P 1• t Q 2 t b • Q 3 t c dt t a H α 1 0 P 1• t Q 2 t b • Q 3 t c dt t H β+γ a H α b C β c H γ
with β + γ < 1 and using α ∈ (0, 1) we have

1 0 Q 3• t P 2 t a • Q 4 t b , c dt t 1 0 Q 3• t P 2 t a • Q 4 t b dt t H α+β c H γ a H α b C β c H γ
which completes the proof since α + β + γ ≥ 0.

-Intertwined paraproducts

The description of solution to elliptic PDEs involving L using paracontrolled calculus necessitate to study how L and P interacte with each other. Following Bailleul, Bernicot and Frey in [START_REF] Bailleul | Spacetime paraproducts for paracontrolled calculus, 3d-PAM and multiplicative Burgers equations[END_REF], we want to define a new paraproduct P intertwined with the paraproduct through

L P f g = P f Lg.
Since L is not invertible, we use L -1 an inverse up to a smooth error term. Hence a more conceivable intertwining relation is

L P f g = P f Lg -e -L (P f Lg) .
Definition. Given any distributions f, g ∈ D (M ), we define P f g as

P f g := L -1 P f Lg
for which we have the explicit formula

P f g = a∈A b ;a2< b 2 Q∈StGC a b Q 1 0 Q 1• t Q 2 t f • Q 3 t g dt t
where

Q 1 t := Q 1 t (tL) -1 and Q 3 t := Q 3 t (tL).
It is immediate that Q 3 belongs to StGC a3+2 . The cancellation property of Q 1 is given by the following lemma. Remark that it is not in standard form anymore, this is where the GC class comes into play. t . This immediatly follows from

Lemma . . Let Q ∈ StGC b 2 . Then Q t := Q t (tL) -1
Q t (tL) -1 = (t |I| 2 V I )(tL) j 2 (tL) -1 P (c) t = (t |I| 2 V I )(tL) j-2 2 P (c) t (Id -e L ).
This lemma immediatly yields the following proposition, that is P has the same structure as P hence the same continuity estimates.

Proposition . . For any distribution f, g ∈ D (M ), P f g is given as a linear combination of terms of the form

1 0 Q 1• t Q 2 t f • Q 3 t g dt t where Q 1 ∈ GC b 2 -2 , Q 2 ∈ StGC [0,b] and Q 3 ∈ StGC b 2 +2
. Thus for any regularity exponent α, β ∈ (-2b, 2b), we have the following continuity results.

If α ≥ 0, then (f, g) → P f g is continuous from C α × C β to C β . If α < 0, then (f, g) → P f g is continuous from C α × C β to C α+β .
We also have the same associated Sobolev estimates.

If α > 0, then (f, g) → P f g is continuous from C α × H β to H β and from H α × C β to H β . If α < 0, then (f, g) → P f g is continuous from C α × H β to H α+β and from H α × C β to H α+β .

-Correctors and commutators

The study of elliptic PDEs with singular product involves resonant term given a function u paracontrolled by a noise dependent function X ∈ C α , that is u = P u X + u with u ∈ C α and u ∈ C 2α a smoother remainder. If α < 1, the product uζ is singular for ζ ∈ C α-2 however we have the formal decomposition

Π(u, ζ) = Π P u X, ζ + Π(u , ζ) = u Π(X, ζ) + C(u , X, ζ) + Π(u , ζ)
with the corrector C introduced by Gubinelli, Imkeller and Perkowski in [START_REF] Gubinelli | Paracontrolled distributions and singular PDEs[END_REF] defined as

C(a 1 , a 2 , b) := Π P a1 a 2 , b -a 1 Π(a 2 , b). If 2 3 < α < 1, then the product Π(u , ζ) is well-defined.
Thus we are able to give a meaning to the product uζ for u paracontrolled by X once we have a proper continuity estimate for C and a meaning to the product Xζ; this is the controlled rough path philosophy. This last task is only a probabilistic one and does not impact the analytical resolution of the equation, this is the renormalisation step. We state here a continuity estimate for C while its proof is given in Proposition B.4 in Appendix B.

Proposition . . Let α 1 ∈ (0, 1) and α 2 , β ∈ R. If α 2 + β < 0 and α 1 + α 2 + β > 0, then (a 1 , a 2 , b) → C(a 1 , a 2 , b) extends in a unique continuous operator from C α1 × C α2 × C β to C α1+α2+β .
We also have the following proposition to work with Sobolev spaces.

Proposition . . Let α 1 ∈ (0, 1) and α 2 , β ∈ R. If α 2 + β < 0 and α 1 + α 2 + β > 0, then (a 1 , a 2 , b) → C(a 1 , a 2 , b) extends in a unique continuous operator from H α1 × C α2 × C β to H α1+α2+β .
Remark : Note that the first paramater α 1 has to be smaller than 1. This is due to the fact that for any function f ∈ C α with α ≥ 0, one has

|f (x) -f (y)| ≤ f C α d(x, y) α∧1
with a factor no greater than 1 even if α is. This means that we are not able to benefit from regularity greater than 1 with only a first order Taylor expansion. To work with a function of regularity α 1 ∈ (1, 2), one have to consider the refined corrector defined in the flat one dimensional case by

C (1) a 1 , a 2 , b (x) := Π P a1 a 2 , b (x) -a 1 (x)Π a 2 , b (x) -a 1 (x)Π P (x-•) a 2 , b (x)
that we interpret as a first order refined corrector for x ∈ T. There is an analog refined corrector on a manifold M , see [START_REF] Bailleul | High order paracontrolled calculus[END_REF]. However, this will not be needed in this work.

We need the corrector C to study ill-defined product, this is the condition α 2 + β < 0. However, we also have to investigate well-defined product to get more accurate descriptions. For this purpose, we introduce the commutator

D(a 1 , a 2 , b) := Π P a1 a 2 , b) -P a1 Π(a 2 , b). Proposition . . Let α 1 ∈ (0, 1) and α 2 , β ≥ 0. Then (a 1 , a 2 , b) → D(a 1 , a 2 , b) extends in a unique continuous operator from C α1 × C α2 × C β to C α1+α2+β and from H α1 × C α2 × C β to H α1+α2+β .
Again, one can bypass the condition α 1 ∈ (0, 1) using refined commutators. Note that in their initial work [START_REF] Gubinelli | Paracontrolled distributions and singular PDEs[END_REF], Gubinelli, Imkeller and Perkowski call C a commutator whereas with the point of view of high order paracontrolled calculus of [START_REF] Bailleul | High order paracontrolled calculus[END_REF], the operator D is closer to be a commutator than C. We need one final commutator that swaps paraproducts defined by

S(a 1 , a 2 , b) := P b P a1 a 2 -P a1 P b a 2 . Proposition . . Let α 1 , α 2 ∈ R and β < 0. Then (a 1 , a 2 , b) → S(a 1 , a 2 , b) extends in a unique continuous operator from C α1 × C α2 × C β to C α1+α2+β and from H α1 × C α2 × C β to H α1+α2+β .

-The Anderson Hamiltonian

In this section, we define and study the Anderson Hamiltonian

H := L + ξ
where -L is the Laplace-Beltrami operator on a compact two-dimensional manifold M without boundary or with a smooth boundary and Dirichlet conditions. To apply the construction of the first section, one needs to have an Hörmander representation for L. This is possible in this case with a number of vector fields possibly greater than the dimension, see for example Section 4.2.1 from Stroock's book [START_REF] Stroock | An Introduction to the Analysis of Paths on a Riemanian Manifold[END_REF]. The random potential ξ is a spatial white noise and belongs almost surely to C α-2 for any α < 1. For a generic function u ∈ L 2 , the product uξ is ill-defined hence one needs to find a proper domain for the operator. A natural method would be to take the closure of the subspace of smooth functions with the domain norm u L 2 + Hu L 2 . However this yields a trivial domain since Hu has the same regularity as the noise, because of the product uξ if u is smooth, thus it does not belong to L 2 . Following the recent study of singular SPDEs, one can construct a random domain D Ξ depending on an enhancement Ξ of the noise obtained through a renormalisation procedure. One can use the paraproduct to decompose the product for u ∈ H α as uξ = P u ξ + P ξ u + Π(u, ξ).

In this expression, the roughest term is P u ξ ∈ C α-2 while P ξ u + Π(u, ξ) formally belongs to H 2α-2 . For a function u in the domain, we want to cancel out the roughest part of the product using the Laplacian term Lu, hence we want

Lu = P u ξ + v with v ∈ H 2α-2
. This suggests the paracontrolled expansion

u = P u X + u with X := L -1 ξ
and u ∈ H 2α . We insist that we want functions in the domain to encode exactly what is needed to have a cancellation between the Laplacian and the product. In particular, H is not treated at all like a perturbation of the Laplacian.

At this point, two natural questions arise. Is the subspace of such paracontrolled functions dense in L 2 and can one make sense of the singular product? 1) For the first question, one can introduce a parameter s > 0, in the spirit of what Gubinelli, Ugurcan and Zachhuber did in [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF], and consider the modified paracontrolled expansion u = P s u X + u s with the truncated paraproduct P s defined below. For s = s(Ξ) small enough, the map Φ s (u) := u -P s u X is invertible as a perturbation of the identity and one can show that the subspace of such paracontrolled functions is indeed dense. The parameter s will also be a very useful tool to investigate the different properties of H. Indeed, the Anderson operator will be given as

Hu = Lu s + F Ξ,s (u)
with F Ξ,s : D(H) ⊂ L 2 → L 2 an explicit operator and as s goes to 0, u s gets closer to u while F Ξ,s diverges. These different representations of H will yield a family of bounds on the eigenvalues λ n (Ξ) n≥1 of H of the form

m -(Ξ, s)λ n -m(Ξ, s) ≤ λ n (Ξ) ≤ m + (Ξ, s)λ n + m(Ξ, s)
with (λ n ) n≥1 the eigenvalues of L. In partiular, m -(Ξ, s) and m + (Ξ, s) converge to 1 while m(Ξ, s) diverges almost surely as s goes to 0. A particular choice for s implies the simpler bounds

λ n -m 1 δ (Ξ) ≤ λ n (Ξ) ≤ (1 + δ)λ n + m 2 δ (Ξ)
for any δ ∈ (0, 1) which is expected but was not present in the initial work of Allez and Chouk [START_REF] Allez | The continuous Anderson Hamiltonian in dimension two[END_REF] or the work of Labbé [START_REF] Labbé | The continuous Anderson Hamiltonian in d ≤ 3[END_REF]. These bounds are sharp enough to imply the Weyl law

lim λ→∞ λ -1 {n ≥ 0; λ n (Ξ) ≤ λ} = Vol(M ) 4π
and give bounds for the tails of the eigenvalues.

2) For the second question, one introduces the corrector C with

Π(u, ξ) = uΠ(X, ξ) + C(u, X, ξ) + Π(u , ξ)
for u paracontrolled by X. One has to define the product Π(X, ξ) independently of the operator, this is the renormalisation step. To do so, we use the Wick product and set

Π(X, ξ) := lim ε→0 Π(X ε , ξ ε ) -E Π(X ε , ξ ε )
with ξ ε a regularisation of the noise. In some sense explained in Proposition 2.8, the operator H is the limit of the renormalised operators

H ε := L + ξ ε -c ε with c ε := E Π(X ε , ξ ε )
a smooth function diverging almost surely as ε goes to 0. Note that on the torus, the noise is invariant by translation and c ε is constant.

The approach sketched above yields an operator H : D(H) ⊂ L 2 → H 2α-2 with D(H) the space of paracontrolled functions. In two dimensions, 2α -2 < 0 hence one needs to refine the definition of the domain to get an unbounded operator in L 2 . To this purpose, Allez and Chouk introduced in [START_REF] Allez | The continuous Anderson Hamiltonian in dimension two[END_REF] the subspace of D(H) of strongly paracontrolled functions still dense in L 2 . This was also used by Gubinelli, Ugurcan and Zachhuber in [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF] and adapted to the dimension 3 using a Hopf-Cole type transformation. We present here a different approach based on a higher order expansion. In particular, the domain of H will consist of functions u such that u = P u X 1 + P u X 2 + u where X 1 ∈ C α , X 2 ∈ C 2α are noise-dependent functions and u ∈ H 2 . Note that since we want to get bounds in Ξ, quantitative estimates are needed and we keep track of the different explicit constants that appear, in particular how small s needs to be with respect to the noise. If one is only interest in qualitative results, details of almost all computations can be skipped.

We shall first construct in Section 2.1 the enhanced noise Ξ from ξ by a renormalisation procedure and prove exponential moments for its norm. The domain D Ξ of H is constructed in Section 2.2 and proved to be dense using a truncated paraproduct P s . We show in particular in Proposition 2.6 that the natural norms of D Ξ are equivalent to the norm operator; this will give the upper bound for the eigenvalues. Section 2.2 is ended with the computation of the Hölder regularity of the elements of the domain. After showing that the operator is closed, we show in Section 2.3 that H is the limit of the operators H ε in some sense which yields the symmetry of H. We then control in Proposition 2.9 the H 1 norm of u from the associated bilinear form applied to u; this will give the lower bound for the eigenvalues. This gives self-adjointness and pure point spectrum using the Babuška-Lax-Milgram theorem and we conclude the section with a bound on the convergence of the eigenvalues of H ε to H. Section 2.4 treats the Schrödinger equation.

As in the work of Allez and Chouk [START_REF] Allez | The continuous Anderson Hamiltonian in dimension two[END_REF], Labbé [START_REF] Labbé | The continuous Anderson Hamiltonian in d ≤ 3[END_REF] and of Gubinelli, Ugurcan and Zachhuber [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF], we construct a dense random subspace of L 2 though a renormalisation step to get a selfadjoint operator with pure point spectrum. Our approach is different since we perform a second order expansion using paracontrolled calculus based on the heat semigroup on the manifold M .

We refine the upper bounds on the eigenvalues obtained in [START_REF] Allez | The continuous Anderson Hamiltonian in dimension two[END_REF] on the torus while also providing lower bounds. We get upper bounds for P(λ n (Ξ) ≤ λ) for λ to +∞ and -∞. For λ to -∞, a bound was first given in [START_REF] Labbé | The continuous Anderson Hamiltonian in d ≤ 3[END_REF] for a bounded domain with different boundary conditions. We have a more explicit dependence on n while a less precise bound with respect to λ. To the best of our knowledge, no bounds for λ to +∞ were known. We also prove that the eigenfunctions of H belong to C 1-while the works [START_REF] Allez | The continuous Anderson Hamiltonian in dimension two[END_REF][START_REF] Labbé | The continuous Anderson Hamiltonian in d ≤ 3[END_REF][START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF] only gave Sobolev regularity. For the Schrödinger equation, we get on a manifold the same result as Gubinelli, Ugurcan and Zachhuber get on the torus, see [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF]. As in their work, our construction of the Hamilton Anderson on M could be used to study other evolution PDEs. All these results are new in our geometrical framework.

-Renormalisation

As explained in the introduction, an element of the domain of H should behave like the linear part X := -L -1 ξ hence the product uξ does not make sense in two dimensions. Using the corrector, we are able to define the product uξ for u paracontrolled by X once the product Xξ is defined. To do so, a naive approach would be to regularize the noise where ξ ε = Φ(εL)ξ is a regularisation of the noise and take ε to 0. The only condition we take is Φ such that (Φ(εL)) ε belongs to the class G, for example Φ(εL) = e εL works. Since the product is ill-defined, the quantity Π(X ε , ξ ε ) diverges as ε goes to 0 with X ε := -L -1 ξ ε . The now usual way is to substract another diverging quantity c ε such that the limit Π(X, ξ) := lim ε→0 Π(X ε , ξ ε ) -c ε exists and take this as the definition of the product. This is the Wick renormalisation and the purpose of the following theorem with the renormalised Anderson Hamiltonian

H ε := L + ξ ε -c ε .
We insist here that c ε is not a priori a constant in the general geometrical setting of a twodimensional manifold.

Theorem . . Let α < 1 and

c ε := E Π(X ε , ξ ε ) .
Then there exists a random distribution Π(X, ξ) that belongs almost surely to C 2α-2 and such that

lim ε→0 E Π(X, ξ) -(Π(X ε , ξ ε ) -c ε ) p C 2α-2 = 0
for any p ≥ 1.

Proof : Since the noise is Gaussian, we only need to control second order moment using hypercontractivity. The resonant term Π(X ε , ξ ε ) is a linear combination of terms of the form

I ε := 1 0 P • t Q 1 t X ε • Q 2 t ξ ε dt t with P ∈ StGC [0,b] and Q 1 , Q 2 ∈ StGC b 2 .
We also define the renormalised quantity

J ε := I ε -E[I ε ]. Let u ∈ (0, 1), x ∈ M and Q ∈ StGC r with r > |2α -2|. The expectation E |Q u I ε (x)| 2 is given by the integral over M 2 × [0, 1] 2 of K QuP • t (x, y)K QuP • s (x, z)E Q 1 t X ε (y)Q 2 t ξ ε (y)Q 1 s X ε (z)Q 2 s ξ ε (z)
against the measure µ(dy)µ(dz)(ts) -1 dtds. Using the Wick formula, we have

E Q 1 t X ε (y)Q 2 t ξ ε (y)Q 1 s X ε (z)Q 2 s ξ ε (z) = E Q 1 t X ε (y)Q 2 t ξ ε (y) E Q 1 s X ε (z)Q 2 s ξ ε (z) + E Q 1 t X ε (y)Q 1 s X ε (z) E Q 2 t ξ ε (y)Q 2 s ξ ε (z) + E Q 1 t X ε (y)Q 2 s ξ ε (z) E Q 1 s X ε (z)Q 2 t ξ ε (y) = (1) + (2) + (3)
and this yields

E |Q u I ε (x)| 2 = I (1) ε (x) + I (2) ε (x) + I (3) ε (x).
The first term corresponds exactly to the extracted diverging quantity since

I (1) ε = E 1 0 Q u P • t Q 1 t X ε • Q 2 t ξ ε dt t 2 = E Q u (I ε ) 2
and we have

E |Q u J ε (x)| 2 = E Q u I ε (x) -E[Q u I ε ](x) 2 = I (2) ε (x) + I (3) ε (x).
Using that (Ψ(εL)) ε belongs to G, ξ is an isometry from L 2 to square-integrable random variables and lemma 1.6, we have

I (2) ε (x) + I (3) ε (x) M 2 [0,1] 2 K QuP • t (x, y)K QuP • s (x, z) G 2ε+t+s (y, •), G 2ε+t+s (z, •) 2 µ(dy)µ(dz)tsdtds M 2 [0,1] 2 K QuP • t (x, y)K QuP • s (x, z)G 2ε+t+s (y, z) 2 µ(dy)µ(dz)tsdtds M 2 [0,1] 2 G u+t (x, y)G u+s (x, z)G 2ε+t+s (y, z) 2 µ(dy)µ(dz)tsdtds M 2 [0,1] 2 (2ε + t + s) -d 2 G u+t (x, y)G u+s (x, z)G 2ε+t+s (y, z)µ(dy)µ(dz)tsdtds [0,1] 2 (2ε + t + s) -d 2 (ε + u + t + s) -d 2 tsdtds (ε + u) 2-d hence the family Π(X ε , ξ ε ) -c ε ε>0 is bounded in C 2α-2 for any α < 1 since d = 2.
These computations also show that the associated linear combination of

J := 1 0 P • t Q 1 t X • Q 2 t ξ -E P • t Q 1 t X • Q 2 t ξ
dt t yields a well-defined random distribution of C 2α-2 for α < 1 that we denote Π(X, ξ). The same type of computations show the convergence and completes the proof.

The enhanced noise is defined as

Ξ := ξ, Π(X, ξ) ∈ X α
where X α := C α-2 × C 2α-2 . One has to keep in mind that the notation Π(X, ξ) is only suggestive.

In particular for almost every ω, one has

Π X, ξ (ω) = Π X(ω), ξ(ω)
since the product is almost surely ill-defined. We also denote the regularized enhanced noise

Ξ ε := ξ ε , Π(X ε , ξ ε ) -c ε with the norm Ξ -Ξ ε X α := ξ -ξ ε C α-2 + Π(X, ξ) -Π(X ε , ξ ε ) + c ε C 2α-2
which goes to 0 as ε goes to 0. Using that the noise is Gaussian and almost surely in C -1-κ for all κ > 0, we have exponential moment for the norm of the enhanced noise.

Proposition . . There exists h > 0 such that

E e h ξ 2 C α-2 +h Π(X,ξ) C 2α-2 < ∞.
Proof : Let t ∈ (0, 1) and Q ∈ StGC r with r > |α -2|. Using the Gaussian hypercontractivity,

we have

E Q t ξ p L p x = M E [|Q t ξ| p (x)] µ(dx) ≤ (p -1) p 2 M E |Q t ξ| 2 (x) p 2 µ(dx)
hence we only need to bound the second moment, which is bounded by

E |Q t ξ| 2 (x) = K Qt (x, •) 2 L 2 1 µ B(x, √ t) . Using that B α-2+ 1 p 2p,2p → B α-2 ∞,∞ , we have E e h ξ 2 C α-2 = p≥0 h p p! E ξ 2p C α-2 ≤ p0 p=0 h p p! E ξ 2p C α-2 + p>p0 h p p! E ξ 2p B α-2+ 1 p 2p,2p p0 p=0 h p p! E ξ 2p0 C α-2 p p 0 + p>p0 h p (2p -1) p p! Vol(M )
for p 0 > 2 1-α hence the result for h small enough. For the bound on Π(X, ξ), the computations are the same without the square since it belongs to the second Wiener chaos hence Gaussian hypercontractivity gives

E |Q t Π(X, ξ)| p (x) ≤ (p -1) p E |Q t Π(X, ξ)| 2 (x) p 2 .

-Domain of the Hamiltonian

We first motivate the definition of the domain. Let α ∈ ( 2 3 , 1) such that ξ belongs almost surely to C α-2 . Let X ∈ C α be a noise-dependent function and consider u = P u X + u a function paracontrolled by X with u ∈ H α and u ∈ H 2α . Then

Hu = Lu + ξu = L P u X + u + P u ξ + P ξ u + Π P u X + u , ξ = P u LX + P u ξ + Lu + P ξ u + u Π(X, ξ) + C(u , X, ξ) + Π(u , ξ) .
Taking u = u and -LX = ξ, the first two terms cancel each other and we get

Hu = Lu + P ξ u + uΠ(X, ξ) + C(u, X, ξ) + Π(u , ξ) ∈ H 2α-2 .
This yields an unbounded operator in L 2 with values in H 2α-2 . Since 2α -2 < 0, Hu does not belong to L 2 hence we do not have an operator from L 2 to itself and this makes harder to study the spectral properties of H. To get around this, Allez and Chouk introduced in [START_REF] Allez | The continuous Anderson Hamiltonian in dimension two[END_REF] the subspace of functions u paracontrolled by L -1 ξ such that Hu does belong to L 2 called strongly paracontrolled functions. This approach was also used by Gubinelli, Ugurcan and Zachhuber in [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF] however we proceed differently and use higher order expansions. Let X 1 := X and X 2 ∈ C 2α be another noisedependent function. Given u 2 ∈ H α and u ∈ H 3α , we consider u = P u X 1 + P u2 X 2 + u and we have

Hu = P u2 LX 2 + uΠ(X 1 , ξ) + C(u, X 1 , ξ) + P u2 Π(X 2 , ξ) + D(u 2 , X 2 , ξ) + P u P ξ X 1 + S(u, X 1 , ξ) + P ξ P u2 X 2 + P ξ u + Lu + Π(u , ξ).
Taking u 2 = u and -LX 2 = Π(X 1 , ξ) + P ξ X 1 cancels the terms of Sobolev regularity 2α -2 and we get

Hu = Π u, Π(X 1 , ξ) + P Π(X1,ξ) u + C(u, X 1 , ξ) + P u Π(X 2 , ξ) + D(u, X 2 , ξ) + S(u, X 1 , ξ) + P ξ P u X 2 + P ξ u + Lu + Π(u , ξ)
hence Hu ∈ H 3α-2 ⊂ L 2 . This motivates the following definition for the domain D Ξ of H with -LX 1 := ξ and -LX 2 := Π(X 1 , ξ) + P ξ X 1 .

Definition. We define the set D Ξ of functions paracontrolled by Ξ as

D Ξ := u ∈ L 2 ; u := u -P u X 1 -P u X 2 ∈ H 2 .
The domain D Ξ is the random subspace of functions u ∈ L 2 paracontrolled by X 1 and X 2 up to a remainder u ∈ H 2 given by the explicit formula

u = Φ(u) := u -P u X 1 -P u X 2 .
With this notation, we have D Ξ = Φ -1 (H 2 ) and since X 1 + X 2 ∈ C α , we actually have D Ξ ⊂ H β for every β < α. However, we have no idea at this point if this domain is trivial or dense in L 2 and an inverse to Φ would be useful. However, it is not necessarily invertible so we introduce a parameter s > 0 and consider

Φ s : D Ξ → H 2 u → u -P s u X 1 -P s u X 2
where P s is defined as

P s f g := a∈A b ;a2< b 2 Q∈StGC a b Q s 0 Q 1• t Q 2 t f • Q 3 t g dt t .
The important property is that while still encoding the important information of the paraproduct P, the truncated paraproduct P s is small as an operator for s small; this is quantified as follows and proved in Proposition B.2 in Appendix B.

Proposition . . Let γ ∈ (0, 1) be a regularity exponent and X ∈ C γ . For any β ∈ [0, γ), we have

u → P s u X L 2 →H β s γ-β 4 γ -β X C γ
Since X 1 and X 2 depends continuously on Ξ, this implies the existence of m > 0 such that

P s u X 1 + P s u X 2 H β ≤ m s α-β 4 α -β Ξ X α (1 + Ξ X α ) u L 2 thus the operator u → P s u (X 1 + X 2 ) is continuous from L 2 to H β for β ∈ [0, α)
and arbitrary small as s goes to 0. Hence we get that Φ s : H β → H β is invertible for s = s(Ξ, β) small enough as a perturbation of the identity. Since P u X i -P s u X i is a smooth function for any s > 0, the domain is still given by

D Ξ = Φ -1 (H 2 ) = (Φ s ) -1 (H 2 )
and we have a decomposition given by Φ s for any u ∈ D Ξ , that is

u = P s u X 1 + P s u X 2 + Φ s (u).
In particular, we emphasize that the domain does not depend on s while the decomposition we consider for element of the domain might. We denote

x := Ξ X α
to keep track of the quantitative dependance with respect to the enhanced noise Ξ and lighten the notation. We use the letter x as a reminder of the noise dependance. For any 0 ≤ β < α, we define

s β (Ξ) := α -β mx(1 + x)
such that for s < s β (Ξ), the operator Φ s : H β → H β is invertible and we denote Γ its inverse. We choose to drop the parameter s in the notation to lighten the computations however the reader should keep in mind that the map Γ depends on s. It is implicitly characterized by the relation

Γu = P s Γu X 1 + P s Γu X 2 + u
for any u ∈ H β . Our choice of P s is motivated by the preservation of the intertwining relation

P s = L -1 • P s • L
with P s defined as P s . The map Γ will be a crucial tool to study the domain D Ξ , in particular to show density in L 2 . Continuity estimates for Φ s and Γ are given in the next proposition. Note that in the following, this bound of the form a -b ≤ c will be used as a ≤ b + c or b ≤ a + c.

Proposition . . Let β ∈ [0, α) and s ∈ (0, 1). We have

Φ s (u) -u H β ≤ m α -β s α-β 4 x(1 + x) u L 2 .
If moreover s < s β (Ξ), this implies

Γu H β ≤ 1 1 -m α-β s α-β 4 x(1 + x) u H β .

Proof :

The bounds on Φ s follows directly from proposotion 2.3. Moreover since

m α -β s α-β 4 x(1 + x) < 1 for s < s β (Ξ), the map Φ s : H β → H β is invertible and we have Γu H β ≤ 1 1 -m α-β s α-β 4 x(1 + x) u H β .
Let us insist that u s H β is always controlled by u H β while s need to be small depending for u H β to be controlled by u s H β . We also define the map Γ ε associated to the regularized noise Ξ ε as

Γ ε u = P s Γεu X (ε) 1 + P s Γεu X (ε) 2 + u with -LX (ε) 1 := ξ ε and -LX (ε) 2 := Π(X (ε) 1 , ξ ε ) -c ε + P ξε X (ε)
1 . It satisfies the same bound as Γ with Ξ ε X α and the following approximation lemma holds. We do not need to explicit the constant, it depends polynomialy on the noise Ξ and diverges as s goes to s β (Ξ).

Lemma . . For any 0 ≤ β < α and 0 < s < s β (Ξ), we have

Id -ΓΓ -1 ε L 2 →H β Ξ,s,β Ξ -Ξ ε X α .
In particular, this implies the norm convergence of Γ ε to Γ with the bound

Γ -Γ ε H β →H β Ξ,s,β Ξ -Ξ ε X α .
and the natural norms of the domain

u 2 DΞ := u 2 L 2 + Φ s (u) 2 H 2
with the following proposition. For s ∈ (0, 1) and δ > 0, we introduce the constant

m 2 δ (Ξ, s) := k s α-2 2 x(1 + x 2 ) + s α-β 4 x 2 (1 + x 3 ) + δ -3 1 + s α 4 x(1 + x) x 4 (1 + x 8 )
where the "2" refers to H 2 and for a constant k > 0 large enough depending only on M and L. In particular, it depends polynomialy on the enhanced noise and diverges as s or δ goes to 0.

Proposition . . Let u ∈ D Ξ and s > 0. For any δ > 0, we have

(1 -δ) u s H 2 ≤ Hu L 2 + m 2 δ (Ξ, s) u L 2 and Hu L 2 ≤ (1 + δ) u s H 2 + m 2 δ (Ξ, s) u L 2 with u s = Φ s (u).
Proof : For any s > 0, we have

Hu = Lu s + P ξ u s + Π(u s , ξ) + R(u) + Ψ s (u).
Then Lu s ∈ L 2 and for β = 1 2 ( 2 3 + α), we have

R(u) L 2 x(1 + x 2 ) u H β Ψ s (u) L 2 s α-2 2 x(1 + x 2 ) u L 2 P ξ u s + Π(u s , ξ) L 2 ξ C α-2 u s H 4 3 .
One can bound the H β norm of u using Proposition 2.4 with

u H β ≤ u s H β + m α -β s α-β 4 x(1 + x) u L 2
and since β < 1, one has

Lu s -Hu L 2 s α-2 2 x(1 + x 2 ) + s α-β 4 x 2 (1 + x 3 ) u L 2 + x(1 + x 2 ) u s H 4 3 .
Since 0 < β < 2, we have for any t > 0

u s H 4 3 t 0 (t L)e -t L u s dt t H 4 3 + e -tL u s H 4 3 t 2 3 u s H 2 + t -4 2 1 + s α 4 x(1 + x) u L 2 . Take t = δ kx(1 + x 2 ) 3 2
with k the constant from the previous inequality and δ > 0. This yields

Lu s -Hu L 2 m 2 δ (Ξ, s) u L 2 + δ u s H 2 .
and completes the proof.

Finally, we can compute the Hölder regularity of the domain. In particular, this will implies the α-Hölder regularity of the eigenfunctions of H.

Proposition. We have

D Ξ ⊂ C 1-κ
for any κ > 0.

Proof :

We have

H ε u ε = Lu s + P ξε u s + Π(u s , ξ) + R ε (u ε ) + Ψ s ε (u ε ) where R ε and Ψ s ε are defined as R and Ψ s with Ξ ε instead of Ξ. For β = 1 2 ( 2 3 + α), we have R(u)-R ε (u ε ) L 2 ≤ R(u -u ε ) L 2 + (R -R ε )(u ε ) L 2 x(1 + x 2 ) u -u ε H β + (1 + x) Ξ -Ξ ε X α u ε H β x(1 + x 2 ) Γ -Γ ε H β →H β + (1 + x) Γ ε H β →H β Ξ -Ξ ε X α u s H 2
and the same reasoning gives

Ψ s (u) -Ψ s ε (u) L 2 s,Ξ Ξ -Ξ ε X α .
Thus one completes the proof with the bound Γ -Γ ε H β →H β from Lemma 2.5.

The symmetry of H immediately follows.

Corollary. The operator H is symmetric.

Proof : Let u, v ∈ D Ξ and consider u := Φ s (u) and v := Φ s (v) for s < s 0 (Ξ). Since H ε is a symmetric operator, we have

Hu, v = lim ε→0 H ε Γ ε u , Γ ε v = lim ε→0 Γ ε u , H ε Γ ε v = u, Hv using that H ε Γ ε converges to HΓ and Γ ε to Γ in norm convergence.
The next proposition states that the quadratic form associated to H is bounded from below by the H 1 norm of u . This weak coercivity property will give below self-adjointness with the Babuška-Lax-Milgram theorem. This was already used in the work [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF] of Gubinelli, Ugurcan and Zachhuuber, where the proof of self-adjointness relies on the reasoning of almost duality encoded in the operator A. For s ∈ (0, 1) and δ > 0, introduce the constant

m 1 δ (Ξ, s) := k x(1 + x 2 ) + s α-β 4 x 2 (1 + x 3 ) + s α-2 2 x(1 + x 2 ) + s α-4 2 x + δ -β 1-β x(1 + x 2 ) + s α-β 4 x 2 (1 + x) β 1-β 1 + s α 4 x(1 + x)
where β = 1 2 ( 2 3 + α) and for a constant k > 0 large enough depending only on M and L while the "1" refers to H 1 . In particular, it depends polynomialy on the enhanced noise and diverges as s or δ goes to 0. Proposition . . Let u ∈ D Ξ and s > 0. For any δ > 0, we have

(1 -δ) ∇u s , ∇u s ≤ u, Hu + m 1 δ (Ξ, s) u 2 L 2 and (1 -δ) ∇u s , ∇u s ≤ u, H ε u + m 1 δ (Ξ, s) u 2 L 2
where u s = Φ s (u).

Proof : For u ∈ D Ξ , we have

Hu = Lu s + P ξ u s + Π(u s , ξ) + R(u) + Ψ s (u)
with u s = Φ s (u) ∈ H where A(u, ξ, u ) = P u ξ, u -u, Π(u , ξ) . For

β := 1 2 ( 2 3 + α), we have u, R(u) u L 2 R(u) L 2 x(1 + x 2 ) u L 2 u H β , u, P ξ u s u H β P ξ u s C 2β-2 x u H β u s H β , P u LX 2 , u s P u LX 2 H 2β-2 u s H β x 2 u L 2 u s H β .
Using Proposition 1.5, we have

A(u, ξ, u s ) ξ C α-2 u H β u s H β x u H β u s H β .
Finally, we have

u, Ψ s (u) u L 2 Ψ s (u) L 2 s α-2 2 x(1 + x 2 ) u 2 L 2 (P u -P s u )ξ, u s (P u -P s u )ξ L 2 u s L 2 s α-4 2 x u L 2 u s L 2 with Proposition B.3 in Appendix B. Since u ∈ D Ξ , we have u H β ≤ u s H β + m α -β s α-β 4 x(1 + x) u L 2
hence there exists k > 0 such that

∇u s , ∇u s ≤ u, Hu + k x(1 + x 2 ) + s α-β 4 x 2 (1 + x 3 ) + s α-2 2 x(1 + x 2 ) + s α-4 2 x u 2 L 2 + k x(1 + x 2 ) + s α-β 4 x 2 (1 + x) u s H β .
Since 0 < β < 1, we have for any t > 0

u s 2 H β t 0 (t L)e -t L u s dt t 2 H β + e -tL u s 2 H β t 1-β u s 2 H 1 + t -β 1 + s α 4 x(1 + x) 2 u 2 L 2 .
Given any δ > 0, we set

t =   δ k x(1 + x 2 ) + s α-β 4 x 2 (1 + x)   1 1-β
where k > 0 the constant from the previous inequality and this yields

(1 -δ) ∇u s , ∇u s ≤ u, Hu + m 1 δ (Ξ, s) u L 2 . The same computations show (1 -δ) ∇u s , ∇u s ≤ u, H ε u + m 1 δ (Ξ ε , s) u 2 L 2 .
Since Ξ ε -Ξ α goes to 0 as ε goes to 0, the result holds uniformly in ε with m 1 δ (Ξ, s).

This implies that H is almost surely bounded below by the random variable -m 1 δ (Ξ, s) for any δ > 0 and s > 0. Using the Babuška-Lax-Milgram theorem, one gets an invertible operator via the solution of

(H + k Ξ )u = v for k Ξ > m 1 δ (Ξ, s) and v ∈ L 2 .
Proposition . . Let δ ∈ (0, 1) and s > 0. Then for any constant k Ξ > m 1 δ (Ξ, s), the operators H + k Ξ and H ε + k Ξ are invertible. Moreover the operators

H + k Ξ -1 : L 2 → D Ξ H ε + k Ξ -1 : L 2 → H 2
are bounded.

Proof :

We want to use the theorem of Babuška-Lax-Milgram, see [START_REF] Babuška | Error-bounds for finite element method[END_REF]. This is a generalization of the Lax-Milgram theorem with a weaker condition of coercivity. Since k Ξ > m 1 δ (Ξ, s), Proposition 2.9 gives

k Ξ -m 1 δ (Ξ, s) u 2 L 2 < (H + k Ξ )u, u for u ∈ D Ξ . Considering the norm u 2 DΞ = u 2 L 2 + u s 2 H 2
on D Ξ , this yields a weakly coercive operator using Proposition 2.6 in the sense that

u DΞ Ξ (H + k Ξ )u L 2 = sup v L 2 =1 (H + k Ξ )u, v
for any u ∈ D Ξ . Moreover, the bilinear map

B : D Ξ × L 2 → R (u, v) → (H + k Ξ )u, v
is continuous since Proposition 2.6 implies

|B(u, v)| ≤ (H + k Ξ )u L 2 v L 2 Ξ u DΞ v L 2
for u ∈ D Ξ and v ∈ L 2 . The last condition we need is that for any v ∈ L 2 \{0}, we have

sup u D Ξ =1 |B(u, v)| > 0.
Let assume that there exists v ∈ L 2 such that B(u, v) = 0 for all u ∈ D Ξ . Then

∀u ∈ D Ξ , u, v DΞ,D * Ξ = 0. hence v = 0 as an element of D * Ξ . By density of D Ξ in L 2 , this implies v = 0 in L 2
hence the property we want. By the theorem of Babuška-Lax-Milgram, for any f ∈ L 2 there exists a unique u ∈ D Ξ such that ∀v ∈ L 2 , B(u, v) = f, v .

Moreover, we have u DΞ Ξ f L 2 hence the result for (H + k Ξ ) -1 . The same argument works for H ε + k Ξ since proposition 2.9 also holds for H ε with bounds uniform in ε.

Using that a closed symmetric operator on a Hilbert space is self-adjoint if it has at least one real value in its resolvent set, this immediatly implies that H and H ε are self-adjoint, see [START_REF] Reed | Methods of modern mathematical physics. II. Fourier Analysis, self-adjointness[END_REF]. Moreover, the resolvant is a compact operator from L 2 to itself since D Ξ ⊂ H β for any β ∈ [0, α) hence the following result.

Corollary . . The operators H and H ε are self-adjoint with discret spectrum λ n (Ξ) n≥1 and λ n (Ξ ε ) n≥1 which are nondecreasing diverging sequences without accumulation points. Moreover, we have

L 2 = n≥1 Ker H -λ n (Ξ)
with each kernel being of finite dimension. We finally have the min-max principle

λ n (Ξ) = inf D sup u∈D; u L 2 =1
Hu, u

where D is any n-dimensional subspace of D Ξ that can also be given as

λ n (Ξ) = sup v1,...,vn-1∈L 2 inf u∈Vect(v1,...,vn-1) ⊥ u L 2 =1
Hu, u .

A natural question now is to estimate the size of the eigenvalues of H and try to get back geometric informations on the manifold M as one can do from the Laplacian. Let λ be an eigenvalue of H and u ∈ D Ξ such that Hu = λu.

Then there exists u ∈ H 2 such that u = Γu thus HΓu = λΓu .

This yields HΓu = λu + λ Γ -Id u hence one can relate the spectrum of H to the one of HΓ and the parameter s measures the error since Γ -Id u = P s Γu X 1 + P s Γu X 2 . And since HΓ is a perturbation of L, one can relate the spectrum of HΓ to the spectrum of L, as stated in the following proposition using the min-max result. We denote by (λ n ) n≥1 the nondecreasing positive sequence of the eigenvalues of L, since it corresponds to the case Ξ = 0. For s ∈ (0, 1) and δ > 0, introduce the constant

m + δ (Ξ, s) := (1 + δ) 1 + m α s α 4 x(1 + x) .
If s < s 0 (Ξ), we also introduce

m - δ (Ξ, s) := (1 -δ) 1 1 -m α s α 4 x(1 + x)
.

In particular, the constants depend polynomialy on the enhanced noise Ξ and converge to 1 as δ and s goes to 0. Moreover, m - δ (Ξ, s) diverges as s goes to s 0 (Ξ). Write a, b ≤ c to mean that we have both a ≤ c and b ≤ c.

Proposition . . Let s ∈ (0, 1) and δ > 0. Given any n ∈ Z + , we have

λ n (Ξ), λ n (Ξ ε ) ≤ m + δ (Ξ, s)λ n + 1 + m α s α 4 x(1 + x) + m 2 δ (Ξ, s).
If moreover s < s 0 (Ξ), we have

λ n (Ξ), λ n (Ξ ε ) ≥ m - δ (Ξ, s)λ n -m 1 δ (Ξ, s).
Proof : Let u 1 , . . . , u n ∈ H 2 be an orthonormal family of eigenfunctions of L associated to λ 1 , . . . , λ n and consider

u i := Γu i ∈ D Ξ for 1 ≤ i ≤ n. Since Γ is invertible, the family (u 1 , . . . , u n ) is free thus the min-max representation of λ n (Ξ) yields λ n (Ξ) ≤ sup u∈Vect(u1,...,un) u L 2 =1
Hu, u .

Given any normalised u ∈ Vect(u 1 , . . . , u n ), we have

Hu, u ≤ Hu L 2 ≤ (1 + δ) u s H 2 + m 2 δ (Ξ, s)
for u s = Φ s (u) using Proposition 2.6. Moreover

u s H 2 ≤ (1 + λ n ) u s L 2 ≤ (1 + λ n ) 1 + m α s α 4 x(1 + x)
hence the upper bound

λ n (Ξ) ≤ m + δ (Ξ, s)λ n + 1 + m α s α 4 x(1 + x) + m 2 δ (Ξ, s).
For the lower bound, we use the min-max representation of λ n (Ξ) under the form

λ n (Ξ) = sup v1,...,vn-1∈L 2 inf u∈Vect(v1,...,vn-1) ⊥ u L 2 =1
Hu, u .

Introducing F := Vect(u m ; m ≥ n),
we have that F ⊥ is a subspace of L 2 of finite dimension n -1 thus there exists a orthogonal family (v 1 , . . . , v n-1 ) such that F ⊥ = Vect(v 1 , . . . , v n-1 ). Since F is a closed subspace of L 2 as an intersection of hyperplans, we have

F = Vect(v 1 , . . . , v n-1 ) ⊥ hence λ n (Ξ) ≥ inf u∈F u L 2 =1
Hu, u .

Let u ∈ F with u L 2 = 1. Using Proposition 2.9, we have

Hu, u ≥ (1 -δ) ∇u s , ∇u s -m 1 δ (Ξ, s) ≥ (1 -δ) u s , Lu s -m 1 δ (Ξ, s) ≥ (1 -δ)λ n u s 2 L 2 -m 1 δ (Ξ, s).
Finally using Proposition 2.4 for s < s 0 (Ξ), we get

Hu, u ≥ 1 -δ 1 -m α s α 4 x(1 + x) λ n -m 1 δ (Ξ, s)
and the proof is complete.

There is a wide range of choices for the constants s ∈ (0, 1) and δ > 0. For example, one can take

s = mx(1 + x) 4 α
for any δ ∈ (0, 1) and get

λ n -m 1 δ (Ξ) ≤ λ n (Ξ) ≤ (1 + δ)λ n + m 2 δ (Ξ)
for explicit constants m 1 δ and m 2 δ , where the lower bound holds since δ < 1 gives s < s 0 (Ξ). This implies the following estimate for the tail of all the eigenvalues. A more precise result of this type was already obtained in [START_REF] Labbé | The continuous Anderson Hamiltonian in d ≤ 3[END_REF] by Labbé in the flat case for λ to -∞ with a = 1 where he also obtained a lower bound on the convergence of the form e -anλ ≤ P(λ n (Ξ) ≤ -λ) ≤ e -bnλ for λ > 0 large enough and a n > b n > 0 two constants. Here we get upper bounds for λ to +∞ and -∞.

Corollary . . For any n ∈ Z + and λ ∈ R, we have

1 -me -h(λ-2λn) 1 12 ≤ P λ n (Ξ) ≤ λ ≤ me -h(λn-λ) 1 5
where m = E e h Ξ X α .

Proof : Fix δ ∈ (0, 1) and let λ ∈ R. Denote m 1 = m 1 δ and m 2 = m 2 δ . We have

P λ n (Ξ) ≤ λ ≤ P λ n -m 1 ≤ λ and P λ n (Ξ) > λ ≤ P (1 + δ)λ n + m 2 > λ thus P m 2 ≤ λ -(1 + δ)λ n ≤ P λ n (Ξ) ≤ λ ≤ P m 1 ≥ -λ + λ n . There exists two constants a 1 , a 2 > 0 such that m i ≤ 1 + Ξ ai X α
for i ∈ {1, 2}, take for example a 1 = 5 and a 2 = 12. Hence P m i ≥ y = P Ξ X α ≥ (y -1) Proof : As stated, Proposition 2.9 yields

1 a i = P e h Ξ X α ≥ e hy 1 a i ≤ e -hy 1 a i E e h Ξ X
u s H 1 Ξ,s u DΞ( √ H + ) .
In fact, the inequality that is proved is

H + u, u -∇u s , ∇u s ≤ k Ξ u L 2 + δ u s H 1
thus one also get the other estimate

u DΞ( √ H + ) Ξ,s u s H 1 .
This yields a version of Brezis-Gallouët inequality for the Anderson Hamiltonian. In some sense, it interpolates the L ∞ -norm between the energy norm and the logarithm of the domain norm. This was already obtained in [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF] by Gubinelli, Ugurcan and Zachhuber on the torus.

Theorem . . For any

v ∈ D Ξ (H + ), we have v L ∞ Ξ v DΞ( √ H + )   1 + log 1 + v DΞ(H + ) v D( √ H + )   .
For any v ∈ H 2 , we have

v L ∞ Ξ H + ε v L 2   1 + log 1 + H + ε v L 2 H + ε v L 2   .
In particular, the second inequality holds uniformly in ε.

Proof : For any t > 0, we have

v L ∞ ≤ t 0 (t L)e -t L v dt t L ∞ + e -tL v L ∞ .
One can bound the integral as

t 0 (t L)e -t L v dt t L ∞ t 0 Lv L 2 dt t v H 2
which converges to u 0 in L 2 . We also have the convergences

u ε (t) → u(t) in L 2 , H + ε u ε (t) → H + u(t) in L 2 , ∂ t u ε (t) → ∂ t u(t) in L 2 for all t ∈ [0, T ].
Remark : From the solution to

i∂ t u = H + u -|u| 2 u,
on the torus, one easily gets the solution to the initial equation

i∂ t v = Hv -|v| 2 v
via the change of variable u(t, •) = e tkΞ v(t, •) since k Ξ is a constant. One could want to do the same in a manifold setting and compare the initial regularized equation

i∂ t u = Lu + ξ ε u -|u| 2 u.
with the renormalised equation

i∂ t v = Lv + ξ ε v -c ε v -|v| 2 v
as Tzvetkov and Visciglia's Theorem 1.1 from [START_REF] Tzvetkov | Two dimensional nonlinear Schrödinger equation with spatial white noise potential and fourth order nonlinearity[END_REF]. It is not clear what the change of variable should be on a manifold since c ε is a function and not a constant. It should still be possible to find an appropriate change of variable even though this requires some work.

A -Approximation operators

We describe in this Appendix technical estimates needed in our continuous setting analog of the discrete Paley-Littlewood decomposition. The following proposition is the analog of the inclusions of p spaces.

Proposition A. . Let p, q 1 , q 2 ∈ [1, ∞] with q 1 ≤ q 2 . For f ∈ L p and α ∈ R, we have

t -α 2 Q t f L p x L q 2 (t -1 dt) t -α 2 Q t f L p x L q 1 (t -1 dt) . Proof : We prove • L ∞ (t -1 dt)
• L q (t -1 dt) for any q ∈ [1, ∞) and the result follows from duality. To get this, we use

Q t = 2 t t 2 Q s t s a+1 P (c) t-s ds s for any Q ∈ StGC a and t ∈ (0, 1] which yields Q t f L p t t 2 Q s f L p dt t t t 2 Q s f q L p dt t 1 q
.

One needs the following bound to keep an accurate track of the constant in different estimates.

Lemma A. . Let r > 0 and α ∈ (-r, r). We have

∞ 0 u 1 + u 2 r u α du u ≤ 2r r 2 -α 2 . t -α |f (t)| dt t ≤ 2r r 2 -α 2 1 0 t -α |f (t)| dt t .
The result then follows for any q ∈ (1, ∞) by interpolation. For the dependance with respect to s, we also interpolate between q = 1 and q = ∞ and conclud with

s 0 tu (t + u) 2 r f (t) dt t ≤ t -β f (t) L ∞ s 0 tu (t + u) 2 r t β dt t ≤ s β-α t -β f (t) L ∞ s 0 tu (t + u) 2 r t α dt t ≤ 2r r 2 -α 2 s β-α u α t -α f (t) L ∞ and 1 0 u -α s 0 tu (t + u) 2 r f (t) dt t du u ≤ s 0 1 0 tu (t + u) 2 r u -α du u |f (t)| dt t ≤ 2r r 2 -α 2 s 0 t -α |f (t)| dt t ≤ 2r r 2 -α 2 s β-α 1 0 t -β |f (t)| dt t . 1 0 Q s Q 1• t P t f • Q 2 t g dt t L 2 x 1 0 ts (t + s) 2 r 2 P t f • Q 2 t g L 2 x dt t g C β 1 0 ts (t + s) 2 r 2 t β 2 P t f L 2 x dt t .
This yields

s -α+β 2 1 0 Q s Q 1• t P t f • Q 2 t g dt t L 2
x L 2 (s -1 ds)

g C β s -α+β 2 1 0 ts (t + s) 2 r 2 t β 2 P t f L 2 x dt t L 2 (s -1 ds) g C β s -α 2 P s f L 2
x L 2 (s -1 ds)

f H α g C β
where we used that α < 0 since P can encode no cancellation and this complete the proof for the third estimate. The proofs for the other estimates on P f g are similar and we only give the details for the resonant term. Let Q ∈ StGC r with r > |α + β| and recall that Π(f, g) is a linear combination of terms

1 0 P • t Q 1 t f • Q 2 t g dt t with Q 1 , Q 2 ∈ StGC b
2 and P ∈ StGC [0,b] . Given s ∈ (0, 1], we have

1 0 Q s P • t Q 1 t f • Q 2 t g dt t L 2 x s 0 Q 1 t f • Q 2 t g L 2 x dt t + 1 s s t r 2 Q 1 t f • Q 2 t g L 2 x dt t
and the result follows again from the lemmas using that α + β > 0.

The dependance of P s with respect to s in given in the following proposition.

Proposition B. . Let s ∈ (0, 1) and a regularity exponent β ∈ (0, 1). Given g ∈ C β , we have

f → P s f g L 2 →H γ s β-γ 4
β -γ g C β for any γ ∈ [0, β).

Proof : Given f ∈ L 2 and γ ∈ [0, β), we want to bound the H γ norm of P s f g hence let Q ∈ StGC r with r > |γ|. Recall that P s f g is a linear combination of terms of the form . Given u ∈ (0, 1], we have

s 0 Q 1• t P t f • Q 2 t g dt t with Q 1 ∈ GC b 2 -2 , Q 2 ∈ StGC
s 0 Q u Q 1• t P t f • Q 2 t g dt t L 2 x s 0 tu (t + u) 2 r 2 P t f • Q 2 t g L 2 x dt t g C β s 0 tu (t + u) 2 r 2 t β 2 P t f L 2 x dt t .
This yields

u -s 2 s 0 Q u Q 1• t P t f • Q 2 t g dt t L 2
x L 2 (u -1 du)

g C β u -γ 2 s 0 tu (t + u) 2 r 2 t β 2 P t f L 2 x dt t L 2 (u -1 du) g C β 4r r 2 -γ 2 s β -γ 2 u -β -β 2 P u f L 2 L 2 (u -1 du) g C β 4r r 2 -γ 2 s β -γ 2 2 k + β -β f H β -β g C β 1 -β s β -γ 2 k + β -β f H β -β
for any β ∈ (γ, β) and P ∈ StGC k using that r ≥ 1. For k ≥ 1, one can take β = β and get

u -γ 2 s 0 Q u Q 1• t P t f • Q 2 t g dt t L 2 x L 2 (u -1 du) s β-γ 2 1 -β g C β f L 2 .
For k = 0, we have

u -γ 2 s 0 Q u Q 1• t P t f • Q 2 t g dt t L 2
x L 2 (u -1 du)

g C β 1 -β s β -γ 2 β -β f L 2 hence taking β = γ+β 2 yields u -γ 2 s 0 Q u Q 1• t P t f • Q 2 t g dt t L 2 x L 2 (u -1 du) s β-γ 4 (1 -β)(β -γ) g C β f L 2 .
Proposition B. . Let s ∈ (0, 1) and a regularity exponent β < 2. Given g ∈ C β , we have

( P f -P s f )g H 2 s β-2 2 f L 2 g C β for any f ∈ L 2 .
Proof : Given f ∈ L 2 , we want to bound the H 2 norm of ( P f -P s f )g hence let Q ∈ StGC r with r > 2. It is a linear combination of terms

1 s Q 1• t P t f • Q 2 t g dt t with Q 1 ∈ GC b 2 -2 , Q 2 ∈ StGC b
2 and P ∈ StGC [0,b] . Given u ∈ (0, 1], we have

1 s Q u Q 1• t P t f • Q 2 t g dt t L 2 x 1 s tu (t + u) 2 r 2 P t f • Q 2 t g L 2 x dt t f L 2 g C β 1 s tu (t + u) 2 r 2 t β 2
dt t using that P t f L 2 f L 2 . This yields

u -1 1 s Q u Q 1• t P t f • Q 2 t g dt t L 2
x L 2 (u -1 du)

g C β u -1 1 s tu (t + u) 2
P 1• t (Q 1 t a • Q 2 t b) dt t ,
while P b a is a linear combination of terms of the form

1 0 Q 3• t P 2 t b • Q 4 t a dt t with Q 1 , Q 2 , Q 4 ∈ StGC b 2 , Q 3 ∈ GC b
2 -2 and P 1 , P 2 ∈ StGC [0,b] . For P 2 ∈ StGC [1,b] , we already have the correct regularity since using that α 1 ∈ (0, 1). We consider P 2 ∈ StGC 0 for the remainder of the proof. For all x ∈ M , we have

1 0 1 0 Q u P 1• t Q 1 t Q 3• s P 2 s a 1 • Q 4 s a 2 • Q 2 t b
P 1• t Q 1 t Q 3• s P 2 s a 1 -a 1 (x) • Q 4 s a 2 • Q 2 t b (x)
K Qu (x, x )P 1• t Q 1 t Q 3• s P 2 s a 1 -a 1 (x ) • Q 4 s a 2 • Q 2 t b (x ) ds s dt t ν(dx ) = K Qu (x, x )K P 1• t (x , x ) Q 1 t Q 3• s P 2 s a 1 -a 1 (x ) • Q 4 s a 2 • Q 2 t b (x )
ds s dt t ν(dx )ν(dx )

+ u 0 K Qu (x, x )K P 1• t (x , x ) a 1 (x ) -a 1 (x ) Q 1 t a 2 • Q 2 t b (x ) dt t ν(dx )ν(dx ) + 1 u K Qu (x, x )K P 1• t (x , x ) a 1 (x ) -a 1 (x ) Q 1 t a 2 • Q 2 t b (x ) dt t ν(dx )ν(dx ) =: A + B + C.
The term A is bounded using cancellations properties. We have , using again that α 1 ∈ (0, 1) and (α 1 + α 2 + β) > 0.

|A| = K QuP 1• t (x, x ) Q 1 t Q 3• s P 2 s a 1 -a 1 (x ) • Q 4 s a 2 • Q 2 t b (x )

=

  b e -tL (b -1)! and -t∂ t P Id. In particular, there exist a polynomial p b of degree (b -1) such that P (b) t = p b (tL) e -tL and p b (0) = 1. The family (P t ) t∈(0,1] regularizes distributions while the family

  (b) t p = p and Q (b)

2

 2

a 2 G

 2 t+s (x, y) with a = min(a, a ).

2 - 2 |I| 2 V

 222 defines a family that belongs to GC b for b large enough. Proof : Since Q ∈ StGC b 2 , there exist I = (i 1 , . . . , i n ), j ∈ N and c ∈ 1, b such that b 2 = |I| + j and Q t = (t I )(tL)

b 2 and

 2 P ∈ StGC[0,b]

2 f

 2 L 2 g C β and the proof is complete. Proposition B. . Let α 1 ∈ (0, 1) and α 2 , β ∈ R. If α 2 + β < 0 and α 1 + α 2 + β > 0, then (a 1 , a 2 , b) → C(a 1 , a 2 , b) extends in a unique bilinear operator from C α1 ×C α2 ×C β to C α1+α2+β and from H α1 × C α2 × C β to H α1+α2+β . Proof : We first consider (a 1 , a 2 , b) ∈ C α1 × C α2 × C β . We want to compute the regularity of C(a 1 , a 2 , b) = Π P a1 a 2 , b -a 1 Π(a 2 , b) using a family Q of StGC r with r > |α 1 + α 2 + β|. Recall that a term Π(a, b) can be written as a linear combination of terms of the form 1 0

C a 1

 1 , a 2 , b (x) = Π P a1 a 2 , b (x) -a 1 (x) • Π(a 2 , b)(x) = Π P a1 a 2 -a 1 (x) • a 2 , b (x) Π P a1-a1(x) a 2 , b (x),since Π is bilinear and a 1 (x) is a scalar and P 1 a 1 = a 1 up to smooth terms. Thus we only have to consider a linear combination of terms of the form

=

  Id up to smooth terms. This gives Q u C(a 1 , a 2 , b) (x) as a linear combination of terms of the form

1 2 ,|B| a 1 α1 a 2 α2 b β x ,x u 0 Ka 1 α1 a 2 α2 b β u 0 t α 1 +α 2 +β 2 dt t a 1 α1 a 2 α2 b β u α 1 +α 2 +β 2

 12002 α1 a 2 α2 b β u α 1 +α 2 +β using that α 1 ∈ (0, 1), P 2 ∈ StGC 0 and (α 1 + α 2 + β) > 0.For the term B, we have Qu (x, x)K P 1• t (x , x )d(x , x ) α1 t

  2 . Thus u, Lu s = P s u X 1 , Lu s + P s u X 2 , Lu s + u s , Lu s = P s u LX 1 , u s + P s u LX 2 , u s + ∇u s , ∇u s and this yields u,Hu = -P s u ξ, u s + P s u LX 2 , u s + ∇u s , ∇u s + u, P ξ u s + Π(u s , ξ) + u, R(u) + Ψ s (u) = -A(u, ξ, u s ) + P s u LX 2 , u s + ∇u s , ∇u s + u, P ξ u s + u, R(u) + Ψ s (u) + (P u -P s u )ξ, u s

α

  Definition. We define the form domain ofH denoted D Ξ ( Proposition . . For s < s 0 (Ξ) and u ∈ L 2 , u ∈ D Ξ (H + ) ⇐⇒ u s = Φ s (u) ∈ H 2 .and we have the boundsu s H 2 Ξ,s H + u L 2 Ξ,s u s H 2 .

					√	H + ) as the closure of the domain under
	the norm
		u DΞ( √	H + ) :=	u, H + u
	Moreover, we have	u ∈ D Ξ (	√	H + ) ⇐⇒ Φ s (u) = u s ∈ H 1
	with the bounds
		u s H 1 Ξ,s u DΞ( √	H + )	Ξ,s u s H 1 .

α-β

Proof : Given any u ∈ H β , we have u = ΓΓ -1 (u) = Γ(u -P s u X 1 -P s u X 2 ). Using proposition 2.4, we get

using the proposition 2.3 and that X (ε) i -X i is i-linear in Ξ ε -Ξ for i ∈ {1, 2}. The second statement follows from

.

This allows to prove density of the domain.

Corollary. The domain D Ξ is dense in H β for any β ∈ [0, α).

Proof : Given f ∈ H 2 , Γ(g ε ) ∈ D Ξ where g ε = Γ -1 ε f ∈ H 2 thus we can conclude with the lemma 2.5 that

The density of H 2 in H β then yields the result.

Taking into account in the previous computation the smooth term e -L coming from the intertwining relation, we are able to define H as an unbounded operator in L 2 with domain D Ξ as follows.

Definition. We define the Anderson Hamiltonian

The parameter s does not appear in the definition of H, it is a tool to study the properties of the operator. Indeed, one has different representations of Hu as

where u s := Φ s (u) and

The different representations of H through the parameter s > 0 will be useful to get different bounds. For example, we can compare the graph norm of H given as

Proof : The Besov embedding in two dimensions implies

Given any u ∈ D Ξ , we get

for any α < 1 and the proof is complete.

-Self-adjointness and spectral properties

We show that H is a closed self-adjoint operator on its dense domain D Ξ ⊂ L 2 . This relies on approximation results and the Babuška-Lax-Milgram theorem. The spectrum is pure point and the eigenvalues verify a min-max principle that allows to get estimates depending on the eigenvalues of L.

Proposition . . The operator H is closed on its domain D Ξ .

Proof :

Proposition 2.6 gives that Φ(u n ) n≥0 is a Cauchy sequence in H 2 hence converges to u ∈ H 2 . Since Φ : L 2 → L 2 is continuous, we have Φ(u) = u hence u ∈ D Ξ . Finally, we have

hence Hu = v and H is closed on D Ξ .

In some sense, the operator H should be the limit of the renormalised H ε as ε goes to 0. Since D(H ε ) = H 2 , one can not compare directly the operators. However given any u ∈ L 2 , we have

where we recall that Γ ε u is implicitly defined for u ∈ L 2 by the relation

Thus for u ∈ D Ξ , the approximation

with u := Φ s (u). The following proposition gives a bound for this quantity which yields the convergence as ε goes to 0 for s is small enough. We do not need to explicit the constant, it depends polynomialy on the enhanced noise Ξ and diverges as s goes to s 0 (Ξ).

Proposition . . Let u ∈ D Ξ and s > 0 small enough. Then

with u s = Φ s (u) and u ε := Γ ε u s . In particular, this implies that H ε Γ ε converges to HΓ in norm as ε goes to 0 as operators from H 2 to L 2 .

using Markov inequality and this yields

where m = E e h Ξ X α .

We proved that H ε converges to H is some sense as ε goes to 0. The following proposition gives the convergence of H ε + k Ξ to H + k Ξ in resolvent sense as ε goes to 0. We do not need to explicit the constant, it depends polynomialy on the enhanced noise Ξ.

Proposition . . Let s ∈ (0, 1) and δ > 0. Then for any constant k Ξ > m 1 δ (Ξ, s) and β ∈ [0, α), we have

We introduce u ε := Γ ε Φ s (u) which converges to u in L 2 and we have

we only have to bound the second term. We have

using Proposition 2.8. In the end, we have

This allows to get a bound on the convergence of λ n (Ξ ε ) to λ n (Ξ) as ε goes to 0.

Corollary . . For all n ∈ N * , we have

In particular, this implies

Proof : We use the min-max principle for (H + k Ξ ) -1 and (H ε + k Ξ ) -1 and denote µ n and µ (ε) n their n-th smallest eigeinvalue with multiplicity. Let D n = Vect(v 1 , . . . , v n ) with v i an eigenfunction associated to µ

n hence with proposition 2.14 we get

Using the same argument with eigeinfunctions associated to (H + k Ξ ) -1 , we get

Thus this gives

and completes the proof with the upper bound on λ n (Ξ).

We conclude this section by giving as corollary the Weyl law for the Anderson Hamiltonian H.

Corollary . . We have

Proof : Proposition 2.12 gives the bounds

for any δ ∈ (0, 1). Using the lower bound, λ n (Ξ) ≤ λ implies

and thus |{n ≥ 1; λ n (Ξ) ≤ λ}| ≤ N λ + m 1 δ (Ξ, s) where N (x) denotes the number of eigenvalues of L smaller than x. Using the upper bound, we get

hence the proof is complete using the result for the Laplace-Beltrami operator.

-Stochastic nonlinear Schrödinger equation

The construction of the Anderson Hamiltonian allows the study of associated evolution equations. This was the motivation for the work [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF] of Gubinelli, Ugurcan and Zachhuber and they studied the nonlinear Schrödinger and wave equations on the torus in two and three dimensions, see the references therein for other approaches. Our work allows to do the same on a two-dimensional manifold. As an example, we give results for the cubic nonlinear Schrödinger equation associated to H. See the work [START_REF] Debussche | The Schrödinger equation with spatial white noise potential[END_REF] of Debussche and Weber for the equation on the torus where they use a Hopf-Cole type transformation. This was extended in [START_REF] Tzvetkov | Two dimensional nonlinear Schrödinger equation with spatial white noise potential and fourth order nonlinearity[END_REF] by Tzvetkov and Visciglia to the fourth order nonlinearity. Define the positive operator

with k Ξ as in Proposition 2.9. Proposition 2.10 yields a characterization of the domain and the form domain which is defined as follows.

and the remainder as

Thus for v D( √

using proposition 2.17. Since every estimates also hold for H + ε with bound uniform in ε, we also get the estimate for H + ε . Applying this result to

yields the general inequality. This inequality can be used for example to study the cubic nonlinear Schrödinger equation with multiplicative noise

with initial condition u 0 ∈ D Ξ . The construction of the operator H immediatly yields the renormalised solution u(t, •) := e -itH u 0 to the linear equation

given any u 0 ∈ D Ξ as done in [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF]. Remark that when one regularizes the question, one also has to consider a suitable sequence of initial data (u (ε) 0 ) ε>0 , it is often refered to as "well-prepared data" in the litterature. This can also be used to solve the associated equation with cubic nonlinearity even if we can not apply the same theorem as Brezis and Gallouët in [START_REF] Brezis | Nonlinear Schrödinger evolution equations[END_REF] since we do not have a control on the cubic term from D Ξ to itself. One could modify the domain taking into account the term Π(X 1 , X 1 ) in X 2 to get a domain stable by multiplication. However since a direct computation as done by Gubinelli, Ugurcan and Zachhuber in [START_REF] Gubinelli | Semilinear evolution equations for the Anderson Hamiltonian in two and three dimensions[END_REF] is enough, it is not necessary. In particular, the proof of the following theorem works exactly as in their work and is left to the reader.

Theorem . . Let T > 0 and u 0 ∈ D Ξ . Then there exists a unique solution

Moreover, u is the L 2 -limit of the solutions

with the initial data u (ε)

The next lemma describes the localisation of the cancellation in our continuous context, including the dependance on s > 0.

Lemma A. . Let r > 0 and α ∈ (-r, r). Given any q ∈ [1, ∞], we have

We also have

for any s > 0 and β ∈ (α, r).

Proof : For q = ∞, we have

which yields the result. For q = 1, we have

Finally, we have the following estimate for integrals.

Lemma A. . Given any α > 0 and q ∈ [1, ∞], we have

Proof : We proceed again by interpolation proving the estimate for q = ∞ and q = 1. Using that α > 0, we have

B -Paracontrolled calculus

We give in this Appendix proofs of estimates needed in paracontrolled calculus. We shall first prove the estimates for the paraproduct P and resonant operator Π in Sobolev spaces. It works as for Hölder spaces with L 2 estimates instead of L ∞ . Proposition B. . Let α, β ∈ (-2b, 2b) be regularity exponent.

Proof : Let f ∈ H α and g ∈ C β with α < 0. We want to compute the regularity H α+β of P f g hence let Q ∈ StGC r with r > |α + β|. Recall that P f g is a linear combination of terms of the form , using that α 1 ∈ (0, 1) and (α 2 + β) < 0. In the end, we have

uniformly in u ∈ (0, 1], so the proof is complete for C. The adaptation of the proof to the case a 1 ∈ H α1 is left to the reader and follows from the estimates of the Appendix A.