arxiv: v1 [math.ap] 21 Mar 2011

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES IN KDV arxiv:113.4135v1 [math.ap] 21 Mar 211 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY Abstract. The Korteweg-de Vries (KdV) equation with periodic boundary conditions is considered. The interaction of a periodic solitary wave (cnoidal wave) with high frequency radiation of finite energy (L 2 -norm) is studied. It is proved that the interaction of low frequency component (cnoidal wave) and high frequency radiation is weak for finite time in the following sense: the radiation approximately satisfies Airy equation. The KdV equation 1. Introduction q t +q xxx +q x q = is one of the most basic dispersive partial differential equations (PDE) with solitary wave solutions. There are two types of solitary waves in KdV posed on the real line: exponentially decaying and spatially periodic solitons. This paper deals exclusively with the periodic case and even more restrictively, we consider KdV with periodic boundary conditions. The periodic traveling waves in KdV (already known to Korteweg and de Vries) are called cnoidal waves, as they may be expressed in terms of the elliptic Jacobi function, see e.g. [13] ( (1) φ c (z) = β 2 +(β 3 β 2 )cn 2 β3 β 1 z;k 12 ) Date: October 3, 218. The authors were partially supported by NSF grants DMS-9865 (B. E.), DMS- 91222 (N. T.), and DMS-87897 (V. Z.). 1

2 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY where z = x ct,β 1 < β 2 < β 3,β 1 +β 2 +β 3 = 3c, k 2 = β 3 β 2 β 3 β 1. In this paper we study the dynamics of periodic solutions of the KdV equation with initial data which is a sum of a cnoidal wave and of a high frequency perturbation with finite energy (L 2 norm). For small perturbations this falls within the theory of Lyapunov stability results. In the KdV case, this problem was first considered by Benjamin [2] who conjectured that cnoidal waves are stable with regard to small perturbations. The conjecture was later proved by McKean [12]. The stability of cnoidal waves on the real line with respect to small perturbations with different periods has also been considered and is currently an active area of research, see [3] and the references therein. We, on the other hand, are interested on the behavior of solutions in the case of high frequency perturbation. We consider the evolution of solitary wave and the high frequency perturbation, and prove that the perturbation evolves almost linearly. Within the context of our previous results in [8, 9] where we proved that the evolution of high energy solutions of KdV is nearlinear, our result in this paper can be considered as a superposition principle for a nonlinear dispersive PDE. It has been long suggested in the physics literature that in the regime below collapse, high frequency solutions evolve almost linearly and interact little with the low frequencies. This mechanism has been also used, perhaps implicitly, to prove low regularity results, see e.g. [4, 6, 7]. In a way, our proof of nearly independent evolution of soliton and high frequency radiation, provides some support of this heuristic, see e.g. [14], pg. 118. This nearly independent evolution is due to a subtle averaging effect in the nonlinear dispersive dynamics, making these results possible. Recently KdV was studied with respect to this averaging effect. In [8, 9], nearly linear dynamics was established for high frequency initial data and in [1] a new elegant proof of well-posedness in H s,s was found using explicitly high frequency averaging effects.

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 3 In contrast to the Lyapunov problem, where infinite time stability is usually established, our result is valid only on finite times, which is related to adiabatic invariance phenomena: high frequency wave oscillations are averaged out to produce effective slow evolution. Since KdV as a model is valid only on a finite time scale, it is meaningful to consider the dynamics of the solutions for finite times. Given that even the classical adiabatical invariance theorem (conservation of action of the pendulum with slowly changing frequency) requires careful analysis, our task becomes even harder because of the infinite dimensionality. Here, we use these ideas to study the interaction of a solitary wave (low frequencies) with radiation (high frequencies). Now we state our main theorem: Theorem 1.1. Let φ(x ct) be a 2π-periodic cnoidal wave solution of KdV. Fix s (,1/2). Consider the real valued solution of KdV on T IR with the initial data q(x,) = φ(x)+g(x) satisfying, for some < s < 1/2, Then, for each t >, we have g 2 1, g H s = ε 1. q(x,t) φ(x ct) (e tl g)(x) 2 Ce Ct ε, where L = 3 x φ x, and C depends only on s and φ. As usual, H s is the completion of L 2 under the norm u H s = û(k)/(1+ k 2 ) s/2 l 2 and we use the notation φ = 1 2π φ(r)dr. Remar.2. The statement of the theorem above can be extended to an arbitrary H 4 solution ρ of KdV in the following sense. Given T, there is a constant C = C(s,T,ρ) such that given ε > and g as in the theorem we have where L = 3 x ρ x. q(x,t) ρ(x,t) (e tl g)(x) 2 Cε, t [,T], This variation follows from the proof of Theorem 1.1 by utilizing Remark 2.8 and Remark 4.1.

4 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY As it is well-known, KdV is a completely integrable system with infinitely many conserved quantities. However, our methods in this paper do not rely on the integrability structure of KdV, and thus they can be applied to other dispersive models. On the other hand we use the fact that the smooth solutions of KdV satisfy momentum conservation: π π and the conservation of energy, π π u(x,t)dx = u 2 (x,t)dx = π π π π u(x, )dx, u 2 (x,)dx. The KdV equation is locally well-posed in L 2 (T), [5]. Due to energy conservation it is globally well-posed and the solution is in C(R;L 2 (T)). Kenig, Ponce, Vega, [11], improved Bourgain s result and showed that the solution of the KdV is locally well-posed in H s (T) for any s > 1 2. Later, Colliander, Keel, Staffilani, Takaoka, Tao, [6], showed that the KdV is globally well-posed in H s (T) for any s 1 2 thus adding a local well-posedness result for the endpoint s = 1 2. Recently T. Kappeler and P. Topalov, [1] extended the latter result and prove that the KdV is globally well-posed in H s (T) for any s 1. Since our statements concerns L 2 functions, from the results listed above, we only use the global well-posedness in L 2, [5]. 1.1. Solitary waves description. Since our goal is to illustrate the phenomena, we restrict ourselves to the special case of stationary periodic waves with prescribed period 2π. The proof can be readily extended to arbitrary period. To make the presentation self-consistent, we directly show that there are such waves, which can also be found by using properties of Jacobi functions (1). Let q = f(x ct), and substitute this in KdV Integrating once, we obtain cf = ff +f. a cf = f2 2 +f,

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 5 which can be written in the potential form (2) f +W f (f) =, where W(f) = f3 6 +cf2 2 af. Under the assumption c 2 + 2a >, this cubic polynomial has one local maximum and one local minimum: f = c 2 +2a c, f + = c 2 +2a c. Taking the 2nd derivative of W at f +, W ff (f + ) = f + +c = c 2 +2a, we obtain the period of small oscillations: T (a,c) = 2π Wf (f + ) = 2π 4 c 2 +2a. Therefore, moving through the family of periodic solutions nested between the minimum and the separatrix, we will see the period assuming all intermediate values between T (a,c) and. Thus, if T (a,c) < 1 W (f + ) < 1 (which can be achieved by taking a or c sufficiently large), by continuity, somewhere between the critical point and the separatrix, there will be a 2π periodic solution. In particular, this can be done by setting c = and taking a sufficiently large. We conclude by noting the well known fact that the cnoidal wave φ is real analytic with exponentially decaying Fourier coefficients. 1.2. Notation. To avoid the use of multiple constants, we write A B to denote that there is an absolute constant C such that A CB. We also write A B to denote both A B and B A. We define the Fourier sequence of a 2π-periodic L 2 function u as u k = 1 2π With this normalization we have u(x) = k u(x)dx, k Z. e ikx u k, and (uv) k = u k v k = m+n=k u n v m.

6 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY 2. Proof of the Main Theorem First we discuss that it suffices to prove the theorem for time-independent cnoidal waves. Consider KdV with periodic boundary conditions (on the circle) q(x+2π) = q(x). Note that if q(x,t) is a solution, then q(x+ct,t)+c is the solution with initial data q(x,)+c. In particular, for a cnoidal wave φ(x ct), the function φ(x)+c is also a cnoidal wave. Applying the statement of the theorem with initial data φ(x) +c+g(x), we obtain at time t, that the solution is of the form q(x+ct,t)+c = φ(x)+c+(e Lct g)(x)+o L 2(εe Ct ), wherel c = x 3 ( 1 2π 2π φ(r)dr+c ) x. Notingthat(e Lct g)(x) = (e Lt g)(x+ ct), we obtain q(x,t) = φ(x ct)+(e Lt g)(x)+o L 2(εe Ct ). Thus, it suffices to prove the theorem for stationary cnoidal wave. Let φ(x) be a time-independent, 2π periodic cnoidal wave. Consider a solution of KdV of the form q(x,t) = φ(x) +u(x,t). Substituting in KdV, we obtain (3) u t +u xxx + φ u x +(Φu) x +uu x =, u(x,) = g(x), where φ = 1 2π φ(x)dx, and Φ = φ φ. Remark 2.1. Our assumption on the H s norm of the initial data g, and the momentum conservation implies that u(x,t) = ĝ() = O(ε). In the proof of our theorem we will restrict ourselves to the case when u(x,t) =. This makes the proof more presentable. Removing this assumption introduces more terms in the differentiation by parts formulas which are smaller then the ones we have. In particular, in Theorem 2.3 below, the formulas for B(v) and R(v) would have additional terms which satisfy the a priori estimates given in Proposition 2.4. Remark 2.2. Note that for a mean-zero L 2 function u, u H s u k / k s l 2, we will use this formula without further comments. For a sequence u k, with u =, we will use u H s notation to denote u k / k s l 2.

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 7 Using the notation u(x,t) = k u k (t)e ikx and Φ(x,t) = k Φ k (t)e ikx, we write (3) on the Fourier side, t u k = ik 2 u k1 u k2 ik Φ k1 u k2 +i(k 3 ak)u k, u k () = ĝ(k), where a = φ. Because of the mean zero assumption on Φ and u, and conservation of momentum, there are no zero harmonics in this equation. Without the mean zero assumption this equation would have an additional term of the form iu kφ k which is of order ε and has fast decay in k. Using the transformations and the identity u k (t) = v k (t)e i(k3 ak)t, Φ k (t) = S k (t)e i(k3 ak)t, ( +k 2 ) 3 k 3 1 k 3 2 a( +k 2 )+a +ak 2 = 3( +k 2 ) k 2, the equation can be written in the form (4) t v k = ik 2 e i3kk 2 t (v k1 +2S k1 )v k2. The following theorem will be proved in Section 3 by distinguishing the resonant and nonresonant sets and using differentiation by parts. Theorem 2.3. The system (4) can be written in the following form: (5) t [v +K(v)+B(v)] k = L (v) k +R(v) k, where we define K(v) = B(v) = L (v) + R(v) =, and for k, we define K(v) k = B(v) k = e 3ikk 2 t S k1 v k2 3 k 2, e 3ikk 2 t v k1 v k2 6 k 2

8 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY 1 9 L (v) k = 2i 3 R(v) k = i v k v k 2 i S k v k v k 2i 3 k 3 k 3 v k + + 2i 3 1 9 + i 9 + i 18 k 2 +k 3 e 3ikk 2 t 3 k 2 ( t S k1 )v k2 +k 2 +k 3 +k 4 =k +k 2 +k 3 +k 4 =k e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) ( +k 2 )(k 2 +k 3 )(k 3 + ) (v +S k1 )v k2 v k3, e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) S k1 S k2 v k3, e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) j k S j v j j v k1 S k2 v k3 e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) v k2 v k3 t S k1 ( +k 2 )(k 2 +k 3 )(k 3 + ) 2i S k 3 k S j v j e it ψ(,k 2,k 3,k 4 ) (k 3 +k 4 )S k1 v k2 (v k3 +2S k3 )v k4 ( +k 2 )( +k 3 +k 4 )(k 2 +k 3 +k 4 ) e it ψ(,k 2,k 3,k 4 ) ( +2k 3 +2k 4 )v k1 v k2 (v k3 +2S k3 )v k4. ( +k 2 )( +k 3 +k 4 )(k 2 +k 3 +k 4 ) Here means the sum does not contain the terms which makes the denominator zero. Proposition 2.4. Assume that v 2 1 and < s < 1/2, then j (6) (7) (8) (9) K(v) H s K(v) 2 v H s, B(v) H s B(v) 2 v 2 H, s L (v) H s v H s, L (v) 2 v 2, R(v) H s R(v) 2 v H s. We will prove this proposition in Section 4. Now, we continue with the proof of the main theorem. First we will prove the near-linear behavior using a modified linear operator (Theorem 2.5), then we will prove that the modified linear evolution is close to the Airy evolution (Theorem 2.6). These two theorems imply Theorem 1.1.

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 9 Theorem 2.5. Let < s < 1/2. Let u be a mean zero solution of (3) with u(,) = g, where Then, for each t >, we have g L 2 1, and g H s ε 1. u(,t) e tl 1 g 2 Ce Ct ε. Here L 1 = L+P, where L = 3 x φ( ) x and P is defined in the Fourier side as (Pu) k = [ e Lt L (e Lt u) ] k = 2i 3 and (Pu) =. Φ k1 Φ k2 u k3, k, Theorem 2.6. Let g, L 1, L be as in the previous theorem. Then, for each t >, we have e tl 1 g e tl g 2 Ce Ct ε. Proof of Theorem 2.5. First we will prove that thenorm assumptions on the initial data remain intact up to times of order log(1/ε). By L 2 conservation in KdV, we have u(,t) L 2 q(,t) L 2 + φ( ) L 2 = q(,) L 2 + φ( ) L 2 1. Now we prove that for some C, (1) u(,t) H s Ce Ct ε. To prove this integrate (5) from to T to obtain (I +K)v(T) = (I +K)v()+B(v)() B(v)(T)+ (L (v)+r(v))dt. Since u k = v k e iψ(k)t and Φ k = S k e iψ(k)t with ψ(k) = k 3 ak, we have (11) (I + K)u(T) = e iψ(k)t (I +K)u()+e iψ(k)t B(u)() e iψ(k)t B(e iψ(k)t u)(t) +e iψ(k)t (L (e iψ(k)t u)+r(e iψ(k)t u)dt,

1 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY where K is the time-independent operator: K(u) k = e 3ik 2 kt iψ( )T iψ(k 2 )T+iψ(k)T Φ k1 u k2 3 k 2 = Φ k1 u k2 3 k 2. Lemma 2.7. For < s 1, (I + K)u H s u H s. Proof. This follows from Fredholm alternative. First note that for u L 2 with mean zero Ku 2 = Φ k1 u k2 3 k 2 2 Φ k 3k l 1 u k k u H s. 2 By the density of L 2 in H s, this inequality holds for each u H s. Therefore, by Rellich s theorem, K is a compact operator on H s. It suffices to show that the kernel of I + K is trivial. Note that if (I + K)u = for some u H s, then by the discussion above, Ku L 2, and hence u L 2. Using the definition of K, we have (I + K)u = u(x) 1 3 Φ 1(x)u 1 (x) =, where f 1 (x) denotes the mean-zero antiderivative of a mean-zero function, (12) f 1 (x) = x f(r)dr + 1 2π rf(r)dr. Let U(x) = u 1 (x), then we obtain equivalent 1st order linear ODE U 1 3 Φ 1(x)U(x) =, which has the general solution U(x) = U()exp( x Φ 1(x)dx), which can be mean-zero only if U(x). This implies u(x). Remark 2.8. Let K t (u) := ρ k1 (t)u k2 3 k 2, where ρ k (t) = ρ(,t)(k), and ρ is an L 2 solution of KdV. Then, the statement of the lemma is valid for K t i) for t [,T] with a constant C T depending on T,

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 11 ii) on R with a constant independent of t if ρ(,t) 2 = ρ(,) 2 < 1/2. To prove these statements first note that (13) K t H s H s 2 ρ(,t) 2. Then, using the resolvent identity, (I + K t ) 1 (I + K s ) 1 = (I + K s ) 1 ( K s K t )(I + K t ) 1, the linearity of Kt in ρ, and (13), we see that the operator (I + K t ) 1 is continuous in time in the operator norm H s H s. Thus, (i) follows form the proof of the lemma and compactness. To see (ii), note that by (13) the operator norm of K in H s is < 1, and invert the operator using Neumann series. Using Lemma 2.7 and Proposition 2.4 in (11), we obtain u(t) H s u() H s + u() 2 H s + u(t) 2 H s + Cε+C u(t) 2 H s +C u(t) H sdt. u(t) H sdt Therefore on [,T ], where T = inf{t : u(t) H s 1 2C }, we have u(t) H s 2Cε+2C u(t) H sdt. This implies by Gronwall that u(t) H s 2Cεe 2CT as claimed. Now, note that K(u) = e iψ(k)t K(e iψ(k)t u). Therefore, using Proposition 2.4 (for K, B, and R) and (1) in (11), we obtain (14) u(t) = e iψ(k)t u()+e iψ(k)t L (e iψ(k)t u)dt+o L 2(εe CT ). Using the definition of L, we have L (e iψ(k)t u) = = 2i 3 = 2i 3 e iψ(k)t e it[ 3(+k 2 )(k 2 +k 3 )(k 3 + ) ψ( ) ψ(k 2 ) ψ(k 3 )] Φ k1 Φ k2 u k3 = e Lt Pu. Φ k1 Φ k2 u k3

12 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY In the last line we used the fact that for k = +k 2 +k 3, ψ(k) = 3( +k 2 )(k 2 +k 3 )(k 3 + ) ψ( ) ψ(k 2 ) ψ(k 3 ), and the definition of L and P. Using this in (14), we have u(t) =e LT u()+e LT e Lt P(u(t))dt+O L 2(εe CT ) =e LT u()+ e L(T t) P(e L 1t u())dt + e L(T t) P[u(t) e L1t u()]dt+o L 2(εe CT ) =e L 1T u()+ e L(T t) P[u(t) e L 1t u()]dt+o L 2(εe CT ). Let h(t) := u(t) e L 1t u() 2. Using the equality above and the bound for L in Proposition 2.4 for the operator P(u(t)) = e Lt L (e Lt u(t)), we obtain h(t) εe CT + The theorem follows from this by Gronwall. h(t)dt. Proof of Theorem 2.6. First we prove that our assumptions on the initial data remain intact for times of order log(1/ε): Lemma 2.9. For s [ 1,1], the operator L 1 defined above satisfies e tl 1 H s H s ect. Proof. First note that we can rewrite P(u) (for mean zero u and Φ) in the following form which is valid for each k Z, P(u) k = 2 3 Φ k1 i Φ k2 u k3 + 2 3 +k 2 +k 3 = Φ k1 i Φ k2 u k3. The constant term makes the right-hand side vanish for k =, which makes P(u) mean-zero in the space side. Using the formula (12) for the function Φ, we write L 1 = L+P in the space side as L 1 u = 3 x u a xu+gu 1 2π u,g,

where Also note that HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 13 G(x) = 2 3 Φ(x) ( x Φ(r)dr + 1 ) rφ(r)dr 2π L 1u = 3 xu+a x u+gu 1 2π G u,1. Note that by duality and interpolation it suffices to prove the assertion of the lemma for s = 1 for L 1 and L 1. We will give the proof for L 1, for L 1 the proof is essentially the same since they have very similar forms. Consider the equation First, we calculate 1 d 2dt u 2 2 = = = u t = L 1 u = 3 x u a xu+gu 1 2π u,g. u t udx u xxx udx a Gu 2 dx 1 2π u,g u x udx+ udx. Gu 2 dx 1 2π u,g udx This implies that d dt u 2 2 2 G u 2 2 + G 2 u 2 2 u 2 2, Similarly, 1 d 2dt u x 2 2 = = u xt u x dx u xxxx u x dx a u xx u x dx+ = Gu 2 x dx 1 u 2 G xx dx 2 ) = O ( u x 2 2 G + u 2 2 G xx = O ( u 2 ) H 1. Combining the two inequalities, we obtain d dt u 2 H u 2 1 H 1, (Gu) x u x dx which finishes the proof.

14 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY We return to the proof of Theorem 2.6. Consider the equation u t = L 1 u, u(,x) = g(x). Repeating the discussion in the beginning of Section 2, and introducing the variables v k and S k as above, we have (15) t v k = 2i e 3it(+k 2 )(k 2 +k 3 )( +k 3 ) S S k2 v k3. 3 We will prove the following proposition, using differentiation by parts, in Section 3. Proposition 2.1. The system (15) is equivalent to the following (16) t ( vk +D(v) k ) = E(v)k, where D(v) k = i 3 E(v) k = 2i 3 S S j v j k j j k i 3 + 2 9 +k 2 +k 3 +k 4 +k 5 =k e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) ( +k 2 )(k 2 +k 3 )(k 3 + ) S S k2 v k3, + 2i S k 3 k S j v j j e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) t (S k1 S k2 )v k3 ( +k 2 )(k 2 +k 3 )(k 3 + ) e it ψ(,k 2,k 3,k 4,k 5 ) S k1 S k2 S k3 S k4 v k5 k 3 ( +k 2 )(k )(k k 2 ). Here ψ is a real valued phase function which is irrelevant for the proof of the Theorem. Proposition 2.11. The following a priori estimates hold D(v) 2 v H 1, E(v) 2 v H 1. Proof. Using +k 3 k 3 (for nonzerointeger values of and +k 3 ), and eliminating the product ( +k 2 )(k 2 +k 3 ) from the denominator, we have D(v) k S k1 S k2 v k 3 k 3.

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 15 This implies by Young s inequality D(v) 2 S l 1 S l 1 v k /k 2 v H 1. The proof for the contribution of the second line in the definition of E(v) is exactly the same using the fast decay of t S k k 3 S k. The L 2 norm of the first line is S 2 j S j v j j + S k /k 2 j js j v j j [ S 2 2 + S k/k 2 ks k 2 ] v k /k 2 v H 1. The second inequality follows from Cauchy Schwarz. To estimate the third line in the definition of E(v), note that k k 2 k 3 = +k 3 +k 4 +k 5 k 3 k 4 k 4 k 5 k 4 for nonzero integer values of the factors. Using this in the sum and eliminating the rest of the factors, we estimate the third line as S k1 S k2 S k3 k 4 S k4 v k 5 k 5. +k 2 +k 3 +k 4 +k 5 =k As above, this implies that the L 2 norm is S 3 l 1 ks k l 1 v k /k 2 v H 1. To complete the proof of Theorem 2.6, integrate (16) from to T: v k (T) = v k () D(v) k (T)+D(v) k ()+ Using the transformation u k = v k e iψ(k)t, we have E(v) k (t)dt. u k (T) = e iψ(k)t u k () e iψ(k)t D(e iψ(k)t u) k (T)+e iψ(k)t D(u) k () + e iψ(k)t E(e iψ(k)t u) k (t)dt. Noting that (e LT u) k = e iψ(k)t u k (), and using the a priori estimates in Proposition 2.11, we have u(t) e LT u() 2 u(t) H 1 + u() H 1 + εe CT. u(t) H 1dt

16 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY In the last line we used Lemma 2.9 and the hypothesis on u() H s u() H 1. 3. Differentiation by Parts In this section we prove Theorem 2.3 and Proposition 2.1. Proof of Theorem 2.3. We need to obtain (5) from (4) by using differentiation by parts. It will be useful to name the terms appearing in the formula (5). We will denote the terms in the first and second line of the definition of B(v) by B 1 (v) and B 2 (v), respectively. We also denote the term in the jth line of the definition of R(v) in the statement of the Theorem by R j (v), j = 1,2,...,6, and we further denote the four summands in R 1 (v) by R 1,m (v), m = 1,2,3,4. Since e 3ikk 2 t = t ( i 3k k 2 e i3kk 2 t ), using differentiation by parts we can rewrite (4) as t v k = t ( 1 2 k 1 2 k e 3ikk 2 t (v k1 +2S k1 )v ) k2 3k k 2 e 3ikk 2 t 3k k 2 t [(v k1 +2S k1 )v k2 ]. Recalling the definition of K(v) and B 1 (v) from equation (5), we can rewrite this equation in the form: t [v k +K(v) k +B 1 (v) k ] = e 3ikk 2 t 6 k 2 t [(v k1 +2S k1 )v k2 ]. Note that since v = and S =, in the sums above and k 2 are not zero. We now handle the term when the derivative hits v k1 v k2. By symmetry and (4), we have = i = i e 3ikk 2 t k 2 t (v k1 v k2 ) = e 3ikk1k2t k= +k 2 k= +µ+λ v k1 ( µ+λ=k 2 e 3itk2µλ (v µ +2S µ )v λ v k1 (v µ +2S µ )v λ e 3it[k(µ+λ)+µλ(µ+λ)]. )

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 17 We note that µ+λ can not be zero since µ+λ = k 2. Using the identity k +µλ = ( +µ+λ) +µλ = ( +µ)( +λ) and by renaming the variables k 2 = µ,k 3 = λ, we have that = i e 3ikk 2 t k 2 +k 3 k 2 t (v k1 v k2 ) = e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) v k1 (v k2 +2S k2 )v k3. Calculating the term when the derivative hits S k1 v k2 similarly, we have (17) t [v k +K(v) k +B 1 (v) k ] = where Note that Y 1 (v) k = i 3 Y 2 (v) k = i 3 Y 3 (v) k = 2i 3 Y 4 (v) k = 2i 3 k 2 +k 3 k 2 +k 3 k 2 +k 3 k 2 +k 3 Y 5 (v) k = 5 Y j (v) k, j=1 e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) v k1 v k2 v k3, e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) S k1 v k2 v k3, e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) v k1 S k2 v k3, e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) S k1 S k2 v k3, e 3ikk 2 t 3 k 2 ( t S k1 )v k2. (18) Y 3 (v) = R 3 (v), Y 4 (v) = L (v)+r 1,4, and Y 5 (v) = R 2 (v). Due to the fast decay of S k and t S k, the terms R 2 (v), R 3 (v), and R 1,4 are small (as stated in Proposition 2.4. The term L (v) is not small but linear in v, and it is handled separately, see the proof of Theorem 2.5. On the other hand one can not directly estimate the terms Y 1 and Y 2 without performing

18 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY another differentiation by parts. To do that we need to check the resonant terms in Y 1 and Y 2 : (19) ( +k 2 )(k 3 + ) =. The set for which (19) holds is the disjoint union of the following 3 sets (taking into account that k 2 +k 3 ) A 1 = { +k 2 = } {k 3 + = } { = k, k 2 = k, k 3 = k}, A 2 = { +k 2 = } {k 3 + } { = j, k 2 = j, k 3 = k, j k }, A 3 = {k 3 + = } { +k 2 }} { = j, k 2 = k, k 3 = j, j k }. We write Y 2 (v) k = Y 2r (v) k +Y 2nr (v) k where the subscript r and nr stands for the resonant and non-resonant terms respectively. We have Y 2r (v) k = i 3 3 and λ=1 S v v k1 k2 k3 = i 3 A λ = R 1,2 (v) k +R 1,3 (v) k, Y 2nr (v) k = i 3 nr S k v k v k k 2i 3 v k j Z j k S j v j j e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) S k1 v k2 v k3. Since the exponent in Y 2nr (v) is not zero we can differentiate by parts one more time and obtain that where (2) M 3 (v) k = 1 9 and M 4 (v) k = 1 9 Y 2nr (v) k = t M 3 (v) k +M 4 (v) k nr nr e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) ( +k 2 )(k 2 +k 3 )(k 3 + ) S v k2 v k3, e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) ( +k 2 )(k 2 +k 3 )(k 3 + ) (v k2 v k3 t S k1 +S k1 v k3 t v k2 +S k1 v k2 t v k3 )

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 19 A calculation as before, by expressing time derivatives using (4), reveals that M 4 (v) k = 1 9 + i 9 +k 2 +k 3 +k 4 =k = R 4 (v) k +R 5 (v) k. e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) v k2 v k3 t S k1, ( +k 2 )(k 2 +k 3 )(k 3 + ) e it ψ(,k 2,k 3,k 4 ) (k 3 +k 4 )S k1 v k2 (v k3 +2S k3 )v k4 ( +k 2 )( +k 3 +k 4 )(k 2 +k 3 +k 4 ) The phase function ψ is irrelevant for our calculations since it is going to be estimated out by taking absolute values inside the sums. For completeness we note that it can be expressed as Hence ( +k 2 +k 3 +k 4 ) 3 k 3 1 k 3 2 k 3 3 k 3 4. (21) Y 2 (v) = R 1,2 (v) k +R 1,3 (v) k + t M 3 (v)+r 4 (v) k +R 5 (v) k. Similarly, (22) Y 1 (v) = R 1,1 + t N 3 (v)+r 6 (v), where (23) N 3 (v) k := 1 9 e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) v k1 v k2 v k3. ( +k 2 )(k 2 +k 3 )(k 3 + ) Using (2) and (23), note that B 2 (v) = M 3 (v) N 3 (v). Therefore, substituting (18), (21), and (22) in (17), we obtain (5). Proof of Proposition 2.1. We write the right hand side of (15) by distinguishing the resonant and nonresonant terms. The resonant set corresponding to the terms with k 2 +k 3 is the same as above, and thus we get the following 3 terms: 2i 3 ( S k k S kv k +v k j, j k S j S j j +S k j, j k S j v ) j. j Note that, by symmetry, the second term is zero. Combining the other terms we obtain the first summand in the definition of E(v). The terms with k 2 +k 3 = gives the second summand.

2 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY For the nonresonant terms we differentiate by parts as above obtaining i ( 3 e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) ) t ( +k 2 )(k 2 +k 3 )(k 3 + ) S S k2 v k3 i 3 e 3it(+k 2 )(k 2 +k 3 )(k 3 + ) ( +k 2 )(k 2 +k 3 )(k 3 + ) t(s k1 S k2 v k3 ). The first line gives D(v), and the second line gives the remaining terms in the definition of E(v) after using the formula for t v k3 and renaming the variables. We start with the term K: 4. Proof of Proposition 2.4 K(v) 2 S k /k l 1 v k /k 2 v H s. Similarly the L 2 norm of the first summand in the definition of B(v) is v k /k l 1 v k /k 2 v k / k s 2 1/ k 1 s l 2 v k /k 2 v 2 H s, since < s < 1/2. Note that the second summand is v k1 + S k1 +k 2 k 2 +k 3 k 3 + v k 2 v k3 v k1 v k2 v k3 k 2 + S k1 v k2 v k3. k 2 k 3 In the second line we use the following inequalities which are valid for the nonzero integral values of the factors: +k 2 k 2 +k 3 k 3 + k 2, +n n. Therefore by Young s inequality, the L 2 norm of the second summand is v k /k 2 l 1 v k 2 + ks k l 1 v k /k l 1 v k /k 2 v 2 H s. The L 2 bound for L (v) follows from Young s inequality: L (v) 2 S k /k l 1 S l 1 v 2 v 2.

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 21 For the H s bound, using the inequality k 3 s s k 2 s +k 2 +k 3 s = k 2 k s (for nonzero integral values of the factors), we obtain Therefore, L (v) k k s S k1 1 s k 2 s S k2 v k 3 k 3 s. L (v) H s S k / k 1 s l 1 k s S k l 1 v k / k s 2 v H s. WenowestimateR(v). Denotethetermsinthejthlineofthedefinitionof R(v) by R j (v), j = 1,2,...,6. The estimate for R 1 (v) follows as in the proof of Proposition 2.11. The estimate for R 2 (v) is the same as the estimate for K(v) with S replaced with t S. For R 3 (v), note that by Young s inequality R 3 (v) 2 v k /k l 1 S 1 v 2 v H s. For R 4 (v), using +k 2 k 2, we have R 4 (v) 2 t S l 1 v k /k l 1 v 2 v H s. For R 5 (v), using + k 2 k 2, k 3 + k 4 + k 3 + k 4 and k 2 +k 3 +k 4 +k 2 +k 3 +k 4 = k, we have R 5 (v) k 1 k +k 2 +k 3 +k 4 =k Therefore, by Young s inequality R 5 (v) 2 1/k l 2 +k 2 +k 3 +k 4 =k k 2 1 S k1 v k 2 k 2 v k 3 +2S k3 v k4. k1 2 S v k 2 k 2 v k 3 +2S k3 v k4 k 2 S k l 1 v k /k l 1 v +2S 2 v 2 v H s. Finally, we consider R 6 (v). Using +2k 3 +2k 4 +2 +k 3 +k 4, and then Cauchy Schwarz, we have [ R 6 (v) k 2 +k 2 +k 3 +k 4 =k +k 2 +k 3 +k 4 =k l v k1 v k2 (v k3 +2S k3 )v k4 ( 1 +k 2 k 2 +k 3 +k 4 + 1 )] 2 +k 3 +k 4 v 2 v 2 k 2 (v k3 +2S k3 ) 2 2s

22 M. B. ERDOĞAN, N. TZIRAKIS, AND V. ZHARNITSKY +k 2 +k 3 +k 4 =k 2s v 2 k 4 +k 2 2 k 2 +k 3 +k 4 2 ( 1 2 + 1 +k 3 +k 4 2 ) Note that the first factor above is v 2 H v + 2S 2 s 2 v 2 2. Therefore it suffices to prove that the sum of the second factor in k is O(1). We write this sum as,k 2,k 3,k 4 +,k 2,k 3,k 4 v 2 k 4 2 2s +k 2 2 k 2 +k 3 +k 4 2 2s v 2 k 4 +k 2 2 +k 3 +k 4 2 k 2 +k 3 +k 4 2. To estimate the first sum, first take the sum in k 3, then in k 2,, and k 4 in the given order. The estimate for the second sum follows also by summing in the order given above and using the inequality 1 +k 2 2 +k 3 +k 4 2 k 2 +k 3 +k 4 2 1 +k 2 2 k 2 2 1 2. k 2,k 3 k 2 Here we used the inequality (for a b > 1) m 1 n 1 m a n 2 m b 1 n 1 n 2 b. Remark 4.1. We note that in the proof of Proposition 2.4 and Proposition 2.11, the strongest conditions we need for S k = e iψ(k)t Φk are t S k l 1, and k 2 S k l 1. It is easy to see, using the equation, that t ρ H 1 for an H 4 solution ρ of KdV. Therefore, the sequence e iψ(k)t ρ k satisfies the conditions above, and thus the assertions of Proposition 2.4 and Proposition 2.11 remains valid when Φ is replaced with ρ(,t) H 4. Acknowledgement. The authors would like to thank Michael Weinstein for a helpful discussion.

HIGH FREQUENCY PERTURBATION OF CNOIDAL WAVES 23 References [1] A. V. Babin, A. A. Ilyin, E. S. Titi, On the regularization mechanism for the periodic Kortewegde Vries equation, Commun. Pure Appl. Math. 64 (211), 591 648. [2] T. B. Benjamin, Lectures on nonlinear wave motion, in Nonlinear Wave Motion, AMS (1974), 3 47. [3] N. Bottman, B. Deconinck, KdV cnoidal waves are spectrally stable, Discrete and contin. dyn. syst. 25 (24), 1163 118. [4] J. Bourgain,Fourier transform restriction phenomena for certain lattice subsets and applications to nonlinear evolution equations. Part I: Schrödinger equations, Geom. Funct. Anal. 3 (1993), 17 156. [5] J. Bourgain, Fourier transform restriction phenomena for certain lattice subsets and applications to nonlinear evolution equations. Part II: The KdV equation, Geom. Funct. Anal. 3 (1993), 29 262. [6] J. Colliander, M. Keel, G. Staffilani, H. Takaoka, T. Tao, Sharp Global Well- Posedness for KdV and Modified KdV on R and T, J. Amer. Math. Soc. 16 (23), 75 749. [7] J. Colliander, M. Keel, G. Staffilani, H. Takaoka, T. Tao, Transfer of energy to high frequencies in the cubic defocusing nonlinear Schrdinger equation, Invent. Math. 181 (21), 39 113. [8] M. B. Erdoğan, N. Tzirakis, V. Zharnitsky, Near-linear dynamics in KdV with periodic boundary conditions, Nonlinearity 23 (21), 1675 1694. [9] M. B. Erdoğan, N. Tzirakis, V. Zharnitsky, Near-linear dynamics of nonlinear dispersive waves, preprint 21. [1] T. Kappeler and P. Topalov, Global wellposedness of KdV in H 1 (T,R), DukeMath. J. 135 (26), 327-36. [11] C. E. Kenig, G. Ponce, L. Vega, A bilinear estimate with applications to the KdV equation,, J. Amer. Math. Soc. 9 (1996), 573 63. [12] H. P. McKean, Stability for the Korteweg-deVries equation, Commun. Pure Appl. Math. 3 (1977), 347 353. [13] J.A. Pava, J.L. Bona, M. Scialom, Stability of cnoidal waves, Adv. in Diff. Equations 11 (26), 1321-1374. [14] T. Tao, Nonlinear dispersive equations (local and global analysis), CBMS 16, 26. Department of Mathematics, University of Illinois, Urbana, IL 6181, U.S.A. E-mail address: berdogan@uiuc.edu, tzirakis@math.uiuc.edu, vzh@uiuc.edu