arXiv:1905.01849v2 [math.AP] 4 Nov 2019 rte hl h rtato a iiigteIsiu Mittag–Leffler p Institut co last working the favourable The visiting providing for was thank. institutions warmly author these we to first whom grateful the Saut, while J.C. by written organised workshop a operators. Toeplitz NTEITGAIIYO H BENJAMIN–ONO THE OF INTEGRABILITY THE ON .Introduction 1. hspprwssatdwiebt uhr eevstn h al Ins Pauli the visiting were authors both while started was paper This phrases. and words Key of spectrum the potentials On equation gap References BO one the and waves C. for Traveling Appendix pair Lax the B. On Appendix Theorem A. of Appendix Proof potentials gap 8. Finite brackets 7. Poisson Gradients and 6. Analyticity coordinates Birkhoff 5. formulae Complex trace and function 4. Generating operator Lax 3. The 2. 2010 Date coe ,2019. 4, October : ahmtc ujc Classification. Subject Mathematics ino h ejmnOoeuto n nteueo generating hierarchy. a Benjamin–Ono of formula- entire use pair the the Lax encodes on the which and relies in equation functional, coordinates operator Benjamin–Ono Lax such the the of of of tion construction study The spectral the functions variable. on periodic time almost the are equation of Benjamin–Ono the of tions o oitgaei yqartr n ec r lorfre oas to referred the also all are consequence, a hence As and coefficients. quadrature Fourier by nonlinear it integrate to low qaini h togs osbesne tamt lblBirkhoff global admits space it the sense: on possible coordinates strongest (pseudo)differ the integrable an in is equation torus, the on considered when tion, Abstract. A RC G PATRICK QAINO H TORUS THE ON EQUATION nti ae epoeta h ejmnOoequa- Benjamin-Ono the that prove we paper this In RR N HMSKAPPELER THOMAS AND ERARD ´ 1 ejmnOoeuto,Lxpi,Hryspace, Hardy pair, Lax equation, Benjamin–Ono n Theorem and L Contents 2 ( T .Teeaecodntswihal- which coordinates are These ). 1 L 71 rmr,4B5secondary. 47B35 primary, 37K15 u 51 2 L 2 nditions. ( T solu- ) ential eare We . iuefor titute r was art 60 57 55 54 44 37 25 21 15 8 2 ON THE INTEGRABILITY OF THE BO EQUATION 2

1. Introduction In this paper we consider the Benjamin-Ono (BO) equation on the torus, (1.1) ∂ u = H∂2u ∂ (u2) , x T := R/2πZ , t R, t x − x ∈ ∈ where u u(t, x) is real valued and H denotes the Hilbert transform, ≡ inx 1 2π inx defined for f = n Z f(n)e , f(n)= 2π 0 f(x)e− dx, by ∈ PHf(x) := i sign(nR)f(n) einx b b− n Z X∈ with sign( n) := 1 for any n 1, whereasb sign(0) := 0. This pseudodifferential± equation± (ΨDE)≥ in one space dimension has been introduced by Benjamin [3] and Ono [26] to model long, one-way inter- nal gravity waves in a two-layer fluid. It has been extensively studied, in particular the wellposedness problem on the torus as well as on the real line. On appropriate Sobolev spaces, the BO equation (1.1) can be written in Hamiltonian form 2π 1 1 1 ∂ u = ∂ ( (u)) , (u) := ( ∂ 1/2u)2 u3 dx t x ∇H H 2π 2 | x| − 3 Z0   1/2 where ∂x is the square root of the Fourier multiplier operator ∂x given by| | | | ∂ f(x)= n f(n)einx . | x| | | n Z X∈ Note that the L2 gradient of canb be computed to be ∂ u u2 − ∇H H | x| − and that ∂x is the Hamiltonian vector field corresponding to the Gardner bracket,∇H defined for any two functionals F,G : L2 C with sufficiently regular L2 gradients by → − 2π 1 F,G := (∂ F ) Gdx . { } 2π x∇ ∇ Z0 The main result of this paper says that the BO equation (1.1) is an integrable ΨDE. In fact, we show that it admits global Birkhoff coor- dinates and hence is an integrable ΨDE in the strongest possible sense. To state our results more precisely, we first need to introduce some no- tation. Denote by L2 the C L2(T, C) of L2 integrable, complex valued functions with− the standard inner product− 2π 1 (1.2) f g := f(x)g(x) dx h | i 2π Z0 and the corresponding norm f := f f 1/2. Furthermore, we denote by L2 the R Hilbert space Lk 2(kT, Rh),| consistingi of elements u L2 r − ∈ ON THE INTEGRABILITY OF THE BO EQUATION 3 which are real valued and let L2 := u L2 u 1 =0 . r,0 { ∈ r |h | i } s s For any subset J Z 0 and any s R, h (J) h (J, C) denotes the weighted ℓ2 sequence⊂ ≥ space ∈ ≡ − s h (J)= (zn)n J C : (zn)n J s < { ∈ ⊂ k ∈ k ∞} where

2s 2 1/2 (zn)n J s := n zn , n := max 1, n . k ∈ k h i | | h i { | |} n J X∈  N Z s s N In case where J = := n : n 1 we write h+ instead of h ( ). {s ∈ R ≥ } s N R s N R In the sequel, we view h+ as the Hilbert space h ( , ) h ( , ) − s ⊕ by identifying a sequence (zn)n N h+ with the pair of sequences s ∈ ∈ s (Re zn)n N, (Im zn)n N in h (N, R) h (N, R). ∈ ∈ ⊕  The main result of this paper then reads as follows. Theorem 1. There exists a homeomorphism

2 1/2 Φ: Lr,0 h+ , u (ζn(u))n 1 → 7→ ≥ so that the following holds: (B1) For any n 1, ζ : L2 C is real analytic. ≥ n r,0 → (B2) The Poisson brackets between the coordinate functions ζn are well defined and for any n, k 1, ≥ ζ , ζ = iδ , ζ ,ζ =0 . { n k} − nk { n k} 1 (B3) On its domain of definition, Φ− is a (real analytic) function, which only depends on the actionsH◦ζ 2, n 1. | n| ≥ The coordinates ζn are referred to as complex Birkhoff coordinates. Remark 1.1. (i) The Birkhoff map Φ is bounded, meaning that for 2 any bounded subset B of Lr,0, the image Φ(B) is bounded. Indeed, this is a direct consequence of the trace formula of Proposition 3.1, saying that for any u L2 , ∈ r,0 ∞ u 2 =2 n ζ 2 . k k | n| n=1 X In analogy to Parseval’s identity in Fourier analysis, we refer to this trace formula as Parseval’s identity for the nonlinear map Φ. (ii) When restricted to submanifolds of finite gap potentials (cf. Def- inition 2.2), the map Φ is a canonical, real analytic diffeomorphism onto corresponding Euclidean spaces – see Theorem 3 in Section 7 for details. ON THE INTEGRABILITY OF THE BO EQUATION 4

In subsequent work [11] we plan to further study the regularity of the Birkhoff map of Theorem 1 and its restrictions to the scale of Sobolev spaces Hs , s 0, where r,0 ≥ 1/2 Hs := u L2 u < , u := n 2s u(n) 2 r,0 { ∈ r,0 |k ks ∞} k ks h i | | n Z X∈  and to make a detailed analysis of the solution map of theb Benjamin– Ono equation, including qualitative properties of solutions of (1.1) such as their long time behaviour. As an immediate application of Theorem 1 in this direction we mention the following result on the solutions of 2 the BO equation for initial data in Lr,0 obtained by quadrature, when equation (1.1) is expressed in Birkhoff coordinates. 2 Theorem 2. For every initial data u(0) in Lr,0, the solution of the BO equation with initial data u(0), t R u(t) L2 , ∈ 7→ ∈ r,0 2 is almost periodic. Its orbit is relatively compact in Lr,0. Remark 1.2. (i) The solutions of the BO equation of Theorem 2 coin- cide with the ones obtained by Molinet in [20] (cf. also [21]). Since the average is conserved by (1.1), the results of Theorem 2 easily extend to 2 solutions with initial data in Lr. (ii) Expressing a given solution t R u(t) L2 of (1.1) in Birkhoff ∈ 7→ ∈ r,0 coordinates one immediately sees that Corollary 1.3 in [31] on the re- currence of solutions extends in the sense that for any initial datum u(0) L2 there exists a sequence 0 < t < t < of times with ∈ r,0 1 2 ··· tn so that limn u(tn) u(0) =0. (iii)→∞ As another immediate→∞ k application− k of Theorem 1, we present a new proof of the result due to Amick&Toland [2], characterizing the trav- eling wave solutions of the BO equation as solutions with initial data given by one gap potentials – see Proposition B.1 in Appendix B for details.

Outline of the proof of Theorem 1: At the heart of the proof is the formulation ∂tLu = [Bu, Lu] of the BO equation, discovered by Nakamura [24] and Bock&Kruskal [4]. It is reviewed at the beginning 2 of Section 2 – see also Appendix A. For any u Lr, Lu is a selfadjoint pseudodifferential operator of order one whose∈ spectrum is conserved along the flow of the BO equation. In contrast to integrable PDEs with a Lax pair formulation such as the Korteweg-de Vries or the nonlinear Schr¨odinger equation, the Lax operator Lu for the BO equation is not a differential operator. When considered on the L2 := h + { ∈ L2 : h(n)=0 n < 0 , the operator L is given by L = i∂ T ∀ } u u − x − u where Tu is the Toeplitz operator Tuh = Π(uh) and Π denotes the Szeg˝oprojectorb Π : L2 L2 . Its spectrum consists of real eigenvalues r → + ON THE INTEGRABILITY OF THE BO EQUATION 5 which when listed in increasing order and with their multiplicities take the form λ0 λ1 λ2 ... with λn as n . In our analysis ≤ ≤ ≤ →∞ 2 →∞2 ix of the spectrum of Lu, the shift operator S : L+ L+, h e h plays an important role. We point out that this operator→ played7→ also a major role in the study of the Szeg˝oequation (cf. [8], [9]). Using commutator relations between Lu and S (cf. Proposition 2.1) we infer that for any u L2 ∈ r γn(u) := λn(u) λn 1(u) 1 0 , n 1 . − − − ≥ ∀ ≥ Since all eigenvalues are simple, γ γ (u), n 1, can be shown to n ≡ n ≥ depend real analytically on u. Note that the eigenvalues λn, n 1, n ≥ can then be expressed in terms of the γk’s by λn = λ0 + n + k=1 γk. Furthermore, we provide an orthonormal basis of eigenfunctions fn, P where each fn depends analytically on u (Definition 2.1). A key ingredient of the proof of Theorem 1 is the generating function 2 λ(u), which we study in Section 3. For any u Lr, λ λ(u) is the Hmeromorphic function, defined by ∈ 7→ H

1 (u) := (L + λId)− 1 1 . Hλ h u | i Note that λ(u) is holomorphic on C λn(u): n 0 and might have simpleH poles at λ = λ (u), n \0. {− By Proposition≥ } 3.1, (u) − n ≥ Hλ can be expressed as a function of λ0 and the γn’s, and hence is a family of integrals of the BO equation, parametrized by the spectral param- eter λ. For our study of the BO equation it plays a role comparable to the one of the discriminant in the normal form theory of the KdV equation. In particular, λ admits a product representation (Proposi- tion 3.1) and an expansionH at λ = whose coefficients constitute the BO hierarchy. The functionals u ∞1 , u 2 as well as the BO Hamil- tonian are elements in this hierarchy.h | i k Ask an application we derive a H 2 2 trace formula for u (Proposition 3.1), implying that for any u Lr,0 k k 1/2 ∈ the sequence (√γn)n 1 is in h+ . Furthermore, we show that for any ≥ n 0, the absolute value 1 fn is a function of λ0 and the γn’s alone (Corollary≥ 3.1). |h | i| Our candidates for the Birkhoff coordinates (ζn)n 1 are then ob- ≥ tained by normalizing ( 1 fn )n 1 so that ζn = √γn for any n 1 (cf. h | i ≥ | | ≥ Section 4) and the corresponding action angle coordinates are (γn)n 1, ≥ (arg 1 fn )n 1. These variables are very reminiscent of the canonical variablesh | i constructed≥ for the Toda system with Dirichlet boundary con- dition in [23]. In fact, the Toda system admits a Lax pair formulation where the Lax operator L is a Jacobi matrix. In the case of Dirichlet boundary conditions, all eigenvalues of L are simple. Canonical vari- ables are then given by these eigenvalues and the first components of appropriately normalized eigenvectors – see [23] for more details. 2 In the sequel, the map Φ : u (ζn(u))n 1, defined on Lr,0 and ≥ 1/2 7→ with values in h+ , is studied in detail. A key observation is that any ON THE INTEGRABILITY OF THE BO EQUATION 6

2 u Lr,0 can be reconstructed from 1 fn , n 0, allowing to prove that∈ Φ is one to one (formula (4.5), Propositionh | i ≥ 4.2). The generating function also plays a major role for proving that (ζn)n 1,(ζn)n 1 satisfy the canonical relations stated in Theorem 1. Since for≥ any n≥ 1, the 2 1 ≥ L gradient of ζn is in H (Proposition 5.1), it follows that the Pois- son− brackets among∇ these functionals are well defined. A key ingredient for proving that they satisfy canonical relations, is the observation that, for λ sufficiently large, the ΨDE ∂tu = ∂x λ(u) admits the Lax pair dLu λ λ ∇H formulation dt = [Bu , Lu] where Bu is a certain skew adjoint operator (Proposition 6.1). As a consequence, the eigenvalues λp commute with , implying that γ ,γ vanishes on L2 for any p, n 1 (Corol- Hλ { p n} r,0 ≥ lary 6.1). In a similar fashion one computes λ, 1, fn and then {H h i} 2 derives that for any p 1, n 0, γp, 1, fn = i 1, fn δpn on Lr,0 (Proposition 6.2). ≥ ≥ { h i} h i As a last key ingredient into the proof of the claimed canonical re- lations for (ζn)n 1, (ζn)n 1 we consider finite gap potentials: we say 2 ≥ ≥ that u Lr is a finite gap potential if the set n 1 : γn(u) > 0 is finite. We∈ then show that for any N 1, the{ restriction≥ of Φ to} the submanifold of special finite gap potentials,≥ defined by UN := u L2 : γ (u) > 0, γ (u)=0 n > N , UN { ∈ r,0 N n ∀ } N 1 is a symplectic diffeomorphism onto C − C∗ (cf. Theorem 3 in Section 7). × In the course of the proof we derive a formula for an arbitrary element ix in N , showing in particular that Πu is of the form Πu(x) = R(e ) U ′ where R is a rational function, R(z)= z Q (z) , and Q(z) a polynomial − Q(z) of degree N with roots outside the unit disc – see (7.9), Remark 7.1. 2 We note that N 0 N is the set of all finite gap potentials in Lr,0 2≥ U and is dense in Lr,0. In Section 8 we make a synopsis of the proof of Theorem 1 andS as an application of Theorem 1 show Theorem 2.

Related work: The Benjamin-Ono equation has been extensively stud- ied. For an excellent account we refer to the recent survey by J.C. Saut [30]. Besides the foundational work of Benjamin [3] and Ono [26], let us point out a few highlights, relevant in the context of the present paper. The first results about (1.1) by PDE methods can be found in [29], [1] and the most recent ones in [20], [21], [13] – we refer to [21] and [30] for numerous other references. Concerning the investigation of the integrability of (1.1), besides the discovery of the Lax pair formulation of this equation already referred toabove (cf. [24], [4]), we mention the pioneering work of Fokas&Ablowitz [7], Coifman&Wickerhauser [5], Kaup&Matsuno [16], and the more re- cent contributions by Wu [32], [33], on the scattering transform for the BO equation, which concerns the integrability of the BO equation ON THE INTEGRABILITY OF THE BO EQUATION 7 when considered on the real line. Surprisingly, the integrability of the BO equation on the torus has not been studied in detail so far. We point out the work by Satsuma&Ishimori [28] on the construction of multi-phase solutions of (1.1) by Hirota’s bilinear method, the one by Amick&Toland [2] on the characterisation of the solitary wave solutions as well as the work by Dobrokhotov&Krichever [6] where multi-phase solutions are constructed by the method of finite zone integration. We refer to these solutions as finite gap solutions and they are treated in detail in Section 7. In their work on the quantum BO equation, Nazarov&Sklyanin [25] introduced a generating function for the classi- cal BO hierarchy of the type introduced in Section 3, which admits a quantum analogue. The approach of Nazarov&Sklyanin has been fur- ther developed by Moll [22] in the context of the classical BO equation. Finally, we mention the work by Tzvetkov&Visciglia on the construc- tion of invariant Gaussian measures and applications to the long time behaviour of solutions of (1.1) (cf. [31] and references therein). Our analysis of the integrability of (1.1) borrows from the one carried out for the Korteweg-de Vries and the defocusing nonlinear Schr¨odinger equations as well as the Szeg˝oequation. As the BO equation, the Korteweg-de Vries and the defocusing nonlinear Schr¨odinger equations are integrable in the sense that they admit global Birkhoff coordinates – see [14], [12] and references therein. Furthermore, it turned out that the Szeg˝oequation, an integrable ΨDE introduced and studied in detail by G´erard and Grellier, is closely related to the BO equation and tools developed for the study of the integrability of this equation were of great use in the present paper – cf. [8], [9], [10] and references therein.

2 2 2 Notation: We have already introduced the Hilbert spaces L , Lr,Lr,0, 2 2 2 and the Hardy space L+ as well as the Szeg˝oprojector Π : L L+. More generally, for any 1 p , Lp Lp(T, C) denotes the standard→ complex Lp space and the≤ L≤∞p norm≡ of an element f Lp is denoted − − ∈p by f Lp . The subspace of real valued functions in L is denoted by p k k p T R Lr L ( , ). Furthermore, we have already introduced the Sobolev ≡ s s spaces Hr,0. More generally, we denote by H the standard Sobolev s T C s space H ( , ) and by Hr the corresponding real subspace. In order to define the L2 gradient of a differentiable functional µ : f µ(f)= − C 2 7→ Reµ(f)+ iImµ(f) , defined on Lr, introduce the complex bilinear form, ∈

2π 1 f,g := f(x)g(x) dx , f,g L2 . h i 2π ∀ ∈ Z0

Note that the function g in the latter integral is not complex conju- gated. The L2 gradient µ(f) of µ at f L2 is then defined as the − ∇ ∈ r ON THE INTEGRABILITY OF THE BO EQUATION 8 element in L2, uniquely determined by µ(f),g = dµ(f)[g] , g L2 . h∇ i ∀ ∈ r 2 2 The L gradient of a functional, defined on a subspace of Lr such as a Sobolev− space, is defined in a similar fashion by extending the bilinear form f,g as a dual paring between the subspace and its dual. In case h i 2 2 the functional µ is defined on a complex neighborhood of Lr in L and analytic there, the L2 gradient of µ is defined in the same way. Finally, we recall that− we have already introduced the sequence spaces hs(J) where J Z. In case s = 0, we write ℓ2(J) for h0(J) and endow it with the⊂ standard inner product

(zn)n J (wn)n J ℓ2 = znwn . h ∈ | ∈ i n J X∈ More generally, for any 1 p < and s R, denote by ℓp,s(J) ℓp,s(J, C) the weighted ℓp ≤sequence∞ space, consisting∈ of sequences ξ ≡= − (ξn)n J C so that ∈ ⊂ ps p 1/p ξ p,s := n ξ < . k kℓ h i | n| ∞ n J X∈  Correspondingly we define the sequence spaces ℓp,s(J, R). The spaces ,s ,s ,s ℓ∞ (J) ℓ∞ (J, C) and ℓ∞ (J, R) are defined in a similar fashion. In ≡ p the sequel we will often use the notation ℓn to denote the n–th element of a sequence in ℓp.

2. The Lax operator Nakamura [24] and Bock&Kruskal [4] discovered that the BO equa- tion admits a Lax pair. In the periodic setup it can be described as 2 follows. For any u Lr, introduce the unbounded operators Lu and 2 ∈ Bu on L+. (2.1) L = D T , B := iD2 +2iT 2iDT . u − u u D(Πu) − u 2 where D is the on L+, defined by D := i∂x, Tu 2 − the Toeplitz operator on L+ with symbol u, T : L2 L2 , h T h = Π(uh) u + → + 7→ u 2 2 2 and Π : L L+ the Szeg˝oprojector. For any u Lr, Tu is an unbounded operator,→ satisfying the estimate ∈

T (h) u h ∞ . k u k≤k kk kL 1 Hence Lu = D Tu is an unbounded operator with domain H+ := 1 2 − 2 2 2 H L+ and Bu is such an operator with domain H+ := H L+. Here∩ Hm Hm(T, C) denotes the standard of∩ order m 1. Furthermore,≡ it follows that T is D compact. Arguing as in ≥ u − ON THE INTEGRABILITY OF THE BO EQUATION 9

Appendix C (cf. Lemma C.1) one verifies that Lu is selfadjoint and Bu skew-adjoint. The Lax pair formulation of the BO equation then reads d (2.2) L = [B , L ] . dt u u u Remark 2.1. In order to make the paper self-contained, we verify the Lax pair formulation of the BO equation in Appendix A. Clearly, the operator Bu is not uniquely determined. Instead of Bu one might take any operator of the form Bu + Pu where Pu is a skew adjoint operator commuting with Lu. In particular, one might choose instead of Bu the operator B˜ := B iL2 . By a straightforward computation one gets u u − u ˜ 2 Bu = i(T ∂x u Tu ) , | | − which is a pseudodifferential operator of order 0. This Lax pair formu- lation will be also a consequence of Proposition 6.1. The Lax pair formulation of the BO equation implies that at least formally, the spectrum of the operator Lu is left invariant by the so- lution map of the BO equation. For this reason, we want to analyze it in more detail. In addition of being selfadjoint, Lu has a compact resolvent. Furthermore note that for u = 0, L0 = D is a nonnegative operator and its spectrum consists of simple eigenvalues, which we list in increasing order, λn = n, n 0. Since Tu is D–compact it then fol- 2 ≥ lows that for any u Lr, Lu is bounded from below and its spectrum consists of eigenvalues∈ which are bounded from below and have finite multiplicities. We list these eigenvalues in increasing order and with their multiplicities so that λ (u) λ (u), n 0. n ≤ n+1 ∀ ≥ Our first result says that all eigenvalues of Lu are simple. More pre- cisely, we have the following Proposition 2.1. For any u L2 and n 0, the eigenvalues λ ∈ r ≥ n ≡ λn(u) satisfy λ λ +1 . n+1 ≥ n As a consequence, all eigenvalues of Lu are simple. Remark 2.2. Since for any u L2, L 1 1 = u 1 , it follows from ∈ r h u | i −h | i the variational characterization of λ0 that λ0 u 1 . Furthermore, if λ = u 1 , then f =1 is a normalized eigenfunction≤ −h | i corresponding 0 −h | i 0 to λ0. Since Lu1 = Πu, one has Πu = u 1 1, implying that u = u 1 . − h | i h | i The proof of Proposition 2.1 is based on properties of the shift oper- ator S and its commutator relations with Lu. More precisely, S is the 2 operator on L+, defined by S : L2 L2 , h Sh(x) = eixh(x) . + → + 7→ ON THE INTEGRABILITY OF THE BO EQUATION 10

Its adjoint with respect to the L2 inner product (1.2) is given by − S∗ = Te−ix , implying that the following identities hold

(2.3) S∗S = Id, SS∗ = Id 1 1 . −h·| i Furthermore, since for any f L2, Π(eixf)= SΠ(f)+ eixf 1 one has for any h H1 , ∈ h | i ∈ + (2.4) T (Sh)= ST h + uSh 1 1 . u u h | i Using that S∗1 = 0, one then gets

(2.5) S∗TuS = Tu and consequently, since DS = SD + S and hence S∗DS = D + Id,

(2.6) S∗LuS = Lu + Id . Moreover, combining DS = SD + S with (2.4), one sees that (2.7) L S = SL + S uS 1 1 . u u −h · | i With these preparations we are now ready to prove Proposition 2.1.

Proof. Apply the max-min formula for λn, 1 λn = max m(F ) , m(F ) := min Luh h : h H+ F ⊥ , h =1 , dim F =n {h | i ∈ ∩ k k } C 2 where, in the above maximum, F describes all subspaces of L+ of dimension n. If F is such a vector subspace, observe− that G(F ) := C 1 1 S(F ) is a vector subspace of dimension n+1, and that G(F )⊥ H+ ⊕ 1 ∩ is precisely S(F ⊥ H ). Therefore, we have ∩ + 1 m(G(F )) = min L Sh Sh : h H F ⊥ , h =1 . {h u | i ∈ + ∩ k k } From (2.6), we infer L Sh, Sh = L h h + h 2 and consequently h u i h u | i k k m(G(F )) = m(F )+1 .

Substituting this identity into the above formula for λn+1, we conclude that λ max m(G(F )) = λ + 1.  n+1 ≥ dim F =n n We shall prove in Section 5 that the eigenvalues λn are real analytic. At this stage, we establish that they are Lipschitz continuous.

Proposition 2.2. For every n 0, the functional λn is uniformly ≥ 2 T Lipschitz continuous on bounded subsets of Lr( ). Proof. For any h H1 , ∈ + 2 2 (2.8) L h, h L h, h = (u v), h u v h 4 |h u i−h v i| |h − | | i|≤k − kk kL and if in addition h = 1, one has by the Sobolev embedding theorem k k 2 1/2 1/2 (2.9) h 4 C (1 + D) h h C (1 + D)h h . k kL ≤ h | i ≤ h | i ON THE INTEGRABILITY OF THE BO EQUATION 11

Applying the max-min formula and (2.8) for v = 0, we obtain 1/2 1 λn(u) max min Dh + C u (1 + D) h h : h F ⊥ H+, h =1 ≤ dim F =n {h k kh | i ∈ ∩ k k } n + C u (1 + n)1/2 . ≤ k k Given a subspace F L2 with dim F = n, select h H1 n ker(L ⊂ + ∈ +∩⊕k=0 u− λ (u)Id) F ⊥ with h = 1. From (2.8) and (2.9) one then infers that k ∩ k k L h h L h h + C u v (1 + D)h h 1/2 . h v | i≤h u | i k − kh | i On the other hand, L h h λ (u), implying that h u | i ≤ n Dh h = uh h + λ (u) C u (1+ D)h h 1/2 + n + C u (1+ n)1/2 , h | i h | i n ≤ k kh | i k k yielding Dh h C˜(n +1+ u 2) . h | i ≤ k k We then conclude by the max-min formula that λ (v) λ (u)+ K u v (n +1+ u 2)1/2 n ≤ n k − k k k where K > 0 is an absolute constant.  In view of Proposition 2.1 we introduce for any n 1 and u L2, ≥ ∈ r (2.10) γn(u) := λn(u) λn 1(u) 1 , − − − and refer to it as the n–th gap length – see Proposition C.1 in Appendix C for an explanation of this terminology. It then follows that for any n 1 ≥ n (2.11) λ (u)= n + λ (u)+ γ (u) n + λ (u) . n 0 k ≥ 0 Xk=1 Next we obtain a useful characterisation for the vanishing of γn. For 2 2 any u Lr and n 0, we denote by fn a L -normalized eigenfunction of L for∈ the eigenvalue≥ λ . Hence f H1 and f = 1. In view of u n n ∈ + k nk Proposition 2.1, fn is determined uniquely up to a phase factor for any n 0. Later we will show that this indeterminacy can be removed in ≥ such a way that the eigenfunctions fn depend real analytically on u. 2 Lemma 2.1. For any u Lr, f0 1 =0. Furthermore, for any n 1 and u L2, the following∈ statementsh | i6 are equivalent: ≥ ∈ r (1) fn 1 =0 , (2)Lu(Sfn 1)= λnSfn 1 , (3) Sfn 1 u =0 , (4)γn =0 . h | i − − h − | i Proof. Assume that, for some n 0, we have f 1 = 0. Then f = ≥ h n| i n Sgn for some gn, and, by (2.7), λ Sg = L Sg = S(L g + g ) Sg u 1 . n n u n u n n −h n| i Applying S∗ to both sides, we infer, in view of (2.3), that L g =(λ 1)g . u n n − n In the case n = 0, the latter identity contradicts the fact that λ0 is the smallest eigenvalue of L and hence f 1 = 0. In the case n 1, the u h 0| i 6 ≥ ON THE INTEGRABILITY OF THE BO EQUATION 12 identity implies that λn 1 = λn 1 and that gn is collinear to fn 1, showing that (1) implies− (2). − − Applying (2.7) to fn 1 immediately yields the equivalence of (2) and (3), and that (3) implies− (4). Finally, assume that γ = 0 for some n 1. Then, again by applying n ≥ (2.7) to fn 1, − LuSfn 1 =(λn 1 + 1)Sfn 1 Sfn 1 u 1= λnSfn 1 Sfn 1 u 1 . − − − −h − | i − −h − | i Taking the inner product with fn of both sides of the latter identity, one gets

uSfn 1 1 1 fn =0 , h − | ih | i which means that either (1) holds, or (3) holds, hence (2) holds, so that Sfn 1 is collinear to fn, implying (1).  −

In a next step we show that the eigenfunctions fn can be normalized in such a way that they depend continuously in u. First we need to make some preliminary considerations. Lemma 2.2. For any u L2 and k,p 0 ∈ r ≥ (λ λ 1) f Sf = f 1 u Sf . p − k − h p| ki −h p| ih | ki Proof. Recall from (2.7) that for any k 0, ≥ (2.12) L (Sf )=(λ + 1)Sf Sf u 1 . u k k k −h k | i Computing the inner product of fp with both sides of (2.12) leads to the claimed formula. 

2 Lemma 2.3. For any u Lr and n 1, fn Sfn 1 =0 . ∈ ≥ h | − i 6 2 Proof. Assume that for some u Lr and n 1, fn Sfn 1 = 0. From Lemma 2.2, we infer ∈ ≥ h | − i

fn 1 u Sfn 1 =0 . h | ih | − i By Lemma 2.1, this implies that Sfn 1 and fn are collinear, contra- − dicting the assumption fn Sfn 1 = 0.  h | − i Lemma 2.1(i) and Lemma 2.3 allow to define an orthonormal basis of eigenfunctions of Lu as follows. Definition 2.1. For any u L2, f is defined to be the L2 normalized ∈ r 0 − eigenfunction of Lu corresponding to the eigenvalue λ0, uniquely deter- mined by the condition 1 f > 0 . h | 0i The L2 normalized eigenfunction f of L , corresponding to the eigen- − 1 u value λ1, is then defined to be the eigenfunction uniquely determined by the condition f Sf > 0 . h 1| 0i ON THE INTEGRABILITY OF THE BO EQUATION 13

Assume that for any 0 k n, the eigenfunctions fk have been de- fined. Then the L2 normalized≤ ≤ eigenfunction of L , corresponding to − u the eigenvalue λn+1, is defined by the condition f Sf > 0 . h n+1| ni 2 In the remaining part of the paper, for any u Lr, fn( ,u), n 0, 2 ∈ · ≥ will always denote the orthonormal basis of L+ introduced in Definition 2.1.

Remark 2.3. It is easy to see that the map u fn( ,u) is continuous 2 1 (k) 2 (k)7→ · (k) from L to H . Indeed, if u u in L , fn := f ( ,u ) satisfies r + → r n · Df (k) 2 = λ (u(k)) f (k) Df (k) + u(k)f (k) Df (k) k n k n h n | n i h n | n i (k) which, by Proposition 2.2 and Sobolev inequalities, implies that fn is 1 bounded in H+ as k . By Rellich’s theorem and Proposition 2.2, → ∞ (k) 1 we infer that any weak * limit point gn of (fn )k 1 in H+ satisfies ≥ 2 (k) 2 Lugn = λn(u)gn, gn = 1 and Dgn = limk Dfn . Further- k k k k →∞ k k more, the conditions 1, f0 0 , fn+1 Sfn 0 are also stable by h 1 i ≥ h | i ≥ weak * convergence in H+, implying that g0 = f0 (Lemma 2.1) and gn = fn for any n 1 (Lemma 2.3). In Section 5, we will prove that u f ( ,u) is a real≥ analytic map from L2 to H1 . 7→ n · r + The complex numbers 1 fn , n 1, appropriately rescaled, are our candidates for the complexh | Birkhoffi ≥ coordinates. Note that for any 2 u Lr, one can express Πu in terms of 1 fn , n 0, as follows. Since ∈ 2 h | i ≥ fn, n 0, is an orthonormal basis of L+, one has Πu = n∞=0 Πu fn fn. When≥ combined with the identity h | i P (2.13) λ 1 f = 1 L f = L 1 f = Πu f = u f nh | ni h | u ni h u | ni −h | ni −h | ni one is led to the following trace formula

∞ (2.14) Πu = λ 1 f f − nh | ni n n=0 X and hence

∞ (2.15) Πu 2 = λ2 1 f 2 k k n |h | ni| n=0 X and

∞ (2.16) u 1 = λ 1 f 2 . h | i − n|h | ni| n=0 X 2 Similarly, the expansion 1 = ∞ 1 f f leads to ∞ 1 f =1 . n=0h | ni n n=0 |h | ni| We record the following property of the eigenfunctions fn. P P Lemma 2.4. For any u L2 with γ = 0 for a given n 1, f = ∈ r n ≥ n Sfn 1. − ON THE INTEGRABILITY OF THE BO EQUATION 14

Proof. By Lemma 2.1(iii), for any n 1 with γn = 0, the function 2 ≥ Sfn 1 is an L normalized eigenfunction of Lu, corresponding to the − − eigenvalue λn and satisfies the condition fn Sfn 1 > 0.  h | − i In analogy with the notion of a finite gap potential, introduced in the context of the Korteweg-de Vries equation and the nonlinear Schr¨odinger equation, we define the corresponding notion in the con- text of the BO equation. Definition 2.2. A potential u L2 is said to be a finite gap potential ∈ r if J(u) = n 1 : γn(u) > 0 N is finite. In case J(u) consists of one element{ only,≥ u is referred} to ⊂ as a one gap potential. Note that for any finite gap potential u, the expansion (2.14), is a finite sum. Finite gap potentials play a similar role as trigonometric polynomials in harmonic analysis. For any given subset J N intro- duce ⊂ G := u L2 : γ (u) > 0 j J; γ (u)=0 j N J . J { ∈ r j ∀ ∈ j ∀ ∈ \ } Elements in the set GJ are referred to as J gap potentials. We refer to Section 7 for a detailed studied of J− gap potentials with J = 1, , N and to Appendix B for a description− of one gap po- tentials.{ ··· } Remark 2.4. One might ask if there exists a potential u L2 with the ∈ r property that γn(u) > 0 for every n 1. We claim that u(x)=2α cos x is such a potential for every α R ≥0 . Indeed, we have for any n 1, ∈ \{ } ≥ ix Sfn 1 u = α Sfn 1 e = α fn 1 1 . h − | i h − | i h − | i Hence, by Lemma 2.1, fn 1 = 0 if and only if fn 1 u = 0. Since f 1 =0, we concludeh that| if 1 =0 for every nh. − | i h 0| i 6 h n| i 6

We finish this section with a brief discussion of the symmetry induced by translation. For any τ R, denote by Qτ the linear isometry on 2 ∈ Lr, given by the translation by τ, Q : L2 L2, , Q u(x) := u(x + τ) . τ r → r τ Lemma 2.5. For any τ R, u L2, and n 0, one has ∈ ∈ r ≥ inτ λn(Qτ u)= λn(u) , fn(x, Qτ u)= e− fn(x + τ,u) . As a consequence, (2.17) 1 f ( , Q u) = einτ 1 f ( ,u) , n 0 . h | n · τ i h | n · i ∀ ≥ R 2 Proof. For any τ , u Lr, and n 0, note that fn fn( ,u) satisfies ∈ ∈ ≥ ≡ ·

LQτ u(Qτ fn)= Qτ (Lufn)= λn(u)Qτ fn . ON THE INTEGRABILITY OF THE BO EQUATION 15

It implies that λ (Q u)= λ (u) and that Q f ( ,u) is an L2 normalized n τ n τ n · − eigenfunction of LQτ u, corresponding to the eigenvalue λn. Further- more, 1 f ( + τ,u) = 1 f ( ,u) > 0 h | 0 · i h | 0 · i implying that f ( , Q u)= f ( + τ,u). Similarly, for any n 1, 0 · τ 0 · ≥ inτ i(n 1)τ iτ e− fn( +τ,u) Se− − fn 1( +τ,u) = fn( +τ,u) e Sfn 1( +τ,u) h · | − · i h · | − · i which equals fn( ,u) Sfn 1( ,u) > 0. Arguing by induction one then h · | − · inτ i concludes that f ( , Q u)= e− f ( + τ,u) as claimed.  n · τ n · 3. Generating function and trace formulae 2 One of the main results of this section are trace formulas for 1 fn , stated in Corollary 3.1. They will be used to define our candidates|h | i| of (complex) Birkhoff coordinates by appropriately scaling 1 fn . A key ingredient for the proof of these trace formulas is the generah | iting function λ. We remark that λ plays a role for the normal form theory ofH the BO equation comparableH to the one of the discriminant in the normal form theory of the KdV equation. 2 For any u Lr, the generating function is the meromorphic function λ (u),∈ defined by 7→ Hλ 1 (3.1) (u) := (L + λId)− 1 1 . Hλ h u | i Note that λ(u) is holomorphic on C λn(u): n 0 and might have simpleH poles at λ = λ (u), n \0. {− Substituting≥ the} expansion − n ≥ 1= ∞ 1 f f into the expression for (u) one obtains n=0h | ni n Hλ 2 P ∞ 1 fn (3.2) (u)= |h | i| . Hλ λ + λ n=0 n X Proposition 3.1. For any u L2, the following identities hold: ∈ r 1 ∞ γ (i) (u)= 1 n , Hλ λ + λ − λ + λ 0 n=1 n Y  ∞ ∞ (ii) u 2 u 1 2 =2 nγ , u 1 = λ γ . k k −h | i n h | i − 0 − n n=1 n=1 X X Hence γ (u)= 1 ℓ1 and both u 1 and u 2 are spectral invariants. n n n h | i k k Proof. (i) The proof makes use of the shift operator S and in particular identity (2.6). Assume that λ is real and λ λ0 +1, so that Lu +λId : H1 L2 is invertible. By (2.6), one has ≥ − + → + (3.3) S∗(Lu + λId)S = Lu +(λ + 1)Id , 1 2 implying that the operator S∗(Lu + λId)S : H+ L+ is invertible. To compute its inverse we need the following → ON THE INTEGRABILITY OF THE BO EQUATION 16

1 2 Lemma 3.1. Assume that A : H+ L+ is positive and selfadjoint so → 1 2 that A and S∗AS are invertible and A− is positive on L+. Then 1 1 1 S∗A− 1 1 (S∗AS)− = S∗A− S h·| iS∗A− 1 . − A 11 1 h − | i 2 1 Proof of Lemma 3.1. Consider f L+ and let h := (S∗AS)− f. Ap- plying S to both sides of the equation∈

S∗ASh = f , we infer from (2.3) that ASh = Sf + ASh 1 1 , h | i hence 1 1 (3.4) Sh = A− Sf + ASh 1 A− 1 . h | i In particular, we have 1 1 A− Sf 1 = ASh 1 A− 1 1 . h | i −h | ih | i 1 1 Since A− is positive, A− 1 1 > 0 and hence we obtain h | i 1 1 A− Sf 1 f S∗A− 1 ASh 1 = h | i = h | i . h | i − A 11 1 − A 11, 1 h − | i h − i Substituting this identity into the equation (3.4) for Sh and applying S∗ to both sides, we obtain the claimed statement. 

Let us now go back to the proof of item(i). Since Lu +λId is selfadjoint and positive we can apply Lemma 3.1 to A := Lu + λId to get, in view of (3.3), 1 1 (L +(λ + 1)Id)− S∗(L + λId)− S u − u 1 S∗(Lu + λId)− 1 1 (3.5) = h·| iS∗(L + λId)− 1. − (L + λId) 11 1 u h u − | i Since the latter operator is of rank 1, it is of and hence so 1 1 is the operator (Lu +(λ +1)Id)− S∗(Lu + λId)− S. To compute its trace we write it as a sum of two operators,− 1 1 (L +(λ + 1)Id)− S∗(L + λId)− S = A + A u − u 1 2 where 1 1 A := (L +(λ + 1)Id)− (L + λId)− , 1 u − u 1 1 A := (L + λId)− S∗(L + λId)− S. 2 u − u In view of Proposition 2.1, A1 is of trace class and

∞ 1 1 TrA = . 1 λ + λ +1 − λ + λ n=0 n n X   ON THE INTEGRABILITY OF THE BO EQUATION 17

Since by (3.5), A1 + A2 has rank 1 and hence is of trace class, the operator A2 is also of trace class. Computing its trace with respect to inx 2 the orthonormal basis (e )n 0 of L+, we obtain ≥ 1 1 1 TrA = Tr (L + λId)− S∗(L + λId)− S = (L + λId)− 1 1 . 2 u − u h u | i On the other hand, by (3.5), Tr(A1 + A2) equals 1 1 2 S∗(Lu + λId)− 1 1 S∗(Lu + λId)− 1 Tr h·| iS∗(L + λId)− 1 = k k − (L + λId) 11 1 u − L + λId) 11 1 h u − | i h u − | i and hence by (2.3)  1 2 (Lu + λId)− 1 1 Tr(A + A )= k k + (L + λId)− 1 1 . 1 2 − (L + λId) 11 1 h u | i h u − | i Altogether we have proved that 1 2 ∞ 1 1 (L + λId)− 1 = k u k λ + λ +1 − λ + λ − (L + λId) 11 1 n=0 n n u − X   h | i or (3.6) 2 ∞ 1 fn |h | i| 2 ∞ 1 1 (λn + λ) d = n=0 = log (u) . X 2 λ λn + λ +1 − λn + λ − ∞ 1 fn dλ H n=0   |h | i| X λ + λ n=0 n X 1 Isolating the term (λ0+λ)− on the left hand side of the latter identity and then rewriting− the remaining sum as a telescopic series yields

d 1 ∞ γ (3.7) log (u)= + n dλ Hλ −λ + λ (λ + λ + 1)(λ + λ ) 0 n=1 n 1 n X − or, written in a slightly more convenient form,

d ∞ γn log (λ0 + λ) λ(u) = . dλ H (λ + λn 1 + 1)(λ + λn) n=1 −  X Integrating both sides of the latter identity from λ to + , we infer ∞ ∞ γ log (λ + λ) (u) = log 1 n , − 0 Hλ − − λ + λ n=1 n  X  implying that the formula of item (i) holds for any λ> λ0 +1 . By analyticity, it is then valid for any λ C λ : n |0 −. | ∈ \ {− n ≥ } (ii) The claimed trace formulas are obtained by expanding λ at λ = . For ε := 1 > 0 with λ> λ +1 , define H ∞ λ | − 0 | 2 1 ∞ 1 fn (3.8) ˜ε := 1 = |h | i| , H εH ε 1+ ελ n=0 n X ON THE INTEGRABILITY OF THE BO EQUATION 18 so that in view of (2.15)-(2.16), ˜ = 1, d ˜ = u 1 , and H0 dε |ε=0Hε h | i d2 (3.9) ˜ =2 Πu 2 = u 2 + u 1 2 dε2 |ε=0Hε k k k k h | i and hence, using these identities, d d log ˜ = ˜ = u 1 dε|ε=0 Hε dε|ε=0Hε h | i and d2 d d2 (3.10) log ˜ = ˜ 2 + ˜ = u 2 . dε2 |ε=0 Hε − dε|ε=0Hε dε2 |ε=0Hε k k Furthermore, identity (3.7) becomes 

d λ ∞ γ (3.11) log ˜ = 0 + n . − dε Hε 1+ ελ (1 + ε(λ + 1))(1 + ελ ) 0 n=1 n 1 n X − Passing to the limit as ε 0 in both sides, we obtain that the series → n∞=1 γn is summable — notice that γn 0 and λn 0 for n large enough — and that ≥ ≥ P ∞ (3.12) u 1 = λ + γ . −h | i 0 n n=1 X Taking the derivatives of both sides of (3.11) with respect to ε, one 2 concludes that d log ˜ equals dε2 Hε 2 ∞ λ0 γn λn 1 +1 λn − 2 + + . (1 + ελ0) (1 + ε(λn 1 + 1))(1 + ελn) 1+ ε(λn 1 + 1) 1+ ελn n=1 − − X  Passing to the limit as ε 0 and using (3.10), one concludes that → ∞ 2 2 u = λ0 + γn(λn + λn 1 + 1) . k k − n=1 X From identity (3.12) and the definition of γ we have for every n 0, n ≥ ∞ (3.13) λ = n u 1 γ . n −h | i − k k=Xn+1 As a consequence,

∞ λn + λn 1 +1=2n 2 u 1 γn 2 γk , − − h | i − − k=Xn+1 so that

∞ ∞ ∞ u 2 = λ2 +2 nγ 2 u 1 γ γ 2 , k k 0 n − h | i n − n n=1 n=1 n=1 X X X  ON THE INTEGRABILITY OF THE BO EQUATION 19 which, in view of (3.12), leads to

∞ u 2 = u 1 2 +2 nγ , k k h | i n n=1 X whence the claimed formula for u 2.  k k Remark 3.1. (i) It follows from the proof of Proposition 3.1 that λ(u) 1 H 1 equals the trace of the rank one operator (Lu+λId)− S∗(Lu+λId)− S. 2 − 2 (ii) Writing u Lr as v + c with c = u 1 and hence v Lr,0, one computes ∈ h | i ∈ 2π 1 1 (u)= (v) c v2dx c3 . H H − 2π − 3 Z0 By the computations in the proof of Proposition 3.1 one has (cf. also (8.1)) 1 d3 d2 (v)= ˜ (v) , v 2 = ˜ (v) , H −6 dε3 |ε=0Hε k k dε2 |ε=0Hε and c = u 1 = d ˜ (u). h | i dε |ε=0Hε The identity of Proposition 3.1(i) can be used to obtain product representations for 1 f 2, which will be used later to appropriately |h | ni| scale 1 fn in order to obtain our candidates for (complex) Birkhoff coordinates.h | i Corollary 3.1. For any u L2 and n 1, ∈ r ≥ ∞ γ 1 f 2 = κ , κ := 1 p , |h | 0i| 0 0 − λ λ p=1 p 0 Y  −  2 1 γp 1 fn = γnκn , κn := 1 . |h | i| λn λ0 − λp λn 1 p=n − ≤Y6  −  In particular, 1, f 2, n 0, are spectral invariants. |h ni| ≥ Remark 3.2. The formulas of Corollary 3.1 provide a new proof of the fact that 1 f = 0 and that, for n 1, 1 f = 0 if and only if h | 0i 6 ≥ h | ni γn =0. Proof. The claimed product representations are obtained by computing with the help of the identity of Proposition 3.1(i) the residue of λ(u) at λ = λ for any n 0. H  − n ≥ As an application of Proposition 3.1 and Corollary 3.1 we analyze the set Iso(u) of an arbitrary potential u L2, defined as ∈ r Iso(u) := v L2 : λ (v)= λ (u) n 0 . { ∈ r n n ∀ ≥ } Note that Iso(u) is closed by Proposition 2.2 and that Iso(u)= v L2 : (v)= (u) λ C λ (u): n 0 . { ∈ r Hλ Hλ ∀ ∈ \{ n ≥ }} ON THE INTEGRABILITY OF THE BO EQUATION 20

First we need to make some preliminary considerations. Recall that 1 by Proposition 3.1, (γn(u))n 1 is in the weighted ℓ sequence space ℓ1,1(N, R) and γ (u) 0 for any≥ n 1. Hence the map− Γ n ≥ ≥ 2 1,1 (3.14) Γ : Lr,0 + ,u (γn(u))n 1 → C 7→ ≥ is well defined. Here 1,1 denotes the positive cone in ℓ1,1(N, R), C+ 1,1 1,1 N R + := (rn)n 1 ℓ ( , ): rn 0 n 1 . C { ≥ ∈ ≥ ∀ ≥ } Proposition 3.2. The map Γ is proper.

(k) 2 Proof. Let (u )k 1 be a sequence in Lr,0 so that ≥ (k) 1,1 Γ(u ) γ =(γn)n 1 + as k . → ≥ ∈ C →∞ By Proposition 3.1(ii) and the assumption,

∞ ∞ u(k) 2 =2 nγ (u(k)) 2 nγ as k , k k n → n →∞ n=1 n=1 X X and thus we may assume, after extracting a subsequence if needed, that (k) 2 u ⇀u weakly in Lr,0. The proof will be complete if we establish that this convergence is strong, or equivalently that

∞ u 2 =2 nγ . k k n n=1 X Since (again by Proposition 3.1(ii))

∞ ∞ (k) (k) L (k) λ (u )= γ (u ) γ , u ≥ 0 − n → − n n=1 n=1 X X we infer that there exists c > λ0 +1 so that for any k 1 and λ c, | − | ≥ ≥ 1 2 L (k) + λId : H L u + −→ + is an isomorphism whose inverse is bounded uniformly in k. Therefore (k) 1 wλ := (Lu(k) + λId)− [1] , 1 is a well defined, bounded sequence of H+. Let us choose an arbitrary countable subset Λ of [c, + ) with a cluster point. By a diagonal ∞ (k) (k) procedure, we extract a subsequence of wλ , again denoted by wλ , so (k) 1 that for every λ Λ, wλ converges weakly in H+ to some element w H1 as k ∈ . By Rellich’s theorem we infer that, weakly in L2 , λ ∈ + →∞ + (k) (L (k) + λId)w ⇀ (L + λId)w as k . u λ u λ →∞ (k) But since (Lu(k) + λId)wλ = 1 for any k 1 it then follows that for every λ Λ, (L + λId)w = 1 and hence≥ by the definition of (u), ∈ u λ Hλ (u(k))= w(k) 1 w 1 = (u) λ Λ. Hλ h λ | i → h λ| i Hλ ∀ ∈ ON THE INTEGRABILITY OF THE BO EQUATION 21

On the other hand, for every n 0, ≥ ∞ ∞ λ (u(k))= n γ (u(k)) n γ as k n − j → − j →∞ j=n+1 j=n+1 X X uniformly with respect to n 0. Hence, setting λ := n ∞ γ ≥ n − j=n+1 j for any n 0 one infers that for any λ [c, ), ≥ ∈ ∞ P 1 ∞ γ (u(k)) 1 ∞ γ 1 n 1 n . λ (u(k))+ λ − λ (u(k))+ λ → λ + λ − λ + λ 0 n=1 n 0 n=1 n Y   Y   By Proposition 3.1(i), it then follows that

1 ∞ γ (u)= 1 n Hλ λ + λ − λ + λ 0 n=1 n Y   for any λ Λ and hence by analyticity, for any λ [c, + ). Thus ∈ 2 ∈ ∞ (3.9) applies and we conclude that u =2 ∞ nγ .  k k n=1 n Our result on Iso(u) then reads as follows.P

2 2 Proposition 3.3. For every u Lr, Iso(u) is a compact subset of Lr. ∈ 2 In fact, Iso(u) u 1 is a compact subset of the sphere in Lr,0 of radius u u 1 , centered−h | i at 0. In particular, for any c R, Iso(c) consists ofk the−h constant| ik potential c only. ∈ Proof. Let v Iso(u). Since by Proposition 3.1(ii) ∈ ∞ v 1 = λ (u) γ (u)= u 1 h | i − 0 − n h | i n=1 X and γ (u u 1 )= γ (u) by the definition of γ , n 1, one has n −h | i n n ≥ 1 v u 1 Γ− (γn(u))n 1 . −h | i ∈ { ≥ } Using that Γ is proper by Proposition 3.2 and that u u 1 is a spectral invariant, it follows that Iso(u u 1 ) = Iso(k −hu) | iku 1 is − h | i 2 −h | i contained in a compact subset of the sphere in Lr,0 of radius u u 1 , centered at 0 . Since Iso(u) is closed one then concludesk that− hIso(| uik) u 1 and in turn Iso(u) are compact.  −h | i 4. Complex Birkhoff coordinates In this section we introduce our candidates of complex Birkhoff co- ordinates, define the corresponding Birkhoff map Φ, and discuss first properties of Φ. For any u L2 and n 0, define ∈ r,0 ≥ 1 fn (4.1) ζn(u) := h | i √κn ON THE INTEGRABILITY OF THE BO EQUATION 22 where we recall that κ = 1 γp and, for any n 1, 0 p 1 λp λ0 ≥ − − ≥   1Q γp κn = 1 . λn λ0 − λp λn 1 p=n − ≤Y6  −  The functionals ζ , n 1, are our candidates for complex Birkhoff n ≥ coordinates of the BO equation. By Corollary 3.1, ζ0 = 1 and for any n 1, ζ 2 = γ . ≥ | n| n Proposition 4.1. The map 2 1/2 (4.2) Φ : Lr,0 h+ , u (ζn(u))n 1 → 7→ ≥ is continuous and proper.

1 1 1/2 Proof. Since by Proposition 3.1, γn = n ℓn and in addition ζn = γn , 1/2 2 | | | | it follows that Φ(u) h+ for any u Lr,0. Furthermore, since for any n 1, we know from∈ Proposition ∈2.2, Remark 2.3 and Proposition 3.1≥(ii), that 1, f and κ depend continuously on u, we infer that h ni n ζn depends continuously on u. Since by Proposition 3.1 (ii), Φ is a 2 bounded map, it then follows that whenever uk u in Lr,0, Φ(uk) 1→/2 converges to Φ(u) weakly in the Hilbert space h+ . Since, again by Proposition 3.1, Φ(v) 2 =2 v 2 , v L2 , k k1/2 k k ∈ r,0 we infer that Φ(uk) 1/2 Φ(u) 1/2 and hence that Φ(uk) Φ(u) 1/2 k k → k k → in h+ . This shows that Φ is continuous. Finally, by Proposition 3.2, the map Γ is proper and so is Φ.  In a next step we want to show that the map Φ is one-to-one. We prove this fact by showing that any potential u L2 can be expressed ∈ r,0 in terms of the complex numbers ζn(u), n 1. First note that the function Πu extends as a holomorphic function≥ to the unit disc z < 1, which by a slight abuse of notation, we denote by Πu(z). It is given| | by its Taylor expansion at z =0,

∞ ∞ ∞ k k k k k Πu(z)= u(k)z = Πu S 1 z = (S∗) Πu 1 z . h | i h | i Xk=0 Xk=0 Xk=0 2 2 Denote by S∗ theb norm of the operator S∗ : L+ L+. Since k k k 1 → z S∗ = z < 1, one has ∞ (zS∗) =(Id zS∗)− and hence | |k k | | k=0 − 1 (4.3) Πu(z)= (Id zS∗)− Πu 1 , z < 1. h P− | i ∀| | Denote by M M(u) the matrix representation of the operator S∗ : 2 2 ≡ L+ L+ in the basis (fn)n 0, → ≥ (4.4) M =(Mnp)n,p 0, Mnp Mnp(u)= S∗fp fn = fp Sfn . ≥ ≡ h | i h | i Then (4.3) together with the formula (2.14) for Πu yields the following ON THE INTEGRABILITY OF THE BO EQUATION 23

Lemma 4.1. For any u L2 and any z C with z < 1, ∈ r ∈ | | 1 (4.5) Πu(z)= (Id zM(u))− X(u) Y (u) 2 , h − | iℓ with X(u) and Y (u) being the column vectors

(4.6) X(u) := (λp 1 fp )p 0 ,Y (u) := ( 1 fn )n 0 . − h | i ≥ h | i ≥ Recall that by (2.3), MM ∗ = Id and hence Id zM is invertible for any z < 1. − | | Proposition 4.2. The map Φ is one-to-one. Proof. First recall that that for any u L2 , 1 f = √κ with κ = ∈ r,0 h | 0i 0 0 1 γp and for any n 1 p 1 λp λ0 ≥ − − ≥   Q 1 γp 1 fn = √κnζn , κn = 1 h | i λn λ0 − λp λn 1 p=n − ≤Y6  −  2 where γ = ζ and λ = n ∞ γ . Therefore the components p | p| n − k=n+1 k 1 fn of Y (u) and the components λp 1 fp of X(u) can be expressed h | i P − h | i in terms of the ζk, k 1. We claim that the coefficients Mnp can also be expressed in this≥ way. Indeed, for any given n 0 one argues as ≥ follows: if ζn+1 = 0, then γn+1 = 0 and hence fn+1 = Sfn, implying that Mnp = δp,n+1.

If ζn+1 = 0, then γn+1 = 0 and 1 fn+1 = 0. We then use that by Lemma 6 2.2 for any p 0,6 h | i 6 ≥ (λ λ 1)M = u Sf f 1 p − n − np −h | nih p| i and hence in particular for p = n+1, γn+1 fn+1 Sfn = u Sfn fn+1 1 . Combining these two identities then yieldsh | i −h | ih | i γ f Sf f 1 (4.7) M = n+1h n+1| ni h p| i . np f 1 λ λ 1 h n+1| i p − n − 2 2 Since fn+1 Sfn > 0 the identity 1 = Sfn = p∞=0 Mnp then reads h | i k k | | P γ2 ∞ f 1 2 (4.8) 1 = f Sf 2 n+1 |h p| i| h n+1| ni f 1 2 (λ λ 1)2 n+1 p=0 p n |h | i| X − − Altogether we thus have shown that all terms in the formula (4.7) for Mnp can be expressed in terms of ζk, k 1, and hence also Πu(z) for z < 1 can be expressed in this way. Since≥ Πu(z), z < 1, completely |determines| u we have shown that Φ is one-to-one. | |  Remark 4.1. By the definition (3.2) of and formula (4.8), Hλ 2 2 γn+1 1= fn+1 Sfn ′ λn 1 . h | i f 1 2 − H− − |h n+1| i|  ON THE INTEGRABILITY OF THE BO EQUATION 24

In view of Proposition 3.1(i) this leads to fn+1 Sfn = √µn+1 with h | i γ 1 p γn+1 − λp λn+1 (4.9) µn+1 := 1  −γ  . − λn+1 λ0 p − 1 p=n+1 1   ≤Y6 − λp λn 1  − −  We finish this section with discussing two symmetry properties of 2 2 the map Φ. For any u Lr, denote by u the element in Lr, given by u (x) := u( x). Note that∈ the fixed points∗ of the involution u u ∗ − 2 7→ ∗ are the even functions in Lr. Proposition 4.3. For any u L2, ∈ r λn(u )= λn(u) , fn(x, u )= fn( x, u) , n 0 . ∗ ∗ − ∀ ≥ As a consequence, for any u L2 , ∈ r,0

ζn(u )= ζn(u) n 1 . ∗ ∀ ≥ Hence u is even if and only if for any n 1, ζ (u) R. ≥ n ∈ Proof. Given any u L2 and n 0, let g (x) := f ( x, u) and, ∈ r,0 ≥ n n − to shorten notation, write λn for λn(u) and fn(x) instead of fn(x, u). Taking the complex conjugate of both sides of the identity Lufn =

λ f and then evaluate them at x one obtains L ∗ g = λ g . Hence n n − u n n n λn = λn(u ) and gn is an eigenfunction of Lu∗ , corresponding to the ∗ eigenvalue λn(u ). One verifies inductively that the eigenfunctions gn, n 0, satisfy the∗ normalisation conditions of Definition 2.1, implying ≥ that gn(x, u)= fn(x, u ). ∗ Since κn(u), n 0, are spectral invariants of Lu one has κn(u )= κn(u) and in turn for≥ any n 1, ∗ ≥ 1 1 ζn(u )= 1 fn( ,u ) = 1 gn = ζn(u) . ∗ ∗ κn(u)h | · i κn(u)h | i It then follows thatp for u L even,pζ (u) R for any n 1. ∈ r,0 n ∈ ≥ Conversely, if for a given u Lr,0, ζn(u) R for any n 1, then Φ(u)=Φ(u ). Since Φ is one-to-one∈ we conclude∈ that u = u≥.  ∗ ∗ 2 The second result concerns potentials in Lr,0 which are 2π/K pe- 2,K − 2 riodic for some integer K 2. Let Lr,0 denote the subspace of Lr,0, ≥ 2,K 2,K consisting of such elements and by L , L+ the corresponding sub- 2 2 2,K (K) spaces of L , L+. For any u Lr,0 , let Lu be the Lax operator (K) 2,K ∈ (K) i∂x Tu , acting on L+ where Tu denotes the Toeplitz oper- − − (K) (K) (K) ator given by Tu (f) = Π (uf) and Π is the Szeg˝oprojector 2,K 2,K (K) L L+ . The spectrum of Lu is given by a sequence of eigen- → (K) values which we list in increasing order, (λn (u))n 0. Following the ≥ ON THE INTEGRABILITY OF THE BO EQUATION 25 arguments of the proof of Proposition 2.1 with the operator S replaced by (K) 2,K 2,K iKx S : L+ L+ , f e f (K) (K→) 7→ one verifies that λn (u) λn 1(u)+ K for any n 1. Since all the ≥ − (K)≥ (K) eigenvalues of Lu are simple and any eigenvalue λn of Lu is also an eigenvalue of Lu it follows by a simple homotopy argument applied to (K) the path t tu, 0 t 1, that λn (u)= λ (u) for any n 0. 7→ ≤ ≤ nK ≥ 2 Proposition 4.4. Assume that u Lr,0 is 2π/K periodic for some integer K 2. Then for any (n, k)∈with n 0 and−1 k K 1, ≥ ≥ ≤ ≤ − ikx λnK+k(u)= λnK(u)+ k , fnK+k(x, u)= e fnK(x, u).

Hence for such pairs (n, k), γnK+k(u)=0 and thus ζnK+k(u)=0. 2 Proof. Let u Lr,0 be 2π/K periodic for some integer K 2. Since ∈ (K) − ≥ for any n 0, λn (u)= λnK(u), and in view of Lemma 2.4 it suffices to prove that≥ for any (n, k) with n 0, 1 k K 1 one has ≥ ≤ ≤ − (4.10) λnK+k(u)= λnK(u)+ k . (K) (K) To verify the latter identity, denote by fn the eigenfunction of Lu (K) corresponding to the eigenvalue λn . (K) (K) (K) (K) Since u and fn are both 2π/K periodic one has Π(ufn ) = Π (ufn ) ikx (K) ikx (K) and hence for any 1 k K 1, e Π(ufn )=Π(e ufn ) imply- ing that ≤ ≤ − i∂ (eikxf (K)) Π(K)(eikxuf (K))=(λ(K) + k)eikxf (K) − x n − n n n By the homotopy argument mentioned above, (4.10) then follows.  Remark 4.2. Conversely, any element u L2 with ζ (u)=0 for ∈ r,0 nK+k any (n, k) with n 0 and 1 k K 1, is 2π/K periodic. Indeed, (K) ≥2,K 1/≤2 ≤ −(K) − denote by Φ : Lr,0 h+ ,u (ζn (u))n 1 the Birkhoff map on ≥ 2,K (K) → 7→ Lr,0 , where (ζn (u))n 1 are the Birkhoff coordinates of u, viewed as 2,K ≥ an element in Lr,0 . Arguing as in the proof of Theorem 1 concerning the case K = 1 (cf Section 8) one infers that Φ(K) is bijective. Hence 2,K (K) there exists v Lr,0 so that Φ (v)=(ζnK(u))n 1. By Proposition 4.4 ∈ ≥ it then follows that Φ(v)=(ζn(u))n 1 and hence by the uniqueness of Φ one concludes that v = u. ≥

5. Analyticity and Gradients

In this section, we establish that the Birkhoff coordinates ζn are real analytic functions, and that their gradients ζn are real analytic maps 1 ∇ with values in H . Hence the Gardner brackets ζn, F are well defined { } 2 real analytic functions for every real analytic functional F on Lr,0. As a first result we establish that the eigenvalues λn are real analytic. Note that for any v L2 ( L2(T, C)), the Lax operator L = i∂ T ∈ ≡ v − x− v ON THE INTEGRABILITY OF THE BO EQUATION 26

2 1 on L+ with domain H+ is no longer selfadjoint, but it is still a closed unbounded operator with compact resolvent. Hence its spectrum con- sists of a sequence of complex eigenvalues, each with finite multiplicity. For any r> 0, λ C and for any ε> 0, u L2, let ∈ ∈ r D (λ)= z C : z λ 0 so ∈ r ≥ ≥ n that for any v Bεn (u) and 0 k n, # spec(Lv) Drk (τk(u)) =1 and ∈ ≤ ≤ ∩  spec(L ) Box (u) D (τ (u)) . v ∩ N,n ⊂ rp(u) p 0 p n ≤[≤

The unique eigenvalue of Lv in Drk(u)(τk(u)), denoted by λk(v), is real analytic on B (u). It then follows that for any 1 k n, γ is an εn ≤ ≤ k analytic functional on Bεn (u) as well.

Proof. For simplicity, we do not indicate the dependence of rk and τk on u. In a first step we prove that there exists δn > 0 so that for any v B (u) ∈ δn spec(L ) Box (u) D (τ ) . v ∩ N,n ⊂ rk k 0 k n ≤[≤ Assume to the contrary that such a δn > 0 does not exist. Then 2 there exists a sequence (vℓ)ℓ 1 L with vℓ u < 1/ℓ and an ≥ ⊂ k − k eigenvalue µℓ spec(Lvℓ ) in BoxN,n(u) so that µℓ τk rk for any 0 k n. Let∈ g be an eigenfunction of L in| H1−corresponding| ≥ to ≤ ≤ ℓ vℓ + µℓ, Lvℓ gℓ = µℓgℓ, with gℓ = 1. Since (µℓ)ℓ 1 is bounded, (gℓ)ℓ 1 is a k1 k ≥ ≥ bounded sequence in H+. Hence by choosing subsequences if needed, we can assume without loss of generality that (µℓ)ℓ 1 converges to a ≥ 1 complex number µ and (gℓ)ℓ 1 converges weakly in H+ to an element 1 ≥ g H+. Then µ BoxN,n(u) 0 k n Drk (τk), and limℓ µℓgℓ = µg ∈ ∈ \ ≤ ≤ 1 →∞ as well as limℓ Lvℓ gℓ = Lug weakly in H+. It then follows that →∞ S Lug = µg and hence that µ is an eigenvalue of Lu, contradicting the fact that specLu BoxN,n(u) 0 k n Drk (τk). Hence there exists ∩ ⊂ ≤ ≤ ε > 0 so that for any v B (u) n ∈ εn S specL Box (u) D (τ ) . v ∩ N,n ⊂ rk k 0 k n ≤[≤ ON THE INTEGRABILITY OF THE BO EQUATION 27

For any 0 k n, denote by C the circle of radius 1/3+ γ (u)/2 ≤ ≤ k k with counterclockwise orientation, centered at τk. Then for any λ Ck and v B (u), L λ is invertible and ∈ ∈ εn v − 2 1 1 C B (u) (L ,H ), (λ, v) (λ L )− k × εn →L + + 7→ − v is analytic and so is the Riesz projector

2 1 1 1 P : B (u) (L ,H ), v (λ L )− dλ . k εn →L + + 7→ 2πi − v CZk Since for any given v B (u), L has a compact resolvent, P (v) is ∈ εn v k an operator of finite rank. Hence its trace, TrPk(v), is finite. Actually, TrPk(v) = 1 since TrPk(v) is constant and TrPk(u) = 1. Hence for any v B (u) and 0 k n, L has precisely one eigenvalue in D (τ ) ∈ εn ≤ ≤ v rk k and this eigenvalue is simple. We denote it by λk(v). By functional 1 1 calculus one has λk(v)= Tr λ(λ Lv)− dλ and hence it follows 2πi Ck − that λ : B (u) C, v λ (v) is analytic.  k εn → 7→ k R As a second step, we prove that the eigenfunctions fn defined by Def- 1 inition 2.1 are real analytic maps with values in H+. To state this result in more detail, we first need to make some preliminary considerations. Denote by L2 L2 (T, C) the Hardy space − ≡ − 0 L2 := h L2 : h(x)= h(k)eikx − { ∈ } k= X−∞ 2b 2 2 and the corresponding Szeg˝oprojector Π− : L L . For any u Lr, → − ∈ denote by Lu− the operator 1 2 1 1 2 Lu− = i∂x Tu− : H L , H := H L − − → − − ∩ − where Tu− denotes the Toeplitz operator 1 2 Tu− : H L , h Π−(uh) . − → − 7→ Using that u is real valued, one verifies that the spectrum of Lu− coin- cides with the one of Lu and that for any n 0, Lu−fn− = λnfn− where ≥ 2 fn− = fn. Hence (fn−)n 0 is an orthonormal basis of L , normalized ≥ − in such a way that 1 f0− > 0 and fn−+1 S−fn− > 0 for any n 0. 2 2 h | i ix h | i ≥ Here S− : L L , h e− h denotes the shift operator to the left. − → − 7→ Furthermore, the Riesz projector Pn−(u) onto span(fn−) is given by

1 1 P −(u)= (λ L−)− dλ n 2πi − u CZn where Cn is the same circle appearing in the definition of the Riesz pro- jector Pn(u), defined in the proof of Lemma 5.1. Finally we introduce for any n 0 the functions ≥ (5.1) f ( ,u) := f ( ,u)+ f −( ,u) /2 n,1 · n · n ·  ON THE INTEGRABILITY OF THE BO EQUATION 28 and

(5.2) f ( ,u) := f ( ,u) f −( ,u) /2i n,2 · n · − n · implying that f ( ,u)= f ( ,u)+if ( ,u). Since for u real valued, n · n,1 · n,2 ·  fn,1(x, u) is the real part of fn(x, u) and fn,2(x, u) its imaginary part one has

1 f ( ,u) = 1, f ( ,u) if ( ,u) = 1, f −( ,u) h | n · i h n,1 · − n,2 · i h n · i and in turn P (u)1 = 1, f −( ,u) f ( ,u) . n h n · i n · Lemma 5.2. For any u L2 and n 0 there exists 0 < δ ε , ∈ r ≥ n ≤ n with εn as in Lemma 5.1, so that for any 0 k n, fk and fk− admit analytic extensions ≤ ≤ 1 1 fk : Bδn (u) H+ , fk− : Bδn (u) H . → → − As a consequence, fk,1 and fk,2, defined by (5.1) - (5.2), admit analytic 1 extensions, f : B (u) H and so does 1, f − : B (u) C and k,j δn → h k i δn → Pk1, meaning that 1 P ( )1 : B (u) H , v P (v)1 = 1, f −( , v) f ( , v) k · δn → + 7→ k h k · i k · is analytic. Furthermore, fk satisfies the normalisation condition 2 2 f , f − = f , f + f , f =1 . h k k i h k,1 k,1i h k,2 k,2i Proof. We argue inductively and begin with f0. By the normalisation of 2 f0, for any v Bεn (u) Lr, 1 f0( , v) > 0. Since for such a v the oper- ∈ ∩ h | 1 · i 1 2 ator Lv is selfadjoint, P0(v)= (λ Lv)− dλ is the L orthogonal 2πi C0 − − projector onto span(f ( , v)) and hence P (v)1 = 1 f ( , v) f ( , v). 0 R 0 0 0 It implies that P (v)1 1· = 1 f ( , v) 2 > 0 and thush | · i · h 0 | i h | 0 · i 1 f ( , v) = P (v)1 1 = P (v)1, 1 . h | 0 · i h 0 | i h 0 i 2 1 Since P0 : Bεn (u) (L+,Hp+) is analytic wep can shrink εn if needed so that Re P (v)1, 1→L> 0 is uniformly bounded away from 0 on B (u). h 0 i εn As a consequence, the principal branch of the root P0(v)1, 1 is analytic on B (u). The analytic extension of f is then definedh byi εn 0 p 1 1 f0( , v) := P0(v)1 H . · P (v)1, 1 ∈ + h 0 i By , p 1 1 L P (v)1 = λ(λ L )− 1dλ = λ (v)P (v)1 v 0 2πi − v 0 0 CZ0 and hence L f ( , v)= λ (v)f ( , v). Using the Riesz projector P −(v) v 0 · 0 0 · 0 instead of P0(v) one sees by the same arguments that the analytic extension of f − is defined for v B (u) by 0 ∈ εn 1 1 f0−( , v) := P0−(v)1 H , · P −(v)1, 1 ∈ − h 0 i p ON THE INTEGRABILITY OF THE BO EQUATION 29 and that L−(v)f −( , v)= λ (v)f −( , v). v 0 · 0 0 · Now assume that after shrinking ε if necessary, f and f −, 0 ℓ k, n ℓ ℓ ≤ ≤ have been extended analytically to Bεn (u) for a given k n 1. More precisely, for any 0 ℓ k ≤ − ≤ ≤ 1 1 fℓ : Bεn (u) H+, fℓ− : Bεn (u) H → → − are analytic maps and for any v B (u), ∈ εn L f ( , v)= λ (v)f ( , v) , L−f −( , v)= λ (v)f −( , v) . v ℓ · ℓ ℓ · v ℓ · ℓ ℓ · 2 In case v Bεn Lr, one has fk+1( , v) Sfk( , v) > 0. Since in such ∈ ∩ h · | · i 1 a case Lv is selfadjoint, the Riesz projector Pk+1(v) = (λ 2πi Ck+1 1 2 − Lv)− dλ is the L orthogonal projector onto span(fk+1( R, v)) and hence − · P (v)Sf ( , v)= Sf ( , v) f ( , v) f ( , v) . k+1 k · h k · | k+1 · i k+1 · It implies that Sf ( , v) f ( , v) = P (v)Sf ( , v) Sf ( , v) . h k · | k+1 · i h k+1 k · | k · i 2 1 Since Pk+1 : Bεn (u) (L+,Hp+) is analytic we can once more shrink →L εn if needed so that Re αk > 0 is uniformly bounded away from 0 on B (u) where α α (v) is defined by εn k ≡ k α (v) := P (v)Sf ( , v) , S−f −( , v) . k h k+1 k · k · i As a consequence, the principal branch of the root √αk is analytic on

Bεn (u). The analytic extension of fk+1 to Bεn (u) is then given by

1 1 fk+1( , v)= Pk+1(v)Sfk( , v) H+ · αk(v) · ∈ and by functional calculus one has L f ( , v) = λ (v)f ( , v). p v k+1 · k+1 k+1 · Using the Riesz projector Pk−+1(v) instead of Pk+1(v) one sees by the arguments used above that the analytic extension of fk−+1 is given for v B (u) by ∈ εn 1 1 fk−+1( , v)= Pk−+1(v)S−fk−( , v) H − · αk−(v) · ∈ and that Lv−(v)fk−+1( , v)p = λk+1(v)fk−+1( , v). Here αk−(v) is defined by · · α−(v) := P − (v)S−f −( , v) , Sf ( , v) . k h k+1 k · k · i We then denote by 0 < δ ε the number obtained after shrinking n ≤ n εn successively so that for any 0 k n , fk extends to an analytic 1 ≤ ≤ map, fk : Bδn (u) H+ and similarly, fk− extends to one on Bδn (u) with values in H1 .→  − 2 We now present a formula for the L gradient of the eigenvalues λn, − 2 R n 0. By Lemma 5.1, for any n 0, λn : u Lr λn(u) is ≥ 2 ≥ ∈ 7→ ∈ real analytic and hence its L gradient λn is well defined. More precisely, we have the following− ∇ ON THE INTEGRABILITY OF THE BO EQUATION 30

2 2 Corollary 5.1. For any u Lr, n 0, and 0 k n, the L gradient λ (u) of λ at u is given∈ by λ (≥u)= f ≤2( ,u≤) and extends− to an ∇ k k ∇ k −| k| · analytic function on Bδn (u) with δn > 0 given as in Lemma 5.2. More precisely, λ : B (u) H1, v f ( , v)2 f ( , v)2 ∇ k δn → 7→ − k,1 · − k,2 · is analytic where fk,1 and fk,2 are given as in Lemma 5.2. 2 Proof. Let u Lr and 0 k n be given. To compute the gradient λ (u), consider∈ for any v≤ L≤2 the one parameter family u := u+εv. ∇ k ∈ r ε It is convenient to introduce the notation fk := fk( ,u), λk := λk(u), and · d d δf := f ( ,u ), δλ := λ (u ). k dε|ε=0 k · ε k dε|ε=0 k ε Taking the derivative with respect to ε of both sides of the identity (L λ (u ))f ( ,u )=0at ε = 0 one obtains uε − k ε k · ε (L λ )δf T f δλ f =0. u − k k − v k − k k Taking the inner product with fk of both sides of the latter identity and using that L λ is selfadjoint one then concludes v − k T f f δλ f f =0. −h v k| ki − kh k| ki 2 Since fk = 1 and Πfk = fk, the definition of the L gradient then impliesk thatk − λ , v = δλ = vf f = f 2 v , h∇ k i k −h k| ki h−| k| | i yielding the claimed formula for λ (u). Since by Lemma 5.2, the real ∇ k and imaginary parts fk,1, fk,2 of fk, extend to analytic functions on 1 Bδn (u) with values in H and by Lemma 5.1, and λk also extends to an analytic functional on Bδn (u), it follows that the gradient λk and its formula extend analytically, ∇ λ : B (u) H1, v f ( , v)2 f ( , v)2 ∇ k δn → 7→ − k,1 · − k,2 · proving the lemma.  As a consequence of Corollary 5.1, we obtain a formula for the gra- dient of γn = λn λn 1 1, − − − 2 2 (5.3) γn = fn + fn 1 . ∇ −| | | − | 2 As a third step, we show that for any n 0, the L gradient 1 fn 2 1 ≥ − ∇h | i is a real analytic map Lr H , u 1, fn( ,u) . Recall that for any u L2 and n 1, the→ eigenfunctions7→ ∇h f ,· 0 ik n, admit an ∈ r ≥ k ≤ ≤ analytic extension fk = fk,1 + ifk,2 to Bδn (u) , i.e. f : B (u) H1, v f ( , v) , j =1, 2, k,j δn → 7→ k,j · are analytic maps (cf. Lemma 5.2). For any v Bδn (u), we denote by df ( , v)[w] the derivative of f at v in direction∈ w L2. k,j · k,j ∈ ON THE INTEGRABILITY OF THE BO EQUATION 31

Lemma 5.3. For any n 0 and v Bδn (u) there exists a constant c > 0 so that for any 0 ≥k n, j =1∈ , 2, and w L2, n ≤ ≤ ∈ df ( , v)[w] c w −1 . k k,j · k ≤ nk kH The constant cn can be chosen locally uniformly with respect to v. As a consequence, 1, f : B H1, j =1, 2, are analytic. ∇h k,ji δn → Proof. The latter assertion easily follows from the first one, since for any 0 k n, j =1, 2, v B (u), and w L2, ≤ ≤ ∈ δn ∈ 1, f ( , v) ,w = 1, df ( , v)[w] |h∇h k,j · i i| |h k,j · i| df ( , v)[w] c w −1 , ≤k k,j · k ≤ nk kH 1 implying that gk,j(v) := 1, fk,j( , v) is bounded in H uniformly ∇h · i 2 with respect to v Bδn (u). Since gk,j : Bδn (u) L is an ana- ∈ inx → C lytic map (cf. Lemma 5.2), so is gk,j e : Bδn (u) for any h | i 1 → n Z. Together with the boundedness of gk,j in H it then follows ∈ 1 that gk,j : Bδn (u) H is analytic (cf. e.g. [12, Appendix A]). To prove the first→ assertion of the lemma, we proceed by induction. First we need to make some preliminary considerations. Without further reference, we will use the notation introduced in the deriva- tion of Lemma 5.1 and Lemma 5.2. According to Lemma 5.2, fk,1 = (fk + fk−)/2 and fk,2 = (fk fk−)/2i where fk−, k 0, are the eigen- 1 2 − ≥ functions of Lv− : H L , h i∂xh Π−(vh) and Π− is the Szeg˝o projector L2 L2−.→ The− first7→ assertion− of the lemma then follows → − by proving corresponding estimates for dfk( , v)[w] and dfk−( , v)[w], 0 k n. By the definition of r , · · ≤ ≤ n 2 1 1 (Dr + 1 Dr + 1 ) Bδn (L+,H+), (λ, v) (λ Lv)− n 2 \ n 4 × →L 7→ − 1 is analytic. Hence the adjoint of (λ Lv)− with respect to the dual 1 1 − pairing of H+− and H+, also gives rise to an analytic map, 1 2 1 (Dr + 1 Dr + 1 ) Bδn (H+− , L+), (λ, v) (λ Lv)− . n 2 \ n 4 × →L 7→ −

Since the circles Ck, 0 k n, are compact, it implies that for any v B there exists M ≤> 0≤ so that for any λ C , 0 k n, ∈ δn n ∈ k ≤ ≤ 1 1 − (5.4) (λ Lv)− L2 H1 Mn , (λ Lv)− H 1 L2 Mn , k − k +→ + ≤ k − k + → + ≤ where Mn can be chosen locally uniformly with respect to v. Further- more, for any 0 k n, the Riesz projector ≤ ≤ 1 1 P (v)= (λ L )− dλ k 2πi − v CZk

2 1 2 gives rise to an analytic map Pk : Bδn (L+,H+). For any h L+, the derivative of P (v)h with respect to→Lv in direction w L2 is∈ given k ∈ ON THE INTEGRABILITY OF THE BO EQUATION 32 by 1 1 1 (dP (v)[w])h = (λ L )− Π(w(λ L )− h)dλ . k −2πi − v − v CZk

We want to estimate the operator norm of dPk(v)[w], viewed as an operator on L2 . Note that for any g, h L2 , + ∈ + 1 1 1 (dP (v)[w])h, g = w(λ L )− h, (λ L )− g dλ . h k i −2πi h − v − v i CZk Taking into account (5.4) and that multiplication by f H1 yields a 1 ∈ on H− , we infer that for any v Bδn , there exists a constant d > 0 so that for any 0 k n ∈ n ≤ ≤ (5.5) dPk(v)[w] L2 L2 dn w H−1 k k +→ + ≤ k k where dn can be chosen locally uniformly with respect to v. We are now ready to present the induction argument with respect to 0 k n for proving the claimed estimates for df ( , v)[w] and ≤ ≤ k · df −( , v)[w]. We start with k = 0. Recall that the formula for the k · analytic extension of f0 to Bδn of Lemma (5.2) reads 1 f0( , v)= P0(v)1, · P (v)1, 1 h 0 i implying that for any w L2 ∈ p 1 1 (dP0(v)[w])1, 1 df0( , v)[w]= (dP0(v)[w])1 h 3 iP0(v)1 . · P0(v)1, 1 − 2 P (v)1, 1 2 h i h 0 i Applying (5.5)p for k = 0, we immediately conclude that there exists An > 0 so that df (v)w 2 A w −1 . k 0 kL ≤ nk kH The corresponding estimate for df0−(v)[w] is proved in an analogous way by using the estimate for dP0−(v)[w] corresponding to the one of dP0(v)[w]. Now assume that the claimed estimates for dfℓ( , v)[w] and dfℓ−( , v)[w] are valid for any 0 ℓ k for a given 0 k · n 1. We claim· that they also hold for k≤+ 1.≤ Let us first consider≤ the≤ one− for df ( , v)[w] k+1 · Again we use the formula for the analytic extension of fk+1 to Bδn derived in Lemma 5.2, 1 fk+1( , v)= Pk+1(v) Sfk( , v) · αk+1(v) · where  p ix α (v)= P (v)(Sf ( , v)) e− f −( , v) k+1 k+1 k · | k · The directional derivative f ( , v)[w] can then be computed as k+1 · 1 1 dαk+1(v)[w] dfk+1( , v)[w]= I 3/2 Pk+1(v) Sfk( , v) · αk+1(v) − 2 αk+1(v) ·  p ON THE INTEGRABILITY OF THE BO EQUATION 33 where I := dP (v)[w](Sf ( , v))+ P (v)(Sdf ( , v)[w]) k+1 k · k+1 k · and dαk+1(v)[w]= II1 + II2 + II3 with ix II := dP (v)[w](Sf ( , v)), e− f −( , v) 1 h k+1 k · k · i ix II := P (v) Sdf ( , v)[w] , e− f −( , v) 2 h k+1 k · k · i ix II3 := Pk+1(v)[Sfk( , v)], e− df −( , v)[w] . h · k · i By increasing An if needed, the estimate (5.5) for k + 1 together with the induction hypothesis then yield

df (v)[w] 2 A w −1 , k k+1 kL ≤ nk kH The corresponding estimate for dfk−+1(v)[w], is proved in an analogous way by using the estimate for dPk−(v)[w] corresponding to the one of dPk(v)[w]. Altogether we have proved the induction step. Going through the arguments of the proof one verifies that the constant An can be chosen locally uniformly with respect to v.  Remark 5.1. At u = 0, one has f (x) = eikx for any k 0 and thus k ≥ 1 fn =0 for any n 1. Let us compute the gradient of 1 fn at u =0 h | i ≥ 2 R d h | i for n 1. For any v Lr and ε , define δfn := dε ε=0fn( , εv) and ≥ ∈ ∈ | · δλ := d λ ( , εv). Arguing as in the proof of Corollary 5.1 one has n dε ε=0 n · for any n| 0 ≥ (D n)δf = T (einx)+ δλ einx . − n v n If n 1, taking the inner product of both sides with 1, we obtain ≥ n δf 1 = einx, v − h n| i h i which leads to 1 inx 1 f = e− . ∇h | ni −n

As a fourth step, we discuss the analytic extension of κn, n 0, given by the formulae in Corollary 3.1. A difficulty here is that these≥ formulae involve infinitely many eigenvalues λn, while the domain of analyticity of λ might shrink as n . Therefore we need an extra n →∞ argument. We shall appeal to the generating functional λ .First we H2 need to make some preliminary considerations. Let u Lr be given. If needed, shrink the radius δ > 0 of the ball B ∈B (u) so that n δn ≡ δn for any v Bδn and 1 k n, λk 1(v) + 1 is in the disk Drn (τn). ∈ ≤ ≤ − 1 For any v Bδn and 0 k n, λ(v) = (Lv + λId)− 1, 1 is holomorphic∈ for λ in D ≤ (τ≤) λH(v) andh possibly has a simplei − rk+1/4 k \{ k } pole at λ = λk(v). Let Pk⊥(v)= Id Pk(v) where Pk(v) denotes the Riesz projector,− introduced in the proof− of Lemma 5.1. Note that 1 ⊥ (v) := (L + λId)− P ⊥(v)1, 1 Hλ,k h v k i ON THE INTEGRABILITY OF THE BO EQUATION 34 is holomorphic on D (τ ). By Lemma 5.2, − rk+1/4 k P (v)1 = 1 , f ( , v) if ( , v) f ( , v) k h k,1 · − k,2 · i k · 1 implying that (v) := (L + λId)− P (v)1, 1 equals Hλ,k h v k i 2 2 1 fk,1( , v), 1 + fk,2( , v), 1 . h · i h · i λk(v)+ λ  Hence the residue resk(v) of λ(v)= λ,k(v)+ λ,k⊥ (v) at λ = λk(v) is given by H H H − (5.6) res (v)= f ( , v), 1 2 + f ( , v), 1 2 . k h k,1 · i h k,2 · i

Furthermore, for any v Bδn , (λk(v)+ λ) λ(v) is holomorphic on D (τ ). Actually, the∈ map H − rk+1/4 k D (τ ) B C, (λ, v) (λ (v)+ λ) (v) − rk+1/4 k × δn → 7→ k Hλ is analytic. Let us now discuss this map in more detail. We begin with the case k = 0. For any (λ, v) in D (τ ) B set − r0+1/4 0 × δn χ (λ, v) := (λ (v)+ λ) (v) , κ (v) := χ ( λ (v), v) . 0 0 Hλ 0 0 − 0 Notice that Proposition 3.1 (i) implies that, for v B L2, the above ∈ δn ∩ r definition of κ0(v) is consistent with the one in Corollary 3.1. With this notation we then get res (v)= f ( , v), 1 2 + f ( , v), 1 2 = κ (v) . 0 h 0,1 · i h 0,2 · i 0 By Corollary 3.1, for any v B L2, one has κ (v) > 0. Hence by ∈ δn ∩ r 0 shrinking δn if necessary, we conclude that Re κ0(v) > 0 is uniformly bounded away from 0 on Bδn . Now let us consider the case 1 k n. Recall that our choice of δ ≤ ≤ n mentioned above assures that λk 1(v)+1 Drk (τk) for any v Bδn . Hence it follows from Lemma 5.1−that the map∈ ∈ C Bδn , v (λk(v)+ λ) λ(v) λ= λ − (v) 1 → 7→ H | − k 1 − is also analytic. Since by Proposition 3.1, for any v B L2, ∈ δn ∩ r (λk(v)+ λ) λ(v) λ= λ − (v) 1 =0 , H | − k 1 −

(λk(v)+ λ) λ(v) λ= λ − (v) 1 vanishes identically on Bδn . Thus the H | − k 1 − function (λk(v)+ λ) λ(v) is of the form (λ + λk 1(v)+1)χk(λ, v) where H − − χ : D (τ ) B C k − rk k × δn → is analytic. Hence we have proved that for any (λ, v) D (τ ) B , ∈ − rk k × δn λk 1(v)+1+ λ (5.7) λ(v)= − χk(λ, v) . H − λk(v)+ λ

Using that λk 1(v)+1= λk(v) γk(v) one then concludes from (5.6) that the residue− res (v) of (v−) at λ = λ (v) is given by k Hλ − k f ( , v), 1 2+ f ( , v), 1 2 = γ (v)κ (v) , κ (v) := χ ( λ (v), v) . h k,1 · i h k,2 · i k k k k − k ON THE INTEGRABILITY OF THE BO EQUATION 35

Again Proposition 3.1 (i) implies that, for v B L2, the above ∈ δn ∩ r definition of κk(v) is consistent with the one in Corollary 3.1. Moreover, by Corollary 3.1, κ (v) > 0 for any v B L2. Hence by shrinking δ k ∈ δn ∩ r n if necessary, we conclude that Re κk(v) > 0 is uniformly bounded away from 0 on Bδn . For later reference we record our findings as follows:

Lemma 5.4. After shrinking δn > 0 if needed, the following holds: for any u L2, v B (u), ∈ r ∈ δn f ( , v), 1 2 + f ( , v), 1 2 = κ (v) h 0,1 · i h 0,2 · i 0 and for 1 k n, ≤ ≤ f ( , v), 1 2 + f ( , v), 1 2 = γ (v)κ (v) . h k,1 · i h k,2 · i k k For any 0 k n, the functional κ : B (u) C is analytic with ≤ ≤ k δn → positive real part and Re κk(v) > 0 is uniformly bounded away from 0. Finally we analyze the L2 gradient of κ . − n 1 2 Lemma 5.5. For any n 1, κn takes values in Hr and κn : Lr 1 ≥ ∇ ∇ → Hr is real analytic. 2 1 Remark 5.2. (i) The map κ0 : Lr Hr is real analytic. Indeed, by ∇ 2 → 2 Lemma 5.4, κ0(u) = 1, f0,1( ,u) + 1, f0,2( ,u) . The claim then follows from Lemma 5.3h . · i h · i (ii) For u =0 one infers from the product representation of κn and the inx facts that fn(x, 0) = e , n 0, and both γp(0) and γp(0), p 1, vanish, that κ (0) = 0 for any≥ n 0. ∇ ≥ ∇ n ≥ 2 Proof. Let n 1 and u Lr. For any 1 k n, κn is an analytic map B (u) ≥C (cf. Lemma∈ 5.4), implying≤ that≤ κ : B (u) L2 δn → ∇ n δn → is such a map. To simplify notation, we write Bδn for Bδn (u). We now 1 prove that κn is in fact an analytic map Bδn H . Let us start from the identity∇ (5.7), →

(5.8) (λk + λ) λ(v)= (λ + λk 1(v)+1)χk(λ, v) , H − − which holds for (λ, v) Dr + 1 (τk) Bδn . Here χk is a holomor- ∈ − k 4 × phic function on the latter domain and is related to κk by κk(v) = χ ( λ (v), v). By the chain rule we then have k − k ∂χ κ (v)= k ( λ (v), v) λ (v)+ χ ( λ (v), v) . ∇ k − ∂λ − k ∇ k ∇ k − k Since λ : B H1 is an analytic map (cf. Corollary 5.1, Lemma ∇ k δn → 5.2 ), it remains to prove that v χk( λk(v), v) is an analytic map B H1. It suffices to prove that7→ ∇ − δn → 1 χk : Dr + 1 (τk) Bδn H ∇ − k 4 × → is analytic. Indeed, let us take the gradient with respect to v of both sides of the identity (5.8). Writing G (λ, v) := (λ (v)+ λ) (v) , k k Hλ ON THE INTEGRABILITY OF THE BO EQUATION 36 we obtain the following identity,

(5.9) Gk(λ, v)= (λ + λk 1(v)+1) χk(λ, v) χk(λ, v) λk 1(v) . ∇ − − ∇ − ∇ − 1 1 We claim that Gk is H valued and Gk : Dr + 1 (τk) Bδn H ∇ ∇ − k 4 × → is analytic. Indeed, Gk satisfies the identity (cf. derivation of Lemma 5.4) 2 2 G (λ, v)= f ( , v), 1 + f ( , v), 1 +(λ + λ) ⊥ (v) k h k,1 · i h k,2 · i k Hλ,k where 1 ⊥ (v)= (L + λId)− P ⊥(v)1, 1 , P ⊥(v)= Id P (v) , Hλ,k h v k i k − k and P (v) denotes the Riesz projector onto span(f ( , v)). For any k k · (λ, v) Dr + 1 (τk) Bδn , we rewrite ∈ − k 4 × 1 1 1 (L + λId)− P ⊥(v)=(L + λId)− (L + λId)− P (v) , v k v − v k using 1 1 1 1 (L + λId)− = (L + λId)− dµ , v 2πi v µ + λ CZk together with

1 1 1 1 (L + λId)− P (v)= (L + λId)− (µId L )− dµ , v k 2πi v − v CZk and the resolvent identity 1 1 1 1 (µ + λ)(L + λId)− (µId L )− =(L + λId)− +(µId L )− . v − v v − v It yields

1 1 1 1 (L + λId)− P ⊥(v)= (µId L )− dµ . v k −2πi µ + λ − v CZk Consequently,

1 1 1 1 ⊥ (v)= (µId L )− [1](µId L ) [1] dµ ∇Hλ,k 2πi − v − v − λ + µ CZk 1 1 implying that ⊥ (v) H and ⊥ : D 1 (τk) Bδn H λ,k ,k rk+ 4 ∇H ∈ ∇H· − 1 × → is analytic. Finally, since fk,j( , v), 1 : Bδn H , j = 1, 2, are analytic (cf. Lemma 5.3), it∇h follows· i → 1 Gk : Dr + 1 (τk) Bδn H ∇ − k 4 × → is analytic. Coming back to (5.9), we infer that the holomorphic map 1 Gk + χk λk 1 : Dr + 1 (τk) Bδn H ∇ ∇ − − k 4 × → vanishes at λ = λk 1(v) 1 for any v Bδn . Therefore there exists a holomorphic map− − − ∈ 1 Fk : Dr + 1 (τk) Bδn H − k 4 × → ON THE INTEGRABILITY OF THE BO EQUATION 37 such that for any (λ, v) Dr + 1 (τk) Bδn ∈ − k 4 × Gk(λ, v)+ χk(λ, v) λk 1(v)= (λ + λk 1(v)+1)Fk(λ, v) . ∇ ∇ − − − Comparing with (5.9), one deduces that χ (λ, v)= F (λ, v) ∇ k k and thus has the claimed property. 

With these preparations made, we now can study functionals ζn. Recall from the definition (4.1) that for any n 1, u L2 , ζ (u) = ≥ ∈ r,0 n 1 fn( ,u) 2 C h | · i . The functionals 1 fn : Lr,0 (cf. Lemma 5.2), and √κn(u) h | i → κ : L2 R (Lemma 5.4) are real analytic. Hence we claim that n r,0 → ζ : L2 C is real analytic. Indeed, let u L2 , n 1, and n r,0 → ∈ r,0 ≥ 1 k n. By Lemma 5.2, 1, f if is analytic on B (u) L2 ≤ ≤ h k,1 − k,2i δn ∩ r,0 and by Lemma 5.4, so is κk. Since by the same lemma, Re κk > 0 is uniformly bounded away from 0 on B (u) L2 it then follows that the δn ∩ r,0 principal branch of the root √κk is a well defined, analytic functional which is also uniformly bounded away from 0 on this ball. We thus have 1 f 1 f proved that ζ = h | k,1i i h | k,2i extends analytically to B (u) L2 . k √κk √κk δn r,0 − ∩ 2 Now we study the gradient of ζn. By the chain rule we have on Lr,0

1 κn 1 (5.10) ζn = ∇3/2 1 fn + 1 fn . ∇ −2 κn h | i √κn ∇h | i Since 1 f : L2 H1 (cf. Lemma 5.3) and κ : L2 H1 (cf. ∇h | ni r,0 → ∇ n r,0 → Lemma 5.5) are real analytic, the following result holds true. Proposition 5.1. For any n 1, the functional ζ : L2 C and the ≥ n r,0 → map ζ : L2 H1 are real analytic. ∇ n r,0 → 2 Remark 5.3. The L gradient of ζn can be computed explicitly at − 1 1 inx u =0. Indeed, since κn(0) = n and 1 fn (0) = n e− (cf. Remark 5.1) one has ∇h | i − 1 inx ζ (0) = e− . ∇ n −√n

6. Poisson brackets The main goal of this section is to compute the Poisson brackets γ ,γ and γ , 1 f . As a first step we show that the flow, corre- { p n} { p h | ni} sponding to the generating function λ, leaves the eigenvalues of Lu H 2 invariant. Recall from (3.1) that for u Lr, the generating function is given by ∈ 1 (u)= (L + λId)− 1 1 . Hλ h u | i Consider with λ R λ (u): n 0 . For any v L2 one has Hλ ∈ \ {− n ≥ } ∈ r 1 1 2 d (u)[v]= (L + λId)− T (L + λId)− 1 1 = v w , Hλ h u v u | i h || λ| i ON THE INTEGRABILITY OF THE BO EQUATION 38 where w w ( ,u) is given by λ ≡ λ · 1 1 w := (L + λId)− 1 H . λ u ∈ + Hence by the definition of the L2 gradient − (u)= w ( ,u) 2 ∇Hλ | λ · | and the Hamiltonian equation associated to is Hλ (6.1) ∂ u = ∂ w ( ,u) 2 . t x| λ · | 2 1 Since the map u Lr wλ( ,u) H+ is Lipschitz on bounded sub- ∈ 7→ · ∈ 2 2 2 sets, so is the Hamiltonian vector field u Lr ∂x wλ( ,u) Lr and hence the corresponding initial value problem∈ 7→ admits| · a local| ∈ in time solution for any initial data u L2 and λ C λ (u): n 0 . Fur- ∈ r ∈ \{ n ≥ } thermore, the Toeplitz operators Twλ and Twλ are bounded operators 2 on L+. The following result says that equation (6.1) has a Lax pair representation. Proposition 6.1. For any λ R λ (u): n 0 , ∈ \ {− n ≥ } dL u = [Bλ, L ] , Bλ := iT T dt u u u wλ wλ along the evolution (6.1).

1 Remark 6.1. Recall that for ε > 0 we introduced ˜ε = 1 and H ε H ε dk ˜ that dεk ε=0 ε, k 0, are referred to as the Hamiltonians in the BO H ≥ 1 ˜ 1 ε hierarchy. Introducing Bε := ε Bu , Proposition 6.1 then leads to a Lax pair formulation for the equations corresponding to the Hamiltonians in the BO hierarchy, dL dk u = [ B˜ , L ] , dt dεk ε=0 ε u dk ˜ where now u evolves according to the flow of dεk ε=0 ε. In particular, 3 H in the case k = 3, for any u L2 , (u) = 1 d ˜ (u) and one r,0 6 dε3 ε=0 ε ∈ 3 H − H can check that the operator 1 d B˜ coincides with B˜ 1 i u 2Id, 6 dε3 ε=0 ε u 2 ˜ 2 − − k k where Bu = i(T ∂x u Tu ) was introduced in Remark 2.1. | | − 2 Recall that by Corollary 5.1, λn = fn for any n 0. Hence 1 ∇ −| | ≥ 2 λn Hr and the Poisson bracket λp,λn is well-defined on Lr for any∇ p∈ 0. Clearly, the same is true for{ γ ,γ} . ≥ { p n} 2 Corollary 6.1. (i) For any p, n 0 and u Lr, λp,λn (u)=0. (ii) For any p, n 1 and u L2,≥ γ ,γ (u∈)=0 { } ≥ ∈ r { p n} Remark 6.2. By the arguments of the proof of Corollary 6.1, , =0 , ,λ =0 {Hλ Hµ} {Hλ n} for any λ,µ C λ (u): k 0 and n 0. ∈ \ {− k ≥ } ≥ ON THE INTEGRABILITY OF THE BO EQUATION 39

Proof of Corollary 6.1. (i) In view of Proposition 6.1, for any n 0 ≥ and λ in R λk(u): k 0 , λn(u) is constant along the evolution (6.1). Hence\ on{−L2, ,λ≥ =} 0. Using (3.2), we therefore have r {Hλ n} ∞ 1 f 2 1 f 2,λ |h | pi| λ ,λ + {|h | pi| n} =0 . −(λ + λ)2 { p n} λ + λ p=0 p p X   By analytic continuation, the latter identity continues to hold for any λ in C λk(u): k 0 . Since the left hand side of this identity is a meromorphic\ {− function≥ in}λ, the principal part of the expansion at any of its pole vanishes, implying that λ ,λ (u)=0 for any p 0.  { p n} ≥ Let us come to the proof of Proposition 6.1. Proof of Proposition 6.1. Assume that u u(t) is a local in time solution of (6.1) with λ R λ (u(0)) : n ≡ 0 . In the sequel of the ∈ \{ n ≥ } proof, we denote also by D the operator i∂x, when viewed as operator 1 − acting on H . Since by the definition of the operator Lu one has

dLu = T∂tu = iTD w 2 , dt − − | λ| the claimed result is equivalent to the identity

2 (6.2) [D Tu, Twλ Twλ ] TD w =0 . − − | λ| We will use Hankel operators with symbols in L+∞ to provide a formula for the commutators of Toeplitz operators. Given a symbol v L∞, ∈ + the Hankel operator Hv is given by H : L2 L2 , h H (h)=Π(vh) . v + → + 7→ v Notice that Hv(h) is antilinear in h, but linear in v.

Lemma 6.1. For any v,w L∞ ∈ + (6.3) [T , T ]= w v H H . v w h·| i − v w Furthermore, for any u L∞, ∈ r (6.4) H = T H + H T Πu w Tuw w Πu w u −h |·i (6.5) H = H T + T H w Πu . Tuw Πu w u w −h |·i Proof of Lemma 6.1. We shall make use of the following elementary identities, valid for any f L2 and h L2 , ∈ ∈ + f = Πf + Πf f, 1 , Π(fh) = Π((Πf)h) . −h i Then, for any h L2 , ∈ + [T , T ]h = Π(vΠ(wh)) Π(wΠ(vh)) = Π(vΠ(wh)) Π(wvh) . v w − − Using that wh = Π(wh)+ Π(wh) wh 1 one then gets −h | i [T , T ]h = Π(vΠ(wh)) + wh 1 v , v w − h | i ON THE INTEGRABILITY OF THE BO EQUATION 40 which yields (6.3). As for (6.4), we have

HTuwh = Π(Π(uw)h)=Π(uwh) . Since uh = Π(uh)+ Π(uh) u h and u is real valued one obtains −h | i H h = Π(wΠ(uh))+Π(wΠ(uh)) u h w Tuw −h | i = Π(wΠ((Πu)h))+Π(wΠ(uh)) Πu h w −h | i which establishes (6.4). The identity (6.5) is proved in a similar way or alternatively, by taking adjoints with respect to the real inner product, induced by .  h ·|·i Let us now return to the proof of Proposition 6.1 and complete the verification of identity (6.2). We write w for wλ for simplicity. Using (6.3), the Leibniz rule, and the antilinearity of Hw, we get [D, T T ] = [D, T T ] + [D, w w H2 ] w w w w h·| i − w = T T + T T + w Dw Dw w H H + H H Dw w w Dw h·| i −h·| i − Dw w w Dw = TD w 2 + w Dw Dw w HDwHw + HwHDw . | | h·| i −h·| i − Regarding [Tu, TwTw], note that [Tu, TwTw] = [Tu, Tw]Tw + Tw[Tu, Tw]. Using u = Πu + Πu u 1 and again (6.3), one gets −h | i [Tu, TwTw] = [TΠu, Tw]Tw + Tw[TΠu, Tw] = H H T T H H + w T (Πu) T (Πu) w . w Πu w − w Πu w h ·| i w −h·| w i Combining these two identities, we get

[D Tu, TwTw] TD w 2 = (HDw TwHΠu)Hw + Hw(HDw HΠuTw) − − | | − − − + w (Dw T (Πu)) ,Dw T (Πu) w . h·| i − w −h· − w i We now have to analyze the terms HDw TwHΠu and HDw HΠuTw in 1 − − more detail. Since w = (Lu + λ)− 1, hence Dw = Tuw +1 λw, and since by (6.4) − T H = H H T + Πu w w Πu Tuw − w u h |·i one sees that

HDw TwHΠu = HTuw+1 λw HTuw + HwTu Πu w − − − −h |·i = H λH + H T Πu w . 1 − w w u −h |·i Using that u is real valued and hence Πu Hwh = h Hw(Πu) one is led to the following identity for the termh | (H i hT |H )H ,i − Dw − w Πu w (H T H )H = H H +λH H H T H + H (Πu) w . − Dw − w Πu w − 1 w w w − w u w h·| w i Similarly, by (6.5) H T = H T H + w Πu Πu w Tuw − u w h | ·i one sees that

HDw HΠuTw = HTuw+1 λw HTuw + TuHw w Πu − − − −h |·i = H λH + T H w Πu . 1 − w u w −h |·i ON THE INTEGRABILITY OF THE BO EQUATION 41

Finally, since Hw is antilinear and λ is real one concludes that H (H H T )= H H λH H + H T H w H (Πu) . w Dw − Πu w w 1 − w w w u w −h·| i w Substituting these identities into the formula for [D Tu, TwTw] TD w 2 , obtained above, one infers that − − | |

[D Tu, TwTw] TD w 2 = H1Hw + HwH1 − − | | − + w Dw T (Πu) H (Πu) Dw T (Πu) H (Πu) w . h·| i − w − w −h·| − w − w i By (6.4),  T (w)= H (1) = T H (1) + H T (1) Πu 1 w u Tu(w) w Πu w u −h | i = T (Πu)+ H (Πu) u 1 w . w w −h | i The identity Dw = T w +1 λw then reads u − Dw T (Πu) H (Πu)=1 λw u 1 w . − w − w − −h | i Hence

[D Tu, TwTw] TD w 2 = H1Hw + HwH1 + w 1 1 w − − | | − h·| i −h·| i which vanishes since w 1= H H and 1 w = H H .  h·| i 1 w h·| i w 1 The second main result of this section calculates the Poisson brackets of 1 fn with the functionals γp. Recall that for any p 1, γp is in H1h(cf.| i Corollary 5.1) and that for any n 0, also 1≥ f ∇is in H1 ≥ ∇h | ni (cf. Lemma 5.3). Hence the Poisson bracket γp, 1 fn is well-defined on L2 for any p 1 and n 0. { h | i} r ≥ ≥ Proposition 6.2. For any p 1, n 0, and u L2 , ≥ ≥ ∈ r,0 γ , 1 f (u)= i 1 f ( ,u) δ . { p h | ni} h | n · i pn Proof. Again we use the Hamiltonian flow of the generating function and its Lax pair formulation provided by Proposition 6.1, Hλ dL u = [Bλ, L ] , Bλ = iT T . dt u u u wλ wλ Given u L2 , n 0, and λ λ +1 , denote by u = u(t) the 0 ∈ r,0 ≥ ≥ | − 0 | trajectory issuing from u at t = 0, and define gλ gλ(t) by 0 n ≡ n dgλ n = Bλgλ , gλ(0) = f (u ) . dt u n n n 0 Since Bλ is skew symmetric, f ( ,u ) = 1 , gλ(t) = 1 and since u k n · 0 k k n k d (L λ )gλ = Bλ(L λ )gλ , (L λ )gλ =0 , dt u − n n u u − n n u − n n |t=0 λ    gn(t) is an eigenfunction of Lu(t) corresponding to the eigenvalue λn, coinciding with the eigenfunction f ( ,u(t)) up to a phase factor. To n · ON THE INTEGRABILITY OF THE BO EQUATION 42

λ determine this phase factor we need to compute the evolution of g0 1 λ λ h | i and of gn Sgn 1 , n 1. Using that h | − i ≥ λ 1 λ 1 λ Twλ 1= wλ 1 1 , gn wλ = (Lu + λ)− gn 1 = gn 1 h | i h | i h | i λn + λh | i one gets for any n 0 ≥ d gλ 1 = Bλgλ 1 = i gλ T T 1 = i w 1 gλ w dth n| i h u n| i h n| wλ wλ i h λ| ih n | λi implying that λ d λ gn 1 (6.6) gn 1 = i λ h | i . dth | i H λn + λ λ Similarly, using in addition that Bu is skew symmetric and that by (2.4) ST = T S Sw 1 = T S wλ wλ −h λ · | i wλ one has for any n 1 ≥ d λ λ λ λ λ λ λ gn Sgn 1 = gn [S, Bu ]gn 1 = i gn Twλ [Twλ ,S]gn 1 . dth | − i h | − i h | − i Using (2.4) once more one gets

λ λ λ Twλ [Twλ ,S]gn 1 = Swλgn 1 1 Twλ 1= wλ Sgn 1 wλ − h − | i h | − i and in this way concludes that λ d λ λ λ λ gn 1 λ gn Sgn 1 = i gn wλ wλ Sgn 1 = i h | i wλ Sgn 1 . dth | − i h | ih | − i λn + λh | − i Note that by (2.7) λ λ λ (λn 1 + λ + 1)Sgn 1 =(Lu + λ)Sgn 1 + Sgn 1 u 1 , − − − h − | i and that λ λ wλ (Lu + λ)Sgn 1 = 1 Sgn 1 =0 , h | − i h | − i yielding λ λ u Sgn 1 λ wλ Sgn 1 = h | − iH . h | − i λn 1 + λ +1 − Since by Lemma 2.2 λ λ λ λ λ λ gn 1 u Sgn 1 =(λn λn 1 1) gn Sgn 1 = γn gn Sgn 1 −h | ih | − i − − − h | − i h | − i we end up with the following identity

d λ λ 1 γn λ λ gn Sgn 1 = i λ gn Sgn 1 dth | − i − H · λn + λ · λn 1 + λ +1 ·h | − i − 1 1 λ λ (6.7) = i λ gn Sgn 1 H · λn + λ − λn 1 + λ +1 h | − i  −  In view of (6.6) and (6.7) we set

λ λ itαp(λ) g (t)= g (t) e− , p 0 p p ≥ e ON THE INTEGRABILITY OF THE BO EQUATION 43

λ where α0(λ) := H and for any n 1, λ0+λ ≥ 1 n 1 1 α (λ) := + , n Hλ · λ + λ λ + λ − λ + λ +1 0 p=1 p p 1 ! X  −  so that d λ d λ λ g0 1 =0 , gn Sg˜n 1 =0 , n 1 , dth | i dth | − i ≥ implying that for any p 0, g (t)= f (u(t)). Taking into account (6.6), e ≥ p e p one sees that for any λ λ +1 and n 1, , 1 f = d 1 f , ≥ |− 0 | ≥ {Hλ h | ni} dt h | ni when evaluated at t = 0, is givene by

, 1 f = i 1 f α Hλ {Hλ h | ni} h | ni n − λ + λ  n  1 1 n 1 1 = i 1 f + h | ni · Hλ · λ + λ − λ + λ λ + λ − λ + λ +1 0 n p=1 p p 1 ! X − or expressed in a shorter form,

n 1 − 1 1 , 1 f = i 1 f . {Hλ h | ni} h | ni · Hλ · λ + λ − λ + λ +1 p=0 p p X   On the other hand, using formula (3.2) for , one computes Hλ ∞ 1 f 2 1 f 2, 1 f , 1 f = |h | pi| λ , 1 f + {|h | pi| h | ni} . {Hλ h | ni} −(λ + λ)2 { p h | ni} λ + λ p=0 p p X   Since the left and right hand side of the latter identity are meromorphic functions in λ, they have the same poles. In particular the principal part of the expansion at any of their poles of order two coincide, hence for any p 0 with 1 f = 0 one obtains ≥ h | pi 6 i 1 f if 0 p n 1 , λ , 1 f = − h | ni ≤ ≤ − { p h | ni} 0 if p n . ( ≥ Since, for p 0, λ = p ∞ γ , it follows that for any k 1, ≥ p − k=p+1 k ≥ (6.8) γ ,P1 f = i 1, f δ { k h | ni} h ni kn at least on the set

2 = u Lr,0 : γj(u) > 0 j 1 . O∞ { ∈ ∀ ≥ } Since for any k 1, γk is a real analytic, not identically vanishing ≥ 2 functional, is a dense Gδ subset of Lr,0, and hence by continuity, O∞ 2  identity (6.8) holds on all of Lr,0. ON THE INTEGRABILITY OF THE BO EQUATION 44

For any n 1, on the dense open subset L2 Z of L2 , consisting of ≥ r,0 \ n r,0 potentials u with γn(u) > 0, the argument ϕn of 1 fn is well defined. 1/2 1/2 h | i Since 1 f = γn κn , Corollary 6.1 and Proposition 6.2 imply that |h | ni| (6.9) γ ,ϕ = δ , p 1 , u L2 Z . { p n} pn ≥ ∀ ∈ r,0 \ n In order to show that the angle variables ϕn, n 1, are conjugate to ≥ 2 the action variables γp, p 1, it remains to establish on Lr,0 (Zp Zn) the complementary identities≥ \ ∪ (6.10) ϕ ,ϕ =0 . { p n} In the subsequent section, we prove these identities for 1 p, n N with N 1 arbitrary on the set ≤ ≤ ≥ = u L2 : γ (u) > 0 1 j N; γ (u)=0 j N +1 . ON { ∈ r,0 j ∀ ≤ ≤ j ∀ ≥ } In fact, we show that is a symplectic manifold of dimension 2N ON with γ1,...,γN and ϕ1,...,ϕN as action and angle variables in the classical sense.

7. Finite gap potentials

In this section we study the set GJ of J gap potentials, introduced in Section 2, with J N finite (cf. Definition− 2.2). In fact, to simplify the exposition we will⊂ consider only the case J = 1,...,N where N 1 is an arbitrary integer. Then { } ≥ 2 G 1,...,N = u Lr : γj(u) > 0 1 j N; γj(u)=0 j > N . { } { ∈ ∀ ≤ ≤ ∀ } Note that 2 G 1,...,N = (c + N ) , N = G 1,...,N Lr,0 , { } O O { } ∩ c R [∈ where c + N are the level sets of the restriction of the average func- 2 O R 2 tional Lr , u u 1 , to G 1,...,N . Furthermore, Lr,0 is endowed with the symplectic→ 7→ form h | i { } 1 (7.1) ω(u, v)= u ∂− v , h | x i so that one has the identity F,G = ω(∂x( F ), ∂x( G)). Note that the average functional is a Casimir{ } for the Gardner∇ bracket.∇ The aim of this section is to characterize elements in and to show that the ON restriction of Φ to N is a real analytic symplectic diffeomorphism onto its image where weO recall that Φ denotes the Birkhoff map, introduced in (7.1). As an important application we will prove that Φ is onto (cf. Corollary 7.1). We refer to Appendix B for a study of the set of one gap potentials G N with N 1 arbitrary. { } ≥ Actually, we consider a slightly larger set than N : for any N 1, define O ≥ (7.2) := u L2 : γ (u) > 0 , γ (u)=0 j > N . UN { ∈ r,0 N j ∀ } ON THE INTEGRABILITY OF THE BO EQUATION 45

It is convenient to define 0 := 0 . Note that for any u N and n N one has by (3.13) andU Lemma{ } 2.4 ∈ U ≥ i(n N)x λn(u)= n , fn(x, u)= e − fN (x, u) .

Furthermore, denote by CN [z] respectively C N [z] the space of poly- ≤ nomials P in the complex variable z of degree N of degree at most  N and by C+ [z] the open subset of polynomials P C [z] with the N ∈ N property that P (z)=0 z > 1 . We identify CN [z] with the  N { } ⊂ {| | } space C C∗ of coefficients of polynomials in CN [z]. The following theorem characterizes× elements in . UN Theorem 3. For any N 1 ≥ Q (eix) (7.3) = u = h + h : h(x)= eix ′ , Q C+ [z] UN − Q(eix) ∈ N  where Q′(z) := ∂zQ(z). Furthermore, N is a connected, real analytic, 2 U symplectic submanifold of Lr,0 of dimension 2N. The restriction ΦN of Φ to N , U N 1 ΦN : N C − C∗, u (ζn(u))1 n N U → × 7→ ≤ ≤ is a real analytic, symplectic diffeomorphism,

N

(ΦN ) ω = i dζn dζn . ∗ ∧ n=1 X N In particular, is a dense open subset of with Φ ( )=(C∗) . ON UN N ON Remark 7.1. (i) Alternatively, potentials u N with N 1 can be written in the form (cf. formula (7.9) below) ∈ U ≥

N 1 N 1 ∞ − ∞ − k ikx k ikx (7.4) u(x)= ( qj )e + ( qj )e− j=0 j=0 Xk=1 X Xk=1 X where 1/q0,..., 1/qN 1 C∗ denote the roots of the polynomial Q(z) of degree N, defined in −(7.6∈) below. Note that the above formula for u can 1 r2 also be expressed in terms of the Poisson kernel Pr(x)= 1 2r cos(− x)+r2 , − N 1 − u(x)= (P (x + α ) 1) , q = r eiαj 0 j N 1 . rj j − j j ∀ ≤ ≤ − j=0 X (ii) Real coordinates on are given by the real and imaginary parts UN (Qj,1)1 j N , (Qj,2)1 j N , of the coefficients (Qj)1 j N of the polyno- ≤ ≤ ≤ ≤ ≤ ≤ mial Q C+ [z] normalized by Q(0) = 1, Q(z)=1+ N (Q + ∈ N j=1 j,1 iQ )zj. The corresponding potential u can then be written as j,2 P ix ix ix Q′(e ) ix P ′(e− ) u(x)= e ix e− ix − Q(e ) − P (e− ) ON THE INTEGRABILITY OF THE BO EQUATION 46 where N P (z)=1+ P zj , P := Q iQ . j j j,1 − j,2 j=1 X Note that this formula for u is real analytic in the coordinates Qj,1, Qj,2, 1 j N. Complex valued coefficients Qj,1, Qj,2, 1 j N then will≤ lead≤ to complex valued potentials. ≤ ≤ Furthermore note that by Proposition 4.1, for any n 1, the real part ≥ ζn,1 and the imaginary part ζn,2 of ζn are real analytic functionals on 2 Lr,0 and so are their restrictions to N . It then follows that ΦN is real U 2(N 1) 2 analytic when viewed as map R − R (0, 0) . UN → × \{ } Proof. We first prove that any element u N has the claimed form. To simplify notation within this proof, we∈ do U not indicate the depen- dence of various quantities on u. Our starting point is the inverse formula (4.5), stated in Lemma 4.1, 1 Πu(z)= (Id zM)− X Y 2 , z C, z < 1, h − | iℓ ∈ | | where X = (λp 1 fp )p 0 and Y = ( 1 fn )n 0. According to the computations− in theh | proofi ≥ of Propositionh4.2| ,i the≥ n–th row of M with n N in the case at hand satisfies ≥ (Mnp)p 0 =(δn+1,p)p 0 n N. ≥ ≥ ∀ ≥ Since the coefficients of X = ( λp 1 fp )p 0 with p N + 1 vanish, 1 − h | i ≥ ≥ (ξn(z))n 0 := (Id zM)− [X] satisfies ≥ − ξ (z) zξ (z)= λ 1 f , ξ (z) zξ (z)=0, n N +1, N − N+1 − N h | N i n − n+1 ∀ ≥ or n N 1 ξ (z)= z − − ξ (z) , n N +1. N+1 n ∀ ≥ 2 Combined with the fact that n 0 ξn(z) < for any z < 1, we infer that ≥ | | ∞ | | P ξ (z)=0 , n N +1 , ξ (z)= λ 1 f . n ≥ N − N h | N i This means that 1 (7.5) Πu(z)= (Id zM )− X Y CN+1 , h − N N | N i where

XN =( λp 1 fp )0 p N ,YN =( 1 fn )0 n N , MN =(Mnp)0 n,p N . − h | i ≤ ≤ h | i ≤ ≤ ≤ ≤ Notice that the last row of MN is 0. With the same arguments used to derive the formula for Πu(z) (cf. (4.3)) one obtains a similar formula for the extension fn(z) of the eigenfunctions fn of Lu to the unit disc. For 0 n N and z < 1 one has ≤ ≤ | | 1 f (z)= (Id zM )− 1 Y CN+1 , n h − N n| N i ON THE INTEGRABILITY OF THE BO EQUATION 47 where 1n := (δpn)0 p N . We thus conclude that Πu and f0,...,fN belong to the C vector≤ ≤ space − P (z) (7.6) N := : P C N [z] , Q(z) := det(Id zMN ), R Q(z) ∈ ≤ −   and that f0,...,fN are N + 1 linearly independent elements in N and thus form a basis. As a consequence, 1/Q(z) belongs to the HardyR space of the unit disc, hence Q(z) cannot vanish at any point z with z 1. Furthermore one has Q(0) = 1. To see that deg(Q)= N first | | ≤ note that the last row of MN is identically zero and thus

(7.7) det(Id zMN ) = det(Id zMN 1), − − − where MN 1 = (Mnp)0 n,p N 1. It is to prove that det(MN 1) = 0 . The formula− (4.7) for the≤ coefficients≤ − M for ζ = 0 and the− defini-6 np n+1 6 tion (4.9) of µn+1 imply that in the case at hand the coefficients Mnp with 0 n, p N 1 read ≤ ≤ − δp,n+1 if ζn+1 =0 , (7.8) Mnp = fp 1 µ γ h | i if ζ =0 . (√ n+1 n+1 (λp λn 1) fn+1 1 n+1 − − h | i 6 Expanding det(MN 1) by the formula for Cauchy determinants, one sees that −

fp 1 det(MN 1) = det √µnγn h | i − | | (λp λn 1 1) fn 1 n J,p J˜  − − − h | i ∈ ∈

f0 1 1 = √µnγn h | i det = 0 ! fN 1 λn 1 +1 λp n J,p J˜ 6 n J |h | i|  − −  ∈ ∈ Y∈ where J := 1 j N γ > 0 , J˜ := 0 J N . Thus we showed j that deg(Q)={ ≤N. ≤ | } { } ∪ \{ } Next we prove that Πu(z)= zQ′(z)/Q(z). Indeed since u 1 = 0, − h | i there exists a polynomial R C N 1[z] so that Πu(z)= zR(z)/Q(z). Furthermore, since L f = λ∈f ≤and− since f ,...,f is a basis of , u n n n 0 N RN Lu leaves N invariant. Let us look at the action of Lu on N in more detail. ForR any h L2 , one has R ∈ + Tuh = (Πu)h + Hh(Πu) , Hh(f) := Π(hf) By elementary properties of Hankel operators — see e.g. [8], Prop. 11—, for every h , the range of the Hankel operator H is included ∈ RN h in N . Hence the linear map h Luh Hh(Π(u)) = i∂xh (Πu)h leavesR invariant. Expressed in7→ an alternative− way, it− means− that RN R(z) z∂ z z − Q(z) 1 leaves C N [z] invariant. Using the latter fact, one verifies in a Q(z) ≤ straightforward way that R(z) = Q′(z), showing that u has indeed − ON THE INTEGRABILITY OF THE BO EQUATION 48 the claimed form Q (z) Π(u)(z)= z ′ . − Q(z) C+ Now let us prove the converse. Assume that Q N [z] with Q(0) = 1. Thus Q is of the form ∈ N 1 − Q(z)= (1 q z) , q C , 0 < q < 1 , 0 j N 1 , − j j ∈ | j| ∀ ≤ ≤ − j=0 Y and in turn N 1 − q z (7.9) Πu(z)= j . 1 q z j=0 j X − Let N 1 − z q f(z) := − j . 1 q z j=0 j Y − Clearly, f is in N . Writing Lu(f) as z∂zf (Πu)f Hf (Πu) one computes R − − N 1 − z q z q z q L (f)= f + j j j = Nf . u z q 1 q z − 1 q z − z q j=0 j j j j X  − − − −  It means that f is an eigenfunction of Lu with eigenvalue N. Further- more, since for any z = 1, | | f(z)= zN Q(z)/Q(z) , k one verifies that for any k 1, S f is orthogonal to N and thus in particular to Πu. By Lemma≥2.1 this implies that for anyR k 1, Skf is ≥ an eigenfunction of Lu with eigenvalue N +k. Taking into account that is of dimension N + 1 and that the eigenvalues λ satisfy (3.13) RN n with Πu 1 = 0 one concludes that λN+k = N + k for any k 0 and h | i N N ≥1 γ = 0 forany k 1. On the other hand, 1 f =( 1) − q =0 N+k ≥ h | i − j=0 j 6 implying that γN = 0. Altogether this shows that u N . 6 ∈ U Q Having established identity (7.3) we know that Π( N ) is a connected, complex manifold of dimension N, parametrized byU the open subset of CN C+ described by the coefficients of the polynomials Q N [z] satisfy- ing Q(0) = 1. Furthermore, since for any u, v L2 , ∈ ∈ r,0 1 1 ω(u, v)= u ∂− v = 2Im D− Πu Πv h | x i − h | i 1 Π( N ) is a K¨ahler manifold with Hermitian form D− Πu Πv and U h 2 | i hence N is a real analytic symplectic submanifold of Lr,0 of real di- mensionU 2N. We claim that, for every 1 n N, γ does not identi- ≤ ≤ n cally vanish on N . Indeed, denote by J the set of such indices n. By the definition ofU ,γ (u) = 0 for any u and hence N J. Since UN N 6 ∈ UN ∈ each γn is a real analytic functional and since N is connected, GJ N is an open dense subset of . From the proofU of Proposition 4.2∩Uwe UN ON THE INTEGRABILITY OF THE BO EQUATION 49 know that the coefficients Mnp of the matrix MN can be expressed in terms of the Birkhoff coordinates (ζj)1 j N . Similarly, this is the case ≤ ≤ for the vectors XN and YN , introduced above. Hence the right hand side of the identity (7.5) locally extends analytically, yielding a real an- J alytic map ΨN on (C∗) so that for any u in the open subset GJ N of , ∩ U UN Ψ Φ (u)=Π(u) . N ◦ N This implies that the rank of the map ΦN at any point in GJ N has to be 2N. Since in view of the definition of J it cannot be bigger∩U than 2♯J, we infer that J = 1, 2,...,N . As a consequence, N = G 1,...,N N is an open dense subset{ of } and a symplectic submanifoldO { of } ∩Uof UN UN dimension 2N. On N , the smooth functions γ1,...,γN , ϕ1,...,ϕN , referred to as actionsO and angles, are well defined. From Corollary 6.1 and identity (6.9), we know that (7.10) γ ,γ =0 , γ ,ϕ = δ , 1 n, p N. { n p} { n p} pn ≤ ≤ This implies that Φ is a local diffeomorphism on and therefore, N ON since ΦN is one to one, a diffeomorphism onto its range which is con- N N tained in (C∗) . In order to prove that the range is (C∗) , we observe N that (C∗) is connected and ΦN proper (cf. Proposition 4.1). Being N a local diffeomorphism, ΦN is open and hence ΦN ( N )=(C∗) . To prove that Φ is symplectic on introduce the two–formO N ON N ω˜ = ω dγ dϕ . − n ∧ n n=1 X In view of the commutation relations (7.10), ∂ ⌟ ω˜ =0 , 1 p N. ϕp ∀ ≤ ≤ Since moreoverω ˜ is a closed two form there exist smooth functions anp so that (7.11)ω ˜ = a (γ ,...,γ ) dγ dγ . np 1 N n ∧ p 1 n

1 iϕn iϕn dζn = e dγn + √γne idϕn 2√γn ON THE INTEGRABILITY OF THE BO EQUATION 50 and thus dζn dζn = idγn dϕn. When expressed in the coordinates ζ , the pull back∧ of the− symplectic∧ form ω to is therefore n ON N ω = i dζ dζ . n ∧ n n=1 X Since is dense in , this identity holds on all of . Using again ON UN UN that ΦN is one to one and proper, we conclude that ΦN is a symplectic N 1 diffeormorphism from N onto C − C∗. By Remark 7.1 (ii), ΦN is a real analytic map andU hence the proof× of Theorem 3 is complete. 

As an important application of Theorem 3, we prove that the Birkhoff map Φ is onto. Recall that u L2 is a finite gap potential if J(u)= ∈ r,0 n 1 γ (u) > 0 is finite. The set of finite gap potentials in L2 is { ≥ | n } r,0 thus given by 0 N 1 N . { } ∪ ≥ U S 2 1/2 Corollary 7.1. The map Φ: Lr,0 h+ ,u (ζn)n 1 is a homeo- −→ 7→ ≥ 2 morphism and the set of finite gap potentials is dense in Lr,0. Further- more, the identities ζ ,ζ =0 , ζ , ζ = iδ , n, m 1 , { n m} { n m} − nm ≥ 2 are valid on all of Lr,0. Proof. We already know that Φ is continuous and proper (cf. Proposi- tion 4.1) and one to one (cf. Proposition 4.2). To prove that Φ is onto, 1/2 let ζ := (ζn)n 1 be any nonzero sequence in h+ N 1 Φ( N ) where N ≥ \ ≥ U U is defined by (7.2). Then there exists an increasing sequence (Nk)k 1 S ≥ with Nk such that ζNk = 0 for any k 1. We approximate ζ by the truncated→∞ sequence 6 ≥

(k) ζn if n Nk ζn = ≤ (0 if n > Nk . By Theorem 3, there exists a unique element u(k) such that ∈ UNk Φ(u(k))= ζ(k). Since Φ is proper and continuous, u(k) has a limit point 2 u in Lr,0 which satisfies Φ(u)= ζ . Thus we have proved that Φ is onto. Since Φ is proper and continuous 1 it then follows that Φ− is continuous as well and in turn N 1 N is ≥ U dense in L2 . It remains to verify the claimed Poisson bracket relations. r,0 S By Proposition 5.1, for any n, m 1, the Poisson brackets ζ ,ζ ≥ { n m} and ζ , ζ are real analytic functionals on L2 . Since by Theorem 3 { n m} r,0 the claimed Poisson bracket relations hold on N 1 N and N 1 N ≥ U ≥ U is dense in L2 it then follows that they hold on all of L2 .  r,0 S r,0 S ON THE INTEGRABILITY OF THE BO EQUATION 51

8. Proof of Theorem 1 and Theorem 2 In this section, we make a detailed synopsis of the proof of Theorem 1. As an application, we derive formulas for the BO frequencies at the end of this section (cf. Proposition 8.1) and prove Theorem 2. Proof of Theorem 1 By Corollary 7.1, Φ is a homeomorphism and the Poisson bracket identities (B2) hold. It remains to express the BO Hamiltonian (u), defined for u H1/2, in terms of the Birkhoff H ∈ r,0 coordinates ζ , n 1. For u H1/2 one computes, using that n ≥ ∈ r,0 Πu = ∞ λ 1 f f , − n=0 nh | ni n 2π P1 1 1 ∞ (u)= D u u u3 dx = L (Πu) Πu = λ3 1 f 2 . H 2h| | | i − 2π 3 h u | i n|h | ni| n=0 Z0 X With the help of ˜ , defined by (3.8) as Hε ∞ 1 f 2 ˜ = |h | ni| , Hε 1+ ελ n=0 n X one can express (u) as follows, H 1 d3 1 d3 (8.1) (u)= ˜ = log ˜ . H −6 dε3 ε=0Hε −6 dε3 ε=0 Hε d ˜ Using the formula for dε log ε derived in Section 3 (cf. (3.11)), one then obtains H

∞ 1 3 1 2 2 (u)= λ0 + γn λn +(λn 1 + 1) + λn(λn 1 + 1) H 3 3 − − n=1 X  or using that by (3.13), λ = n s with s := ∞ γ , n 0, n − n+1 n+1 k=n+1 k ≥ 1 ∞ 2 P 1 (8.2) (u)= s3 + γ n s + γ n s + γ2 . H −3 1 n − n n − n 3 n n=1 X    This proves Theorem 1.  Actually, formula (8.2) of the BO Hamiltonian (u) can be further simplified as follows. H Proposition 8.1. For any u H1/2, ∈ r,0 ∞ ∞ ∞ (8.3) (u)= n2γ γ 2 . H n − k n=1 n=1 k=n X X X  ∂ As a consequence, the BO frequencies ωn = H , n 1, read ∂γn ≥ ∞ (8.4) ω = n2 2 min(k, n)γ . n − k Xk=1 ON THE INTEGRABILITY OF THE BO EQUATION 52

Thus for any u L2 , the solution u of the BO equation (1.1) with 0 ∈ r,0 u(0) = u0, when expressed in Birkhoff coordinates, is given by (8.5) ζ (u(t)) = ζ (u ) eiωn(u0)t , t R. n n 0 ∈ Remark 8.1. (i) Note that the right hand side of (8.4) is well defined for any u L2 . Hence the BO frequencies continuously extend to L2 . ∈ r,0 r,0 (ii) Note that by the trace formula of Proposition 3.1 (ii), the formula (8.4) for ωn can be written as

∞ ω = n2 u 2 +2 (k n)γ . n −k k − k k=Xn+1 The sum 2 ∞ (k n)γ satisfies the asymptotics k=n+1 − k P ∞ 2 (k n)γ = o(1) as n − k →∞ k=Xn+1 2 2 and hence in particular, limn (ωn n )= u . →∞ − −k k Proof. We view (u) as a function of the actions γn, n 1, given by formula (8.2) andH write the claimed identity (8.3) as ≥

∞ ∞ (8.6) (u)= n2γ s2 . H n − n n=1 n=1 X X Note that both sides are real analytic functionals on the space ℓ1,2(N, R), defined as the of real sequences (ξn)n 1 satisfying ≥ ∞ n2 ξ < + . | n| ∞ n=1 X ∂ By (8.2), the partial derivative H , r 1, is given by ∂γr ≥ ∂ H = s2 +(r s )2 γ (2(n s )+ γ ) , ∂γ − 1 − r − n − n n r n

Appendix A. On the Lax pair for the BO equation The purpose of this appendix is to derive the Lax pair formulation of the Benjamin-Ono equation, described in Section 1. The correspond- ing computations for the Benjamin-Ono equation on the line can be found in [32, Appendix A]. In order to be comprehensive, we include its derivation. 2 Recall that for any u Lr, we introduced the (unbounded) operators 2 ∈ Lu and Bu on L+, L = i∂ T , B := i∂2 +2T 2∂ T u − x − u u − x ∂x(Πu) − x u 2 2 where Tu : L+ L+ denotes the Toeplitz operator given by Tuf = →2 2 Π(uf) and Π : L L+ the Szeg˝oprojector. We claim that for any smooth function u(→t, x) on R T, and any h L2 × ∈ + d (A.1) L + [L , B ] h = Π ∂ u +2u∂ u H∂2u h dt u u u − t x − x   To verify this identity, we write   [L , B ]= i[∂ , B ] [T , B ]= Q + Q + Q u u − x u − u u 0 1 2 where for any 0 k 2, Qk is an operator, homogenous of degree k in u. One computes≤ Q≤ = [∂ , ∂2] = 0 and 0 − x x Q = 2i[∂ , T ∂ T ]+ i[T , ∂2] , 1 − x ∂x(Πu) − x u u x Q = 2[T , T ∂ T ] . 2 − u ∂x(Πu) − x u Let us first discuss Q . For any h L2 , one has 2 ∈ +

Tu T∂x(Πu)h = Π ∂x(Πu)uh , T∂x(Πu) Tuh = Π ∂x(Πu)Π(uh)). Since Π(I Π) = 0 one then obtains −    [T , T ]h = Π Π(∂ u)(I Π)(uh) = Π ∂ u)(I Π)(uh) . u ∂x(Πu) x − x − On the other hand, [Tu, ∂xTu]h = Π ∂xuΠ(uh) . Combining the two results then leads to −  Q h = 2Π (∂ u)uh = Π (2u∂ u)h . 2 − x − x Concerning the operator Q1, note that 

2i[∂ , T ]= 2iT 2 , 2i[∂ , ∂ T ]=2i∂ T , − x ∂x(Πu) − ∂x(Πu) x x u x ∂xu and 2 i[T , ∂ ]= iT 2 2iT ∂ . u x − ∂xu − ∂xu x Combining the three results then leads to 2 2 Q h = 2iT 2 h + iT 2 h = 2iΠ Π(∂ u)h + iΠ (∂ u)h 1 − ∂x(Πu) ∂xu − x x 2 2 2 or Q1h = Π i(2Π I)(∂xu)h . Since i(2Π I)(∂xu) = H∂xu we conclude that− − − −  2 Q1h = Π (H∂xu)h .  ON THE INTEGRABILITY OF THE BO EQUATION 55

Finally, one has ∂tLuh = Π(∂tuh). Altogether, we thus have proved that − d L + [L , B ] (h)= Π ∂ u +2u∂ u H∂2u h . dt u u u − t x − x    

Appendix B. Traveling waves and one gap potentials In this appendix we identify traveling wave solutions of the Benjamin Ono equation as one gap potentials, which we describe in more detail, obtaining a new proof of the result of Amick–Toland in [2]. As we already observed, we may restrict our analysis to solutions with average 0. Proposition B.1. The nonzero traveling wave solutions of the Benjamin– Ono equation coincide with the solutions with initial data given by one gap potentials.

Proof. A traveling wave solution is of the form u(t, x)= u0(x + ct) for some c R, or, with our notation, u(t) = Qctu0. In view of identity (2.17), we∈ have ζ (Q v)= eiτnζ (v) , v L2 . n τ n ∈ r,0 Comparing with the general expression (8.5) of the BO solution, we infer icnt iωn(u0)t e ζn(u0)= e ζn(u0) . Hence for any n 1 with ζ (u ) = 0, one has cn = ω (u ). In ≥ n 0 6 n 0 view of formula (8.4), the mapping n ωn/n is strictly increasing. Consequently, there cannot be more than7→ one n with ζ (u ) =0. If u is n 0 6 0 not identically 0, this precisely means that u0 is a one gap potential. 

2 Let N 1 be given and introduce N := G N Lr,0. Without further reference,≥ we use the notation introducedO{ } { in} the∩ main body of the paper. Denote by Φ N the restriction of the Birkhoff map Φ to { } N . It is given by O{ } 1 Φ N : N C∗, u ζN (u)= 1 fN { } O{ } → 7→ √κN h | i where κN =1/(N + γN ) and hence 1 (B.1) 1 fN = ζN . h | i √N + γN By Theorem 3 in Section 7,

Πu(z)= zQ′(z)/Q(z) − ON THE INTEGRABILITY OF THE BO EQUATION 56 where by (7.7), Q(z) = det(Id zMN 1) and MN 1 is the N N − − − × matrix (Mnp)0 n,p N 1 with coefficients Mnp given by ≤ ≤ −

δp,n+1 if ζn+1 =0 , Mnp = fp,1 µn+1γn+1 h i if ζn+1 =0 . (√ (λp λn 1) fn+1 1 − − h | i 6 Since u G N and hence 1 fp = 0 for any p =0, N one computes ∈ { } h | i 6 √µN γN f 1 Q(z)=1+ zN h 0| i N f 1 h N | i where by (4.9) γ N µ =1 N = . N − λ λ N + γ N − 0 N To compute 1 f , use again that 1 f = 0 for k = 0, N to infer h | 0i h | ki 6 1 = 1 f0 f0 + 1 fN fN . By the definition of f0 one has 1 f0 > 0 and henceh | i h | i h | i

2 1/2 √N (B.2) 1 f0 = (1 1 fN ) = . h | i − |h | i| √N + γN

iϕN Since ζN = √γN e and fN 1 = 1 fN , one obtains, using (B.1) and (B.2), h | i h | i

N √µN γN √N N 1 Q(z)=1+ z =1+ z ζN . N ζN √N + γN Combining the results obtained leads to the following formula wzN 1 Πu(z)= zQ′(z)/Q(z)= N , w := 1 f = ζ − 1 wzN −h | N i −√N + γ N − N or by substituting eix for z, weiNx we iNx u(x)= N + N − . 1 weiNx 1 we iNx − − − Note that 0 < w = √γN < 1, showing that the potential u(x + iy) | | √N+γN is analytic in a strip y

Appendix C. On the spectrum of Lu 2 By Proposition 2.1, for any u Lr, the spectrum of Lu is discrete and the eigenvalues λ (u), n 0,∈ of L satisfy λ (u) < λ (u) < n ≥ u 0 1 ··· with γn(u) = λn(u) λn 1(u) 1 0 for any n 1. The purpose of this appendix is to provide− − a spectral− ≥ interpretation≥ of the real numbers γn(u), n 1, by considering the operator Lu on the real line R. Denote by≥L2(R) L2(R, C) the standard L2 space with correspond- ≡ − ing inner product f,g = ∞ f(x)g(x)dx and for any f L2(R) by h i ∈ (f) its Fourier transform, −∞ F R ∞ ixη (f)(η)= f(x)e− dx . F Z −∞ 2 R 2 R Furthermore, let L+( ) be the closed subspace of L ( ) given by L2 (R)= f L2(R): (f)(η)=0 η < 0 , + { ∈ F ∀ } R 2 R 2 R 1 R Π : L ( ) L+( ) the corresponding orthogonal projector, H ( ) 1 R C → 1 1 R 1 R ≡ H ( , ) the standard H Sobolev space and H+( ) := H ( ) 2 R 2 − R 2 R ∩ L+( ). For any u Lr, denote by Lu the operator on L+( ) with 1 R ∈ domain H+( ) given by LRf = i∂ f ΠR(uf) , f H1 (R) . u − x − ∀ ∈ + R To see that Lu is selfadjoint, we first need to establish the following Lemma C.1. For any ǫ> 0, u L2, and f H1(R), ∈ r ∈ ∞ 2 2 2 1 2 uf dx u ǫ f ′ 2 R + 1+ f 2 R . | | ≤k k k kL ( ) ǫ k kL ( ) Z   −∞   Proof. For any u L2 and f H1(R) one has ∈ r ∈ 2π + ∞ ∞ u(x)f(x) 2 dx = u(x) 2 f(x +2nπ) 2 dx . | | | | | | Z Z0 n= −∞ X−∞ Let g(t; x) := f(x +2πt). It then suffices to show that for any 0 x 2π, ≤ ≤ + ∞ 2 2 1 2 g(n; x) ǫ f ′ 2 R + 1+ f 2 R . | | ≤ k kL ( ) ǫ k kL ( ) n=   X−∞ For any given 0 x 2π, n Z, and t [n, n + 1], ≤ ≤ ∈ ∈ t g(t; x) 2 g(n; x) 2 = 2Re g(s; x)∂ g(s; x) ds | | −| | s Z n  ON THE INTEGRABILITY OF THE BO EQUATION 58 and hence n+1 1 g(t; x) 2 g(n; x) 2 ε ∂ g(s; x) 2 + g(s; x) 2 ds . | | −| | ≤ | s | ε| | Zn  Integrating in t [n, n +1], one gets ∈ n+1 1 g(n; x) 2 ε ∂ g(t; x) 2 + 1+ g(t; x) 2 dt | | ≤ | t | ε | | Z n    and then summing over n Z yields the claimed inequality.  ∈ 2 R 1 R Clearly, i∂x is a selfadjoint operator on L+( ) with domain H+( ). By the Kato-Rellich− theorem (cf. [27]) one then infers from Lemma C.1 R that Lu is selfadjoint as well. 2 R R Proposition C.1. For any u Lr, the spectrum spec(Lu ) of Lu is absolutely continuous and consists∈ of a union of bands,

R ∞ spec(Lu )= [λn(u),λn(u) + 1] . n=0 [ Hence, for any n 1, γn(u) is the length of the nth gap in the spectrum R ≥ of Lu . To prove Proposition C.1, we first need to make some preliminary considerations. Introduce the following L2 space − 1 L2 T [0, 1]) L2 T [0, 1], C; dxdξ × ≡ × 4π2 2 2 T 1  which is clearly isometric to L [0, 1], L ( ); 2π dξ . Furthermore for any f L2(R), denote by ψ[f] L2(T [0, 1]) the function ∈ ∈ ×  ψ[f](x, ξ)= einx (f)(ξ + n) . F n Z X∈ Since 2π 1 ∞ 1 (C.1) f(x) 2 dx = ψ[f](x, ξ) 2 dξ dx , | | 4π2 | | Z Z0 Z0 −∞   1 1 (C.2) f(x)= eixξ ψ[f](x, ξ) dξ , 2π Z0 it follows that the linear map L2(R) L2 T [0, 1] , f ψ[f] is → × 7→ unitary. Restricting it to the Sobolev space H1(R) H1(R, C) one 1 R 2 1 T 1 ≡ sees that H ( ) is isometric to L ([0, 1],H ( ); 2π dξ). Finally, we claim that for any 0 ξ < 1, ≤ (C.3) ψ[ΠRf]( ,ξ) = Π ψ[f]( ,ξ) . · ·  ON THE INTEGRABILITY OF THE BO EQUATION 59

R inx R Indeed, for any 0 ξ < 1, ψ[Π f](x, ξ)= n Z e (Π f)(ξ + n) is given by ≤ ∈ F P ∞ ψ[ΠRf](x, ξ)= einx (f)(ξ + n)=Π(ψ[f]( ,ξ))(x) . F · n=0 X R inx Note that for ξ = 1, one has ψ[Π f](x, 1) = n∞= 1 e (f)(1 + n) ix − F which equals Π(ψ[f]( , 1))(x)+ e− (f)(0) . · F P Lemma C.2. For any f H1 (R) and 0 ξ < 1, ∈ + ≤ ψ[LRf]( ,ξ)=(ξ + L )(ψ[f]( ,ξ)) . u · u · Proof. Let f H1 (R) and 0 ξ < 1. In view of the definition of ψ[f], ∈ + ≤ ψ[Df](x, ξ)=(ξ + Dx)ψ[f](x, ξ) . Since u is 2π periodic, one has ψ[uf](x, ξ)= u(x)ψ[f](x, ξ) and hence one infers from (C.3) that ψ[ΠR(uf)]( ,ξ) = Π ψ[uf]( ,ξ) = Π uψ[f]( ,ξ) = T ψ[f]( ,ξ) . · · · u ·  Combining these computations one obtains the claimed formula. 2 2 R Proof of Proposition C.1. Let u Lr. For any f L+( ) introduce u ∈ ∈ the sequence (ψn[f](ξ))n 0 where for any 0 ξ 1, ≥ ≤ ≤ ψu[f](ξ) := ψ[f]( ,ξ) f ( ,u) . n h · | n · i Then 2 R 2 Z 2 1 u L+( ) ℓ n 0, L [0, 1], dξ , f (ψn[f])n 0 → ≥ 2π 7→ ≥ is a , and, for every f H1 (R), Lemma C.2 implies ∈ + (C.4) ψu[LRf](ξ)=(ξ + λ (u))ψu[f](ξ) , 0 ξ < 1 . n u n n ≤ 1 R Consequently, the spectral measure µf of a function f H+( ), acting on any given continuous function ϕ : R R with compact∈ support, ϕ(λ) dµ (λ)= ϕ(LR)f f , can be computed→ as R f h u | i 1 R ∞ ∞ 1 ϕ(λ) dµ (λ) = ϕ(ξ + λ (u)) ψu[f](ξ) 2 dξ f 2π n | n | Z n=0 Z0 −∞ X λn(u)+1 1 ∞ = ϕ(λ) gf (λ) 2 dλ , 2π | n | n=0 Z X λn(u) f u where for any λn(u) λ<λn(u)+1, gn(λ) is given by ψn[f](λ λn(u)) . ≤ R − This proves that the spectrum of Lu is absolutely continuous and is contained in the union ∞ [λn(u),λn(u) + 1] . n=0 [ ON THE INTEGRABILITY OF THE BO EQUATION 60

R In order to prove that spec(Lu ) coincides with the latter union it is enough to show that any measurable subset E of [0, 1] with the property that, for some n, (C.5) 1 (ξ) ψu[f](ξ)=0 , 0 ξ 1, f L2(R) , E n ∀ ≤ ≤ ∈ has Lebesgue measure 0. Note that (C.5) says that 2π 1 1 (ξ)ψ[f](x, ξ) f (x, u) dxdξ =0 , f L2(R) , E n ∀ ∈ Z0 Z0 or, since ψ is unitary, 1 ψ− [f 1 ]=0 . n ⊗ E In view of (C.2), one then concludes that 1 f (x) 1 (ξ)eixξ dξ =0 , x R . n E ∀ ∈ Z0

Since fn is continuous and is not identically 0, the Fourier transform of 1E then vanishes on some nonempty open set of R. Since the Fourier transform of 1E is an entire function, it then must vanish identically, implying that E has measure 0. This completes the proof of Proposition C.1. 

References [1] L. Abelouhab, J. Bona, M. Felland, J.-C. Saut, Non local models for nonlinear dispersive waves, Physica D, Nonlinear Phenomena, 40(1989), 360– 392 [2] C. Amick, J. Toland, Uniqueness and related analytic properties for the Benjamin-Ono equation a nonlinear Neumann problem in the plane, Acta Math., 167(1991), 107–126 [3] T. Benjamin, Internal waves of permanent form in fluids of great depth, J. Fluid Mech., 29(1967), 559–592 [4] T. Bock, M. Kruskal, A two-parameter Miura transformation of the Benjamin-Ono equation,Phys. Lett. A, 74(1979), 173–176 [5] R. Coifman, M. Wickerhauser, The scattering transform for the Benjamin-Ono equation, Inverse Problems 6(1990), 825–861 [6] S. Dobrokhotov, I. Krichever, Multi-phase solutions of the Benjamin- Ono equation and their averaging, Math. Notes 49(1991), no. 5–6, 583–594 [7] A. Fokas, M. Ablowitz, Inverse scattering transform for the Benjamin-Ono equation: A pivot to multidimensional problems, Stud. Appl. Math. 68(1983), 1–10 [8] P.Gerard,´ S. Grellier, The cubic Szeg˝oequation, Ann. Sci. Ec.´ Norm. Sup´er. (4) 43(2010), no. 5, 761–810 [9] P.Gerard,´ S. Grellier, Invariant tori for the cubic Szeg˝oequation, Invent. Math. 187(2012), no. 3, 707–754 [10] P.Gerard,´ S. Grellier, The cubic Szeg˝oequation and Hankel operators, Ast´erisque No. 389 (2017) ON THE INTEGRABILITY OF THE BO EQUATION 61

[11] P.Gerard,´ T. Kappeler, P. Topalov, in preparation [12] B. Grebert,´ T. Kappeler, The defocusing NLS and its normal form, EMS Series of Lectures in Mathematics, EMS, Z¨urich, 2014 [13] M. Ifrim, D. Tataru, Well-posedness and dispersive decay of small data solu- tions for the Benjamin-Ono equation, Ann. Sci. Ec.´ Norm. Sup´er. (4) 52(2019), no.2, 297-335. [14] T. Kappeler, J. Poschel¨ , KdV&KAM, vol. 45, Ergeb. der Math. und ihrer Grenzgeb., Springer, 2003 [15] T. Kappeler, P. Topalov, Global well-posedness of KdV in H−1(T, R), Duke Journal of Mathematics, 135(2006), no. 2, 327–360 [16] D. Kaup, Y. Matsuno, The inverse scattering for the Benjamin-Ono equa- tion, Stud. Appl. Math. 101(1998), 73–98 [17] B. Levitan, V. Zhikov, Almost periodic functions and differential equations, Cambridge University Press, 1982 [18] H. McKean, E. Trubowitz, Hill’s operator and hyperelliptic function theory in the presence of infinitely many branch points, Comm. Pure Appl. Math., 29(1976), no. 2, 143–226 [19] H. McKean, K. Vaninsky, Action-angle variables for the cubic Schr¨odinger equation, Comm. Pure Appl. Math. 50(1997), no. 6, 489–562 [20] L. Molinet, Global well-posedness in L2 for the periodic Benjamin-Ono equa- tion, American J. Math.130 (3) (2008), 2793–2798 [21] L. Molinet, D. Pilod, The Cauchy problem for the Benjamin-Ono equation in L2 revisited, Anal. and PDE 5 (2) (2012), 365–395 [22] A. Moll, Finite gap conditions and small dispersion asymptotics for the clas- sical Benjamin–Ono equation, arXiv:1901.04089 [23] J. Moser, Finitely many mass points on the line under the influence of an exponential potential – an integrable system, in Dynamical Systems, Theory and Applications, Springer Verlag, New York, 1975, 467–497 [24] A. Nakamura, Backlund transform and conservation laws of the Benjamin- Ono equation, J. Phys. Soc. Japan 47(1979), 1335–1340 [25] M. Nazarov, E. Sklyanin, Integrable hierarchy of the quantum Benjamin– Ono equation, SIGMA Symmetry Integrability Geom. Methods Appl. 9(2013), paper 078, 14 pp [26] H. Ono, Algebraic solitary waves in stratified fluids, J. Physical Soc. Japan 39(1975), 1082–1091 [27] M. Reed, B. Simon, Methods of modern mathematical physics, vol. 2: Fourier analysis, self-adjointness, Academic Press, 1975 [28] J. Satsuma, Y. Ishimori, Periodic wave and rational soliton solutions of the Benjamin-Ono equation, J. Phys. Soc. Japan 46(1979), 681–687 [29] J.-C. Saut, Sur quelques g´en´eralisations de l’´equation de Korteweg-de Vries, J. Math. Pures Appl. 58(1979), 21–61 [30] J.-C. Saut, Benjamin-Ono and Intermediate Long Wave equations: modeling, IST, and PDE, ArXiv:1811.08652, 2018 [31] N. Tzvetkov, N. Visciglia, Invariant measures for BO, Intern. Math. Re- search Notices, vol 2014, no 17, p 4679–4714 [32] Y. Wu, Simplicity and finiteness of discrete spectrum of the Benjamin-Ono scattering operator, SIAM J. Math. Anal. 48(2016), no. 2, 1348–1367 [33] Y. Wu, Jost solutions and the direct scattering problem of the Benjamin-Ono equation, preprint, arXiv.1704.01692.v1, to appear in SIMA. ON THE INTEGRABILITY OF THE BO EQUATION 62

Laboratoire de Mathematiques´ d’Orsay, Univ. Paris-Sud, CNRS, Universite´ Paris–Saclay, 91405 Orsay, France E-mail address: [email protected] Institut fur¨ Mathematik, Universitat¨ Zurich,¨ Winterthurerstrasse 190, 8057 Zurich, Switzerland E-mail address: [email protected]