Asymptotic Stability of high-dimensional Zakharov-Kuznetsov solitons
We prove that solitons (or solitary waves) of the Zakharov-Kuznetsov (ZK) equation, a physically relevant high dimensional generalization of the Korteweg-de Vries (KdV) equation appearing in Plasma Physics, and having mixed KdV and nonlinear Schrödinger (NLS) dynamics, are strongly asymptotically stable in the energy space. We also prove that the sum of well-arranged solitons is stable in the same space. Orbital stability of ZK solitons is well-known since the work of de Bouard . Our proofs follow the ideas by Martel  and Martel and Merle , applied for generalized KdV equations in one dimension. In particular, we extend to the high dimensional case several monotonicity properties for suitable half-portions of mass and energy; we also prove a new Liouville type property that characterizes ZK solitons, and a key Virial identity for the linear and nonlinear part of the ZK dynamics, obtained independently of the mixed KdV-NLS dynamics. This last Virial identity relies on a simple sign condition, which is numerically tested for the two and three dimensional cases, with no additional spectral assumptions required. Possible extensions to higher dimensions and different nonlinearities could be obtained after a suitable local well-posedness theory in the energy space, and the verification of a corresponding sign condition.
- 1 Introduction
- 2 Linear Liouville property
- 3 Nonlinear Liouville Property
- 4 Proof of the asymptotic stability result
- 5 Stability of the sum of -solitons
- A Numerical Estimates for the Spectral Property
- B Proof of (5.19)
- C Linear waves versus Asymptotic stability in the energy space
- D Proof of Theorem 1.5
We are interested in the the Zakharov-Kuznetsov (ZK) equation
where is a real-valued function, , , denotes the laplacian. The ZK equation is a particular case of the generalized Zakharov-Kuznetsov (gZK) equation
where is such that if and if . We observe that when the spatial dimension is equal to , equation (1.2) becomes the well-known generalized Korteweg- de Vries (gKdV) equation.
The ZK equation was introduced by Zakharov and Kuznetsov in  to describe the propagation of ionic-acoustic waves in uniformly magnetized plasma in the two dimensional and three dimensional cases. The derivation of ZK from the Euler-Poisson system with magnetic field in the long wave limit was carried out by Lannes, Linares and Saut in . The ZK equation was also derived by Han-Kwan  from the Vlasov-Poisson system in a combined cold ions and long wave limit. Moreover, the following quantities are conserved by the flow of ZK,
The well-posedness theory for ZK and gZK has been extensively studied in the recent years. In the two dimensional case, Faminskii proved that the Cauchy problem associated to the ZK equation is globally well-posed in the energy space . The local well-posedness result was pushed down to for by Linares and Pastor  and to by Grünrock and Herr  and Molinet and the third author . The best result for the ZK equation in the three dimensional case was obtained last year by Ribaud and Vento . They proved local well-posedness in for . Those solutions were extended globally in time in . Note however that it is still an open problem to obtain well-posedness in and for the ZK equation. Finally, we also refer to [27, 28, 14, 45, 17] for more well-posedness results for the gZK equation with and to [42, 6, 7] for unique continuation results concerning ZK.
Note that if solves (1.2) with initial data , then is also a solution to (1.2) with initial data for any . Hence, , so that the scale-invariant Sobolev space for the gZK equation is , where . In particular, the gZK equation is -critical (or simply critical) if . In the sequel, we will say that the problem is subcritical if and supercritical if .
1.1. The elliptic problem
For , equation (1.2) admits special solutions of the form
where and satisfies
Observe that .
Assume that if and if . Then there exists a unique positive radially symmetric solution to (1.6) in , which is called a ground state. In addition, , for all and there exists such that
The solutions of (1.2) of the form (1.5) with are called solitary waves or solitons. They were proved by de Bouard in  to be orbitally stable in if and unstable for . In other words, the solitary waves associated to (1.2) are orbitally stable in the subcritical case and unstable in the supercritical case.
In the following, for any , we will denote by the operator which linearizes (1.6) around , i.e.,
In the case , we also denote .
Next, we gather some well-known facts about the operator (see Weinstein ).
Assume that if and if .
Then, the following assertions are true.
(i) is a self-adjoint operator and
(iii) has a unique single negative eigenvalue (with ) associated to a positive radially symmetric eigenfunction . Without loss of generality, we choose such that . Moreover, there exists such that , for all .
(iv) Let us define
1.2. Statement of the results
As already mentioned, de Bouard proved in  that the solitary waves of the gZK equation are stable in the subcritical case in the following sense.
Theorem 1.1 (Asymptotic stability).
Assume . Let . For any , there exists such that if and is a solution of (1.1) satisfying
then the following holds true.
There exist with , for some positive constant independent of , and such that
It will be clear from the proof that the convergence in (1.15) can also be obtained in regions of the form
Note that the maximal angle of improvement must be strictly less than on each side of the vertical line (see Figure 1). We also refer to Lemma 4.3 for more details and a relation with the nonlinear dynamics of the equation. Moreover, we expect the range of for which asymptotic stability occurs in to be sharp. Indeed, it will be shown in Appendix C that linear plane waves of ZK exist if and only if the velocity group vector has a negative -component and forms an angle with (as in Figure 1) satisfying .
The angle is also related to the linear part of ZK. In , Carbery, Kenig and Ziesler proved that
where is the Fourier multiplier associated to the symbol . This Strichartz estimate was used in  to improve the well-posed results for ZK at low regularity. Note that the multiplier cancels out along the cone . We also refer to Apendix C for an interesting relation between the angle and the linear plane waves of ZK.
Our proof does not rely on the structure of the nonlinearity of (1.1) (i.e. ) neither on the dimension . Actually, our main theorem could be extended to (1.1) in dimension or to the following generalization of gZK
where is a real number under the following conditions:
The spectral condition holds true. (Note that makes sense since is radial and orthogonal to , and we choose orthogonal to .)
This spectral condition was shown in the appendix to be true in dimension for , where is a real number satisfying .
On the other hand, in dimension , it is shown in the appendix that . Note however that in this case, one could try to verify the more general property: the operator restricted to the space is positive definite.
Theorem 1.2 (Nonlinear Liouville property around ).
Assume . Let . There exists such that if and is a solution of (1.1) satisfying for some function and some positive constant
then, there exist (close to ) and such that
We also prove a rigidity theorem for the solutions of the linearized gZK (or (1.18)) equation in spatial dimension around which are uniformly localized in the direction .
(Linear Liouville property around ) Assume . There exists such that for all , the following holds true. Let and be a solution to
where is defined in (1.8). Moreover, assume that there exists a constant such that
Then, there exists such that
Recall that the first result of asymptotic stability of solitons for generalized KdV equations was proved by Pego and Weinstein  in weighted spaces (see also  for some refinements on the weights). In , Martel and Merle have given the first asymptotic result for the solitons of gKdV in the energy space . They improved their result in  and generalized it to a larger class of nonlinearities than the pure power case in .
Their proof relies on a Liouville type theorem for -compact solutions around a soliton (similar to Theorem 1.2 in one dimension). Then, it is proved that a solution near a soliton converges (up to subsequence) to a limit object, whose emanating solution satisfies a good decay property. Due to the rigidity result, this limit object has to be a soliton.
It is worth noting that this technique of proof was also adapted to prove asymptotic stability in the energy space for other one dimensional models such as the BBM equation  and the BO equation .
We also refer to [38, 1, 40] for stability results for KdV and mKdV in and to [4, 16] for asymptotic results for the Gross-Pitaevskii equation in one dimension. For other results on asymptotic stability for nonlinear Schrödinger and wave equations, see [47, 5, 24, 21] and references therein.
About the proofs. Comparison with previous results. When proving Theorems 1.1, 1.2 and 1.3, we generalize the ideas of Martel and Merle [31, 34, 35] and Martel  to a multidimensional model. However, compared with these previous results, the higher dimensional case describing the ZK dynamics presents new challenges, that we explain in the following lines.
First of all, as far as we know, our results represent the first two dimensional model where asymptotic stability is proved, in the energy space, and with no nonstandard spectral assumptions on the linearized dynamical operator. As we expressed before, we only need to check the numerical condition
Obtaining a direct proof of this result seems far from any reasonable approach because the soliton , and therefore, the function , have no closed and explicit forms. This is the first difference with respect to the one dimensional case: we work with a solitary wave that is not explicit at all.
We will see through the proofs that ZK behaves as a KdV equation in the direction, and as a nonlinear Schrödinger (NLS) equation in the variable. In particular, we are able to prove monotonicity properties (see Lemma 3.3) along the direction and along a slighted perturbed cone around the direction (Lemma 4.3). This last result is, to our knowledge, new in the literature and makes use of the geometrical properties of the nonlinear ZK dynamics around a solitary wave. Remark 1.1 and the asymptotic stability result inside the set (see (1.17))
are deep consequences of these geometrical properties. Recall that such rich foliations are not present in the one dimensional case. We also complement our results by a simple linear analysis leading to the same formal conclusions, carried out in Appendix C.
Another barrier that appears in the higher dimensional case is the lack of control on the solution if we only assume bounds. We need such a control to ensure pointwise exponential decay around solitons at infinity for a compact part of the solution. In the one dimensional case, the proof of this fact is direct from the Sobolev embedding. However, since is not contained in in , we must prove new monotonicity properties at the level (cf. Lemma 3.6), which are obtained by proving new energy estimates.
No monotonicity property seems to hold for the direction, mainly because of the conjectured existence of trains of small solitons moving to the right in but without restrictions on the coordinate. From the point of view associated to the variable, such solutions represent movement of mass along the direction without a privileged dynamics. In particular, no asymptotic stability result is expected for a half-plane involving the variable only (see Fig. 2). This is the standard situation in many models like KP-I and NLS equations. However, here we are able to prove the asymptotic stability of ZK solitons because the KdV dynamics is exactly enough to control the movement of mass along the direction.
The second ingredient in the proofs of Theorems 1.1, 1.2 and 1.3 is a new virial identity in higher dimensions (cf. (2.36)), which holds only for the half space and for and slightly larger. Compared with the previous works by Martel and Merle, the additional dimensions make things harder because they induce transversal variations that seem to destroy any virial-type inequality. In order to overcome this difficulty, we use a different orthogonality condition for the function employed in the virial (see Lemma 2.3):
We emphasize that this condition is somehow natural and necessary if we want to get full control of the perturbations appearing from the variations of the virial terms. Without using this modified condition, any form of two-dimensional virial identity is no longer true.
Here is when the nonstandard spectral condition (1.26) appears: under the orthogonality condition (1.27), the virial identity holds provided (1.26) is satisfied. We prove that (1.26) holds for the case and , as expressed by some numerical computations obtained in Appendix A. This condition in fact generalizes the Martel one in  and seems to be the natural one for the case, as described when proving the nonlinear stability result (see (3.60)-(3.61) for instance). It is worth noting that this condition has already been used by Kenig and Martel in the Benjamin-Ono context , for different reasons. For powers of the nonlinearity which are definitely larger than , or just the three dimensional case for , we have a strong instability effect at the level of the previous spectral theory, probably associated to the dynamics around the soliton in the variable, and the virial identity seems no longer to hold.111In the one dimensional case, this instability condition does not appear, see Martel . Once again, a good understanding of the dynamics for powers close the the critical case or supercritical as in  needs a deep extension of (2.36) by incorporating now the dynamics in the variable, which could be very complicated, in view of some results by del Pino et al. . The extension of the ideas introduced by Martel  to any power of seems a very interesting problem.
Finally, we mention that another crucial application of the monotonicity formula on perturbed cones, needed in the higher dimensional case, is given in Lemma 4.4. Here, a new compact region of the plane is introduced, outside of which we prove exponential decay. This set is constructed in order to prove the strong convergence of sequences of bounded solutions, thanks to the use of the Sobolev embedding theorem.
One can also ask for the nonlinear dynamics in the remaining part of the plane, namely the region , see (1.17). We believe that in addition to radiation, one can find small solitons moving to the right in a very slow fashion. No finite energy solitary waves with speed along the direction are present, as shows the following (general) definition and result. As usual, we define the symbol by using its corresponding Fourier representation .
We say that , is the profile of a solitary wave of speeds if
is solution of (1.2).
Note that such a must satisfy the equation in
Assume that for some . Then (1.28) has no finite energy solutions.
We prove this result in Appendix D, using adapted Pohozahev identities.
Finally, as a consequence of the monotonicity properties associated to the linear part of the dynamics, in particular, using Lemma 2.1, we are able to prove the stability of the sum of essentially non-colliding solitons.
Let be an integer and . Consider solitons with scalings and centers , where . We say that these solitons are -decoupled if
that is, the solitons centers remains separated by a distance of at least for positive times. (See Fig. 3 below.)
-decoupled solitons can be characterized by a condition on the initial data only, at least up to a constant in : indeed, one can check that if, for all , we have either:
then the solitons are -decoupled.
Theorem 1.7 (Stability of the sum of decoupled solitons).
Assume . Consider a set of solitons of the form
where each is a fixed positive scaling, for all , and . Assume that the solitons are -decoupled, in the sense of Definition 1.6. Then there are , and depending on the previous parameters such that, for all , and for every , the following holds. Suppose that satisfies
Then there are fixed and defined for all such that , solution of (1.1) with initial data satisfies
The proof of this result is obtained by adapting the ideas by Martel, Merle and Tsai  for the generalized, one dimensional KdV case. Note that we do not need strictly well-prepared initial data as in . Instead, from (1.29) we just need sufficient well-separated solitons in the variable (no particular order at the beginning), and in the case where solitons have the same coordinate, we ask for well-ordered solitons to avoid multi-collisions.
The rest of the paper is organized as follows. The linear and nonlinear Liouville properties (Theorems 1.3 and 1.2) are proved respectively in Section 2 and 3. The nonlinear Liouville property is used to show Theorem 1.1 in Section 4. Section 5 is devoted to the proof of Theorem 1.7. Finally, in Appendix A, we present some numerical computations which establish the negativity of a scalar product in the case . Recall that this condition is a crucial element in the proofs of the rigidity results (Theorems 1.2 and 1.3). We also make an interesting observation about the plane wave solutions of the linear part of ZK in Appendix C and give the proof of Theorem 1.5 in Appendix D.
2. Linear Liouville property
Let denote a positive number such that . We define by
so that and . Note also that
The implicit constant appearing in (2.3) depends only on . In particular it does not depend on and .
To deal with the last term appearing on the right-hand side of (2.5), let us define
We claim that
where is the positive constant appearing in (1.7). Let to be fixed later. We consider the three following cases.
Case: . Then , so that
Case: . Then, we get that
since for .
Case: . In this case
On the other hand, it follows from (2.4) that