On the convexity of the KdV Hamiltonian

T. Kappeler ,∗ A. Maspero ,† J.C. Molnar ,‡ P. Topalov §

September 29, 2018

Abstract We prove that the nonlinear part H∗ of the KdV Hamiltonian Hkdv,

when expressed in action variables I = (In)n>1, extends to a real analytic 2 2 function on the positive quadrant `+(N) of ` (N) and is strictly concave near 0. As a consequence, the differential of H∗ defines a local diffeomor- phism near 0 of `2 ( ). C N Furthermore we prove that the Fourier-Lebesgue spaces Ls,p with −1/2 6 s 6 0 and 2 6 p < ∞, admit global KdV-Birkhoff coordinates. 2 In particular, it means that `+(N) is the space of action variables of the underlying phase space L−1/2,4 and that the KdV equation is globally in time C0-well-posed on L−1/2,4.

Keywords. convexity, KdV Hamiltonian, Korteweg-de Vries equation, integrable PDEs

AMS Subject classification. 37K10 (primary) 35Q53, 34B24 (sec- ondary)

1 Introduction

The Korteweg-de Vries (KdV) equation

3 ∂tu = −∂xu + 6u∂xu (1.1) arXiv:1502.05857v1 [math.AP] 20 Feb 2015 on the circle T = R/Z is well known to be a Hamiltonian PDE. When considered s s on the Sobolev spaces H ≡ H (T, R), s > 1, it can be written in the form

kdv ∂tu = ∂x∂uH ,

∗Supported in part by the Swiss National Science Foundation †Supported in part by the Swiss National Science Foundation ‡Supported in part by the Swiss National Science Foundation §Supported in part by NSF DMS-0901443.

1 kdv 2 where ∂uH denotes the L -gradient of the KdV-Hamiltonian Z 1   kdv 1 2 3 H (u) := (∂xu) + u dx 0 2

kdv and ∂x∂uH is the Hamiltonian vector field defined by the Poisson bracket Z 1 {F,G} := ∂uF ∂x∂uG dx. (1.2) 0 The main purpose of this paper is to study convexity properties of the KdV R 1 Hamiltonian. Note that the mean u0 := 0 u(x) dx is a Casimir of the bracket (1.2) and that the level sets

s s Hc := {u ∈ H : u0 = c}, c ∈ R,

s are symplectic leaves. We concentrate on the leaf H0 only as our results can be easily extended to any other leaf. To state them we first need to describe in some detail the property of the KdV equation of being an integrable PDE. For any s ∈ R and 1 6 p < ∞, let

s,p ` := {z = (zn)n∈ ⊂ : z0 = 0, kzk < ∞}, 0,C Z C s,p where !1/p X ps p kzks,p := hni |zn| , hni := 1 + |n|, (1.3) n∈Z s,p s,p and denote by ` the real subspace {(zn)n∈ ∈ ` : z−n = zn ∀ n 1} of 0 Z 0,C > `s,p . The spaces `s,p( ) = `s,p( , ) and `s,p( ) = `s,p( , ) are defined in an 0,C N N R C N N C s analogous way. Note that the Sobolev spaces H0 can then be described by

Hs = {u ∈ S0 :(u ) ∈ `s,2}, 0 C n n∈Z 0

i2nπx 2 where un = hu, e iI , n ∈ Z, and h·, ·iI denotes the L -inner product on 2 R 1 L (T, C), hf, giI = 0 f(x)g(x) dx extended by duality to a pairing of the S = {f ∈ C∞( , ): f(x+1) = f(x) ∀ x ∈ } and its dual S0 . C R C R C −1 s It is shown in [12, 11] that on H0 and by restriction, on H0 , s > −1, the KdV equation admits canonical real analytic coordinates xn, yn, n > 1, so that the initial value problem of the KdV equation, when expressed in these coordinates, can be solved by quadrature (cf. [15]). The coordinates xn, yn are canonical R 1 in the sense that {xn, yn} = 0 ∂uxn∂x∂uyn dx = 1 for any n > 1 whereas all other brackets between coordinate functions vanish. For our purposes it is convenient to introduce the complex coordinates

xn − iyn xn + iyn zn = √ , z−n = √ , ∀ n > 1, z0 = 0, 2 2 referred to as (complex) Birkhoff coordinates.

September 29, 2018 2 In more detail, the above result says the following: there exists a complex −1 −1 neighborhood W of H0 in H0 (T, C) so that Φ(u) := (zn(u))n∈Z is an analytic map from W to `−1/2,2 with Φ(0) = 0 satisfying the property that for any 0,C s s > −1, the restriction of Φ to H0 is a canonical real analytic diffeomorphism s+1/2,2 onto `0 . The actions I = (In)n>1 defined by In = znz−n > 0, n > 1, then 2s+1,1 2s+1,1 take values on the positive quadrant `+ (N) of ` (N, R),

2s+1,1 2s+1,1 `+ (N) = {(rn)n>1 ∈ ` (N, R): rn > 0 ∀ n > 1}.

kdv 1 In addition, Φ has the property that the KdV Hamiltonian H : H0 → R, when expressed in the coordinates (zn)n∈Z, is a function of the actions (In)n>1 alone. kdv kdv 3,1 By [12] this function, denoted by H , is a real analytic map, H : `+ (N) → R and according to [12, Corollary 15.2] admits an expansion at I = 0 of the form

X 3 3 X 2 8n π In − 3 In + ··· (1.4) n>1 n>1 where the dots stand for higher order terms in I. In particular, at I = 0

kdv 3 3 kdv ∂In H I=0 = 8π n , ∂In ∂Im H I=0 = −6δnm, ∀ n, m > 1. (1.5)

See Section2 for further details on these results. The main aim of this paper is to prove a conjecture of Korotyaev and Kuksin [17], motivated by the pertur- bation theory of the KdV equation. It says that the nonlinear part of the KdV Hamiltonian

∗ kdv X 3 3 H (I) := H (I) − 8π n In (1.6) n>1

2 0,2 extends analytically to a function on `+(N) ≡ `+ (N) and that it is strictly kdv kdv concave near 0. Note that the KdV frequencies, defined as ωn = ∂In H , ∗ ∗ differ from ωn := ∂In H by a constant,

kdv 3 3 ∗ 3,1 ωn (I) = 8π n + ωn(I) ∀ n > 1, ∀ I ∈ `+ (N). Furthermore, we recall that in [15] it was shown that the KdV frequencies kdv −1,1 −1,1 (ω )n 1 analytically extend to a neighborhood V0 of ` ( ) in ` ( ) so n > + N C N that

kdv 3 3 ωn − 8n π = O(1) as n → ∞ locally uniformly on V0. (1.7)

Theorem 1.1 The nonlinear part H∗ of the KdV Hamiltonian admits a real analytic extension to a complex neighborhood V ⊂ `2 ( ) of the positive quadrant C N 2 2 ∗ ∗ `+(N) of ` (N). The map of the renormalized frequencies ω = ∂I H , restricted to V0 ∩ V equals

2 kdv 3 3 V ∩ V0 → ` (N) ,I 7→ ωn (I) − 8n π . C n>1

September 29, 2018 3 ∗ It vanishes at I = 0 and its differential d0ω at I = 0 is given by −6 Id`2 ( ). C N In particular near I = 0 in `2 ( ), ω∗ defines a local diffeomorphism and H∗ is C N strictly concave in the sense that for I in a (sufficiently small) neighborhood of 2 0 in `+(N)

2 ∗ X 2 2 d H (J, J) − J ∀ J = (Jn) ∈ ` ( ). (1.8) I 6 n n>1 N n>1 Remark 1.1. Due to a mismatch of norms, it is difficult to put identity (1.5) in 3,1 the setup with the space `+ (N) to use for studying perturbations of the KdV 2 equation. The space `+(N) is the appropriate one to express the strict concavity of H∗ – see [8]. In combination with results from [17] Theorem 1.1 implies that 2 ∗ 2 dI H 6 0 on all of `+(N).

The second main result of this paper says that the Birkhoff coordinates can be analytically extended to a class of Fourier-Lebesgue spaces. To state our result precisely, denote for s ∈ and 1 p < ∞ by Ls,p ⊂ S0 the Fourier- R 6 0,C C Lebesgue space of distributions on T given by

s,p 0 s,p L := {u ∈ S :(un)n∈ ∈ ` } 0,C C Z 0,C s,p s,p s,p s,p and by L ⊂ L the real subspace L = {u ∈ L : u−n = un ∀ n 1}. 0 0,C 0 0,C > Note that Ls,2 = Hs and by the Sobolev embedding theorem, Ls,p ,→ L−1,2 for 0 0 0,C 0,C any −1/2 6 s 6 0 and 1 6 p < ∞.

Theorem 1.2 For any −1/2 6 s 6 0 and 2 6 p < ∞, the restriction of s,p s+1/2,p s,p s+1/2,p Φ to L0 takes values in `0 and Φ: L0 → `0 is a real analytic s,p diffeomorphism, and therefore provides global Birkhoff coordinates on L0 .

−1/2,4 0,4 Remark 1.2. Note that for any u ∈ L0 , Φ(u) = (zn)n∈Z ∈ `0 and hence 2 the corresponding sequence of action variables (In = znz−n)n>1 is in `+(N). 2 Thus the space of action variables `+(N) considered in Theorem 1.1 stems from −1/2,4 the underlying phase space L0 which admits global Birkhoff coordinates. −1/2,4 Although our principal interest lies in the space L0 , Theorem 1.2 is s,p stated for the family of spaces L0 with −1/2 6 s 6 0, 2 6 p < ∞ since the proof of the more general result does not require any additional effort. With some more work, Theorem 1.2 can be proved to hold for a larger family of spaces s,p L0 . For any complex Borel measure µ on T, the distribution µ − µ(T) is an example of an element of Ls,p for any pair −1/2 s < 0, 2 p < ∞ satisfying 0,C 6 6 |s|p > 1. In particular, µ − µ( ) ∈ L−1/2,4. T 0,C Theorem 1.2 has an immediate application to the initial value problem of the KdV equation. As already mentioned above it was shown in [15] that the kdv KdV frequencies (ωn (I))n>1 analytically extend to a complex neighborhood

September 29, 2018 4 −1,1 of the positive quadrant `+ (N), satisfying (1.7). In Birkhoff coordinates on −1 H0 , the KdV equation (1.1) takes the form

kdv kdv z˙n = iωn zn , z˙−n = −iωn z−n ∀ n > 1. We refer to [15] for a discussion of the notion of solutions of the initial value −1 0 problem of KdV in spaces of distributions such as H0 and the one of C - well-posedness. Theorem 1.2 leads to the following extension of the results of [15]:

Theorem 1.3 The initial value problem of the Korteweg-de Vries equation (1.1) 0 s,p is globally in time C -well-posed on the Fourier-Lebesgue spaces L0 with −1/2 6 s 6 0 and 2 6 p < ∞.

−1 Other applications concern the characterization of potentials in H0 which s,p are in L0 for −1/2 6 s 6 0, 2 6 p < ∞, in terms of the decay of the sequence of gap lengths of the spectrum as well as a spectral description of s,p sets in L0 – see Corollary 5.6 and Corollary 5.7.

Method of proof: Inspired from the work in [17] we derive a formula for the nonlinear part of the KdV Hamiltonian Hkdv(u) for finite gap potentials, using 2 the asymptotics of the discriminant ∆(λ, u) of the operator −∂x + u at λ = ∞ kdv P∞ 3 3 and the residue calculus. It leads to a formula for H (u) − n=1 8π n In for P smooth potentials in terms of the weighted sum (8nπ)Rn of functionals n>1 Rn, each of which is an expression involving only the discriminant ∆(λ, u) – 2 or equivalently the periodic spectrum of −∂x + u. Using that by [11], ∆(λ, u) −1 can be analytically extended to H0 one shows that the same holds true for 2 each functional Rn. The asymptotics of the periodic eigenvalues of −∂x + u, P related to Theorem 1.2, then imply that (8nπ)Rn converges absolutely n>1 −1/2,4 for u in a complex neighborhood of L0 , yielding the claimed statement for H∗. The starting point of the proof of Theorem 1.2 are the Birkhoff coordinates −1 constructed in [11] on H0 . The main ingredients are asymptotic estimates of the periodic and Dirichlet eigenvalues in the case where u is in Ls,p , −1/2 0,C 6 s 6 0, 2 6 p < ∞.

Related works: Using the conformal mapping theory together with the heights 2 associated with the periodic eigenvalues of −∂x + u, Korotyaev and Kuksin ∗ 2 proved in [17] that H can be continuously extended to `+(N), leaving open the question if the KdV frequencies, which are the dynamically relevant quantities, 2 can be extended to `+(N) as well. Using different techniques, we succeed in answering this question. Actually we prove much more. Theorem 1.1 and Theorem 1.2 together show that H∗ and hence the KdV frequencies can be 2 2 indeed extended analytically to `+(N) and that `+(N) is the space of action −1/2,4 variables with underlying phase space L0 .

September 29, 2018 5 Using results on the KdV frequencies established in [2], it was also shown ∗ 2 in [17] that H is concave on all of `+(N) and hence Theorem 1.1 implies that 2 ∗ 2 ∗ dI H 6 0 on `+(N) as well. The problem whether H is strictly concave on all 2 of `+(N) is still open – see [2,3, 12, 18] as well as the paper [7], concerned with the study of the frequencies of the Toda lattice. Finally, we note that the initial value problem of the KdV equation has been intensively studied – see [1] for a comprehensive account.

Organization of the paper: In Section2 we introduce additional notation and review known results used throughout the paper. In Section3 we prove asymp- totic estimates for the periodic and Dirichlet eigenvalues of −∂2 +u for u ∈ Ls,p x 0,C with −1/2 6 s 6 0 and 2 6 p < ∞. They are used in the proofs of Theorem 1.1 and Theorem 1.2 provided in Section4 and Section5 respectively.

2 Preliminaries

In [11] it is shown that the Birkhoff map of the KdV equation on T = R/Z, defined and studied in detail in [12], can be analytically extended to a real analytic diffeomorphism

−1 −1/2,2 Φ: H0 → ` . A key ingredient of the construction of the Birkhoff map and its extension to −1 H0 are spectral quantities of the differential operator

2 L(q) = −∂x + q, (2.1) which appears in the formulation of the KdV equation. In this section we review properties of Φ as well as spectral properties of L(q) used throughout the paper. Let q be a complex potential in H−1 := H−1( / , ). In order to treat 0,C 0 R Z C periodic and antiperiodic boundary conditions at the same time, we consider 2 −1 the differential operator L(q) = −∂x + q, on H (R/2Z, C) with domain of definition H1(R/2Z, C). See AppendixB for a more detailed discussion. The of L(q), while classical for q ∈ L2 , has been only fairly recently 0,C extended to the case q ∈ H−1 – see e.g. [5, 10, 14, 16, 21] and the references 0,C therein. The spectrum of L(q), called the periodic spectrum of q and denoted by spec L(q), is discrete and the eigenvalues, when counted with their multiplicities and ordered lexicographically – first by their real part and second by their imaginary part – satisfy

+ − + ± 2 2 2 λ0 (q) 4 λ1 (q) 4 λ1 (q) 4 ··· , λn (q) = n π + n`n. (2.2)

Furthermore, we define the gap lengths γn(q) and the mid points τn(q) by λ+(q) + λ−(q) γ (q) := λ+(q) − λ−(q), τ (q) := n n , n 1. n n n n 2 >

September 29, 2018 6 −1 2 For q ∈ H we also consider the operator Ldir(q) defined as the operator −∂ +q 0,C x −1 1 on Hdir ([0, 1], C) with domain of definition Hdir([0, 1], C). See AppendixB as well as [5, 10, 14, 16, 21] for a more detailed discussion. The spectrum of Ldir(q) is called the Dirichelt spectrum of q. It is also discrete and given by a sequence of eigenvalues (µn)n>1, counted with multiplicities, which when ordered lexicographically satisfies

2 2 2 µ1 4 µ2 4 µ2 4 ··· , µn = n π + n`n. (2.3)

s,p s,p Finally, define ` ( ) ≡ ` ( , ) to be the space of sequences z = (zn)n 1 ⊂ C N N C > C 1/p P sp p s,p s,p with kzk := hni |zn| and denote by ` ( ) ≡ ` ( , ) the real s,p n>1 N N R s,p subspace of real sequences z = (zn)n 1 ⊂ in ` ( ). > R C N In [11] the following asymptotic estimates are proved:

Theorem 2.1 ([11]) There exists a complex neighborhood W of H−1 in H−1 0 0,C so that the following holds:

−1,2 (i) For any q ∈ W, the sequence (γn(q))n 1 is in ` ( ) and the map W → > C N −1,2 ` ( ), q 7→ (γn(q))n 1 is locally bounded. C N > (ii) The mappings

2 2 2 2 q 7→ (n π − τn(q))n>1 and q 7→ (n π − µn(q))n>1

are analytic from W into `−1,2( ). C N By [11] the discriminant ∆(λ, q) of L(q) admits an analytic extension to C × W and we have the product representation Y (λ+ − λ)(λ− − λ) ∆2(λ, q) − 4 = −4(λ − λ+) m m . (2.4) 0 m4π4 m>1 The λ-derivative ∆• of the discriminant ∆ is analytic on C×W, too, and admits the product representation Y λ• − λ ∆•(λ) = − m (2.5) m2π2 m>1

• where (λm)m>1 ⊂ C are ordered lexicographically

• • • 2 2 2 λ1 4 λ2 4 ··· , λn = n π + n`n . −1 For each potential q ∈ H0 there exists a complex neighborhood Wq of q in H−1 such that the closed intervals 0,C + − + G0 = {λ0 + t : −∞ < t 6 0},Gn = [λn , λn ], n > 1, are disjoint from each other for every p ∈ Wq. By [11] there exist mutually disjoint neighborhoods Un (C, n > 0, called isolating neighborhoods, which satisfy:

September 29, 2018 7 • (i) Gn, λn, and µn are contained in the interior of Un for every p ∈ Wq,

(ii) there exists a constant c > 1 such that for m 6= n,

−1 2 2 2 2 c |m − n | 6 dist(Un,Um) 6 c|m − n |,

2 2 (iii) Un = {λ ∈ C : |λ − n π | 6 n/2} for n sufficiently large.

In the sequel, for any q ∈ W, Wq denotes a neighborhood of q such that a common set of isolating neighborhoods for all p ∈ Wq exists. We shrink W, if −1 necessary, such that W is contained in the union of all Wq with q ∈ H0 . Following [11], for q ∈ W, one can define action variables for the KdV equation by

1 Z λ∆•(λ) In = dλ, n 1. (2.6) pc 2 > π Γn ∆ (λ) − 4

Here Γn denotes any sufficiently close counter clockwise oriented circuit around pc 2 S Gn, and the canonical root ∆ (λ) − 4 is defined on C\ n>0 Gn with γ0 = ∞, γn6=0 where, for q real, the sign of the root is determined by

pc 2 + − i ∆ (λ) − 4 > 0, λ0 < λ < λ1 , and for q ∈ W it is defined by continuous extension. The Dirichlet eigenvalues and the discriminant can be used to construct the angles θk(q), k > 1, which are conjugated to the actions In(q), n > 1, (cf. [11, § 3.4]). For any given k > 1 the action Ik is a real analytic function on W (cf. [11, Proposition 3.3]), whereas the angle θk is defined modulo 2π on W \ Zk and is a real analytic function on W \ Zk when considered modulo π. Here

Zk := {q ∈ W : γk(q) = 0} , k > 1. (2.7)

−1 It was shown in [11, Proposition 4.3] that Zk ∩H0 is a real analytic submanifold −1 of H0 of codimension two. By [11, § 6] the following commutator relations hold for any m, n > 1

{Im,In} = 0, on W, (2.8)

{Im, θn} = δnm, on W \ Zn, (2.9)

{θm, θn} = 0, on W \ (Zm ∪ Zn). (2.10)

−1 For any q ∈ H0 \ Zk with k > 1 define

p+ −iθk(q) p+ iθk(q) zk(q) := Ik(q) e , z−k(q) := Ik(q) e . (2.11)

September 29, 2018 8 −1 It is shown in [11, § 5] that the mappings H0 \Zk → C, q 7→ z±k(q), analytically extend to a common neighborhood W of H−1 in H−1. The Birkhoff map is then 0 0,C defined as follows

−1/2,2 Φ: W → ` , q 7→ Φ(q) := (zk(q)) (2.12) 0,C k∈Z with z0(q) = 0.

Theorem 2.2 ([11]) The mapping Φ: W → `−1/2,2 is analytic and satisfies 0,C the following properties: R 1 (i) Φ is canonical in the sense that {zn, z−n} = 0 ∂uzn∂x∂uz−n dx = i for all n > 1, whereas all other brackets between coordinate functions vanish.

s s+1/2,2 (ii) For any s −1, the restriction Φ ≡ Φ| s is a map Φ| s : H → ` > s H0 H0 0 0 which is a bianalytic diffeomorphism.

(iii) The KdV Hamiltonian, when expressed in the new variables Hkdv ◦ Φ−1, 3/2,2 is defined on `0 and depends on the action variables alone. In fact, it 3,1 is a real analytic function of the actions on the positive quadrant `+ (N), introduced in the introduction.

We will also need the following result (cf. [11, § 3]).

Theorem 2.3 ([11]) After shrinking, if necessary, the complex neighborhood W of H−1 in H−1 constructed above satisfies: 0 0,C 2 −1 (i) The quotient In/γn, defined on H0 \ Zn, extends analytically to W for any n. Moreover, for any ε > 0 and any q ∈ W there exists n0 > 1 and an open neighborhood Wq of q in W so that

In 8nπ 2 − 1 6 ε, ∀ n > n0 ∀ p ∈ Wq. (2.13) γn

(ii) The Birkhoff coordinates (zn)n∈Z are analytic as maps from W into C and fulfill locally uniformly in W and uniformly for n > 1, the estimate |γ | + |µ − τ | |z | = O n √ n n . ±n n

(iii) For any q ∈ W and n > 1 one has In(q) = 0 if and only if γn(q) = 0. In particular, Φ(0) = 0.

3 Asymptotic estimates for the periodic and Dirichlet eigen- values

In this section we prove asymptotic estimates of the quantities |γn| and |µn −τn| appearing in the estimates of the Birkhoff coordinates of Theorem 2.3 (i) - (ii)

September 29, 2018 9 in the case where q ∈ Ls,p with −1/2 s 0 and 2 p < ∞. They are needed 0,C 6 6 6 in the proofs of Theorem 1.1 and Theorem 1.2. −1 Recall from (2.2) that for q ∈ H , the sequence of gap lengths (γn(q))n 1 0,C > is in `−1,2( ). For q ∈ Ls,p , this sequence has the following stronger decay: C N 0,C Theorem 3.1 For any q ∈ Ls,p with −1/2 s 0 and 2 p < ∞, one has 0,C 6 6 6 s,p (γn(q))n 1 ∈ ` ( ) and the map > C N s,p s,p L → ` ( ), q 7→ (γn(q))n 1, 0,C C N > is locally bounded. o

−1 −1,2 Recall from (2.2)–(2.3) that for any q ∈ H , (τn(q) − µn(q))n 1 is in ` ( ). 0,C > C N For q in Ls,p , the decay of this sequence is again stronger: 0,C Theorem 3.2 For any q ∈ Ls,p with −1/2 s 0 and 2 p < ∞, one has 0,C 6 6 6 s,p (τn − µn(q))n 1 ∈ ` ( ) and the map > C N s,p s,p L → ` ( ), q 7→ (τn(q) − µn(q))n 1, 0,C C N > is locally bounded. o In the remainder of this section we prove Theorem 3.1 and Theorem 3.2. For this purpose we adapt the methods developed in [9, 20,4] for potentials q ∈ L2 to singular potentials, using results from [10] and [11].

3.1 Setup

As in the introduction we extend the L2-inner product of L2(R/2Z) by duality to S0 ( /2 ) × S ( /2 ). For s ∈ and 1 p < ∞ denote by Ls,p the space C R Z C R Z R 6 ?,C of periodic, complex valued distributions f ∈ S0 ( /2 ) so that the sequence of C R Z s,p their Fourier coefficients fn = hf, eni is in the space ` = {z = (zn)n∈ ⊂ : C Z C kzks,p < ∞}. We write !1/p X sp p kfks,p := k(fn)n∈Zks,p = hni |fn| . n∈Z In the sequel we will identify a potential q ∈ Ls,p with the corresponding 0,C P s,p element qnen in L where qn is the nth Fourier coefficient of the poten- n∈Z ?,C tial obtained from q by viewing it as a distribution on R/2Z instead of R/Z, i.e., q2n = hq, e2ni, whereas q2n+1 = hq, e2n+1i = 0 and q0 = hq, 1i = 0. We denote by V the operator of multiplication by q with domain H1( /2 ). See C R Z AppendixB for a detailed discussion of this operator as well as the operator L(q) introduced in (2.1). When expressed in its Fourier series, the image V f of f = P 1 P P  fnen ∈ H ( /2 ) is the distribution V f = qn−mfm en ∈ n∈Z C R Z n∈Z m∈Z H−1( /2 ). To prove the asymptotic estimates of the gap lengths stated in The- C R Z orem 3.1 we need to study the eigenvalue equation L(q)f = λf for sufficiently

September 29, 2018 10 large periodic eigenvalues λ. For q ∈ H−1, the domain of L(q) is H1( /2 ) and 0,C C R Z hence the eigenfunction f is an element of this space. It is shown in AppendixB that for q ∈ Ls,p one has f ∈ Ls+2,p and ∂2f, V f ∈ Ls,p . Note that for q = 0 0,C ?,C x ?,C + − 2 2 and any n > 1, λn (0) = λn (0) = n π , and the eigenspace corresponding to + − the double eigenvalue λn (0) = λn (0) is spanned by en and e−n. Viewing for n ± ± large L(q) − λn (q) as a perturbation of L(0) − λn (0), we are led to decompose s,p s,p L into the direct sum L = Pn ⊕ Qn, ?,C ?,C

Pn = span{en, e−n}, Qn = span{ek : k 6= ±n}.

The projections onto Pn and Qn are denoted by Pn and Qn, respectively. It is convenient to write the eigenvalue equation Lf = λf in the form Aλf = V f, 2 where Aλf = ∂xf + λf and V denotes the operator of multiplication with q. Since Aλ is a Fourier multiplier, we write f = u + v, where u = Pnf and v = Qnf, and decompose the equation Aλf = V f into the two equations

Aλu = PnV (u + v),Aλv = QnV (u + v), (3.1)

± 2 2 2 referred to as P - and Q-equation. As λn (q) = n π + n`n, we conclude that for ± n sufficiently large, λn (q) ∈ Sn where Sn denotes the closed vertical strip

2 2 Sn := {λ ∈ C : |Rλ − n π | 6 12n}, n > 1. S Note that {λ ∈ : Rλ 0} ⊂ Sn. Given any n 1, u ∈ Pn, and C > n>1 > λ ∈ Sn, we derive in a first step from the Q-equation an equation for w := V v which for n sufficiently large can be solved as a function of u and λ. In a second step, for λ a periodic eigenvalue in Sn, we solve the P equation for u after having substituted in it the expression of V v. The solution of the Q-equation is then easily determined. Towards the first step note that for any λ ∈ Sn, s+2,p Aλ : Qn ∩ L → Qn is boundedly invertible as for any k 6= n, ?,C 2 2 2 2 2 2 min |λ − k π | > min |Rλ − k π | > |n − k | > 1. λ∈Sn λ∈Sn In order to derive from the Q-equation an equation for V v, we apply to it the −1 operator VAλ to get

−1 V v = VAλ QnV (u + v) = TnV (u + v), where

−1 s,p s,p Tn ≡ Tn(λ) := VA Qn : L → L . λ ?,C ?,C It leads to the following equation for w := V v

(Id −Tn(λ))w = Tn(λ)V u. (3.2)

September 29, 2018 11 To show that Id −Tn(λ) is invertible, we introduce for any l ∈ Z the shifted norm of f ∈ Ls,p , ?,C !1/p X p p kfks,p;l := kfelks,p = hn + li |fn| , n∈Z and denote by kTnks,p;l the corresponding operator norm. Furthermore, we s,p denote by RN f, N 1 the tail of the Fourier series of f ∈ L , > ?,C X RN f = fnen. |n|>N Lemma 3.1 Let −1/2 6 s 6 0, 2 6 p < ∞, and n > 1 be given. For any s,p q ∈ L and λ ∈ Sn, 0,C s,p s,p Tn(λ): L → L ?,C ?,C is a bounded linear operator satisfying the estimate   kTn(λ)ks,p;±n 6 cs,p εs,p(n)kqks,p + kRn/2qks,p , where cs,p is a constant depending only on s and p and εs,p is defined here as  1 , 2 p < 1 ,  n1−|s| 6 |s|  1−1/p ε (n) := (loghni) 1 s,p 1−|s| , p = |s| , o  n  1 1 n1−2|s|+1/p , |s| < p < ∞.

−1 s,p s+2,p Proof. Note that A : L → L is bounded for any λ ∈ Sn and hence λ ?,C ?,C s,p s,p −1 for any f ∈ L , q ∈ L , and λ ∈ Sn, the multiplication of A Qnf with q, ?,C 0,C λ defined by   X X qm−kfk T (λ)f = VA−1Q f = e n λ n  λ − k2π2  m m∈Z |k|6=n is a distribution in S0 ( /2 ). The norm kT fk satisfies C R Z n s,p;n  p p X s X |qm−k||fk| kT fk hm + ni . n s,p;n 6  |n2 − k2|  m∈Z |k|6=n We write the latter expression in the form  p |s| |s| X X hk + ni hm − ki |qm−k| |fk| p   6 2 (I + II), |n + k||n − k|hm + ni|s| hm − ki|s| hk + ni|s| m∈Z |k|6=n P where I and II are defined by splitting the sum |k|6=n into a sum over {|k|= 6 n : |m − k| > n/2} and one over {|k| 6= n : |m − k| < n/2}. Using that hm − ki 6 hm + nihn + ki we conclude  p 2|s| X X hk + ni 1 |qm−k| |fk| I 6  1{|m−k| n/2}  . > |n + k| |n − k| hm − ki|s| hk + ni|s| m∈Z |k|6=n

September 29, 2018 12 2|s| As −1/2 6 s 6 0 and hk + ni /|n + k| 6 (1 + |n + k|)/|n + k| 6 2, one then P gets by Hölder’s inequality applied to the sum |k|6=n  p/p0 p0 X 2 p p I 6  p0  kRn/2qks,pkfks,p;n. |k|6=n |n − k| To estimate II, we show that hm − ki 1 1. (3.3) hn − kihm + ni {|m−k|

Indeed, if |n − k| > n/2, this is obvious as |m − k| < n/2. In case |n − k| 6 n/2 one has that |n + k| > 2n − |n − k| > 3n/2 and hence

|m + n| > ||n + k| − |m − k|| > n, which implies (3.3). With (3.3) and Hölder’s inequality it then follows that  p |s| |s| X X 2 2 |qm−k| |fk| II 6  1{|m−k| 0 which only depends on s and p. Altogether we have shown that for some constant cs,p > Cs,p,   kTnfks,p;n 6 cs,p εs,p(n)kqks,p + kRn/2qks,p kfks,p;n. Going through the arguments of the proof one sees that the same kind of esti- v mates also lead to the claimed bound for kTnfks,p;−n. Lemma 3.1 can be used to solve, for n sufficiently large, the equation (3.2) as well as the Q-equation (3.1) in terms of any given u ∈ Pn and λ ∈ Sn. Corollary 3.2 For any q ∈ Ls,p with −1/2 s 0 and 2 p < ∞, there 0,C 6 6 6 exists ns,p = ns,p(q) > 1 so that for any n > ns,p and λ ∈ Sn, Tn(λ) is a 1/2 con- s,p traction on L with respect to the norms shifted by ±n, kTn(λ)k 1/2. ?,C s,p;±n 6 The threshold ns,p(q) can be chosen locally uniformly in q. As a consequence, equation (3.2) and (3.1) can be uniquely solved for any given u ∈ Pn, λ ∈ Sn,

s,p −1 wu,λ = Kn(λ)Tn(λ)V u ∈ L ,Kn ≡ Kn(λ) := (Id −Tn(λ)) , (3.4) ?,C −1 −1 −1 s+2,p vu,λ = A QnV u + A Qnwu,λ = A QnKnV u ∈ L ∩ Qn. (3.5) λ λ λ ?,C

September 29, 2018 13 In particular, one has wu,λ = V vu,λ. We increase ns,p, if necessary, so that in addition

± 2 2 |λn | 6 (ns,p − 1) π + ns,p/2, ∀ n < ns,p, ± 2 2 |λn − n π | 6 n/2, ∀ n > ns,p, (3.6) ± λn are 1-periodic [1-antiperiodic] if n even [odd] ∀ n > ns,p. o Remark 3.3. By the same approach, one can study the inhomogeneous equation

(L − λ)f = g, g ∈ Ls,p , ?,C for λ ∈ Sn and n > ns,p. Writing f = u + v and g = Png + Qng, the Q-equation becomes

Aλv = QnV (u + v) − Qng = QnV v + Qn(V u − g) leading for any given u ∈ Pn to the unique solution w of the equation corre- sponding to (3.2)

s,p w = V v = KnTn(V u − g) ∈ L , ?,C s+2,p and, in turn, to the unique solution v ∈ L ∩ Qn of the Q-equation ?,C −1 −1 −1 v = Aλ Qn(V u − g) + Aλ QnKnTn(V u − g) = Aλ QnKn(V u − g).

3.2 Reduction

In a next step we study the P -equation Aλu = PnV (u + v) of (3.1). For n > ns,p(q), u ∈ Pn, and λ ∈ Sn, substitute in it the solution wu,λ of (3.2), obtained in (3.4),

Aλu = PnV u + Pnwu,λ = Pn(Id +KnTn)V u.

Using that Id +KnTn = Kn one then obtains Aλu = PnKnV u or Bnu = 0, where

Bn ≡ Bn(λ): Pn → Pn, u 7→ (Aλ − PnKn(λ)V )u. (3.7)

Lemma 3.4 Assume that q ∈ Ls,p with −1/2 s 0 and 2 p < ∞. Then 0,C 6 6 6 for any n > ns,p with ns,p given by Corollary 3.2, λ ∈ Sn is an eigenvalue of L(q) if and only if det(Bn(λ)) = 0. o

Proof. Assume that λ ∈ Sn is an eigenvalue of L = L(q). By Lemma B.3 there s+2,p exists 0 6= f ∈ L so that Lf = λf. Decomposing f = u + v ∈ Pn ⊕ Qn ?,C it follows by the considerations above and the assumption n > ns,p that u 6= 0 and Bn(λ)u = 0. Conversely, assume that det(Bn(λ)) = 0 for some λ ∈ Sn. Then there exists 0 6= u ∈ Pn so that Bn(λ)u = 0. Since n > ns,p, there exists s+2,w,p wu,λ and vu,λ as in (3.4) and (3.5), respectively. Then v ≡ vu,λ ∈ L ∩ Qn ?,C

September 29, 2018 14 solves the Q-equation by Corollary 3.2. To see that u solves the P -equation note that Bn(λ)u = 0 implies that

Aλu = PnKnV u = Pn(Id +KnTn)V u.

As by (3.4), wu,λ = KnTnV u, and as wu,λ = V v, one sees that indeed

Aλu = PnV u + PnV v. v

Remark 3.5. Solutions of the inhomogeneous equation (L−λ)f = g for g ∈ Ls,p , ?,C λ ∈ Sn, and n > ns,p can be obtained by substituting into the P -equation

Aλu = PnV u + PnV v − Png the expression for V v obtained in Remark 3.3, V v = KnTn(V u − g), to get

Aλu = PnV u + PnKnTnV u − Png − PnKnTng

= Pn(Id +KnTn)V u − Pn(Id +KnTn)g.

Using that Id +KnTn = Kn one concludes that

Bn(λ)u = −PnKn(λ)g. (3.8)

Conversely, for any solution u of (3.8), f = u + v, with v being the element in Ls+2,p given in Remark 3.3, satisfies (L − λ)f = g and f ∈ Ls+2,p. ?,C ?,C

We denote the matrix representation of a linear operator F : Pn → Pn with respect to the orthonormal basis en, e−n of Pn also by F ,  hF e , e i hF e , e i  F = n n −n n . hF en, e−ni hF e−n, e−ni In particular,

 2 2    λ − n π 0 an bn Aλ = 2 2 ,PnKnV = , 0 λ − n π b−n a−n where for any λ ∈ Sn and n > ns,p the coefficients of PnKnV are given by

an ≡ an(λ) := hKnV en, eni, a−n ≡ a−n(λ) := hKnV e−n, e−ni,

bn ≡ bn(λ) := hKnV e−n, eni, b−n ≡ b−n(λ) := hKnV en, e−ni.

By a straightforward verification it follows from the expression of an in terms of the representation of Kn and V in Fourier space that for any n > ns,p

an = hKnV e−n, e−ni = a−n.

Hence the diagonal of Bn is a complex multiple of the identity and

 2 2  λ − n π − an −bn, Bn = 2 2 . (3.9) −b−n λ − n π − an

September 29, 2018 15 Lemma 3.6 Suppose q ∈ Ls,p with −1/2 s 0 and 2 p < ∞. Then for 0,C 6 6 6 any n > ns,p, with ns,p as in Corollary 3.2, the coefficients an(λ) and b±n(λ) are analytic functions on the strip Sn and for any λ ∈ Sn

|an(λ)| 6 2kTn(λ)ks,p;±nkqks,p, s hni |bn(λ) − q2n| 6 2kTn(λ)ks,p;nkqks,p, s hni |b−n(λ) − q−2n| 6 2kTn(λ)ks,p;−nkqks,p. o

Proof. Since kTn(λ)ks,p;±n 6 1/2 for n > ns,p and λ ∈ Sn, the series expan- sions of an(λ) and b±n(λ) converge uniformly on Sn to functions analytic in λ.

Moreover, we obtain from the identity Kn = Id +TnKn

bn = hV e−n, eni + hTnKnV e−n, eni = q2n + hTnKnV e−n, eni.

Furthermore, for any f ∈ Ls,p we compute ?,C s s hni |hf, eni| = hni |hfen, e2ni| 6 kfks,p;n.

Consequently, using that kTnks,p;n 6 1/2 and hence kKnks,p;n 6 2, one gets

s hni |bn − q2n| 6 kTnKnV e−nks,p;n 6 2kTnks,p;nkV e−nks,p;n

= 2kTnks,p;nkqks,p.

The estimates for |b−n − q−2n| and |an| are obtained in a similar fashion. v

2 2 The preceding lemma implies that for any n > ns,p, det Bn(λ) = (λ−n π − 2 2 2 2 an) − bnb−n is an analytic function of λ ∈ Sn which is close to (λ − n π ) for n > ns,p sufficiently large. Lemma 3.7 Suppose q ∈ Ls,p with −1/2 s 0 and 2 p < ∞. Choose 0,C 6 6 6 ns,p = ns,p(q) > 1 as in Corollary 3.2 and increase it if necessary so that √ kqks,p 6 ns,p/8. (3.10)

Then for any n > ns,p, det(Bn(λ)) has exactly two roots ξn,1 and ξn,2 in Sn counted with multiplicity. They are contained in

2 2 1/2 Dn := {λ : |λ − n π | < 4n kqks,p} ⊂ Sn and satisfy

√ 1/2 |ξn,1 − ξn,2| 6 6 sup |bn(λ)b−n(λ)| . o (3.11) λ∈Sn

Proof. Since for any n > ns,p and λ ∈ Sn, kTn(λ)ks,p;±n 6 1/2, one concludes from the preceding lemma that |an(λ)| 6 kqks,p and, with |b±n(λ)| 6 |q±2n| + |b±n(λ) − q±2n|, that

|s| |b±n(λ)| 6 2hni kqks,p.

September 29, 2018 16 Therefore, for any λ, µ ∈ Sn,

1/2  |s| |an(µ)| + |bn(λ)b−n(λ)| 6 1 + 2hni kqks,p (3.12) 1/2 2 2 < 4n kqks,p = inf |λ − n π |. λ∈∂Dn

It then follows that det Bn(λ) has no root in Sn \ Dn. Indeed, assume that 2 2 1/2 ξ ∈ Sn is a root, then |ξ − n π − an(ξ)| = |bn(ξ)b−n(ξ)| and hence 2 2 1/2 1/2 |ξ − n π | 6 |an(ξ)| + |bn(ξ)b−n(ξ)| < 4n kqks,p implying that ξ ∈ Dn. In addition, (3.12) implies that by Rouché’s theorem 2 2 2 2 the two analytic functions λ − n π and λ − n π − an(λ), defined on the strip

Sn have the same number of roots in Dn counted with multiplicities. As a 2 2 2 consequence (λ − n π − an(λ)) has a double root in Dn. Finally, (3.12) also implies that

1/2 2 2 sup |bn(λ)b−n(λ)| < inf |λ − n π | − sup |an(λ)| λ∈∂D λ∈Sn n λ∈Sn 2 2 6 inf |λ − n π − an(λ)| λ∈∂Dn 2 2 2 and hence again by Rouché’s theorem, the analytic functions (λ−n π −an(λ)) 2 2 2 and (λ − n π − an(λ)) − bn(λ)b−n(λ) have the same number of roots in Dn. 2 2 2 Altogether we thus have established that det(Bn(λ)) = (λ − n π − an(λ)) − bn(λ)b−n(λ) has precisely two roots ξn,1, ξn,2 in Dn.

To estimate the distance of the roots, write det Bn(λ) as a product g+(λ)g−(λ) 2 2 p where g±(λ) = λ − n π − an(λ) ∓ ϕn(λ) and ϕn(λ) = bn(λ)b−n(λ) with an arbitrary choice of the sign of the root for any λ. Each root ξ of det(Bn) is either a root of g+ or g− and thus satisfies 2 2 ξ ∈ {n π + an(ξ) ± ϕn(ξ)}. As a consequence,

|ξn,1 − ξn,2| |an(ξn,1) − an(ξn,2)| + max|ϕn(ξn,1) ± ϕn(ξn,2)| 6 ± (3.13) 6 sup |∂λan(λ)||ξn,1 − ξn,2| + 2 sup |ϕn(λ)|. λ∈Dn λ∈Dn

1/2 1/2 Since kqks,p 6 ns,p /8 6 n /8 and hence

dist(Dn, ∂Sn) > 12n − n/2 > 6n one concludes from Cauchy’s estimate kqk supλ∈Sn |an(λ)| s,p 1 sup |∂λan(λ)| 6 6 6 . λ∈Dn dist(Dn, ∂Sn) 6n 6 Therefore, by (3.13), 2 |ξn,1 − ξn,2| 6 6 sup |bn(λ)b−n(λ)| λ∈Dn as claimed. v

September 29, 2018 17 3.3 Proof of Theorem 3.1

Let q ∈ Ls,p with −1/2 s 0 and 2 p < ∞. The eigenvalues of L(q), when 0,C 6 6 6 listed with lexicographic ordering, satisfy

+ − + ± 2 2 2 λ0 4 λ1 4 λ1 4 ··· , and λn = n π + n`n. ± ± Hence for n > ns,p with ns,p as in (3.6), λn ∈ Sn and λn ∈/ Sk for any k 6= n. It − + then follows from Lemma 3.4 and Lemma 3.7 that {ξn,1,, ξn,2} = {λn , λn } and hence

+ − |γn| = |λn − λn | = |ξn,1 − ξn,2|, ∀ n > ns,p. In view of estimate (3.11) of Lemma 3.7, we will need the following estimates on bn(λ) − q2n and b−n(λ) − q−2n for the proof of Theorem 3.1. Lemma 3.8 Let q ∈ Ls,p with −1/2 s 0 and 2 p < ∞. Then the 0,C 6 6 6 following holds:

s,p (i) For any f ∈ L and λ ∈ Sn with n ns,p and ns,p = ns,p(q) as in ?,C > Lemma 3.7

s hni |hTnf, e±ni| 6 Cs,pεs,p(n)kqks,pkfks,p;±n

where Cs,p is independent of q, n, and λ, and here  1 , 2 p < 2 ,  n 6 |s|  1−|s| ε (n) = (loghni) 2 s,p n , p = |s| ,   1 2 n1−|s|+2/p , |s| < p < ∞.

(ii) For any N > ns,p with ns,p = ns,p(q) as in Lemma 3.7   2p X sp p p kqks,p hni sup |bn(λ) − q2n| + sup |b−n(λ) − q−2n| 6 Cs,p p∗ , λ∈Sn λ∈Sn N n>N ∗ ∗ ∗ + where p > 1 is given by p = p − 1 if 2 6 p < 2/|s|, p = (p − 1) − 0 ∗ if p = 2/|s| and p = p/2 + 1 if 2/|s| < p < ∞, and the constant Cs,p is independent of q and N. o

Proof. (i) As the estimates of hTnf, eni and hTnf, e−ni can be proved in a similar −1 way we concentrate on hTnf, eni. Since by definition Tn = VAλ Qn,

X qn−mfm hT f, e i = . n n λ − m2π2 |m|6=n

Using that hn + mi/|n + m|, hn − mi/|n − m| 6 2 for |m| 6= n one gets s hni |hTnf, eni| 6 In(λ), where

X 1 1 |qn−m| |fm| In(λ) = 2 2 1−|s| |s| |s| |s| |m|6=n |n − m | hni hn − mi hn + mi

September 29, 2018 18 2 2 In the case p = 2, use that 1/|m − n | 6 1/n for all |m| 6= n to conclude by Cauchy–Schwarz that 1 I (λ) kqk kfk . n 6 n s,2 s,2;n In the case 2 < p < ∞, apply Hölder’s inequality for p/2 > 1 to get with p00 = p/(p − 2) > 1,

 1/p00 1 X 1 In(λ) 6  00  · n|s| 2 2 (1−|s|)p |m|6=n |m − n |  !p/2 !p/22/p X |qn−m| |fm| · .  |s| |s|  |m|6=n hn − mi hn + mi

By applying Cauchy–Schwarz to the latter sum one gets  1/p00 1 X 1 In(λ) 6  00  kqks,pkfks,p;n. n|s| 2 2 (1−|s|)p |m|6=n |m − n |

0 Using Lemma A.1 one sees that for some constant Cs,p > 0,  1/p00 1 X 1 0  00  6 Cs,pεs,p(n), n|s| 2 2 (1−|s|)p |m|6=n |m − n | yielding the claimed statement. (ii) Recall that for any n > ns,p and λ ∈ Sn

(n) (n) (n) bn − q2n = hTnu , eni, u ≡ u (λ) = Kn(λ)V e−n.

Using that

(n) ku ks,p;n 6 kKnks,p;nkV e−nks,;n 6 2kqks,p, it follows from item (i) that for any λ ∈ Sn

s 2 hni |bn − q2n| 6 2Cs,pεs,p(n)kqks,p.

Hence for any N > ns,p and after increasing Cs,p, if necessary,

X sp p 2p X p hni sup |bn(λ) − q2n| 6 Cs,pkqks,p (εs,p(n)) . λ∈Sn n>N n>N Taking into account that 2 6 p < ∞ and −1/2 6 s 6 0 one verifies that P p p∗ εs,p(n) = O(1/N ). In the same way, one shows that n>N X sp p 2p X p v hni sup |b−n(λ) − q−2n| 6 Cs,pkqks,p (εs,p(n)) . λ∈Sn n>N n>N

September 29, 2018 19 Proof of Theorem 3.1. By Lemma 3.7,

p p p/2 |γn| = |ξn,1 − ξn,2| 6 sup |6bn(λ)b−n(λ)| λ∈Sn p   3 p p 6 sup |bn(λ)| + sup |b−n(λ)| . 2 λ∈Sn λ∈Sn

It then follows from Lemma 3.8 that with N := ns,p

p  X sp p 6 X sp p p hni |γ | hni |q | + |q | n 6 2 2n −2n n>N n>N  p p + sup |bn(λ) − q2n| + sup |b−n(λ) − q−2n| λ∈Sn λ∈Sn 6p 6p 1 kR kp + C kqk2p . 6 2 2N s,p 2 s,p N s,p s,p In particular, (γn(q))n 1 ∈ ` ( ). As ns,p can be chosen locally uniformly in > C N s,p s,p s,p q ∈ L , the map L → ` ( ), q 7→ (γn(q))n 1 is locally bounded. v 0,C 0,C C N >

3.4 Jordan blocks of L(q)

To treat the Dirichlet problem we need to analyze L(q) further. We again consider potentials in the space Ls,p with −1/2 s 0 and 2 p < ∞. In case 0,C 6 6 6 q ∈ Ls,p is not real valued, the operator L(q) might have complex eigenvalues 0,C and it might happen that the geometric multiplicity of an eigenvalue is less than its algebraic multiplicity. For n > ns,p with ns,p as in Lemma 3.7, let  + − + − Null(L − λn ) ⊕ Null(L − λn ), λn 6= λn , En = + 2 + − Null(L − λn ) , λn = λn .

We need to estimate the coefficients of L(q) when represented with respect En + − to an appropriate orthonormal basis of En. In the case where λn = λn the matrix representation will be in Jordan normal form. By Lemma B.3, s+2,p 2 2 + 2 En ⊂ L ,→ L := L ([0, 2], ). Denote by f ∈ En an L -normalized ?,C C n + 2 eigenfunction corresponding to λn and by ϕn an L -normalized element in En + 2 so that {fn , ϕn} forms an L -orthonormal basis of En. Then the following lemma holds.

Lemma 3.9 Let q ∈ Ls,p with −1/2 s 0 and 2 p < ∞. Then there 0,C 6 6 6 0 0 exists ns,p > ns,p so that for any n > ns,p,

+ + (L − λn )ϕn = −γnϕn + ηnfn , where ηn ∈ C satisfies the estimate

+ + |ηn| 6 16(|γn| + |bn(λn )| + |b−n(λn )|). The threshold n0 can be chosen locally uniformly in q ∈ Ls,p . s,p 0,C o

September 29, 2018 20 + Proof. We begin by verifying the claimed formula for (L − λn )ϕn in the case + − − 2 where λn 6= λn . Let fn be an L -normalized eigenfunction corresponding to − − 2 2 λn . As fn ∈ En there exist a, b ∈ C with |a| + |b| = 1 and b 6= 0 so that 1 a f − = af + + bϕ or ϕ = f − − f +. n n n n b n b n Hence 1 a Lϕ = λ−f − − λ+f +. n b n n b n n − Substituting the expression for fn into the latter identity then leads to a (L − λ+)ϕ = (λ− − λ+)ϕ + (λ− − λ+)f + = −γ ϕ + η f + n n n n n b n n n n n n n + where ηn = −γna/b. In the case λn is a double eigenvalue of geometric multi- plicity two, ϕn is an eigenfunction of L and one has ηn = 0. Finally, in the case + + λn is a double eigenvalue of geometric multiplicity one, (L − λn )ϕn is in the + eigenspace En ⊂ En as claimed. + + To prove the claimed estimate for ηn, we view (L−λn )ϕn = −γnϕn+ηnfn as + a linear equation with inhomogeneous term g = −γnϕn+ηnfn . By identity (3.8) one has

+ BnPnϕn = γnPnKnϕn − ηnPnKnfn ,

+ + 2 where Kn ≡ Kn(λn ) and Bn ≡ Bn(λn ). To estimate ηn, take the L -inner + product of the latter identity with Pnfn to get

+ + ηnhPnKnfn ,Pnfn i = γnI − II, (3.14) where

+ + I = hPnKnϕn,Pnfn i,II = hBnPnϕn,Pnfn i.

+ + We begin by estimating hPnKnfn ,Pnfn i. Using that Kn = Id +TnKn one gets

+ + + 2 + + hPnKnfn ,Pnfn i = kPnfn kL2 + hTnKnfn ,Pnfn i, and by Cauchy-Schwarz  1/2 + + X + 2 + |hTnKnfn ,Pnfn i| 6 |hTnKnfn , emi| kPnfn kL2 . m∈{±n} + + Note that kPnfn kL2 6 kfn kL2 = 1 and further by Lemma 3.8 (i) + |s| + |hTnKnfn , e±ni| 6 hni Cs,pεs,p(n)kqks,pkKnfn ks,p;±n. 2 s,p + + By Corollary 3.2, kKnk 2 and as L [0, 2] ,→ L , kf k kf k = s,p;n 6 C ?,C n s,p;±n 6 n 0,2;±n |s| 1/p 1, whereas by the definition of εs,p(n) in Lemma 3.8 (i), n εs,p(n) 6 Cs,p/n . 0 Hence there exists ns,p > ns,p so that 1 |hT K f +,P f +i| . n n n n n 6 8

September 29, 2018 21 0 + By increasing ns,p if necessary, Lemma 3.10 below assures that kPnfn kL2 > 1/2. Thus the left hand side of (3.14) can be estimated as follows   + + 1 1 1 0 ηnhPnKnf ,Pnf i |ηn| − = |ηn|, ∀ n n . (3.15) n n > 4 8 8 > s,p

+ Next let us estimate the term I = hPnKnϕn,Pnfn i in (3.14). Using again Kn = Id +TnKn one sees that

+ + I = hPnϕn,Pnfn i + hTnKnϕn,Pnfn i.

+ + Clearly, |hPnϕn,Pnfn i| 6 kϕnkL2 kfn kL2 6 1 and arguing as above for the second term, one then concludes that

0 |I| 6 1 + 1/8, ∀ n > ns,p. (3.16)

+ Finally it remains to estimate II = hBnPnϕn,Pnfn i. Using again kϕnkL2 = + kfn kL2 = 1, we conclude from the matrix representation (3.9) of Bn that

+ + + 2 2 |hBnPnϕn,Pnfn i| 6 kBnkkϕnkL2 kfn kL2 6 |λn − n π − an| + |bn| + |b−n|.

+ Since det Bn(λn ) = 0, one has 1 |λ+ − n2π2 − a | = |b b |1/2 (|b | + |b |), n n n −n 6 2 n −n 0 and hence it follows that for all n > ns,p that

|II| 6 2(|bn| + |b−n|). (3.17)

Combining (3.15)-(3.17) leads to the claimed estimate for ηn. v

It remains to prove the estimate of Pn used in the proof of Lemma 3.9. To 2 this end, we introduce the Riesz projector Pn,q : L → En given by (see also AppendixB) Z 1 −1 Pn,q = (λ − L(q)) dλ. 2πi |λ−n2π2|=n Lemma 3.10 Let q ∈ Ls,p with −1/2 s 0 and 2 p < ∞. Then there 0,C 6 6 6 s+2,p exists n˜s,p ns,p so that for any eigenfunction f ∈ L of L(q) corresponding > ?,C to an eigenvalue λ ∈ Sn with n > n˜s,p, 1 kP fk kfk . n L2 > 2 L2 o Proof. In Lemma B.4 we show that

kPn,q − PnkL2→L∞ = o(1), s,p locally uniformly in q ∈ L and uniformly as n → ∞. Clearly, Pn,qf = f, 0,C hence

 v kPnfkL2 > kPn,qfkL2 − k(Pn,q − Pn)fkL2 > 1 + o(1) kfkL2 .

September 29, 2018 22 3.5 Proof of Theorem 3.2

We begin with a brief outline of the proof of Theorem 3.2. Let q ∈ Ls,p with 0,C −1/2 s 0 and 2 p < ∞. Since for any q ∈ H−1 the Dirichlet eigenvalues, 6 6 6 0,C when listed in lexicographical ordering and with their algebraic multiplicities, 2 2 2 µ1 4 µ2 4 ··· , satisfy the asymptotics µn = n π + n`n, they are simple for −1 n ndir, where ndir 1 can be chosen locally uniformly for q ∈ H . For any > > 0,C 2 n > ndir let gn be an L -normalized eigenfunction corresponding to µn. Then g ∈ H1 := {g ∈ H1([0, 1], ): g(0) = g(1) = 0}. n dir,C C Now let q ∈ Ls,p with −1/2 s 0 and 2 p < ∞. Increase n0 of 0,C 6 6 6 s,p 0 Lemma 3.9, if necessary, so that ns,p > ndir and denote by En the two dimen- sional subspace introduced in Section 3.4. We will choose an L2-normalized func- tion G˜ in E so that its restriction G to I = [0, 1] is in H1 and close to g . n n n dir,C n + + We then show that µn − λn can be estimated in terms of h(Ldir − λn )Gn,GniI . As by Lemma 3.9

+ ˜ + + (L − λn )Gn = O(|γn| + |bn(λn )| + |b−n(λn )|),

+ the claimed estimates for µn−τn = µn−λn +γn/2 then follow from the estimates of γn of Theorem 3.1 and the ones of bn − q2n, b−n − q−2n of Lemma 3.8 (ii). ˜ + 2 The function Gn is defined as follows. Let fn , ϕn be the L -orthonormal basis of E chosen in Section 3.4. As E ⊂ H1( /2 ), its elements are con- n n C R Z + tinuous functions by the Sobolev embedding theorem. If fn (0) = 0, then + + fn (1) = 0 as fn is an eigenfunction of the 1-periodic/antiperiodic eigenvalue + ˜ + + ˜ λn of L(q) and we set Gn = fn . If fn (0) 6= 0, then we define Gn(x) = + +  rn ϕn(0)fn (x) − fn (0)ϕn(x) , where rn > 0 is chosen in such a way that R 1 ˜ 2 ˜ ˜ ˜ 0 |Gn(x)| dx = 1. Then Gn(0) = Gn(1) = 0 and since Gn is an element of En its restriction G := G˜ is in H1 . n n I dir,C Denote by Πn,q the Riesz projection, introduced in AppendixB, Z 1 −1 Πn,q := (λ − Ldir(q)) dλ. 2πi |λ−n2π2|=n

It has span(gn) as its range and hence Πn,qGn = νngn for some νn ∈ C.

Lemma 3.11 Let q ∈ Ls,p with −1/2 s 0 and 2 p < ∞. Then there 0,C 6 6 6 00 0 0 00 exists ns,p > ns,p with ns,p as in Lemma 3.9 so that for any n > ns,p (i) the following identity holds

+ +  νn(µn − λn )gn = βn ηnΠn,q(fn I ) − γnΠn,q(ϕn I ) , (3.18)

where βn ∈ C with |βn| 6 1 and ηn is the off-diagonal coefficient in the + + matrix representation of (L − λ ) with respect to the basis {f , ϕn}, n En n introduced in Lemma 3.9, and

September 29, 2018 23 (ii) 1/2 6 |νn| 6 3/2. n00 can be chosen locally uniformly for q ∈ Ls,p . s,p 0,C o

Proof. (i) Write Gn = νngn + hn, where hn = (Id −Πn,q)Gn. Then + + + (Ldir − λn )Gn = νn(µn − λn )gn + (Ldir − λn )hn. ˜ ˜ ˜ + On the other hand, Gn = Gn I , where Gn ∈ En is given by Gn = αnfn + βnϕn with α , β ∈ satisfying |α |2 + |β |2 = 1 and G ∈ H1 . Hence by n n C n n n dir,C 0 Lemma B.1 and Lemma 3.9, for n > ns,p, + + ˜ + (Ldir − λn )Gn = (L − λn )Gn I = βn(ηnfn − γnϕn) I .

Combining the two identities and using that Πn,qhn = 0 and that Πn,q commutes + with (Ldir − λn ), one obtains, after projecting onto span(gn), identity (3.18). (ii) Taking the inner product of Πn,qGn = νngn with gn one gets

νn = νnhgn, gniI = hΠn,qGn, gniI . √ Let sn(x) = 2 sin(nπx) and denote by Πn = Πn,0 the orthogonal projection 2 onto span{sn}. Recall that Pn,q : L → En is the Riesz projection onto En. In Lemma B.4 we show that

kΠn,q − ΠnkL2(I)→L∞(I), kPn,q − PnkL2→L∞ = o(1), (3.19) s,p locally uniformly in q ∈ L and uniformly as n → ∞. Thus using ΠnGn = 0,C ˜ 2 2 Πn(PnGn) I and recalling that kGnkL2(I) = kgnkL2(I) = 1 we obtain ˜ vn = hΠnGn, gniI + h(Πn,q − Πn)Gn, gniI = hPnGn, ΠngniI + o(1). Moreover, it follows from (3.19) that uniformly in 0 6 x 6 1

iφn Πngn(x) = e sn(x) + o(1), n → ∞, with some real φn. Similarly, again by (3.19), uniformly in 0 6 x 6 2

PnG˜n(x) = anen(x) + bne−n(x) + o(1), n → ∞, where, since kGnkL2 = 1 and Gn(0) = 0, the coefficients an and bn can be chosen so that 2 2 |an| + |bn| = 1, an + bn = 0.

iψn That is PnG˜n(x) = e sn(x) + o(1) with some real ψn and hence

˜ iψn−iφn iψn−iφn hPnGn, ΠngniI = e hsn, sniI + o(1) = e + o(1), n → ∞. From this we conclude

|νn| = 1 + o(1), n → ∞.

00 00 0 Therefore, 1/2 6 |νn| 6 3/2 for all n > ns,p provided ns,p > ns,p is sufficiently large. 00 Going through the arguments of the proof one verifies that ns,p can be chosen locally uniformly in q. v

September 29, 2018 24 Lemma 3.11 allows to complete the proof of Theorem 3.2.

Proof of Theorem 3.2. Take the inner product of (3.18) with gn and use that |νn| > 1/2 by Lemma 3.11 to conclude that

1 + +  |µn − λ | |βn| |ηn|hΠn,q(f ), gni + |γn||hΠn,q(ϕn )i , gn| . (3.20) 2 n 6 n I I I I Recall that |β | 1 and note that for any f, g ∈ L2 (I) n 6 C

|hΠn,qf, giI | 6 |hΠnf, giI | + |h(Πn,q − Πn)f, giI | 6 (1 + o(1))kfkL2(I)kgkL2(I).

+ Since kfn kL2(I) = kϕnkL2(I) = 1 and kgnkL2(I) = 1,(3.20) implies that

+ |µn − λn | 6 (2 + o(1))(|ηn| + |γn|) yielding with Lemma 3.9 the estimate

+ + |µn − τn| 6 (3 + o(1))|γn| + (32 + o(1))(|γn| + |bn(λn )| + |b−n(λn )|).

s,p By Theorem 3.1 and Lemma 3.8 (ii) it then follows that (τn − µn)n 1 ∈ ` ( ). > C N Going through the arguments of the proof one verifies that the map Ls,p → 0,C s,p ` ( ), q 7→ (τn − µn)n 1 is locally bounded. v C N >

4 Proof of Theorem 1.1

The proof of Theorem 1.1 relies on a formula for the nonlinear part of the KdV Hamiltonian which we derive in Subsection 4.2. The formula is based on the integral

Z λ ∆• √ F (λ) = c dλ, + 2 λ0 ∆ − 4 whose properties are studied in Subsection 4.1. The proof of Theorem 1.1 is then completed in Subsection 4.3. In the sequel we use without further reference the notation introduced in Section1 and Section2.

4.1 Properties of F

Let q be in W. It is convenient to define the standard root q s + − ςn(λ) = (λn − λ)(λn − λ), λ ∈ C \ Gn, n > 1, by the condition

p+ 2 2 − + ςn(λ) = (τn − λ) 1 − γn/4(τn − λ) , τn = (λn + λn )/2. (4.1) √ Here + denotes the principal branch of the square root on the complex plane minus the ray (−∞, 0]. The standard root is analytic in λ on C\Gn and in (λ, p)

September 29, 2018 25 pc 2 on (C \ Un) × Wq. The canonical root ∆ (λ) − 4, introduced in Section2, can then be written in terms of standard roots as follows q + Y ςm(λ) pc ∆2(λ) − 4 := −2i λ − λ+ (4.2) 0 m2π2 m>1 S S and is analytic in λ on \ Gn and in (λ, p) on ( \ Un) × Wq. C γn6=0 C n>0 A path in the complex plane is said to be admissible for q if, except possibly at its endpoints, it does not intersect any non collapsed gap Gn(q).

Lemma 4.1 For each q ∈ W the following holds:

∆•(λ) S (i) √ is analytic in (λ, p) on ( \ Un) × Wq and analytic in λ c ∆2(λ)−4 C n>0 S on C \ γn6=0 Gn. (We recall that γ0 = ∞.) n>0

λ+ • − + R n √∆ (λ) (ii) For any n 1 and any admissible path from λn to λn in Un, − c dλ = > λn ∆2(λ)−4 R ∆•(λ) 0. As a consequence, for any closed circuit Γn in Un around Gn, √ dλ = Γn c ∆2(λ)−4 0.

∆•(λ) 1 (iii) For q = 0, √ = +√ . c ∆2(λ)−4 2i λ o

Proof. Standard – see e.g. [19, Lemma 1]. v

S Next we define for any q ∈ W and λ ∈ C \ γn6=0 Gn the improper integral n>0

Z λ ∆•(z) ∆•(λ) 1 1 Y λ• − λ F (λ) := dz, = m , pc 2 pc 2 q λ+ ∆ (z) − 4 ∆ (λ) − 4 2i + + ςm(λ) 0 λ − λ0 m>1 (4.3) computed along an arbitrary admissible path. The improper integral F (λ) exists q + + as in the product representation (4.3) the factor 1/ λ − λ0 is integrable on C \ G0. Furthermore, by Lemma 4.1 it is independent of the chosen admissible ± path and hence well defined. Moreover, F (λ) continuously extends to Gn for + − any n > 0, where Gn [Gn ] is the left [right] hand side of Gn. Lemma 4.2 For any q ∈ W the following holds:

S (i) F is analytic in (λ, p) on (C \ n 0 Un) × Wq and F (λ) ≡ F (λ, q) is S > analytic in λ on C \ γn6=0 Gn. (We recall that γ0 = ∞.) n>0 + + − (ii) F (λ0 ) = 0 and F (λn ) = F (λn ) = −inπ for any n > 1.

(iii) Locally uniformly on Wq

sup |F (λ) + inπ| = O(γn/n), n → ∞. + − λ∈Gn ∪Gn

September 29, 2018 26 √ + (iv) For q = 0, F (λ) is analytic on C \ (−∞, 0] and given by F (λ) = −i λ. Proof. Items (i), (ii), and (iv) are well known. See e.g. [19, Lemma 2] for a proof. R λ ∆•(λ) (iii) To obtain the claimed estimate for F (λ) + inπ = − √ dλ on λn c ∆2(λ)−4 + − ± Gn ∪ Gn it suffices to consider p ∈ Wq with γn(p) 6= 0. Parameterizing Gn by ± λs = τn + (s ± i0)γn/2, −1 6 s 6 1, one gets for −1 6 t 6 1 λ± t Z t ∆• Z λ• − τ − sγ /2 ± √ n √n n F (λt ) + inπ = c dz = ∓ + χn(λs) ds, − 2 2 λn ∆ − 4 −1 i 1 − s where

• 1 1 Y λm − λ χn(λ) := q (4.4) 2i + + ςm(λ) λ − λ0 16m6=n

The map χn is analytic on Un × Wq and nχn|Gn = O(1) as n → ∞ locally uniformly on Wq – see [11, Lemma 9.7]. Moreover, by [11, Proposition 2.18] • one can choose n0 > 1 locally uniformly on Wq such that |λn − τn| 6 |γn| for n > n0. Altogether this shows

sup |F (λ) + inπ| = O(|γn|/n) ± λ∈Gn locally uniformly on Wq. v

Occasionally we write for n > 0 Z λ ∆• : √ Fn(λ) = F (λ) + inπ = c dz (4.5) + 2 λn ∆ − 4

• to denote the primitive of √ ∆ normalized by the condition F (λ+) = 0. c ∆2−4 n n 2 Lemma 4.3 For any q ∈ W and any n > 0, Fn is analytic on Un and hence S on C \ n6=m>0 Gm. γm6=0 2 Proof. Let q ∈ W and n > 0 be arbitrary. It remains to show that Fn admits an analytic extension from Un \ Gn to all of Un. In view of (4.3) we can write

• • ∆ (λ) 1 1 Y λm − λ = χ0(λ), χ0(λ) = , pc 2 q ∆ (λ) − 4 + + 2i ςm(λ) λ − λ0 m>1 and for any k > 1, • • • ∆ (λ) λk − λ 1 1 Y λm − λ = χk(λ), χk(λ) = . pc 2 q ∆ (λ) − 4 ςk(λ) 2i + + ςm(λ) λ − λ0 16m6=k q + + For any k > 0, χk is analytic on Uk. Moreover, the root λ − λ0 [ςk] admits opposite signs on opposite sides of G0 [Gk]. For any given λ ∈ Gn denote by ± ± λ the corresponding elements on Gn . Thus for all λ ∈ Gn + − Fn(λ ) = −Fn(λ ).

September 29, 2018 27 2 Consequently, Fn is continuous on all of Un. Since it is already analytic on Un \ Gn, it is analytic on Un. v

4.2 Asymptotics of F R 1 Assume that q is a real valued finite gap potential with 0 q(x) dx = 0. Then q is real analytic and there exists A ⊂ N finite so that γn > 0 if and only if n ∈ A. + − + 2 2 Hence there exists n0 > 1 so that λn = λn and |λn − π n | 6 1 for any n > n0 + ± 2 2 and |λ0 |, |λn | 6 n0π + 1 for any n 6 n0. By Lemma 4.3 it then follows that 2 4 2 2 F (λ) and hence F (λ) are analytic functions in the domain |λ| > n0π + 1. In particular, F 2(λ) admits a Laurent expansion,

2 X k 2 2 F (λ) = akλ , |λ| > n0π + 1. (4.6) k∈Z From the product representation (4.3) one sees that F 2(λ) = O(λ) and hence ak = 0 for any k > 2. It is well known that the coefficients ak can be explicitly computed – see e.g. [19] for a self-contained proof. R 1 Lemma 4.4 For any real valued finite gap potential q with 0 q(x) dx = 0, kdv   Z 1 4 2 H 1 1 1 2 F (λ) = λ − H0 − + O 2 , H0 = q(x) dx. o (4.7) 4 λ λ 2 0 As an application of the asymptotics (4.7) of F 4(λ) we derive a formula for the KdV Hamiltonian Hkdv for finite gap potentials. To this end we introduce for any q ∈ W and n > 1 Z 1 3 Rn ≡ Rn(q) = Fn dλ, Fn(λ) = F (λ) + inπ. (4.8) π Γn S According to Lemma 4.2 (i), Fn(λ) is analytic on C \ γm6=0 Gm, hence Rn m>0 vanishes if γn = 0. Proposition 4.5 For any real finite gap potential q,

kdv X 3 3 X H (q) = 8n π In + 8nπRn , (4.9) n∈A n∈A where A ≡ Aq = {n > 1 : γn 6= 0} is finite. o

4 S Proof. By Lemma 4.3, F (λ) is analytic on C \ n∈A Gn, hence in view of Lemma 4.4 and Cauchy’s integral formula 1 Z Hkdv F 4(λ) dλ = − , (4.10) 2πi Cr 4 where Cr denotes a counter clockwise oriented circle around the origin enclosing all open gaps. By contour deformation we thus obtain 1 Z X 1 Z Hkdv = −4 F 4(λ) dλ = −4 F 4(λ) dλ. (4.11) 2πi 2πi Cr n∈A Γn

September 29, 2018 28 2 Moreover, since Fn is analytic on Un, one concludes that in the expression 4 4 F (λ) = (Fn(λ) − inπ) , the even powers of Fn will not contribute to the integral R F 4(λ) dλ. As a consequence Γn Z Z 1 4 1 3 3  F (λ) dλ = − 4(inπ)Fn (λ) + 4(inπ) Fn(λ) dλ. 2πi Γn 2πi Γn

After integration by parts, formula (2.6) of the action In becomes 1 Z 1 Z In = − F (λ) dλ = − Fn(λ) dλ. π Γn π Γn One then concludes with (4.8) that Z 1 4 3 3 −4 F (λ) dλ = (8n π )In + (8nπ)Rn (4.12) 2πi Γn yielding in view of (4.11) the claimed formula. v

4.3 Nonlinear part of Hkdv

Recall that Rn, n > 1, has been defined in (4.8) as follows Z 1 3 Rn = Fn dλ, Fn(λ) = F (λ) + inπ. π Γn

When expressed in Birkhoff coordinates, we denote Rn by Rn. The following lemma describes the properties of Rn and Rn used in the sequel.

Lemma 4.6 (i) For any n > 1, Rn is analytic on W and satisfies Rn = 4 3 O(γn/n ) locally uniformly on W. Furthermore, if γn = 0, then Rn = 0, and for q ∈ W real valued, Rn 6 0. In particular, for q = 0, Rn = 0 for all n > 1.

(ii) For any n > 1, Rn is an analytic function of the actions alone on a neigh- borhood V (independent of n) of the positive cone `−1,1( ) in `−1,1( ). If + N C N −1,1 In = 0, then Rn = 0, and on `+ (N) we have Rn 6 0. o

Proof. (i) By Lemma 4.2 (i), Rn is analytic on W, and by Lemma 4.2 (iii),

Fn = O(γn/n), ± Gn 4 3 locally uniformly in q, hence Rn = O(γn/n ) locally uniformly on W. If γn = 0, then F (λ) is analytic on Un by Lemma 4.2 (i) hence Rn = 0. If q is real valued, − + then one has for λn < λ < λn Z λ ∆• Z λ (−1)n∆• (−1)n∆(λ) √ √ Fn(λ) = c dλ = ± + dλ = ± arccosh , ± + 2 ± + 2 2 Gn λn ∆ − 4 Gn λn ∆ − 4 and consequently by shrinking the contour Γn in the definition of Rn to Gn λ+ 3 2 Z n  (−1)n∆(λ) Rn(q) = − arccosh dλ 6 0. π − 2 λn

September 29, 2018 29 (ii) Since Fn and hence Rn is defined in terms of the discriminant, Rn depends only on the periodic spectrum of q. Therefore, when expressed as a function of the Birkhoff coordinates, it is a function of the actions alone. Arguing as in the proof of [6, Theorem 20.3] it follows that each Rn, n > 1, is an analytic −1,1 function on a neighborhood V (independent of n) of the positive cone `+ (N) −1,1 −1,1 in ` ( ). Since every sequence (In)n 1 ∈ ` ( ) is the sequence of actions C N > + N −1 −1,1 v for some real valued potential q ∈ H0 , we conclude that Rn 6 0 on `+ .

Lemma 4.6 can be applied to prove the following result which is the key ingredient in the proof of Theorem 1.1. P Proposition 4.7 (i) The series 8nπRn converges absolutely locally n>1 uniformly on W ∩ L−1/2,4 and hence is analytic. 0,C (ii) The series

∗ X H := 8nπRn n>1

is an analytic function of the action variables defined on a complex neigh- borhood of the positive quadrant `2 ( ) in `2 ( ). On `2 ( ) it is non- + N C N + N positive. o

Proof. (i) By Lemma 4.6 (i), every Rn with n > 1 is analytic on W. As L−1/2,4 ,→ H−1, its restriction to W ∩ L−1/2,4 is analytic as well. In addition, 0,C 0,C 0,C 4 2 also by Lemma 4.6, 8nπRn = O(γn/n ) locally uniformly on W, and by The- 1/2 4 −1/2,4 orem 3.1, γn = n ` locally uniformly on W ∩ L . It then follows that n 0,C P −1/2,4 8nπRn is analytic on W ∩ L – cf. [6, Theorem A.4]. n>1 0,C (ii) By Lemma 4.6 (ii), each Rn, n > 1, is analytic on V. According to (2.13) the action variable In is related to the gap length γn as γ2 n = O(I ) n n 2 locally uniformly on V. Thus, 8nπRn = O(In) locally uniformly on V. Clearly, the restriction of each R , n 1, to V = V ∩ `2 is analytic, too. Therefore, n > 2 C P 8nπRn converges absolutely and locally uniformly to an analytic function n>1 ∗ 2 ∗ on V2, which we denote by H . On `+(N) the function H is real valued with values in (−∞, 0]. v

From the property that for q ∈ W ∩ H1 , I = n3`1 , Proposition 4.7 leads 0,C n n to the following

Corollary 4.8 On W ∩ H1 0,C

kdv X 3 3 X H = 8n π In + 8nπRn (4.13) n>1 n>1

September 29, 2018 30 where both sums converge absolutely locally uniformly on W∩H1 . In addition, 0,C kdv ∗ 3 3 3,1 expressing (4.13) in Birkhoff coordinates gives ωn = ∂In H +8n π on `+ (N). o Proof. Let us first consider the case where the potential is real valued. Recall 1 1 that real valued finite gap potentials are dense in H0 and that for any q ∈ H0 and + − n > 1 with λn = λn one has In = 0 and Rn = 0. The claimed formula (4.13) for 1 q ∈ H0 then follows from (4.9) by approximating q by a sequence of finite gap potentials. As both sums in (4.13) are real analytic on W ∩ H1 , the identity 0,C holds on W ∩ H1 . v 0,C The proof of Theorem 1.1 now easily follows from Proposition 4.7. Proof of Theorem 1.1. By Proposition 4.7, H∗ is an analytic function on a com- plex neighborhood V ⊂ `2 ( ) of the positive quadrant `2 ( ). Its differential C N + N ω∗ := d H∗ is an analytic map with values in `2 ( ), I C N ω∗ : V → `2 ( ),I = (I ) 7→ (ω∗ = ∂ H∗) . C N n n>1 n In n>1 ∗ kdv 3 3 In view of Corollary 4.8 and (1.7) one has ωn = ωn − 8n π on V0 ∩ V for kdv 3,1 any n > 1. The expansion (1.4) of H : `+ (N) → R at I = 0 implies that H∗ : `2 ( ) → has an expansion at I = 0 of the form −3 P I2 + ··· where + N R n>1 n ∗ the dots stand for higher order terms in I. In particular, ω (I) = (−6In +

··· )n>1 vanishes at I = 0 and its differential at I = 0 is given by ∗ d0ω = −6 Id`2 ( ) . C N It implies that ω∗ defines a local diffeomorphism near I = 0 and that the Hessian 2 ∗ v dI H satisfies the estimate (1.8).

5 Proof of Theorem 1.2

In [11] and [13], the restrictions of the Birkhoff map

−1 −1/2,2 Φ: H0 → `0 , q 7→ (zn(q))n∈Z, z0(q) = 0, s to the Sobolev spaces H0 , −1 6 s 6 0, are studied. It turns out that the argu- ments developed in these two papers can also be adapted to prove Theorem 1.2. As a first step we prove the following version of [11, Theorem 5.3].

Lemma 5.1 For any −1/2 6 s 6 0 and 2 6 p < ∞

s,p s+1/2,p Φs,p ≡ Φ : L0 → `0 , q 7→ (zn(q))n∈Z, s,p L0 s,p is real analytic and extends analytically to an open neighborhood Ws,p of L0 s,p in L . Its Jacobian d0Φs,p at q = 0 is the weighted Fourier transform 0,C ! s,p s+1/2,p 1 d0Φs,p : L → ` , f 7→ hf, e2ni 0 0 p2π max(|n|, 1) I n∈Z

September 29, 2018 31 with inverse given by

−1 s+1/2,p s,p X p (d0Φs,p) : `0 → L0 , (wn)n∈Z 7→ 2π|n|wne2n. n∈Z

In particular, Φs,p is a local diffeomorphism at q = 0. o

s,p −1 Proof. Since L ,→ H , the coordinate functions zn defined in (2.11) are ana- 0,C 0,C lytic functions on the complex neighborhood W ⊂ H−1 of H−1 of Theorem 2.3. 0,C 0 Furthermore, |γ (q)| + |µ (q) − τ (q)| z (q) = O n √n n ±n n locally uniformly on W and uniformly in n > 1. By the asymptotics of the periodic and Dirichlet eigenvalues of Theorems 3.1- 3.2, Φs,p maps the complex s,p s,p s+1/2,p neighborhood Ws,p := W ∩ L of L into the space ` and is locally 0,C 0 0,C bounded. Arguing as in the proof of [12, Theorem 8.5] one sees that Φs,p is analytic. The formulas for d0Φs,p and its inverse follow from [12, Theorem 9.7] by continuity. v

As a next step we show that Φs,p is a local diffeomorphism.

s,p Lemma 5.2 For any q ∈ L0 with −1/2 6 s 6 0 and 2 6 p < ∞,

s,p s+1/2,p dqΦs,p : L0 → `0 is a linear isomorphism.

Proof. We follow the line of arguments used in [11, Proposition 7.1]. By The- −1 −1/2,2 −1 orem 2.3, dqΦ: H0 → `0 is a linear isomorphism for any q ∈ H0 . Fur- s,p s+1/2,p thermore, by Lemma 5.1, d0Φs,p : L0 → `0 is a linear isomorphism. −1 −1 −1 It is convenient to introduce the maps Ω := (d0Φ) ◦ Φ: H0 → H0 and −1 s,p s,p s,p Ωs,p := (d0Φs,p) ◦ Φs,p : L0 → L0 . It is to prove that, for any q ∈ L0 , −1 −1 dqΩs,p is a linear isomorphism. By Theorem 2.2, dqΩ: H0 → H0 is a linear isomorphism and as dqΩs,p = dqΩ s,p it then follows that dqΩs,p is one-to-one. L0 To show that dqΩs,p is onto, we prove that dqΩs,p is a Fredholm operator of 2 index zero. To this end, denote by ∂qzn the L -gradient of zn : W → C. Then 0 1 |s|,p ∼ s,p 0 0 ∂qzn ∈ H and hence ∂qzn ∈ L = (L ) where p is conjugate to p, 0,C 0,C 0,C 1 1 s,p p + p0 = 1. By Lemma 5.1, for any q ∈ L0 , dqΩs,p has the form

X s,p dqΩs,p(f) = hf, dn(q)iI e2n, ∀ f ∈ L0 , n6=0

p s,p where dn(q) = 2|n|π∂qz−n. We claim that for any q ∈ L0 , T (q) = (dqΩs,p)− s,p Id is a on L0 . Since q 7→ T (q) is analytic and by [12, 0 s,p Theorem B.16] the finite gap potentials are dense in H0 and hence in L0 , it suffices to show that T (q) is compact for any finite gap potential. Recall that

September 29, 2018 32 −1 q ∈ H0 is a finite gap potential if Aq = {n > 1 : γn(q) 6= 0} is finite. Using s,p P that for any f ∈ L0 , f = n6=0hf, e2niI e2n, one gets X T f = hf, dn − e2niI e2n. n6=0

s,p s,p For any N > 1 define the operator TN : L0 → L0 given by X TN f = hf, dn − e2niI e2n. 06=|n|6N

Then TN is a linear operator of finite rank and thus compact. We now show s,p that kT − TN ks,p → 0 in the operator norm as N → ∞. For any f ∈ L0 ,

|hf, dn − e2niI | 6 kfks,pkdn − e2nk|s|,p0 and hence

 1/p p X |hf, dn − e2ni | k(T − T )fk = I N s,p  |s|p  |n|>N |n| 1/p  p  X kdn − e2nk|s|,p0 kfk . 6 s,p |s|p  |n|>N |n|

p 1/p  kdn−e2nk  Since by Lemma 5.3 below P |s|,p0 → 0 as N → ∞, it then |n|>N |n||s|p follows that kT − TN ks,p → 0 in the operator norm as N → ∞. Hence T (q) is compact for any finite gap potential which proves the claim. v

The following lemma, used to prove Lemma 5.2, extends Lemma 7.3 in [11].

Lemma 5.3 For any real valued finite gap potential q with q0 = 0, p X kdn(q) − e2nk|s|,p0 < ∞ |s|p |n|6=0 |n| for any −1/2 6 s 6 0 and any 2 6 p < ∞.

Proof. The case (s, p) = (0, 2) is contained in [12, Lemma 11.12] so it remains to consider the case (s, p) 6= (0, 2). We adapt the proof of [11, Lemma 7.3] to the case at hand. Let q be a finite gap potential with q0 = 0 and introduce + N = {n > 1 : γn = 0; λn 6= 0}. Then N \N is finite and it suffices to show that p X kdn − e2nk|s|,p0 < ∞ |s|p |n|∈N |n| where dn = dn(q). By [12, Theorem 9.7] one has for n ∈ N ,

p + p − dn = 2|n|π∂qz−n = cn∂qzn , d−n = 2|n|π∂qzn = cn∂qzn .

September 29, 2018 33 Here cn is a complex number independent of x and

± 2 2 ∂qzn = hn − gn ± 2ignhn, where gn denotes the real valued eigenfunction of the µn = + 2 R 1 2 0 λn of −∂x + q, normalized by 0 gn dx = 1, gn(0) > 0, and hn denotes the real + − 2 valued eigenfunction of the periodic eigenvalue λn = λn which is L -orthogonal R 1 2 to gn and satisfies 0 hn dx = 1 and hn(0) > 0. Following the arguments of 1 3 0 the proof of Lemma 7.3 in [11] consider the operator Q = − 2 ∂x + 2q∂x + q . 00 One verifies that the product Y = yiyj of any two solutions of −y + qy = λy R 1 is smooth and satisfies QY = 2λ∂xY . Since 0 ∂qzn dx = 0 for any n ∈ N 00 + and ∂qzn is a linear combination of products of solutions of −y + qy = λn y it follows that for any n ∈ N

1 −1 d±n = + D0 Qd±n. 2λn + −1 Recall that λn 6= 0 for n ∈ N . Here D0 denotes the inverse of the restriction s,p s−1,p s,p 2+s,p D0 : L0 → L0 of ∂x to L0 . For f ∈ L0 one has for any n ∈ N

1 −1 1 −1 1 ∗ hd±n, fiI = + hD0 Qd±n, fiI = + hd±n, QD0 fiI = − + hd±n, f iI 2λn 2λn 4λn ∗ 00 0 −1 s,p where f := f − 4qf − 2q D0 f ∈ L0 . This yields

00 ∗ ∗ 00 hd±n, f iI = hd±n, f iI − hd±n, f − f iI + 0 −1 = −4λn hd±n, fiI − hd±n, −4qf − 2q D0 fiI .

00 2 2 Combining this with he±2n, f iI = −4n π he±2n, fiI one obtains for n ∈ N

00 2 2 + + hd±n − e±2n, f iI = 4(n π − λn )he±2n, fiI − 4λn hd±n − e±2n, fiI 0 −1 + he±2n, 4qf + 2q D0 fiI 0 −1 + hd±n − e±2n, 4qf + 2q D0 fiI . (5.1) We now estimate the terms on the right hand side of the latter identity sepa- 2+s,p 1 rately. By standard estimates one has for f ∈ L0 ,→ H0  1  he±2n, fiI = O kfkH1 , n 0  1  λ+ = n2π2 + O , n n   0 −1 1 he±2n, 4qf + 2q D0 fiI = O kfkH1 . n 0

0 −1 1 1 For the latter one we used that qf and q D0 f are in H0 since H0 is an algebra. Finally, by [12, Lemma 11.12] one has for any n ∈ N log n  hd±n − e±2n, fiI = O kfkH1 , n5/2 0

September 29, 2018 34 0 −1 1 and, using again that qf and q D0 f are in H0 ,   0 −1 log n hd±n − e±2n, 4qf + 2q D0 fiI = O kfkH1 . n5/2 0

s,p 00 P 2 2 Note that any element in L0 is of the form f = n6=0 4n π f2ne2n where P 2+s,p f = n6=0 f2ne2n is in L0 and

2 2 p p 00 p X (4π n ) |f2n| 2p X (2−|s|)p p 2p p kf k s,p (2π) hni |f | (2π) kfk L 6 |s|p 6 2n 6 2−|s| 0 L0 ,p n6=0 hni n6=0 Hence, combining the above estimates we conclude from (5.1) that loghni kd − e k sup |hd − e , f 00i | C ±n ±2n |s|,p0 6 ±n ±2n I 6 n1/2 kfk2+s,p61 and, in turn,

p p X kdn − e2nk|s|,p0 X (loghni) C < ∞ |s|p 6 (|s|+1/2)p |n|∈N |n| |n|∈N |n| as (|s| + 1/2)p > 1 for any −1/2 6 s 6 0 and 2 6 p < ∞ with (s, p) 6= (0, 2). v

In a third step, following arguments used in [13], we prove that Φs,p is onto.

s,p s+1/2,p Lemma 5.4 For any −1/2 6 s 6 0 and 2 6 p < ∞, Φs,p : L0 → `0 is onto.

(0) (0) s+1/2,p Proof. Fix −1/2 6 s 6 0, 2 6 p < ∞, and z = (zn )n6=0 ∈ `0 . Since −1 −1/2,2 (0) −1 (0) (0) Φ: H0 → `0 is onto, there exits q ∈ H0 such that z = Φ(q ). Assume that

(0) −1 s,p q ∈ H0 \ L0 . (5.2)

We show that (5.2) leads to a contradiction. As Φs,p(0) = 0 and by Lemma 5.2 the differential of Φs,p at 0 is a linear isomorphism, by the inverse function theo- s,p rem, there exists an open neighborhood U of 0 in L0 and an open neighborhood s+1/2,p V of 0 in `0 so that

Φs,p U : U → V is a diffeomorphism. Without loss of generality we can assume that V is a ball of radius 2ε, centered at 0,

s+1/2,p V = B2ε(0) ⊂ `0 . (5.3)

To obtain the desired contradiction we construct a finite sequence q(1), . . . , q(N) (N) (k) (k−1) 0 with the property that q ∈ U and q − q ∈ H0 for any 1 6 k 6 N. To this end consider the action angle variables (Ik, θk), k > 1, constructed in [11].

September 29, 2018 35 −1 2 For any k > 1, Ik is defined on H0 and satisfies Ik = |zk| whereas the angle −1 variable θk is defined on H0 \ Zk with values in R/2πZ and is real analytic when considered mod π. Here Zk denotes the set Zk = {q ∈ W : γk(q) = 0}, 2 −1 introduced in (2.7). The L -gradient ∂qθk is a real analytic map from H0 \ Zk 1 with values in H0 whereas the corresponding Hamiltonian vector field, −1 0 Yk = ∂x∂qθk : H0 \ Zk → H0 0 is a real analytic map with values in H0 . It defines a dynamical system q˙ = Yk(q) −1 on H0 \ Zk. We now use the flows of these vector fields to construct the sequence q(1), . . . , q(N) recursively as follows. The potential q(0) is given by (n−1) (5.2). For any n > 1 assume that q has already been constructed. If 2s+1 (n−1) 2 2n/p (n) (n−1) 2s+1 (n−1) hni In(q ) < ε /2 , then set q ≡ q . If hni In(q ) > 2 2n/p (n−1) −1 ε /2 , then q ∈ H0 \ Zn and hence the vector field Yn is well defined in a neighborhood of q(n−1). The element q(n) is then chosen to be an element (n−1) on the solution curve of the vector field Yn passing through q so that 2s+1 (n) 2 2n/p (n) hni In(q ) < ε /2 . The existence of such a potential q follows from Lemma 5.5 (i) below. Moreover, by the commutator relations (cf. section2)

Yn(Ik) = {Ik, θn} = δnk, the vector field Yn preserves the value of the action variable Ik with k 6= n. In particular, we have 1 hki2s+1I (q(n)) < ε2 , 1 k n, k 22k/p 6 6 and

(n) (0) Ik(q ) = Ik(q ), k > n. One then obtains p/2 (n) p X (s+1/2)p (n) p X 2s+1 (n)  kΦ(q )ks+1/2,p = hki |zk(q )| = 2 hki Ik(q ) k6=0 k>1 p/2 X 1 X 2s+1  2εp + 2 hki I (q(0)) . 6 2k k 16k6n k>n

(0) s+1/2,p P 2k+1 (0) p/2 Since Φ(q ) ∈ ` , one has (hki Ik(q )) < ∞. Choose N 1 0 k>1 > so that p/2 X  2s+1 (0)  p hki Ik(q ) < ε . k>N (N) p p p As a consequence, since p > 2, kΦ(q )ks+1/2,p < 4ε 6 (2ε) . Since Φ U : U → (N) V = B2ε(0) is a diffeomorphism and Φ(q ) ∈ V it then follows that

(N) s,p q ∈ U ⊂ L0 . (5.4)

(0) −1 s,p On the other hand, it follows from the assumption (5.2) that q ∈ H0 \ L0 (N) −1 s,p v and from Lemma 5.5 (ii) that q ∈ H0 \ L0 contradicting (5.4).

September 29, 2018 36 In the proof of Lemma 5.4 we have used [13, Lemma 1]. We state it for the convenience of the reader.

−1 Lemma 5.5 For any q ∈ H0 \ Zk with k > 1 the initial value problem q˙ = 1 −1 Yk(q), q(0) = q0 has a unique solution t 7→ q(t) in C ((−Ik(q0), ∞),H0 ). The solution has the following additional properties

(i) lim Ik(q(t)) = 0, t&−Ik(q0)

0 (ii) q(t) − q0 ∈ H0 for any −Ik(q0) < t < ∞. o

Proof of Theorem 1.2. The claimed results follow from Lemma 5.1, Lemma 5.2, and Lemma 5.4. v

A first application of Theorem 1.2 is in inverse spectral theory.

−1 Corollary 5.6 Let q be in H0 and −1/2 6 s 6 0, 2 6 p < ∞. If its sequence s,p s,p of gap lengths (γn(q))n>1 is in ` (N), then q ∈ L0 . o

Proof. Let z = (zn)n∈Z = Φ(q). By Theorem 2.3 there exists n0 > 1 depending 2 2 on q such that for any n > n0, γn/2 6 8nπIn(q) 6 2γn, where In(q) denotes 2 the nth action variable of q. Since In = |zn| it then follows that

s+1/2,p s,p (zn)n∈Z ∈ `0 ⇔ (γn)n>1 ∈ ` (N). The claim then follows from Theorem 1.2 saying that

s+1/2,p s,p v (zn)n∈Z ∈ `0 ⇔ q ∈ L0 .

A second application of Theorem 1.2 concerns the isospectral sets defined as −1 follows: For any q ∈ H0 let

−1 Iso(q) := {q˜ ∈ H0 : spec(L(˜q)) = spec(L(q))}.

−1/2,2 For z ∈ `0 let

−1/2,2 2 2 Tz = {w ∈ `0 : |wk| = |zk| ∀ k ∈ Z}.

s+1/2,p s+1/2,p Note that for any z ∈ `0 , Tz ⊂ `0 .

s,p Corollary 5.7 For any q ∈ L0 with −1/2 6 s 6 0 and 2 6 p < ∞, Iso(q) ⊂ s,p s,p L0 and Φ(Iso(q)) = TΦ(q). In particular, Iso(q) is a compact subset of L0 .

−1 Proof. By [11, Theorem 1.2] one has Φ(Iso(q)) = TΦ(q) for any q ∈ H0 . By def- s,p s+1/2,p inition, for any q ∈ L0 , TΦ(q) is a compact subset of `0 . As a consequence, −1 s,p s,p v Iso(q) = Φ (TΦ(q)) ⊂ L0 and Iso(q) is compact in L0 .

September 29, 2018 37 A Auxiliary estimates

In this appendix we collect auxiliary estimates needed in Section3.

Lemma A.1 For any 1/2 < σ < ∞ there exists a constant Cσ > 0 so that for P 1 2σ−1 loghni any n > 1, |m|6=n |m2−n2|σ is bounded by Cσ/n if 1/2 < σ < 1, Cσ n σ if σ = 1, and Cσ/n if σ > 1. o

Proof. Note that for n > 1, n−1 ∞ X 1 1 X 1 X 1 = + 2 + 2 . |m2 − n2|σ n2σ (n2 − m2)σ (m2 − n2)σ |m|6=n m=1 m=n+1

2 2 1 2 2 As n − (n − k) > 3 ((n + k) − n ) for any 0 6 k 6 n it follows that ∞ X 1 X 1 5 · 3σ . |m2 − n2|σ 6 (m2 − n2)σ |m|6=n m=n+1 Furthermore,

∞ 2n ∞ X 1 1 X 1 1 X 1 + + (m2 − n2)σ 6 nσ nσ (m − n)σ (m − n)2σ m=n+1 m=n+2 m=2n+1 1 1 Z 2n 1 Z ∞ 1 6 σ + σ σ dx + 2σ dx. n n n+1 (x − n) 2n (x − n) One computes, Z ∞ 1 1 1 2σ dx = 2σ−1 . 2n (x − n) 2σ − 1 n whereas for σ > 1/2 with σ 6= 1

Z 2n 1 1  1  σ dx 6 σ−1 + 1 n+1 (x − n) |σ − 1| n

R 2n 1 v and n+1 x−n dx = log n. Combining these estimates yields the stated ones.

The next result concerns estimates for the convolution of sequences.

t,2 Lemma A.2 (i) Let −1 t < −1/2. For a = (am)m∈ ∈ ` and b = 6 Z C 1,2 P (bm)m∈ ∈ ` , the convolution a∗b = ( an−mbm)n∈ is well defined Z C m∈Z Z and

ka ∗ bkt,2 6 Ctkakt,2kbk1,2.

(ii) Let −1/2 6 s 6 0, 2 6 p < ∞, and −s − 3/2 < t < 0. For any s,p t+2,p a = (am)m∈ ∈ ` and b = (bm)m∈ ∈ ` , Z C Z C

ka ∗ bks,p 6 Cs,tkaks,pkbkt+2,p. o

September 29, 2018 38 Proof. (i) First note that for any a ∈ `t,2, b ∈ `1,2 C C !2 X 1 X 2 2 |an−m||bm| 6 2 I + 2 II, hni2|t| n∈Z m∈Z where !2 X 1 X I = 1{|n−m| 2|m|}|an−m||bm| hni2|t| 6 n∈Z m∈Z and !2 X 1 X II = 1{|n−m|>2|m|}|an−m||bm| . hni2|t| n∈Z m∈Z To estimate I, use that (hn − mi|t|/hmi)1 2 to conclude by the {|n−m|62|m|} 6 Cauchy-Schwarz inequality that

X 4 2 2 I 6 kakt,2kbk1,2. hni2|t| n∈Z Towards II, note that for |n − m| > 2|m|, |n| > |n − m| − |m| > |m|, and hence |n − m| 6 |n| + |m| 6 2|n| yielding the estimate hn − mi|t| 1{|n−m|>2|m|} 6 2. hni|t| Hence by the Cauchy-Schwarz inequality ! !2 X X 4 X |an−m| X 4 2 2 II 6 hmi|bm| 6 kakt,2kbk1,2. hmi2 hn − mi|t| hmi2 n∈Z m∈Z m∈Z m∈Z This proves item (i). Towards (ii), write

|s| !p p X X hn − mi |an−m| |s| ka ∗ bks,p 6 hmi |bm| . hni|s|hmi|s| hn − mi|s| n∈Z m∈Z Use that hn − mi 6 hnihmi and Young’s inequality to conclude ! X |s| ka ∗ bks,p 6 hmi |bm| kaks,p. m∈Z By the Cauchy-Schwarz inequality we get !1/2 !1/2 X |s| X 2(2−|t|) 2 X 1 hmi |bm| 6 hmi |bm| . hmi(2−|t|−|s|)2 m∈Z m m Since −3/2 − s < t 6 0, it follows that 2 − |t| − |s| > 1/2 and therefore !1/2 X 1 Cs,t := < ∞ hmi(2−|t|−|s|)2 m∈Z implying the claimed result. v

September 29, 2018 39 Finally, we prove the following result on embeddings of sequence spaces.

Lemma A.3 Suppose −1/2 s 0 and 2 p < ∞, then `s,p ,→ `t,2 for any 6 6 6 C C t < s − 1/2 + 1/p. In particular, `s,p ,→ `−1+1/2p,2 for any −1/2 s 0. C C 6 6 o Proof. For p = 2 the claim clearly holds. So it remains to consider the case 2 < p < ∞. Note that for s and p in the given ranges, s − 1/2 + 1/p 6 0. Thus s,p t < 0 and |t| − |s| > 1/2 − 1/p > 0. Let (am)m∈ ∈ ` . By Hölder’s inequality Z C with p/2, p00/2 where 1/p + 1/p00 = 1/2

00 2 2 p !2/p !2/p X |an| X |an| 1 X |an| X 1 6 6 00 . hni2|t| hni2|s| hni2(|t|−|s|) hni|s|p hni(|t|−|s|)p n∈Z n∈Z n∈Z n∈Z

00 −1 Since (|t| − |s|)p > (1/2 − 1/p)(1/2 − 1/p) = 1, one has kakt,2 6 Ckaks,p with C = Ct,s,p < ∞. v

B Schrödinger Operators

In this appendix we review definitions and properties of Schrödinger operators 2 −∂x + q with a singular potential q used in Section3 – see e.g. [10]. Boundary conditions. Denote by H1[0, 1] = H1([0, 1], ) the C C of functions f : [0, 1] → C which together with their distributional derivative ∂xf are in L2 [0, 1]. On H1[0, 1] we define the following three boundary conditions C C (bc),

(per+) f(1) = f(0); (per−) f(1) = −f(0); (dir) f(1) = f(0) = 0.

The corresponding subspaces of H1[0, 1] are defined by C H1 = {f ∈ H1[0, 1] : f satisfies (bc)}, bc C

−1 1 0 1 and their duals are denoted by Hbc := (Hbc) . Note that Hper+ can be canon- ically identified with the Sobolev space H1(R/Z, C) of 1-periodic functions 2 f : R → C which together with their distributional derivative are in Lloc(R, C). 1 1 Analogously, Hper− can be identified with the subspace of H (R/2Z, C) consist- 2 ing of functions f : R → C with f, ∂xf ∈ Lloc(R, C) satisfying f(x + 1) = −f(x) for all x ∈ R. In the sequel we will not distinguish these pairs of spaces. Fur- 1 1 1 thermore, note that Hdir is a subspace of Hper+ as well as of Hper−. Denote by 2 R 1 h·, ·ibc the extension of the L -inner product hf, giI = 0 fg dx to a sesquilinear −1 1 pairing of Hbc and Hbc. Finally, we record that the multiplication

1 1 1 Hbc × Hbc → Hper+, (f, g) 7→ fg, (B.1)

1 1 and the complex conjugation Hbc → Hbc, f 7→ f are bounded operators.

September 29, 2018 40 −1 Multiplication operators. For q ∈ Hper+ define the operator Vbc of multipli- 1 −1 1 cation by q, Vbc : Hbc → Hbc as follows: for any f ∈ Hbc, Vbcf is the element −1 in Hbc given by

1 hVbcf, gibc := hq, fgiper+, g ∈ Hbc.

In view of (B.1), Vbc is a well defined bounded linear operator.

−1 1 Lemma B.1 Assume that q ∈ Hper+. Then for any g ∈ Hdir

(Vper±g) 1 = Vdirg. o Hdir

1 1 Proof. Since any h ∈ Hdir is also in Hper+, the definitions of Vper+ and Vdir imply

hVper+g, hiper+ = hq, ghiper+ = hVdirg, hidir,

v which gives (Vper+g) 1 = Vdirg. Similarly, one sees that Vper−g 1 = Vdirg. Hdir Hdir

1 1 It is convenient to introduce also the space Hper+ ⊕ Hper− and define the multiplication operator V of multiplication by q

1 1 −1 −1 V : Hper+ ⊕ Hper− → Hper+ ⊕ Hper−, (f, g) 7→ (Vper+f, Vper−g).

1 1 1 We note that Hper+ ⊕ Hper− can be canonically identified with H (R/2Z, C),

1 1 1 + − H (R/2Z, C) → Hper+ ⊕ Hper−, f 7→ (f , f ),

+ 1 − 1 where f (x) = 2 (f(x) + f(x + 1)) and f (x) = 2 (f(x) − f(x + 1)). Its dual is denoted by H−1(R/2Z, C). 1 1 1 Fourier basis. The spaces Hper±, H (R/2Z, C) and Hdir and their duals admit the following standard Fourier basis. 1 −1 1 −1 Basis for Hper+, Hper+. Any element f ∈ Hper+ [Hper+] can be represented P 1,2 −1,2 as f = fmem where (fm)m∈ ∈ ` [` ] and m∈Z Z C C

f2m = hf, e2miper+, f2m+1 = 0, ∀ m ∈ Z. P −1 Furthermore, for any q = qmem ∈ H , m∈Z per+ ! X X −1 Vper+f = qn−mfm en ∈ Hper+. n∈Z m∈Z P −1,2 Note that by Lemma A.2, ( qn−mfm)n∈ is in ` . m∈Z Z C 1 −1 1 −1 Basis for Hper−, Hper−. Any element f ∈ Hper− [Hper−] can be represented P 1,2 −1,2 as f = fmem where (fm)m∈ ∈ ` [` ] and m∈Z Z C C

f2m+1 = hf, e2m+1iper−, f2m = 0, ∀ m ∈ Z.

September 29, 2018 41 P −1 Similarly, for any q = qmem ∈ H , m∈Z per+ ! X X −1 Vper−f = qn−mfm en ∈ Hper−. n∈Z m∈Z Basis for H1(R/2Z, C), H−1(R/2Z, C). Any element f ∈ H1(R/2Z, C) [H−1(R/2Z, C)] P can be represented as f = fmem where fm = hf, emi. Here h·, ·i denotes m∈Z the normalized L2-inner product on [0, 2] extended to a sesquilinear pairing between H1(R/2Z, R) and its dual. In particular, for f ∈ H1(R/2Z, C), and m ∈ Z, Z 2 1 −i2mπx hf, emi = f(x)e dx. 2 0 P −1 1 For any q = qmem ∈ H ,→ H ( /2 , ) m∈Z per+ R Z C ! X X −1 V f = qn−mfm en ∈ H (R/2Z, C). n∈Z m∈Z 1 −1 1 Basis for Hdir, Hdir : Any element f ∈ Hdir can be represented as √ √ X X sin f(x) = hf, 2smiI 2sm(x) = fm sm(x), sm(x) = sin(mπx), m>1 m∈Z sin R 1 −1 where fm = 0 f(x)sm(x) dx. For any element g ∈ Hdir one gets by duality X sin sin g = gm sm, gm = hg, smidir. m∈Z One verifies that gsin = −gsin for all m ∈ and P hmi−2|gsin|2 < ∞. For −m m Z m∈Z m any q ∈ H−1 with kqk < ∞ and −1 < t < −1/2, we need to expand for a 0,C t,2 1 −1 P sin given f ∈ H , Vdirf ∈ H in its sine series (Vdirf) sm where by the dir dir m∈Z m definition of Vdir

sin X sin (Vdirf)m = hVdirf, smidir = hq, fsmiper+ = fn hq, snsmiper+. n∈Z sin sin Using that f−n = −fn for any n ∈ Z and that 1 s (x)s (x) = (cos((m − n)πx) − cos((m + n)πx)) m n 2 it follows that

X sin X sin fn hq, snsmiper+ = fn hq, cos((m − n)πx)iper+. m−n even m−n even 1 Note that hq, cos((m − n)πx)iper+ is well defined as cos((m − n)πx) ∈ Hper+ if m − n is even. If m − n is odd, we decompose the difference of the cosines in 1 Hper+ as follows

cos((m − n)πx) − cos((m + n)πx) = (cos((m − n)πx) − cos(πx)) − (cos((m + n)πx) − cos(πx))

September 29, 2018 42 sin sin and then obtain, using again that f−n = −fn for all n ∈ Z, X sin X sin fn hq, snsmiper+ = fn hq, cos((m − n)πx) − cos(πx)iper+. m−n odd m−n odd Altogether we have shown that ! X X cos sin Vdirf = qm−nfn sm, m∈Z n∈Z where  cos hq, cos(kπx)iper+, if k ∈ Z even, qk = (B.2) hq, cos(kπx) − cos(πx)iper+, if k ∈ Z odd.

Taking into account that kqkt,2 < ∞ with −1 < t < −1/2, one argues as in [10, Proposition 3.4], using duality and interpolation, that !1/2 X 2t cos 2 hmi |qm | 6 Ctkqkt,2. (B.3) m∈Z Schrödinger operators with singular potentials. For any q ∈ H−1 denote by 0,C 2 −1 L(q) the −∂x + V acting on H (R/2Z, C) with domain 1 1 1 1 H (R/2Z, C). As H (R/2Z, C) = Hper+ ⊕ Hper− and V = Vper+ ⊕ Vper− the 1 operator L(q) leaves the spaces Hper± invariant and L(q) = Lper+(q)⊕Lper−(q) 2 with Lper±(q) = −∂x + Vper±. Hence the spectrum spec(L(q)) of L(q), also referred to as spectrum of q, is the union spec(Lper+(q)) ∪ spec(Lper−(q)) of the spectra spec(Lper±(q)) of Lper±(q). It is known to be discrete and to consist of complex eigenvalues which, when counted with multiplicities and ordered lexicographically, satisfy

+ − + ± 2 2 2 λ0 4 λ1 4 λ1 4 ··· , λn = n π + n`n, – see e.g. [14]. For any q ∈ H−1 there exists N 1 so that 0,C > 2 2 + 2 2 |λ − n π | 6 n/2, n > N, |λn | 6 (N − 1) π + N/2, n < N, where N can be chosen locally uniformly in q. Since for q = 0 and n > 0, + + + ∆(λ2n(0), 0) = 2 and ∆(λ2n+1(0), 0) = −2, all λ2n(0) are 1-periodic and all + λ2n+1(0) are 1-antiperiodic eigenvalues of q = 0. By considering the compact interval [0, q] = {tq : 0 t 1} ⊂ H−1 it then follows after increasing N, if 6 6 0,C necessary, that for any n > N

+ − λn (q), λn (q) ∈ spec(Lper+(q)), [spec(Lper−(q))] if n even [odd]. (B.4)

For any q ∈ H−1 and n N the following Riesz projectors are thus well 0,C > defined on H−1(R/2Z, C) Z 1 −1 Pn,q := (λ − L(q)) dλ. 2πi |λ−n2π2|=n

September 29, 2018 43 1 1 For n > N even [odd], the range of Pn,q is contained in Hper+ [Hper−]. If q = 0, we write Pn for Pn,0. −1 2 Similarly, for any q ∈ H denote by Ldir(q) the unbounded operator −∂ + 0,C x −1 1 Vdir acting on Hdir with domain Hdir. Its spectrum spec(Ldir(q)) is known to be discrete and to consist of complex eigenvalues which, when counted with multiplicities and ordered lexicographically, satisfy

2 2 2 µ1 4 µ2 4 ··· , µn = n π + n`n, – see e.g. [14]. By increasing N as chosen above, if necessary, we can thus assume that

2 2 2 2 |µn − n π | < n/2, n > N, |µn| 6 (N − 1) π + N/2, n 6 N.

In particular, for any n > N, µn is simple and the corresponding Riesz projector Z 1 −1 Πn,q := (λ − Ldir(q)) dλ 2πi |λ−n2π2|=n

−1 is well defined on Hdir . If q = 0, we write Πn for Πn,0. Regularity of solutions. In Section3 we consider solutions f of the equation (L(q) − λ)f = g in Ls,p and need to know their regularity. ?,C Lemma B.2 For any q ∈ Ls,p with −1/2 s 0 and 2 p < ∞ the following 0,C 6 6 6 holds: For any g ∈ Ls,p and any λ ∈ , a solution f ∈ H1( /2 , ) of the ?,C C R Z C inhomogeneous equation (L(q) − λ)f = g is an element in Ls+2,p. ?,C o

2 s,p Proof. Write (L−λ)f = g as Aλf = V f −g where Aλ = ∂ +λ. Since q ∈ L , x 0,C Lemma A.3 implies that q and g are in Lr,2 with r = s − 1/2 + 1/2p. By ?,C r,2 r,2 Lemma A.2 (i), V f ∈ L and hence Aλf = V f − g ∈ L implying that ?,C ?,C f ∈ Lr+2,2. Since −s − 3/2 −1 < r 0, Lemma A.2 (ii) applies. Therefore, ?,C 6 6 s,p s+2,p V f ∈ L and using the equation Aλf = V f −g once more one gets f ∈ L ?,C ?,C as claimed. v

s,p For any q ∈ L with −1/2 s 0, 2 p < ∞, and n ns,p as in 0,C 6 6 6 > Lemma 3.7, introduce  + − + − Null(L(q) − λn ) ⊕ Null(L(q) − λn ), λn 6= λn , En ≡ En(q) := + 2 + − Null(L(q) − λn ) , λn = λn .

1 1 Then En is a two-dimensional subspace of Hper+ [Hper−] if n is even [odd]. The following result shows that elements in En are more regular.

Lemma B.3 For any q ∈ Ls,p with −1/2 s 0, 2 p < ∞ and for any 0,C 6 6 6 s+2,p 1 s+2,p 1 n ns,p En(q) ⊂ L ∩ H [L ∩ H ] if n is even [odd]. > ?,C per+ ?,C per− o

September 29, 2018 44 Proof. By Lemma B.2 with g = 0, any eigenfunction f of an eigenvalue λ of L(q) is in Ls+2,p. Hence if λ+ 6= λ− or if λ+ = λ− and has geometric multiplicity ?,C n n n n s+2,p + − two, then En ⊂ L . Finally, if λ = λ is a double eigenvalue of geometric ?,C n n + multiplicity 1 and g is an eigenfunction corresponding to λn , there exists an 1 + element f ∈ H (R/2Z, C) so that (L − λn )f = g. As g is an eigenfunction it is in Ls+2,p and by applying Lemma B.2 once more, it follows that f ∈ ?,C s+2,p s+2,p L . Clearly, En = span(g, f) and hence En ⊂ L also in this case. ?,C ?,C ± By the consideration above one knows that λn are 1-periodic [1-antiperiodic] s+2,p 1 s+2,p 1 eigenvalues of q if n is even [odd]. Hence En ⊂ L ∩ H [L ∩ H ] ?,C per+ ?,C per− if n is even [odd] as claimed. v

Estimates for projectors. The projectors Pn,q [Πn,q] with n > N are defined −1 −1 1 1 on H (R/2Z, C) [Hdir ] and have range in H (R/2Z, C) [Hdir]. The following 2 2 result concerns estimates for the restriction of Pn,q [Πn,q] to L = L ([0, 2], C) 2 2 [L (I) = L ([0, 1], C)] needed in Section3 for getting the asymptotics of µn −τn. Lemma B.4 Assume that q ∈ H−1 with kqk < ∞ and −1 < t < −1/2. 0,C t,2 0 Then there exist constants Ct > 0 (only depending on t) and N > N (with N as above) so that (log n)2 (i) kP − P k C kqk n,q n L2→L∞ 6 t n1−|t| t,2 (log n)2 (ii) kΠ − Π k C kqk n,q n L2(I)→L∞(I) 6 t n1−|t| t,2 The constant N 0 can be chosen locally uniformly in q. o Remark B.5. We will apply Lemma B.4 for potentials q ∈ Ls,p with −1/2 0,C 6 s 6 0 and 2 6 p < ∞ using the fact that by the Sobolev embedding theorem there exists −1 < t < −1/2 so that Ls,p ,→ Lt,2 . 0,C 0,C Proof. The proof follows the line of arguments of [10]. Let us first prove (i). 2 2 Suppose |λ−m π | > hmi for a given m > 0. Then we may write with L ≡ L(q) 1/2 1/2 λ − L = Dλ (Iλ − Sλ)Dλ , 2 where the operators Dλ, Iλ and Sλ are defined with respect to the L -orthonormal 2 basis (em)m∈Z of L ([0, 2], C) by (l, m ∈ Z) λ − l2π2 D1/2(l, m) = |λ − l2π2|1/2δ ,I (l, m) = δ , λ lm λ |λ − l2π2| lm ql−m Sλ(l, m) = . |λ − l2π2|1/2|λ − m2π2|1/2 0 0 2 2 Then there exists N with N > N so that any λ ∈ C with |λ − n π | = 0 −1/2 n and n > N is in the of L. Hence Iλ − Sλ = Dλ (λ − −1/2 2 2 L)Dλ : L ([0, 2], C) → L ([0, 2], C) is a boundedly invertible operator and Z 1 −1/2 −1 −1 −1/2 Pn,q − Pn = Dλ (Iλ − Sλ) SλIλ Dλ dλ. 2πi |λ−n2π2|=n

September 29, 2018 45 2 2 2 0 For any f ∈ L ([0, 2], C) and λ ∈ C with |λ − n π | = n, n > N , one gets by Cauchy-Schwarz and Lemma A.1,  1/2 −1/2 X |fm| 2 X 1 kDλ fkL∞ 6 6  +  kfkL2 |λ − m2π2|1/2 n |n2 − m2| m∈Z |m|6=n (log n)1/2 4 kfk . 6 n1/2 L2 Similarly, one has

2 !1/2 −1/2 X |fm| 1 kD fk kfk , λ L2 6 |λ − m2π2| 6 n1/2 L2 m∈Z 2 2 and by [10, Lemma 1.2] using that kqkt,2 < ∞ with −1 < t < −1/2, |λ−n π | = 0 n and n > N , log n kS k C kqk . (B.5) λ L2→L2 6 n1−|t| t,2 0 By increasing N , if necessary, it follows that kSλkL2→L2 6 1/2 and hence 1/2 0 kPn,q − PnkL2→L∞ 6 Ct sup (log n) kSλkL2→L2 , ∀ n > N . |λ−n2π2|=n Combined with (B.5), this yields the estimate of (i). (ii) The estimate for the Dirichlet projectors are proved in a similar way: 2 2 Suppose that |λ − m π | > hmi for a given m > 0. Then with Ldir ≡ Ldir(q) ˜ 1/2 ˜ ˜ ˜ 1/2 λ − Ldir = Dλ (Iλ − Sλ)Dλ where the operators D˜ λ, I˜λ, S˜λ are defined in terms of their coefficients with respect to the sine basis sn, n > 1. For any l, m ∈ Z, λ − l2π2 D˜ 1/2(l, m) = |λ − l2π2|1/2δ , I˜ (l, m) = δ , λ lm λ |λ − l2π2| lm cos ql−m S˜λ(l, m) = , |λ − l2π2|1/2|λ − m2π2|1/2

cos 0 with ql−m given by (B.2). By increasing N , if necessary, we can assure that 2 2 0 any λ ∈ C with |λ − n π | = n and n > N is in the resolvent set of Ldir. Hence ˜ ˜ ˜ −1/2 ˜ −1/2 2 2 (Iλ − Sλ) = Dλ (λ − Ldir)Dλ : L ([0, 1], C) → L ([0, 1], C) is a boundedly invertible operator and Z 1 ˜ −1/2 ˜ ˜ −1 ˜ ˜−1 ˜ −1/2 Πn,q − Πn = Dλ (Iλ − Sλ) SλIλ Dλ dλ. 2πi |λ−n2π2|=n Using the same arguments as in the proof of (i) together with (B.3), the estimate (ii) then follows. Going through the arguments of the proofs of (i) and (ii) one sees that N 0 can be chosen locally uniformly in q and Ct independently of q. v

September 29, 2018 46 References

[1] Dispersive Wiki. http://wiki.math.toronto.edu/DispersiveWiki/

[2] Bikbaev, R.F., Kuksin, S.B.: On the parametrization of finite-gap solutions by frequency and wavenumber vectors and a theorem of I. Krichever. Lett. Math. Phys. 28(2), 115–122, 1993.

[3] Bobenko, A.I., Kuksin, S.B.: Finite-gap periodic solutions of the KdV equation are nondegenerate. Phys. Lett. A 161(3), 274–276, 1991.

[4] Djakov, P., Mityagin, B.: Instability zones of periodic 1-dimensional Schrödinger and Dirac operators. Russian Math. Surveys 61, 663–766, 2006.

[5] Djakov, P., Mityagin, B.: Spectral gaps of Schrödinger operators with periodic singular potentials. Dyn. Partial Differ. Equ. 6(2), 95–165, 2009.

[6] Grébert, B., Kappeler, T.: The defocusing NLS equation and its normal form. European Mathematical Society (EMS), Zürich, 2014.

[7] Henrici, A., Kappeler, T.: Nekhoroshev theorem for the periodic Toda lattice. Chaos 19(3), 033,120–033,113, 2009.

[8] Huang, G., Kuksin, S.B.: The KdV equation under periodic boundary conditions and its perturbations. Nonlinearity 27(9), R61–R88, 2014.

[9] Kappeler, T., Mityagin, B.: Estimates for periodic and Dirichlet eigenvalues of the Schrödinger operator. SIAM J. Math. Anal. 33(1), 113–152, 2001.

[10] Kappeler, T., Möhr, C.: Estimates for periodic and Dirichlet eigenvalues of the Schrödinger operator with singular potentials. J. Funct. Anal. 186(1), 62–91, 2001.

[11] Kappeler, T., Möhr, C., Topalov, P.: Birkhoff coordinates for KdV on phase spaces of distributions. Selecta Math. (N.S.) 11(1), 37–98, 2005.

[12] Kappeler, T., Pöschel, J.: KdV & KAM. Springer, Berlin, 2003.

[13] Kappeler, T., Serier, F., Topalov, P.: On the symplectic phase space of KdV. Proc. Amer. Math. Soc. 136(5), 1691–1698, 2008.

−1 [14] Kappeler, T., Topalov, P.: Riccati representation for elements in H (T) and its applications. Pliska Stud. Math. Bulgar. 15, 171–188, 2003.

− [15] Kappeler, T., Topalov, P.: Global wellposedness of KdV in H 1(T, R). Duke Math. J. 135(2), 327–360, 2006.

[16] Korotyaev, E.: Characterization of the spectrum of Schrödinger operators with periodic distributions. Int. Math. Res. Not. 2003(37), 2019–2031, 2003.

[17] Korotyaev, E., Kuksin, S.B.: KdV Hamiltonian as function of actions. arXiv , 2011.

September 29, 2018 47 [18] Krichever, I.M.: Perturbation Theory in Periodic Problems for Two-Dimensional Integrable Systems. Sov. Sci. Rev. C. Math. Phys. 9, 1–103, 1992.

[19] Molnar, J.C.: On Two-Sided Estimates for the Nonlinear Fourier Transform of KdV. arXiv:1502.04550 , 2015.

[20] Pöschel, J.: Hill’s potentials in weighted Sobolev spaces and their spectral gaps. Math. Ann. 349(2), 433–458, 2011.

[21] Savchuk, A.M., Shkalikov, A.A.: Sturm-Liouville operators with distribution po- tentials. Tr. Mosk. Math. Obs. 64, 159–212, 2003.

September 29, 2018 48