arxiv: v3 [math.ap] 20 Jun 2017

Transcription

1 A modulation technique for the blow-up profile of the vector-valued semilinear wave equation arxiv: v3 [math.ap] 2 Jun 217 Asma AZAIEZ Université de Cergy-Pontoise, AGM, CNRS UMR 888, 9532, Cergy-Pontoise, France. Hatem ZAAG Université Paris 13, Sorbonne Paris Cité LAGA, CNRS UMR 7539, F-9343, Villetaneuse, France. June 21, 217 Abstract We consider a vector-valued blow-up solution with values in R m for the semilinear wave equation with power nonlinearity in one space dimension this is a system of PDEs. We first characterize all the solutions of the associated stationary problem as an m-parameter family. Then, we show that the solution in self-similar variables approaches some particular stationary one in the energy norm, in the non-characteristic cases. Our analysis is not just a simple adaptation of the already handled real or complex case. In particular, there is a new structure of the set a stationary solutions. Keywords: Wave equation, blow-up profile, stationary solution, modulation theory, vector valued PDE. AMS classification: 35L5, 35L81, 35B44, 35B4, 34K21. 1 Introduction We consider the vector-valued semilinear wave equation { 2 t u = 2 x u+ u p u, u = u and u t = u 1, Paper accepted for publication in Bull. Sci. math. 1 where here and all over the paper. is the euclidian norm in R m, ut : x R ux,t R m, m 2, p > 1, u Hloc,u 1 and u 1 L 2 loc,u with v 2 = sup vx 2 dx and L 2 loc,u a R x a <1 v 2 = v 2 + v 2 Hloc,u 1 L 2 loc,u Lloc,u 2 The Cauchy problem for equation 1 in the space Hloc,u 1 L2 loc,u follows from the finite speed of propagation and the wellposedness in H 1 L 2. See for instance Ginibre, Soffer and Velo [9], Ginibre and Velo [1], Lindblad and Sogge [14] for the local in time wellposedness in H 1 L 2. Existence of blow-up solutions follows from ODE techniques or the energy-based This author is supported by the ERC Advanced Grant no , BLOWDISOL. This author is supported by the ERC Advanced Grant no , BLOWDISOL. and by ANR project ANAÉ ref. ANR-13-BS

2 blow-up criterion of [13]. More blow-up results can be found in Caffarelli and Friedman [6], Alinhac [1] and [2], Kichenassamy and Littman [12], [11] Shatah and Struwe [25]. The real case in one space dimension has been understood completely, in a series of papers by Merle and Zaag [18], [19], [21] and [22] and in Côte and Zaag [7] see also the note [2]. Recently, the authors give an extension to higher dimensions in [24] and [23], where the blow-up behavior is given, together with some stability results. For other types of nonlinearities, we mention the recent contribution of Azaiez, Masmoudi and Zaag in [5], where we study the semilinear wave equation with exponential nonlinearity, in particular we give the blow-up rate with some estimations. In [4], we consider the complex-valued solution of 1 or R 2 -valued solution, characterize all stationary solutions and give a trapping result. The main obstruction in extending those results to the vector case m 3 was the question of classification of all self similar solutions of 1 in the energy space. In this paper we solve that problem and show that the real valued and complex valued classification also hold in the vector-valued case m 3 see Proposition 2 below, with an adequate choice in S m. This is in fact our main contribution in this paper, and it allows us to generalize the results of the complex case to the vector valued case m 3. In this paper, we aim at proving similar results for the general case ux,t R m, for m 3. Let us first introduce some notations before stating our results. If u is a blow-up solution of 1, we define see for example Alinhac [1] a continuous curve Γ as the graph of a function x Tx such that the domain of definition of u or the maximal influence domain of u is D u = {x,t t < Tx}. From the finite speed of propagation, T is a 1-Lipschitz function. The time T = inf x R Tx and the graph Γ are called respectively the blow-up time and the blow-up graph of u. Let us introduce the following non-degeneracy condition for Γ. If we introduce for all x R, t Tx and δ >, the cone C x,t,δ = {ξ,τ x,t τ t δ ξ x }, then our non-degeneracy condition is the following: x is a non-characteristic point if δ = δx,1 such that u is defined on C x,tx,δ. 2 If condition 2 is not true, then we call x a characteristic point. Already when u is real-valued, we know from [21] and [7] that there exist blow-up solutions with characteristic points. Given some x R, we introduce the following self-similar change of variables: w x y,s = Tx t pux,t, 2 y = x x Tx t, s = logtx t. 3 This change of variables transforms the backward light cone with vertex x,tx into the infinitecylinder y,s,1 [ logtx,+. Thefunctionw x wewritew forsimplicity satisfies the following equation for all y < 1 and s logtx : s 2 2p+1 w = Lw p 2w + w p w p+3 p sw 2y ys w 4 where Lw = 1 ρ yρ1 y 2 y w and ρy = 1 y 2 2 p. 5 2

3 This equation will be studied in the space H = {q Hloc 1 L2 loc,1,rm q 2 H q q y 2 + q 2 2 ρ dy < + }, 6 which is the energy space for w. Note that H = H L 2 ρ where H = {r Hloc 1,1,Rm r 2 H r 2 1 y 2 + r 2 ρ dy < + }. 7 In some places in our proof and when this is natural, the notation H, H and L 2 ρ may stand for real-valued spaces. Let us define 1 Ew, s w = 2 sw yw 2 1 y 2 + p+1 p 2 w 2 1 p+1 w p+1 ρdy. 8 By the argument of Antonini and Merle [3], which works straightforwardly in the vector-valued case, we see that E is a Lyapunov functional for equation Blow-up rate Only in this subsection, the space dimension will be extended to any N 1. We assume in addition that p is conformal or sub-conformal: 1 < p p c 1+ 4 N. We recall that for the real case of equation 1, Merle and Zaag determined in [15] and [16] the blow-up rate for 1 in the region {x,t t < T} in a first step. Then in [17], they extended their result to the whole domain of definition {x,t t < Tx}. In fact, the proof of [15], [16] and [17] is valid for vector-valued solutions, since the energy structure see8, which is the main ingredient of the proof, is preserved. This is the growth estimate near the blow-up surface for solutions of equation 1. Proposition 1. Growth estimate near the blow-up surface for solutions of equation 1 If u is a solution of 1 with blow-up surface Γ : {x Tx}, and if x R N is noncharacteristic in the sense 2 then, i Uniform bounds on w For all s log Tx 4 : w x s H 1 B + s w x s L 2 B K. ii Uniform bounds on u For all t [ 3 4 Tx,Tx : Tx t 2 ut p L 2 Bx,Tx t Tx t +Tx t 2 p +1 t ut L 2 Bx,Tx t Tx t N/2 + ut L2 Bx,Tx t Tx t N/2 K, where the constant K depends only on N, p, and on an upper bound on Tx,1/Tx, δ x and the initial data in H 1 loc,u L2 loc,u. 3

4 1.2 Blow-up profile This result is our main novelty. In the following, we characterize the set of stationary solutions for vector-valued solutions. Proposition 2. Characterization of all stationary solutions of equation 4 in H. i Consider w H a stationary solution of 4. Then, either w or there exist d,1 and Ω S m such that wy = Ωκd,y where d,y,1 2, κd,y = κ 1 d 2 1 p 1+dy 2 p and κ = 1 2p+1 p. 9 p 2 ii It holds that where E is given by 8. E, = and d,1, Ω S m, Eκd,.Ω, = Eκ, > 1 Thanks to the existence of the Lyapunov functional Ew, s w defined in 8, we show that when x is non-characteristic, then w x approaches the set of non-zero stationary solutions: Proposition 3. Approaching the set of non-zero stationary solutions near a noncharacteristic point Consider u a solution of 1 with blow-up curve Γ : {x Tx}. If x R is non-characteristic, then: A.i inf {Ω S m, d <1} w x.,s κd,.ω H 1,1 + s w x L 2,1 as s. A.ii Ew x s, s w x s Eκ, as s. We write the fundamental theorem of our paper: Theorem 4. Trapping near the set of non-zero stationary solutions of 4 There exist positive ǫ, µ and C such that if w C[s,,H for some s R is a solution of equation 4 such that s s,ews, s ws Eκ,, 11 and ws κd s ws,.ω H ǫ 12 for some d,1,ω S m and ǫ,ǫ ], then there exists d,1 and Ω S m such that argtanhd argtanhd + Ω Ω C ǫ and for all s s : ws κd,.ω s ws C ǫ e µ s s. 13 H Combining Proposition 3 and Theorem 4, we derive the existence of a blow-up profile near non-characteristic points in the following: 4

5 Theorem 5. Blow-up profile near a non-characteristic point If u a solution of 1 with blow-up curve Γ : {x Tx} and x R is non-characteristic in the sense 2, then there exist d x,1, Ω x S m and s x logtx such that for all s s x, 13 holds with ǫ = ǫ, where C and ǫ are given in Theorem 4. Moreover, w x s κd x Ω x H 1,1 + s w x s L 2,1 as s. Remark: From the Sobolev embedding, we know that the convergence takes place also in L, in the sense that w x s κd x Ω x L,1 as s. In this paper, we give the proofs of Proposition 2 and Theorem 4, which present the novelties of this work comparing with the handled real and complex cases, since Propositions 1, 3 and Theorem 5, can be generalized from the real case treated in [18] without any difficulty. Let us remark that our paper is not a simple adaptation of the complex case. In fact, the vector-valued structure of our solution implies a new characterization of the set of stationary solutions in R m see Proposition 2 above. In addition, in order to apply the modulation theory, we need more parameters, and for that, a suitable m m rotation matrix will be defined see 6 and 61 below; see the beginning of the proof of Proposition 3.3 page 19 below, and we have to treat delicately the terms coming from the rotation matrix. This paper is organized as follows: - In Section 2, we give the proof of Proposition 2. - In Section 3, we give the proof of Theorem 4. 2 Characterization of the set of stationary solutions In this section, we prove Proposition 2 which characterizes all H solutions of 1 ρ ρ1 y2 w 2p+1 p 2w + w p w =, 14 the stationary version of 4. Note that since and κ Ω are trivial solutions to equation 4 for any Ω S m, we see from a Lorentz transformation see Lemma 2.6 page 54 in [18] that T d e iθ κ = κd,y is also a stationary solution to 4. Let us introduce the set S {,κd,.ω, d < 1,Ω S m }. 15 Now, we prove Proposition 2 which states that there are no more solutions of 14 in H outside the set S. We first prove ii, since its proof is short. ii Since we clearly have from the definition 8 that E, =, we will compute EΩκd,.,. From 8 and the proof of the real case treated in page 59 in [18], we see that Eκd,.Ω, = Eκd,., = Eκ, >. Thus, 1 follows. iconsiderw H anr m non-zerosolutionof14. Letusprovethattherearesomed,1 and Ω S m such that w = κd,.ω. For this purpose, define ξ = 1 2 log 1+y 1 y that is y = tanhξ and wξ = wy1 y 2 1 p. 16 5

6 As in the real case, we see from straightforward calculations that w is a H 1 R solution to 2 ξ w + w p w 4 w =, ξ R. 17 p 2 Our aim is to prove the existence of Ω S m and ξ R such that wξ = Ω kξ+ξ where kξ = κ. coshpξ 2 Since w H 1 R C 1 2R, we see that w is a strong C 2 solution of equation 17. Since w, there exists ξ R such that wξ. By invariance of 17 by translation, we may suppose that ξ =. Let G = {ξ R wξ }, 18 a nonempty open set by continuity. Note that G contains some non empty interval I containing. We introduce ρ and Ω by From equation 17, we see that ρ = w, Ω = w w, whenever ξ G. ρ Ω+2ρ Ω +ρω +ρ p Ω 4 p 2ρΩ =. 19 Now, since Ω = 1, we immediately see that Ω.Ω = and Ω.Ω+ Ω 2 =. Let Hξ = Ω 2. Projecting equation 19 according to Ω and Ω we see that ξ G, { ρ ξ ρξhξ c ρξ+ρξ p =, c = 4 p 2 4ρ ξhξ+ρξh ξ = 2 Integrating the second equation on the interval I G, we see that for all ξ I, Hξ = Hρ 4. Plugging this in the first equation, we get ρξ 4 Now let ξ I, ρ ξ µ ρξ 3 c ρξ+ρ p ξ = where µ = Hρ G = } {ξ G, ξ I ξ, Hξ = Hρ4 ρξ 4, 22 where I ξ = [,ξ if ξ or I ξ = ξ,] if ξ. Note that I G. Now, we give the following: Lemma 2.1. There exists ǫ > such that ξ G, ξ I ξ, < ǫ wξ 1 ǫ. Proof. The proof is the same as in the complex-case, see page 5898 in [4]. But for the reader s convenience and for the sake of self-containedness, we recall it here. Take ξ G. By definition 22 of G, we see that equation 21 is satisfied for all ξ I ξ. Multiplying ρ ξ µ ρξ 3 c ρξ+ρ p ξ = by ρ and integrating between and ξ, we get: ξ I ξ, Eξ = E, where Eξ = 1 2 ρ ξ µ 2ρξ 2 c 2 ρ2 ξ + ρp+1 ξ p+1,

7 or equivalently, ξ I ξ, Fρξ = 1 2 ρ ξ 2 where Fr = µ 2r 2 + c 2 r2 rp+1 p+1 +E. Since Fr as r or r, there exists ǫ = ǫ µ,e > such that ǫ ρξ 1 ǫ, which yields to the conclusion of the Claim 2.1. We claim the following: Lemma 2.2. It holds that G = R. Proof. Note first that by construction, G is a nonempty interval note that I G where I is defined right before 18. We have only to prove that supg = +, since the fact that infg = can be deduced by replacing wξ by w ξ. By contradiction, suppose that supg = a < +. First of all, by Lemma 2.1, we have for all ξ [,a, < ǫ wξ 1 ǫ. By continuity, this holds also for ξ = a, hence, wa, and a G. Furthermore, by definition of G and continuity, we see that ξ [,a],hξ = Hρ4 ρξ Therefore, we see that a G. By continuity, we can write for all ξ a δ,a +δ, where δ > is small enough, { ρ ξ ρξhξ c ρξ+ρξ p =, c = 4 p 2 4ρ ξhξ+ρξh ξ =. From the second equation and 23 applied with ξ = a, we see that Hξ = Haρa4 = ρξ 4 Hρ 4. Therefore, itfollowsthata,a+δ G, whichcontradictsthefactthata = supg. ρξ 4 Note from Lemma 2.2 that 2 and 21 holds for all ξ R. We claim that H =. Indeed, if not, then by 21, we have µ, and since G = R, we see from Lemma 2.1 that for all ξ R, wξ ǫ, therefore w / L 2 R, which contradicts the fact that w H 1 R. Thus, H =, and µ =. By uniqueness of solutions to the second equation of 2, we see that Hξ = for all ξ R, so Ωξ = Ω, and { w = ρω w = ρ Ω. Let W be the maximal real-valued solution of W c W + W p W = W = ρ W = ρ. By uniqueness of the Cauchy problem of equation 17, we have for all ξ R, wξ = WξΩ, and as w H 1 R, W is also in H 1 R. It is then classical that there exists ξ such that for all ξ R, Wξ = kξ+ξ remember that ρ >, hence we only select positive solutions here. In addition, for Ω = Ω, wξ = kξ+ξ Ω. Thus, for d = tanhξ,1 and y = tanhξ, we get [1 ] 1 [ wξ = κ 1 tanhξ +ξ 2] 1 p tanhξ 2 p +tanhξ Ω = κ Ω 1+tanhξtanhξ ] 1 [ y +d 2 p = κ [1 ] 1 d 2 1 y p Ω = κ Ω = κd,y1 y 2 pω 1. 1+dy 1+dy 7

8 By 16, we see that wy = κd,yω. This concludes the proof of Proposition 2. 3 Outline of the proof of Theorem 4 The proof of Theorem 4 is not a simple adaptation of the complex-case to the vector-valued case, in fact, it involves a delicate modulation. In this section, we will outline the proof, insisting on the novelties, and only recalling the features which are the same as in the real-valued complexvalued cases. This section is organized as follows: - In Subsection 3.1, we linearize equation 4 around κd,ye 1 where e 1 = 1,,..., and figureout that, with respect to the complex-valued case, our linear operator is just a superposition of one copy of the real part operator, with m copies of the imaginary part operator. - In Subsection 3.2, we recall from [18] the spectral properties of the real-part operator. - In Subsection 3.3, we recall from [4] the spectral properties of the imaginary-part operator. - In Subsection 3.4, assuming that Ω = e 1 possible thanks to rotation invariance of 4, we introduce a modulation technique adapted to the vector-valued case. This part makes the originality of our work with respect to the complex-valued case. - In Subsection 3.5, we write down the equations satsified by the modulation parameters along with the PDE satisfied by qy,s and its components. - In Subsection 3.6, we conclude the proof of Theorem The linearized operator around a non-zero stationary solution We study the properties of the linearized operator of equation 4 around the stationary solution κd,y 9. Let us introduce q = q 1,q 2 R m R m for all s [s,, for a given s R, by wy,s κd,ye1 q1 y,s = s wy,s q 2 y,s Let us introduce the coordinates of q 1 and q 2 by q 1 = q 1,1,q 1,2,...,q 1,m, q 2 = q 2,1,q 2,2,...,q 2,m. We see from equation 4, that q satisfies the following equation for all s s : s q1 q 2 q1 = L d +, 25 q 2 f d q 1 where q1 q 2 L d = q 2 Lq 1 + ψd,yq 1,1 e 1 + m ψd,yq j=2 1,j e j p+3 p q, 2 2y y q 2 ψd,y = pκd,y p 2p+1 p 2 26 ψd,y = κd,y p 2p+1 p 2 27 f d q 1 = f d,1 q 1 e 1 + f d,j q 1 e j, j=2 8

9 where f d,1 q 1 = κd,ye 1 +q 1 p κd,y+q 1,1 κd,y p pκ p d,yq 1,1. 28 f d,j q 1 = κd,ye 1 +q 1 p q 1,j κ p d,yq 1,j. 29 Projecting 25 on the first coordinate, we get for all s s : where L d is given by: q1,1 L d q 2,1 s = q1,1 q 2,1 = L q1,1 d +, 3 q 2,1 f d,1 q 1 q 2,1 Lq 1,1 + ψd,yq 1,1 p+3 p q 2,1 2y y q 2,1, 31 Now, projecting equation 25 on the j-th coordinate with j = 2,..,m, we see that q1,j = s q L q1,j d +, 32 2,j q 2,j f d,j q 1 where q1,j L d = q 2,j q 2,j Lq 1,j + ψd,yq 1,j p+3 p q 2,j 2y y q 2,j, 33 Remark: Our linearized operator L d is in fact diagonal in the sens that q1 L d = q L q1,1 d e 2 q 1 + 2,1 j=2 q1,j L d e q j. 2,j We mention that for j = 1, equation 3 is the same as the equation satisfied by the real part of the solution in the complex case see Section 3 page 5899 in [4], whereas for j = 2,..,m, equation 32 is the same as the equation satisfied by the imaginary part of the solution operator in the complex case. Thus, the reader will have no difficulty in adapting the remaining part of the proof to the vector-valued case. Thus, the dynamical system formulation we performed when m = 2 can be adapted straightforwardly to the case m 3. Note from 6 that we have q H = [φq,q] 1 2 < +, where the inner product φ is defined by q1 r1 φq,r = φ, = q 2 r 2 q 1.r 1 +q 1.r 1 1 y2 +q 2.r 2 ρ dy. where q 1.r 1 = m j=1 q 1,j.r 1,j is the standard inner product in R m, with similar expressions for q 1.r 1 and q 2.r 2. Using integration by parts and the definition of L 5, we have the following: φq,r = q 1 Lr 1 +r 1 +q 2 r 2 ρdy. 34 In the following two sections, we recall from [18] and [4] the spectral properties of L d and L d. 9

10 3.2 Spectral theory of the operator L d From Section 4 in [18], we know that L d has two nonnegative eigenvalues λ = 1 and λ = with eigenfunctions F d 1 y = 1 d2 p p 1+dy p+1 p 1+dy p+1 p Note that for some C > and any λ {,1}, we have and F d y = 1 d2 1 p y+d 1+dy p+1 p. 35 d < 1, 1 F λ d C H C and d Fd λ H C 1 d2. 36 Also, we know that L d the conjugate operator of L d with respect to φ is given by L d r1 R d r 2 = r 2 Lr 1 +r 1 + p+3 p r 2 +2yr 2 8 r 2 p 1 y 2 for any r 1,r 2 DL 2, where r = R d r 2 is the unique solution of Lr+r = Lr 2 + ψd,yr 2. Here, the domain DL of L defined in 5 is the set of all r L 2 ρ such that Lr L2 ρ. Furthermore, L d has two nonnegative eigenvalues λ = and λ = 1 with eigenfunctions W λ d such that W 1,2 d y = c 1 y 2 1 d 1 p y +d1 d 1 p 1, Wd 1+dy p+1,2 y = c, p 1+dy p+1 p with = 2 c λ p +λ y 2 1 y 21 λ ρydy, and W λ,1 d is the unique solution of the equation Lr+r = λ p+3 r 2 2yr 2 p + 8 r 2 p1 y 2 with r 2 = W λ,2 d. We also have for λ =,1 W d λ H +1 d 2 d Wd λ H C, d < Note that we have the following relations for λ = or λ = 1 φ W λ d, F λ d = 1 and φ W λ d, F d 1 λ =. 38 Let us introduce for λ {,1} the projectors π λ r, and π d r for any r H by π λ d r = φ W λ d,r, 39 r = π d r F d y+ π d 1r F d 1y+ π d r, 4 1 In section 4 of [18], we had non explicit normalizing constants c λ = c λ d. In Lemma 2.4 in [24], the authors compute the explicit dependence of c λ d. 1

11 and the space H d {r H π d 1r = π d r = }. Introducing the bilinear form ϕ d q,r = ψd,yq 1 r 1 +q 1 r 1 1 y2 +q 2 r 2 ρdy, 41 where ψd,y is defined in 26, we recall from Proposition 4.7 page 9 in [18] that there exists C > such that for all d < 1, for all r H d, 1 C r 2 H ϕ dr,r C r 2 H. 42 Furthermore, if r H, then 1 r H π d C r + πd 1 r + ϕ d r,r C r H where r = π d r. 43 In the following section we recall from [4] the spectral properties of L d. 3.3 Spectral theory of the operator L d From Section 3 in [4], we know that L d has one nonnegative eigenvalue λ = with eigenfunction F d κd,y y =. 44 Note that for some C > we have d < 1, 1 F d C H C and d Fd H C 1 d2. 45 We know also that the operator L d conjugate of L d with respect to φ is given by L d r1 R d r 2 = r 2 Lr 1 +r 1 + p+3 p r 2 +2yr 2 8 r 2 p 1 y 2 46 for any r 1,r 2 DL 2, where r = R d r 2 is the unique solution of Lr+r = Lr 2 + ψd,yr Furthermore, L d have one nonnegative eigenvalue λ = with eigenfunction W d such that W,2 d y = c κd,y and 1 c = 4κ2 ρy p 1 y2dy 48 and W d,1 is the unique solution of the equation with r 2 = W d,2. We also have for λ =,1 Lr+r = p+3 p r 2 2yr p r 2 1 y 2 49 W d H +1 d 2 d Wd λ H C, d <

12 Moreover, we have Let us introduce the projectors π dr and πd r for any r H by and the space Introducing the bilinear form φ W, d F d = π d r = φ W d,r, 52 r = π d r F d y+ πd r. 53 H d {r H π d r = }. 54 ϕ d q,r = ψd,yq 1 r 1 +q 1r 11 y 2 +q 2 r 2 ρdy, 55 where ψd,y is defined in 27, we recall from Proposition 3.7 page 596 in [4] that there exists C > such that for all d < 1, for all r H d, 3.4 A modulation technique 1 C r 2 H ϕ d r,r C r 2 H. 56 We start the proof of Theorem 4 here. Let us consider w C[s,,H for some s R a solution of equation 4 such that s s,ews, s ws Eκ, and ws κd s ws,.ω H ǫ 57 for some d,1, Ω S m and ǫ > to be chosen small enough. Our aim is to show the convergence of ws, s ws as s to some κd,ω,, for some d,ω close to d,ω. As one can see from 57, w, s w is close to a one representative of the family of the non-zero stationary solution S {κd,y,ω, d < 1,Ω S m }. From the continuity of w, s w from [s, to H, ws, s ws will stay close to asoliton from S, at least for a short time after s. In fact, we can do better, and impose some orthogonality conditions, killing the zero directions of the linearized operator of equation 4 see the operator L d defined in 25. From the invariance of equation 4 under rotations in R m, we may assume that Ω = e

13 We recall that at this level of the study in the complex case i.e. for m = 2, we were able to modulate w, s w as follows [ wy,s = e iθs κds,y + s wy,s for some well chosen ds,1 and θs R, such that π ds q1,1 s = π ds q 2,1 s ] q1 y,s. 59 q 2 y,s q1,2 s = q 2,2 s where π d and πd are defined in 39 and 52 and q = q 1,q 2 is small in H. From 59, we see that we have a rotation in the complex plane, which has to be generalized to the vector-valued case. In order to do so, we introduce for i = 2,...,m sinθ i 1. R i sinθ i Note that R i is an m m orthonormal matrix which rotates the e 1,e i -plane by an angle θ i and leaves all other directions invariant. We introduce R θ by R θ R 2 R 3 R m, 61 where θ = θ 2,θ 3,,θ m. Clearly, R θ is an m m orthonormal matrix. We also define A j by In the appendix, we show a different expression for A j : A j = R R θ θ. 62 θ j A j = R θ R θ. 63 θ j In fact, this formalism is borrowed from Filippas and Merle [8] who introduced the modulation technique for the vector-valued heat equation t u = u+ u p u. We are ready to give our modulation technique result well adapted to the vector-valued case: Proposition 3.1. Modulation of w with respect to κd,.ω, where Ω R m There exists ǫ > and K 1 > such that for all ǫ ǫ if v H, d,1 and ˆθ = ˆθ 2,..., ˆθ m R m are such that i = 2,...,m,cos ˆθ i 3 [ ] 4 and ˆq κˆd,.e1 H ǫ where v = Rˆθ + ˆq, 13

14 then, there exist d,1, ˆθ = ˆθ 2,..., ˆθ m R m such that π d q1,1 =, and π d q1,j =, j = 2,..m, 64 q 2,1 where q = q 1,q 2 is defined by: log 1+d log 1 d q 2,j 1+ ˆd + θ 1 ˆd ˆθ C ˆq H K 1 ǫ, i = 2,...,m, 1 2 and q H K 1 ǫ. In order to prove this proposition, we need the following estimates on the matrix A j given in 62 and 63: Lemma 3.2 Orthogonality and continuity results related to the matrix A i 62. i For any i {2,...,m}, A i e 1 = ii For any i {2,...,m}, z R m, we have m j=i+1 cosθ j e i A i z z. Proof. The proof is straightforward though a bit technical. For that reason, we give it in Appendix A Now, we are ready to prove Proposition 3.1. Proof of Proposition 3.1. The proof is similar to the complex-valued case. However, since our notations are somehow complicated, we give details for the reader s convenience. First, we recall that θ = θ 2,θ 3,...,θ m R m. From 39 and 52, we see that the condition 64 becomes Φv,d,θ = where Φ CH,1 R m,r m is defined by where V = V1 V 2 Φv,d,θ Φ 2 v,d,θ Φv,d,θ =. = φ Φ m v,d,θ R m R m is given by V = R θ v. V1,1 φ κd,., W d V 2,1 V1,j V 2,j, W d j=2...m 65 We claim that we can apply the implicit function theorem to Φ near the point ˆv, ˆd, ˆθ with ˆv = Rˆθκˆd,.e 1,. Three facts have to be checked: 1-First, note that ˆv = R ˆv, hence ˆθ ΦRˆθκˆd,.e 1,, ˆd, ˆθ =. 2-Then, we compute from 65, for all u H, U1,1 D v Φv,d,θu = φ, U W d, 2,1 14

15 and for all j = 2...m, we have so we have from 37 and 5 D v Φj v,d,θu = φ U1,j U 2,j, W d, D v Φv,d,θ C and D v Φj v,d,θ C Let J Φ, Φ j,j=2..m the jacobian matrix of Φ with respect to d,θ, and D its determinant so d φ θ2 φ θm φ d φ2 θ2 φ2 θm φ2 J d φm θ2 φm θm φm Then, we compute from 65: d Φ = φ d κd,.,, W d +φ V1,1 V 2,1 κd,., d Wd, 68 and for i,j = 2,..m V1,j d Φj = φ, V d Wd, 69 2,j V1,1 θi Φ = φ θi, V W d = φ < e 1, R θ v 1 >, W, d 7 2,1 < e 1, R θ v 2 > θi Φj = φ < e j, R θ v 1 >, d Wd < e j, R. 71 θ v 2 > Now, we assume that θ + log 1+d log 1 d ˆd 1 ˆd + κˆd,.e1 v H Rˆθ ǫ 1 73 for some small ǫ 1 >. In the following, we estimate each of the derivatives whose expressions where given above. - Since d κd,y = 2κ p1 d 2 F d, by definiftion 35 and 9, it follows from the orthogonality condition 38 that Therefore, from 68, we write φ d κd,.,, W d = d Φ = 2κ p1 d 2 +φ V1,1 V 2,1 15 2κ p1 d 2. κd,., d Wd. 74

16 Since Since, V1,1 V 2,1 < e1,r = θ v 1 > < e 1,R θ v, we write 2 > V1,1 κd,. H R V 2,1 θ v κd,.e1 H 75 R θ R ˆθ v H + R ˆθ v κˆd,.e1 H + κˆd,. κd,. H. θ,θ R, R θ R θ + R θ R θ C θ θ, see 67 below for R θ, and use an adhoc change of variables for R θ, recalling the following continuity result from estimate 174 page 11 in [18]: κd 1,. κd 2,. H C 1+d1 1+d2, 76 1 d 1 1 d 2 we see from the Cauchy-Schwartz inequality, 75, 37 and 73 that V1,1 κd,. H κd,.e1 V H Cǫ V 2, Since V1,j = V 2,j we write 2κ d Φ Cǫ1 p1 d < ej,r θ v 1 > < e j,r θ v 2 >, V1,j H R V θ v κd,.e1 H Cǫ 2,j 1 79 by the same argument as for 75. Using the Cauchy-Schwarz inequality together with 5, we see from 69 that d Φj Cǫ 1 1 d2. 8 From 63, we see that R θ = A i R θ. Therefore using ii of Lemma 3.2 and the fact that the rotation R θ does not change the norm in H, we write R θ κd,.ai e v 1 H κd,.e1 v Rθ H Cǫ 1, 81 by the same argument as for 75. Using the Cauchy-Schwarz identity together with 37, we see from 7 that θi Φ Cǫ By the same argument as for 82, we obtain from 71 θi Φj Cǫ1 if i j

17 Now, if i = j, noting from 63 that R θ v = A i R θ v = A iv, applying the operator A i to 77, then taking the scalar product with e i, we see from Lemma 3.2 that < e i, R θ v 1 > κd i,. < e i, R θ v 2 > Since we know from 44 and 51 that it follows from 71 that θi Φi m k=i+1 φκd,., W d = 1, m k=i+1 Collecting 78, 8, 82, 83 and 84 we see that Since cosθ k H Cǫ 1. cosθ k Cǫ1. 84 D 2κ p1 d 2 cosθ 3cosθ cosθ m m 2 Cǫ 1 1 d by hypothesis, we have the non-degeneracy of Φ voir 65 near the point ˆv, ˆd, ˆθ with ˆv = Rˆθκˆd,.e 1,. Applying the implicit function theorem, we conclude the proof of Proposition Dynamics of q, d and θ Let us apply Proposition 3.1 with v = w, s ws, ˆd = d and ˆθ =. Clearly, from 57 and 58, we have ˆq H ǫ. Assuming that ǫ ǫ defined in Proposition 3.1, we see that the proposition applies, and from the continuity of w, s wfrom[s, toh,wehaveamaximal s > s, suchthatws, s wscanbemodulated in the sense that [ wy,s κds,ye1 = R s wy,s θs + ] q1 y,s, 85 q 2 y,s where the parameters ds,1 and θs = θ 2 s,...,θ m s are such that for all s [s, s π ds q1,1 s =, and π ds q1,j s q 2,1 s =, j = 2,..m, 86 q 2,j s and i = 2,...,m s 1 2 and qs H ǫ 2K K 1 ǫ, 87 where K 1 > is defined in Proposition 3.1 and K 1 > 1 is a constant that will be fixed below in

18 Two cases then arise: - Case 1: s = + ; - Case 2: s < + ; in this case, we have an equality case in 87, i.e. s = 1 2 for some i = 2,...,m, or q s H = 2K K 1 ǫ. At this stage, we see that controlling the solution ws, s ws H is equivalent to controlling q H, d,1 and θs R m. Before giving the dynamics of this parameters, we need to introduce some notations. From 86, we will expand q and q respectively according to the spectrum of the linear operators L d and L d as in 4 and 53: where α 1,1 = π ds 1 q1,1 y,s q,1,1 y,s = α q 2,1 y,s 1,1 Fd 1 y+ q,2,1 y,s q1,j y,s q,1,j y,s j {1,...,m}, = q 2,j y,s q,2,j y,s q1,1 q 2,1 α,j = π ds, α,1 = π ds q1,j q 2,j q1,1 q 2,1 =, α,1 s = ϕ d =, α,j s = ϕ d q,1,j q,2,j q,1,1, q,2,1, q,1,j q,2,j q,1,1 q,2, and q,1,1 q1,1 q,2,1 j {1,...,m}, = π d q,1,j q,2,j q 2,1 = π d From 88, 89, and 56, we see that for all s s, 1 q,1,1 α,1 s C q H C α,1 s,2,1 1 q1,1 α 1,1 s +α,1 s C q H C α 1,1 s +α,1 s 92 2,1 1 q1,j α,j s H C α,j s C q 2,j q1,j for some C >. In the following proposition, we derive from 13 and 14 differential inequalities satisfied by α 1,1 s, α,1 s, α,j s, θ i s and ds. Introducing where q 2,j R s = F d q 1 ρdy, 93 F ds q 1 y,s = κd, e 1 +q 1 p+1 κd, p+1 κd, p q 1,1 p p+1 p+1 2 κd, p q1,1 2 κd, p 2 we claim the following q1,j 2, j=

19 Proposition 3.3. Dynamics of the parameters For ǫ small enough and for all s [s, s, we have: i Control of the modulation parameter i=2 θ i + d 1 d 2 C q 2 H. 95 ii Projection of equation 13 on the different eigenspaces of L d and L d α 1,1s α 1,1 s C q 2 H. 96 R α2 1,1 +α 2,j 4 q,2,1 2 +q 2 ρ p,2,j 1 y 2dy +C qs 3 H, 97 for j {2,...,m} and R s defined in 93, satisfying R s C qs 1+ p H where p = minp,2 > iii An additional relation For j {2,...,m}, we have: iv Energy barrier d 1 q 1,1 q 2,1 ρ 4 5ᾱ2,1 +C ds d 1 q 1,j q 1,j ρ 4 ds 5 α2,j +C q 2,2,1 q 2 2,j ρ 1 y 2 +C qs 2 H. 99 ρ 1 y 2 +C qs 2 H. 1 m α 1,1 s C α,1 s+c 1 α,j s. 11 Proof. The proof follows the general framework developed by Merle and Zaag in the real case see Proposition 5.2 in [18], then adapted to the complex-valued case in [4] 5see Proposition 4.2 page 5915 in [4]. However, new ideas are needed, mainly because we have m rotation parameters in the modulation technique see Proposition 3.1 above, rather than only one in the complex-valued case. For that reason, in the following, we give details only for the new terms, referring the reader to the earlier literature for the old terms. Let us first write an equation satisfied by q defined in 85. We put the equation 4 satisfied by w in vectorial form: j=2 s w 1 = w 2 s w 2 = Lw 1 2p+1 p 2 w 1 + w 1 p w 1 p+3 p w 2 2y y w We replace all the terms of 12 by their expressions from 85. Precisely, for the terms of the right hand side of 12 we have: s w 1 = R θ d d κe 1 + s q 1 + m i=1 θ i R θ θ j κ d e 1 +q 1, s w 2 = R θ s q 1 + m i=1 θ i R θ θ j q 2. For the terms on the left hand side of 12 we have: w 2 = R θ q 2, Lw 1 = R θ Lκ d e 1 +Lq 1, w 1 = κ d e 1 +q 1, y w 2 = R θ y q 2. 19

20 Then, multiplying by R θ, using the fact that κd,, is a stationnary solution and dissociating the first and jth component of these equations, we get for all s [s, s, for all j {2,...,m}: s s q1,1 q 2,1 q1,j q 2,j = L q1,1 ds q 2,1 = L ds q1,j q 2,j + + d d κd,y s f ds,1 q 1 f ds,j q 1 i=2 θ i s ai,1,j a i,2,j i=2 θ i s ai,1,1, 13 a i,2,1, 14 where L ds, L ds,f ds,1 and f ds,j are defined in 31, 33, 28 and 29, and a i by ai,1 a i,2 = Ai κ d e 1 +q 1, 15 A i q 2 with a i,1 = a i,1,1,a i,1,2,...,a i,1,m R m and a i,2 = a i,2,1,a i,2,2,...,a i,2,m R m. Let i {2,,m}, Projecting equation 13 with the projector π λ d 39 for λ = and λ = 1, we write π d λ s q1,1 q 2,1 = π λ d L q1,1 ds + π d q λ d s π d d κd,y 2,1 f ds,1 q 1 λ θ is π λ d ai,1,1. 16 i=2 a i,2,1 Note that, expect the last term, all the terms of 16 can be controled exactly like the real case using 87 for details see page 15 in [18]. So, we recall that we have: π d λ s q1,1 q 2,1 π d λ L d α λ,1 C 1 d 2 d q H, 17 q1,1 q 2,1 = λα λ,1, 18 π λ d C f ds,1 q 1 q 2 H, 19 π λ d d κd,y 2κ 2κ = p1 d 2 πd λ Fd = p1 d 2 δ λ,. 11 Now, we focus on the study of the last term of 16. From the definition of a i 15 and i of Lemma 3.2, we have: ai,1,1 < e1,a = κ i e 1 > < e1,a d + i q 1 > < e1,a = i q 1 >. < e 1,A i q 2 > < e 1,A i q 2 > a i,2,1 Applying the projector π λ d 39, we get π λ d ai,1,1 = π λ d < e1,a i q 1 > < e1,a C i q 1 > < e 1,A i q 2 > < e 1,A i q 2 > H a i,2,1 C < e 1,A i q 1 > H + < e 1,A i q 2 > L 2 ρ

21 Using ii of Lemma 3.2, we have: < e 1,A i q 2 > 2 L 2 ρ = < e 1,A i q 2 > 2 ρdy A i q 2 2 ρdy and by the same way, using ii of Lemma 3.2 and the definition of H 7, we have < e 1,A i q 1 > 2 H = From 111, 112 and 113, we have < e 1,A i q 1 > 2 ρdy + q 2 2 ρdy, 112 < e 1,A i y q 1 > 2 1 y 2 ρdy q y 2 y q 1 2 ρdy. 113 π d λ ai,1,1 a i,2,1 C q H. 114 Using 17, 17, 18, 19, 11, 114, and the fact that α,1 α,1 see 9, we get for λ =,1: 2κ p1 d 2 d C 1 d 2 d q H +C q 2 H +C q H ᾱ 1 s ᾱ 1s C 1 d 2 d q H +C q 2 H +C q H Now, projecting equation 14 with the projector π d π d q1,j s = π d q L q1,j ds + π d 2,j q 2,j f ds,j q 1 θ i 115 j=2 j=2 θ i , where j {2,,m}, we get: θ i ai,1,j s πd i=2 a i,2,j. 117 From thecomplex-valued casewerecall that wehave fordetails seepage5917 in[4], together with Lemma 3.2: π d q1,j s C q 2,j 1 d 2 d q H, 118 π L d q1,j d =, 119 q 2,j π d C q f H ds q 1 12 π d κd = Thus, only the last term in 117 remains to be treated in the following. From the definition of a i 15, we recall that ai,1,j < ej,a = i e 1 > κ d < ej,a + i q 1 >. 122 < e j,a i q 2 > By i of Lemma 3.2: θ is i=2 a i,2,j ai,1,j a i,2,j = θ js m cosθ l κ d l=j θ is < ej,a i q 1 >, 123 < e j,a i q 2 > i=2

22 where by convention m l=m+1 Applying the projection π d θ i s πd i=2 ai,1,j a i,2,j where, we use the fact that cosθ l = 1 if j = m. to 123 and using 121, we see that θ j s m π d l=j+1 cosθ l m i=2 θ i π s d < ej,a i q 1 > < e j,a i q 2 > m C q H θ is, 124 i=2 < ej,a i q 1 > C < e j,a i q 2 > q H, 125 which follows by the same techniques as in 111, 112 and 113. Using 117,118, 119, 12 and 124, and recalling from 87 that we get for any j {2,...,m}: m l=j+1 cosθ l 1 2 m j, θ js C 1 d 2 d q H +C q 2 H +C q H Using 126 together with 115, we see that θ j + s d 1 d 2 C d 1 d 2 q H +C q 2 H +C q H j=2 Thus, using again 87 and taking ǫ small enough, we get θ js + d 1 d 2 C q 2 H, j=2 which yields 95. Then, using 116 together with 95 gives 96. θ i. 126 i=2 θ i, For estimations , the study in the complex case Subsection 4.3 page 5914 in [4] can be adapted without any difficulty to the vector-valued case. For the reader convenience, we detail for example the energy barrier 11: Using the definition of qy,s 59, we can make an expansion of Ews, s ws 8 for q in H and get after from straightforward computations: Ews, s ws = Eκ,+ 1 2 [ ϕ d q1,1 q 2,1, q1,1 q 2,1 + ϕ d where ϕ d, ϕ d and F d q 1 are defined in 41, 55 and 94. i=2 q1,j q 2,j, i=2 q1,j q 2,j ] F d q 1 ρdy

23 Using the argument in the real case see page 113 in [18] we see that for some C,C 1 > we have: q1,1 q1,1 ϕ d, C α 2 1,1 C 1 α 2, q 2,1 q 2,1 From 89, 91 and 56, we see by definition that q1,j ϕ d, q 2,j q1,j q 2,j = α 2,j. 129 Since we have from 93, 98, 87, 129 and 92: m F d q 1 ρdy C qs p+1 H Cǫ p α 2 1,1 +α2,1 + α 2,j, 13 we get Using 11, 127, 128 and 13, we see that taking ǫ small enough so that Cǫ p C 1 4, Ews, s ws Eκ, which yields Exponential decay of the different components i=2 C 2 + C 1 α 2,1 4 C α2 1, C m 1 α 2,j 4. i=2 Our aim is to show that qs H and that both θ and d converge as s. An important issue will be to show that the unstable mode α 1,1, which satisfies equation 92 never dominates. This is true thanks to item iv in Proposition 3.3. If we introduce λs = 1 1+ds 2 log,as = α 1,1 s 2 andbs = α,1 s ds α,j s 2 +R s 131 note that ds = tanhλs, then we see from 98, 92 and 87 that for all s [s, s R s = bs α,1 s 2 + m j=2 α,js 2 C ǫ p α 1,1 s 2 +α,1 s 2 + m j=2 α,js 2, hence 99 1 α,1s α,j s a b 11 1 α,1s α,j s a 132 j=2 for ǫ small enough. Therefore, using Proposition 3.3, estimate 87, 92 and the fact that λ s = d s 1 ds 2, we derive the following: Claim 3.4. Relations between a, b, λ, θ, q 1,1q 2,1 ρ and q 1,jq 2,j ρ There exist positive ǫ 4, K 4 and K 5 such that if ǫ ǫ 4, then we have for all s [s, s] and j = 2,...,m: i Size of the solution and 132 holds. 1 K 4 as+bs qs 2 H K 4as+bs K 2 4 ǫ2, 133 θ s + λ s K 4 as+bs K4 qs 2 2 H, 134 q 1,1 q 1,1 ρ K 4as+bs, 135 q 1,j q 1,j ρ K 4bs, j=2 j=2

24 ii Equations d ds 3 2 a K 4ǫb a 5 2 a K 4ǫb, 137 b 8 q,2,1 2 +q 2 ρ p,2,j 1 y 2dy +K 4ǫa+b, 138 q 1,1 q 2,1 +q 1,j q 2,j ρ 3 5 b+k 4 q,2,1 2 +q2,j 2 ρ 1 y 2 +K 4a. iii Energy barrier we have as K 5 bs. End of the Proof of Theorem 4. Now, we are ready to finish the proof of Theorem 4 just started at the beginning of Section 3.4. Let us define s 2 [s, s] as the first s [s, s] such that as bs 5K where K 4 is introduced in Corollary 3.4, or s 2 = s if 139 is never satisfied on [s, s]. We claim the following: Claim 3.5. There exist positive ǫ 6, µ 6, K 6 and f C 1 [s,s 2 ] such that if ǫ ǫ 6, then for all s [s,s 2 ]: i 1 2 fs bs 2fs and f s 2µ 6 fs, ii qs H K 6 qs H e µ 6s s K 6 K 1 ǫ e µ 6s s. Proof. The proof of Claim 5.6 page 115 in [18] remains valid where fs is given by fs = bs+η 6 q 1,1 q 2,1 + q 1,j q 2,j ρ, j=2 where η 6 > is fixed small independent of ǫ. Claim 3.6. i There exists ǫ 7 > such that for all σ >, there exists K 7 σ > such that if ǫ ǫ 7, then s [s 2,mins 2 +σ, s], qs H K 7 qs H e µ 6s s K 7 K 1 ǫ e µ 6s s and θ i s C K 7K 1 ǫ 2 2µ 6 where µ 6 has been introduced in Claim 3.5. ii There exists ǫ 8 > such that if ǫ ǫ 8, then s s 2, s], bs as 5K 4 e s s K 5 14 where K 4 and K 5 have been introduced in Corollary 3.4. Proof. The proof is the same as the proof of Claim 5.7 page 117 in [18]. 24

25 Now, in order to conclude the proof of Theorem 4, we fix σ > such that 5K σ K 5 1 2K 5, where K 4 and K 5 are introduced in Claim 3.4. Then, we fix the value of K = max2,k 6,K 7 σ, 141 and the constants are defined in Claims 3.5 and 3.6. Then, we fix ǫ i ǫ = min 1,ǫ 1, for i {4,6,7,8} 2K K 1 and the constants are defined in Claims 3.4, 3.5 and 3.6. Now, if ǫ ǫ, then Claim 3.4, Claim 3.5 and Claim 3.6 apply. We claim that for all s [s, s], qs H K qs H e µ 6s s K K 1 ǫ e µ 6s s = ǫ 2 e µ 6s s. 142 Indeed, if s [s,mins 2 + σ, s], then, this comes from ii of Claim 3.5 or i of Claim 3.6 and the definition of K 141. Now, if s 2 + σ < s and s [s 2 + σ, s], then we have from 14 and the definition of σ, bs as 2K 5 on the one hand. On the other hand, from iii in Claim 3.4, we have as K 5 bs, hence, as = bs = and from 133, qy,s, hence 142 is satisfied trivially. In particular, we have for all s [s, s], q H ǫ 2 and 1 C ǫ2 3 µ 2 4, hence, by definition 6 of s given right before 85, this means that s =. From i of Claim 3.6 and 134, we have s s, qs H ǫ 2 e µ 6s s and θ s + λ s K4 2 ǫ 2 4 e 2µ 6s s, 143 where, θs = θ 2 s,...,θ m s. Hence, there is θ R m, λ in R such that θs θ, λs λ as s and s s, λ λs C 1 ǫ 2 e 2µ 6s s = C 2 ǫ 2 e 2µ 6s s s s, θ θs C 1 ǫ 2 e 2µ 6s s = C 2 ǫ 2 e 2µ 6s s for some positive C 1 and C 2. Taking s = s here, we see that λ λ + θ C ǫ, where Ω = R θ e 1. If d = tanhλ, then we see that d d C 3 1 d 2 ǫ. Using the definition of q 85, 143, 144 and 145 we write ws κd, Ω s ws H ws κds, Ω + κds, κd s ws, Ω H H + κd, H R θs e 1 R θ e qs H +C λ λs +C θ θs C 4 ǫ e µ 6s s, 146 where, we used the fact that θ R m O m is a Lipschitz function see 67 to be convinced and λ R κd, H is also Lipschitz, where d = tanhλ see 76.This concludes the proof of Theorem 4 in the case where Ω = e 1 see 58. From rotataion invariance of equation 4, this yields the conclusion of Theorem 4 in the general case. 25

26 A A some technical estimates In this section, we give the proof of estimate 63 and Lemma 3.2. Proof of estimate 63: Using 61, we see that From 6, we see that R θ θ j = R 2...R j R j θ j R j+1...r m. 147 R θ θ j = Π j R j θ j + π 2 = R jθ j + π 2 Π j, 148 where Π j is the orthogonal projection on the plane spanned by e 1 and e j, and the rotation R j α is given by considering the matrix of R j defined in 6, and changing θ j into α. Since R j R j θ j + π 2 = R j π 2, it follows from 61 and 148 that R R θ θ = Rm...Rj+1 θ R j π j 2 Π jr j+1...r m. 149 By the same argument, we drive that R θ θ j R θ has the same expression, thus, 63 holds from 62 and 149. Now, we give the proof of Lemma 3.2. Proof of Lemma 3.2: i We first give the expression of the m m matrix R θ defined 61. Indeed, using 6 and 61, we have: ϕ 2,m sinθ 2 sinθ j ϕ 2,l sinθ m ϕ 2,m sinθ 2 ϕ 3,m cosθ R θ,k,l. R θ. cosθ j. sinθ k ϕ k+1,m sinθ m ϕ m+1,m cosθ m where for k 1, l 2: sinθ l ϕ 2,l if k = 1 sinθ R θ,k,l = k sinθ l ϕ k+1,l if 2 k l cosθ k if k = l if k l

27 with l cosθ ϕ k,l = n if k l n=k 1 if k l In fact, we will prove the following identities, which imply item i: A For all i,j {2,...,m}, such that i j, we have B For all i {2,...,m} < e j,a i e 1 >=. C For all i {2,...,m}, we have < e 1,A i e 1 >=. < e i,a i e 1 >= ϕ i+1,m, where A i and ϕ i+1,m are given in 62 and 152. Proof of A. Let i,j {2,...,m}, such that i j. The idea is to compute < R θ e j, R θ e 1 > instead of < e j,a i e 1 >. In fact, using the conservation of the inner product after a rotation and the fact that A i = R θ R θ by 62, we have: < e j,a i e 1 >=< R θ e j,r θ A i e 1 >=< R θ e j, R θ e 1 >. 153 In the following, we distinguish two cases: - Case 1: i j, - Case 2: i j +1. We first handle Case 1. Case 1: i j. Using 15 and its derivative with respect to θ i, we write: and Therefore, R θ e j = R θ e j k=1,...,m = R θ e 1 = R θ e 1 k=1,...,m = < R θ e j, R θ e 1 >= k=1 sinθ j ϕ 2,j, if k = 1 sinθ k sinθ j ϕ k+1,j if 2 k j cosθ j if k = j if k j +1 R θ,k,j R θ k,1 ϕ sinθ 2,m i if k = 1 ϕ sinθ i sinθ k+1,m k if 2 k i ϕ i+1,m, if k = i if k i+1. ϕ 2,m i ϕ k+1,m = sinθ i sinθ j ϕ 2,j + sinθ k sinθ i sinθ k sinθ j ϕ k+1,j ϕ i+1,m sinθ i sinθ j ϕ i+1,j ϕ 2,m i = sinθ i sinθ j ϕ 2,j + sinθ 2 ϕ k+1,m k ϕ k+1,j ϕ i+1,m ϕ i+1,j

28 In order to transform the sum term in the previous identity, we make in the following a finite induction where the parameter q decreases from i to 1: Lemma A.1. We have: q {1,...,i }, i sinθk 2 ϕ k+1,m ϕ k+1,j ϕ i+1,m ϕ i+1,j = 157 q sinθk 2 ϕ k+1,m ϕ k+1,j ϕ q+1,mϕ q+1,j. Remark: If q = 1, the sum in the right hand side is naturally zero. Proof. See below. Applying this Lemma, we conclude the proof of A in Case 1 i. e. when i j. Indeed, from 156 and Lemma A.1 with q = 1 we write < R θ e j, R θ ϕ2,m e 1 >= sinθ i sinθ j ϕ 2,j ϕ 2,m ϕ 2,j =. It remains now to prove Lemma A.1. Proof of Lemma A.1. First, we give the following: Claim A.2. We have ϕ i,j = ϕ i+1,j. Proof. Since i j, we have two cases: - If i j 2: trivial. - If i = j : ϕ i,j = ϕ i,i = and ϕ i+1,j = ϕ i+1,i = 1, and the result follows. Now, we are ready to start the proof of Lemma A.1. Let us prove the result using an induction with a decreasing index. For q = i 1, 157 is satisfied using Claim A.2. Assume now that 157 is true for q = i,...,2 and let us prove it for q. Using 157 with q, we write i q q q sinθk 2 ϕ k+1,m ϕ k+1,j ϕ i+1,m ϕ i+1,j = sinθk 2 ϕ k+1,m ϕ k+1,j +sinθ 2 ϕ q+1,m q ϕ q+1,j ϕ q+1,mϕ q+1,j = sinθk 2 ϕ k+1,m ϕ k+1,j cosθ 2 ϕ q+1,m q ϕ q+1,j = sinθk 2 ϕ k+1,m ϕ k+1,j ϕ q,mϕ q,j. Thus, 157 is satisfied for q. This concludes the proof of Lemma A.1 and identity A when i j. 28

29 Now, we handle Case 2. - Case 2: i j +1. Using 154 and 155, we write: < R θ e j, R θ e 1 >= k=1 R θ,k,j R θ k,1 ϕ 2,m j = sinθ i sinθ j ϕ 2,j + sin 2 ϕ k+1,m ϕ j+1,m θ k ϕ k+1,j cosθ j. 158 In order to transform the sum term in the previous identity, we make in the following a finite induction where the parameter q decreases from j to 1: Lemma A.3. We have: q {1,...,j }, j sin 2 θ k ϕ k+1,m ϕ k+1,j cosθ j ϕ j+1,m = 159 q sin 2 θ k ϕ k+1,m ϕ k+1,j ϕ q+1,j ϕ q+1,m. Remark: If q = 1, the sum in the right hand side is naturally zero. Proof. See below. Applying this Lemma, we conclude the proof of A in Case 2 i. e. when i j+1. Indeed, from 158 and Lemma A.3 with q = 1 we write < R θ e j, R θ ϕ2,m e 1 >= sinθ i sinθ j ϕ 2,j ϕ 2,m ϕ 2,j =. It remains now to prove Lemma A.3. Proof of Lemma A.3. We prove the result using an induction with a decreasing index. For q = j, 159 is satisfied. Assume now that 159 is true for q = j,...,2 and let us prove it for q. Using 159 with q, we write j q q q sin 2 θ k ϕ k+1,m ϕ k+1,j cosθ j ϕ j+1,m = sin 2 θ k ϕ k+1,m ϕ k+1,j ϕ q+1,j ϕ q+1,m = sin 2 θ k ϕ k+1,m ϕ k+1,j +sin 2 θ q ϕ q+1,m ϕ q+1,j ϕ q+1,j ϕ q+1,m = sin 2 θ k ϕ k+1,m ϕ k+1,j ϕ q,m ϕ q,j. Thus, 159 is satisfied for q. This concludes the proof of Lemma A.3. 29

30 Proof of B: As for 153 we have: From 15, we have: < e 1,A i e 1 >=< R θ e 1, R θ e 1 >. Therefore, using 155, we have: < R θ e 1, R θ e 1 >= sinθ i R θ e 1 = ϕ 2,m,sinθ 2 ϕ 3,m,sinθ i ϕ i+1,m,,sinθ m. ϕ 2,m 2 + i sin 2 θ k ϕ k+1,m 2 + sinθ i ϕ i+1,m 2. In order to transform the sum term in the previous identity, we make in the following a finite induction: Lemma A.4. We have: q {2,...,i }, ϕ i 2,m 2 + l=2 sin 2 θ l ϕ l+1,m 2 = ϕ i q,m 2 + Remark: If q = i, the sum in the right hand side is naturally zero. Using Lemma A.4 with q = i we get < R θ e 1, R θ e 1 >= sinθ i ϕ i+1,m 2 + sinθ i ϕ i+1,m 2 =, which yields the result. In order to conclude B we give the proof of Lemma A.4. l=q sin 2 θ l ϕ l+1,m Proof of Lemma A.4. We proceed by induction for q {2,...,i}. For q = 2, 16 is satisfied. Assume that 16 is true for q = 2,...,i and prove it for q+1. Using 16 with q, we write ϕ 2,m 2 + i l=2 sin 2 θ l ϕ l+1,m 2 = ϕ i q,m 2 + l=q = cos 2 θ q ϕ q+1,m 2 +sin 2 θ q ϕ q+1,m 2 + = ϕ q+1,m 2 + i l=q+1 sin 2 θ l ϕ l+1,m 2. sin 2 θ l ϕ l+1,m 2 i l=q+1 sin 2 θ l ϕ l+1,m 2 Thus 16 is satisfied for q +1. This concludes the proof of Lemma A.4 and identity B. Proof of C: Consider i {2,...,m}. As for 153 we have: < e i,a i e 1 >=< R θ e i, R θ e 1 >. Using 151 and 155 i < e i,a i e 1 >= sin 2 θ i ϕ i+1,m ϕ 2 2,i + sin 2 θ k ϕ 2 k+1,i +cos 2 θ i ϕ i+1,m. 161 In order to transform the sum term in the previous identity, we make in the following a finite induction: 3

31 Lemma A.5. We have: q {2,...,i}, i ϕ 2 2,i + sin 2 θ l ϕ 2 l+1,i = i ϕ2 q,i + sin 2 θ l ϕ 2 l+1,i. 162 l=2 Remark: If q = i, the sum in the right hand side is naturally zero. From 161 and 162 with q = i we get < R θ e i, R θ e 1 >= sin 2 θ i ϕ i+1,m +cos 2 θ i ϕ i+1,m = ϕ i+1,m. which yields the result. In order to conclude C we give the proof of Lemma A.5. Proof of Lemma A.5. We proceed by induction for q {2,...,i}. For q = 2, 162 is satisfied. Assume now that 162 is true for q = 2,...,i and prove it for q+1. Using 162 with q, we write l=q i i ϕ 2 2,i + sin 2 θ l ϕ 2 l+1,i = ϕ 2 q,i + sin 2 θ l ϕ 2 l+1,i l=2 l=q i = cos 2 θ q ϕ 2 q+1,i +sin2 θ q ϕ 2 q+1,i + i = ϕ 2 q+1,i + l=q+1 sin 2 θ l ϕ 2 l+1,i. l=q+1 sin 2 θ l ϕ 2 l+1,i Thus 162 is satisfied for q +1. This concludes the proof of Lemma A.5. iiwe recall from 61 that we have so by 62, A j is given explicitly by R θ θ j = R 2 R j R j θ j R j+1 R m, A j = R m R m R j R j θ j R j+1 R m. From a straightforward geometrical observation, we can see that the rotation conserves the euclidien norm in R m. For R j θ j, it can be seen as a composition of a projection on the plane e 1,e j and a rotation with angle θ j + π 2, which decreases the norm. This concludes the proof of Lemma 3.2. References [1] S. Alinhac. Blowup for nonlinear hyperbolic equations. Progress in Nonlinear Differential Equations and their Applications, 17. Birkhäuser Boston Inc., Boston, MA, [2] S. Alinhac. A minicourse on global existence and blowup of classical solutions to multidimensional quasilinear wave equations. In Journées Équations aux Dérivées Partielles Forges-les-Eaux, 22, pages Exp. No. I, 33. Univ. Nantes, Nantes,

32 [3] C. Antonini and F. Merle. Optimal bounds on positive blow-up solutions for a semilinear wave equation. Internat. Math. Res. Notices, 21: , 21. [4] A. Azaiez. Blow-up profile for the complex-valued semilinear wave equation. Trans. Amer. Math. Soc., 3678: , 215. [5] A. Azaiez, N. Masmoudi, and H. Zaag. Blow-up rate for a semilinear wave equation with exponential nonlinearity in one space dimension. Accepted in the proceedings of the conference in honor of Pr. Abbas Bahri., 215. [6] L. A. Caffarelli and A. Friedman. The blow-up boundary for nonlinear wave equations. Trans. Amer. Math. Soc., 2971: , [7] R. Côte and H. Zaag. Construction of a multisoliton blowup solution to the semilinear wave equation in one space dimension. Comm. Pure Appl. Math., 661: , 213. [8] S. Filippas and F. Merle. Modulation theory for the blowup of vector-valued nonlinear heat equations. J. Differential Equations, 1161: , [9] J. Ginibre, A. Soffer, and G. Velo. The global Cauchy problem for the critical nonlinear wave equation. J. Funct. Anal., 111:96 13, [1] J. Ginibre and G. Velo. Regularity of solutions of critical and subcritical nonlinear wave equations. Nonlinear Anal., 221:1 19, [11] S. Kichenassamy and W. Littman. Blow-up surfaces for nonlinear wave equations. I. Comm. Partial Differential Equations, 183-4: , [12] S. Kichenassamy and W. Littman. Blow-up surfaces for nonlinear wave equations. II. Comm. Partial Differential Equations, 1811: , [13] H. A. Levine. Instability and nonexistence of global solutions to nonlinear wave equations of the form Pu tt = Au+Fu. Trans. Amer. Math. Soc., 192:1 21, [14] H. Lindblad and C. D. Sogge. On existence and scattering with minimal regularity for semilinear wave equations. J. Funct. Anal., 132: , [15] F. Merle and H. Zaag. Determination of the blow-up rate for the semilinear wave equation. Amer. J. Math., 1255: , 23. [16] F. Merle and H. Zaag. Determination of the blow-up rate for a critical semilinear wave equation. Math. Ann., 3312: , 25. [17] F. Merle and H. Zaag. On growth rate near the blowup surface for semilinear wave equations. Int. Math. Res. Not., 19: , 25. [18] F. Merle and H. Zaag. Existence and universality of the blow-up profile for the semilinear wave equation in one space dimension. J. Funct. Anal., 2531:43 121, 27. [19] F. Merle and H. Zaag. Openness of the set of non-characteristic points and regularity of the blow-up curve for the 1 D semilinear wave equation. Comm. Math. Phys., 2821:55 86, 28. [2] F. Merle and H. Zaag. Points caractéristiques à l explosion pour une équation semilinéaire des ondes. In Séminaire X-EDP. École Polytech., Palaiseau,