Bifurcations of Solitary Waves

The paper provides a brief review of the recent results devoted to bifurcations of solitary waves. The main attention is paid to the universality of soliton behavior and stability of solitons while approaching supercritical bifurcations. Near the transition point from supercritical to subcritical bi...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2008
Автори: Kuznetsov, E.A., Agafontsev, D.S., Dias, F.
Формат: Стаття
Мова:English
Опубліковано: Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України 2008
Назва видання:Журнал математической физики, анализа, геометрии
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/106521
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Bifurcations of Solitary Waves / E.A. Kuznetsov, D.S. Agafontsev, F. Dias // Журнал математической физики, анализа, геометрии. — 2008. — Т. 4, № 4. — С. 529-550. — Бібліогр.: 42 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-106521
record_format dspace
spelling irk-123456789-1065212016-09-30T03:03:05Z Bifurcations of Solitary Waves Kuznetsov, E.A. Agafontsev, D.S. Dias, F. The paper provides a brief review of the recent results devoted to bifurcations of solitary waves. The main attention is paid to the universality of soliton behavior and stability of solitons while approaching supercritical bifurcations. Near the transition point from supercritical to subcritical bifurcations, the stability of two families of solitons is studied in the frame-work of the generalized nonlinear Schrodinger equation. It is shown that one-dimensional solitons corresponding to the family of supercritical bifurcations are stable in the Lyapunov sense. The solitons from the subcritical bifurcation branch are unstable. The development of this instability results in the collapse of solitons. Near the time of collapse, the pulse amplitude and its width exhibit a self-similar behavior with a small asymmetry in the pulse tails due to self-steepening. 2008 Article Bifurcations of Solitary Waves / E.A. Kuznetsov, D.S. Agafontsev, F. Dias // Журнал математической физики, анализа, геометрии. — 2008. — Т. 4, № 4. — С. 529-550. — Бібліогр.: 42 назв. — англ. 1812-9471 http://dspace.nbuv.gov.ua/handle/123456789/106521 en Журнал математической физики, анализа, геометрии Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description The paper provides a brief review of the recent results devoted to bifurcations of solitary waves. The main attention is paid to the universality of soliton behavior and stability of solitons while approaching supercritical bifurcations. Near the transition point from supercritical to subcritical bifurcations, the stability of two families of solitons is studied in the frame-work of the generalized nonlinear Schrodinger equation. It is shown that one-dimensional solitons corresponding to the family of supercritical bifurcations are stable in the Lyapunov sense. The solitons from the subcritical bifurcation branch are unstable. The development of this instability results in the collapse of solitons. Near the time of collapse, the pulse amplitude and its width exhibit a self-similar behavior with a small asymmetry in the pulse tails due to self-steepening.
format Article
author Kuznetsov, E.A.
Agafontsev, D.S.
Dias, F.
spellingShingle Kuznetsov, E.A.
Agafontsev, D.S.
Dias, F.
Bifurcations of Solitary Waves
Журнал математической физики, анализа, геометрии
author_facet Kuznetsov, E.A.
Agafontsev, D.S.
Dias, F.
author_sort Kuznetsov, E.A.
title Bifurcations of Solitary Waves
title_short Bifurcations of Solitary Waves
title_full Bifurcations of Solitary Waves
title_fullStr Bifurcations of Solitary Waves
title_full_unstemmed Bifurcations of Solitary Waves
title_sort bifurcations of solitary waves
publisher Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
publishDate 2008
url http://dspace.nbuv.gov.ua/handle/123456789/106521
citation_txt Bifurcations of Solitary Waves / E.A. Kuznetsov, D.S. Agafontsev, F. Dias // Журнал математической физики, анализа, геометрии. — 2008. — Т. 4, № 4. — С. 529-550. — Бібліогр.: 42 назв. — англ.
series Журнал математической физики, анализа, геометрии
work_keys_str_mv AT kuznetsovea bifurcationsofsolitarywaves
AT agafontsevds bifurcationsofsolitarywaves
AT diasf bifurcationsofsolitarywaves
first_indexed 2025-07-07T18:36:11Z
last_indexed 2025-07-07T18:36:11Z
_version_ 1837014316132335616
fulltext Journal of Mathemati al Physi s, Analysis, Geometry2008, vol. 4, No. 4, pp. 529�550Bifur ations of Solitary WavesE.A. KuznetsovP.N. Lebedev Physi al Institute53 Leninsky Ave., Mos ow, 119991, RussiaE-mail:kuznetso�itp.a .ruD.S. AgafontsevL.D. Landau Institute for Theoreti al Physi s2 Kosygin Str., Mos ow, 119334, RussiaE-mail:dmitry�itp.a .ruF. DiasCMLA, ENS Ca han, CNRS, PRES UniverSud61 Av. President Wilson, F-94230 Ca han, Fran eRe eived June 25, 2008The paper provides a brief review of the re ent results devoted to bi-fur ations of solitary waves. The main attention is paid to the universalityof soliton behavior and stability of solitons while approa hing super riti- al bifur ations. Near the transition point from super riti al to sub riti albifur ations, the stability of two families of solitons is studied in the frame-work of the generalized nonlinear S hr�odinger equation. It is shown thatone-dimensional solitons orresponding to the family of super riti al bifur- ations are stable in the Lyapunov sense. The solitons from the sub riti albifur ation bran h are unstable. The development of this instability resultsin the ollapse of solitons. Near the time of ollapse, the pulse amplitudeand its width exhibit a self-similar behavior with a small asymmetry in thepulse tails due to self-steepening.Key words: stability, riti al regimes, wave ollapse.Mathemati s Subje t Classi� ation 2000: 37K50, 70K50.1. Introdu tionA ording to the usual de�nition, solitons are nonlinear lo alized obje ts pro-pagating uniformly with a onstant velo ity (see, for example, [1, 2℄). The solitonvelo ity V represents the main soliton hara teristi s whi h often de�nes the soli-ton shape, in parti ular its amplitude and width. E.A. Kuznetsov, D.S. Agafontsev, and F. Dias, 2008 E.A. Kuznetsov, D.S. Agafontsev, and F. DiasIt is well known that if the velo ity V of a moving obje t is su h that theequation !k = k �V (1)has a nontrivial solution where ! = !k is the dispersion law for linear wavesand k is the wave ve tor, then this obje t will lose energy due to Cherenkovradiation. This also pertains, to a large extent, to solitons as lo alized stationaryentities. They annot exist if the resonan e ondition (1) is satis�ed. Hen efollows the �rst, and simplest, sele tion rule for solitons: the soliton velo ity mustbe either less than the minimum phase velo ity of linear waves or greater thanthe maximum phase velo ity. The boundary separating the region of existen eof solitons from the resonan e region (1) determines the riti al soliton velo ityV r. As it is easily seen, this velo ity oin ides with the group velo ity of linearwaves at the tou hing point where the straight line ! = kV is tangent to thedispersion urve ! = !k (in the multidimensional ase � the point of tangen y ofthe plane ! = k �V to the dispersion surfa e). If tou hing o urs from below, thenthe riti al velo ity determines the maximum soliton velo ity for this parameterrange and, onversely, for tou hing from above V r oin ides with the minimumphase velo ity. Two regimes are possible in rossing this boundary orrespondingto super riti al or sub riti al bifur ations (soft or rigid ex itation regimes).While approa hing the super riti al bifur ation point from below or abovethe soliton amplitude vanishes smoothly a ording to the same � Landau � law(/ jV � V rj1=2) as for the phase transitions of the se ond kind (see, for instan e,[3℄). The behavior of solitons in this ase is ompletely universal, both for theiramplitudes and their shapes. As V ! V r solitons transform into os illating wavetrains with the arrying frequen y orresponding to the extremal phase velo ityof linear waves V r. The shape of the wave train envelope oin ides with thatfor the soliton of the standard � ubi � nonlinear S hr�odinger equation (NLS).The soliton width happens to be proportional to jV � V rj�1=2.Bifur ations of solitons were �rst observed for gravity- apillary waves in nu-meri al simulations by Longuet�Higgins [4℄ and explained later in [5�9℄. Thenthe bifur ation � a transition from periodi solutions to a soliton solution � wasstudied in [5℄ and [6℄ using normal forms. The stationary NLS for gravity- apillarywave solitons was derived in [8℄. In [10℄ it was shown that this me hanism an beextended to opti al solitons. In fa t, this paper provided the �rst demonstrationof the universality of soliton behavior near a super riti al bifur ation for waves ofarbitrary nature. It is worth noting that the universal hara ter of solitons allowsnot only to �nd their shapes, but also to study their stability. This analysis, asstated in [11℄, shows that near super riti al bifur ation the solitons are stable onlyin the one-dimensional ase.The question of whether the bifur ation is super riti al or sub riti al dependson the hara ter of nonlinear intera tion. The super riti al bifur ation o urs for530 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 Bifur ations of Solitary Wavesa fo using nonlinearity when the produ t !00T < 0, where !00 = �2!=�k2 is these ond derivative of the frequen y with respe t to the wave number, taken at thetou hing point k = k0, and T is the value of the matrix element Tk1k2k3k4 of thefour-wave intera tion for ki = k0. If !00T > 0, whi h orresponds to a defo usingnonlinearity, then there are no solitons � lo alized solutions � with amplitudevanishing gradually as V ! V r. In the theory of phase transitions this orre-sponds to a �rst-order phase transition, and in the theory of turbulen e, usingLandau's terminology [12℄, it orresponds to a rigid regime of ex itation. Thetransition through the riti al velo ity is a ompanied by a jump in the solitonamplitude. The magnitude of the jump is determined by the next higher-orderterms in the expansion of the Hamiltonian. Like for the �rst-order phase transi-tions, the universality of soliton behavior is no longer guaranteed in this situation.When the amplitude jump at this transition is small, it is enough to keep a �nitenumber of next order terms in the Hamiltonian expansion to des ribe su h a bifur- ation. In the phase transitions this orresponds to a �rst-order phase transition lose to a se ond-order transition, whi h o urs, for example, near a tri- riti alpoint. As shown in [13℄, this situation arises for one-dimensional internal-wavesolitons propagating along the interfa e between two ideal �uids with di�erentdensities in the presen e of both gravity and apillarity. A ording to [13℄ thematrix element T in this ase vanishes for density ratio �1=�2 = (21 � 8p5)=11.This type of bifur ations an also be met for gravity water waves with �nite depthwhen the matrix element T = 0 at � r = k0h � 1:363. In nonlinear opti s, asshown in [11℄, a de rease of T (Kerr onstant) an be provided by the intera tionof light pulses with a ousti waves (Mandelstamm�Brillouin s attering). If thejump in soliton amplitude is of order one, then we need to keep all the remainderterms in the Hamiltonian expansion.In this paper we give a brief review of the re ent results devoted to thissubje t. The main attention will be paid to the universality of soliton behaviorand stability of solitons while approa hing the super riti al bifur ation point.The paper is organized as follows. The next se tion is devoted to stationarysolitons for arbitrary nonlinear wave media and their properties near the super- riti al bifur ation. Se tion 2 deals with the stability of solitons based on theLyapunov theorem and the Hamiltonian approa h. It is shown by means of inte-gral estimates of Sobolev type in their multipli ative variant (Gagliardo�Nirenberginequalities) that only one-dimensional solitons are Lyapunov stable. It is worthnoting that, in ontrast to the method of normal forms, whi h is extensively usedin [5, 6, 13, 14℄ for studying bifur ations of solitons, the Hamiltonian approa h isfundamental for investigating soliton stability. In the method of normal forms, theintrodu tion of envelopes is not unique. Consequently, the Hamiltonian equationsof motion lose their initial Hamiltonian stru ture after their averaging. In thisse tion it is shown that near the bifur ation point the multi-dimensional solitonsJournal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 531 E.A. Kuznetsov, D.S. Agafontsev, and F. Diasare unstable due to the modulational instability. In the last se tion we onsiderwhi h nonlinear e�e ts must be taken into a ount near the transition from super- riti al to sub riti al bifur ations and how they hange the shape of solitons andtheir stability. 2. Super riti al Bifur ationsLet us onsider a purely onservative nonlinear wave medium whi h an bedes ribed by the HamiltonianH = Z !kjakj2dk+Hint; (2)where !k is the dispersion law of small-amplitude waves, ak are normal amplitudesof the waves, and the Hamiltonian Hint des ribes the nonlinear intera tion of thewaves.The equations of motion of the medium an be written in terms of the ampli-tudes ak in the standard manner�ak�t + i!kak = �iÆHintÆa�k ; (3)so that in the absen e of an intera tion the system onsists of a olle tion ofnonintera ting os illators (waves):ak(t) = ak(0)e�i!kt:Equation (3) des ribes the dynami s in the wave number spa e. To go ba k tothe physi al spa e one needs to perform the inverse Fourier transform (x; t) = 1(2�)d=2 Z ak(t)eik�rdk: (4)Originally, the fun tion (x; t) is related to the hara teristi s of the medium(�u tuations of the density and velo ity of the medium, ele tri and magneti �elds, and so on) by a linear transformation (see, for example, [15℄). It is impor-tant that if (x; t) is a periodi fun tion of the oordinates, then its spe trumak(t) onsists of a sum of Æ-fun tions. For lo alized distributions (x; t) ! 0 asjxj ! 1. The Fourier amplitude ak(t), being a lo alized fun tion of k, does not ontain Æ-fun tion singularities.Let us now onsider the solution of (3) in the form of a soliton propagatingwith the onstant velo ity V: (x; t) = (x�Vt):532 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 Bifur ations of Solitary WavesIn this ase the entire dependen e of ak on time t is ontained in the os illatingexponent: ak(t) = ke�ik�Vt;where by virtue of (3) the amplitude k will satisfy the equation(!k � k �V) k = �ÆHÆ �k � fk: (5)The di�eren e !k � k �V appearing in this equation will be positive for all k ifthe soliton velo ity is less than the minimum phase velo ityjVj < min(!k=k): (6)Conversely, the di�eren e will be negative for all k if the soliton velo ity is greaterthan the maximum phase velo ityjVj > max(!k=k): (7)We will show that a soliton solution is possible if the ondition (6) (or (7)) issatis�ed. Let us assume the opposite to be true � let the onditions (6) and (7)be violated, i.e., the equation (1) has a solution. For simpli ity, we will assumethat it is unique: k = k0. Then, sin e xÆ(x) = 0, the homogeneous linear equation(!k � k �V)Ck = 0possesses a nontrivial solution in the form of a mono hromati waveCk = AÆ(k � k0):In this ase (5) an be written (by virtue of the Fredholm alternative) k = AÆ(k� k0) + fk!k � k �V with fk0 = 0: (8)This equation, in ontrast to (5), ontains a free parameter � the omplexamplitude A. It an be solved, for example, by iterations, taking AÆ(k � k0) asthe zeroth term. It is important that be ause of the nonlinearity as a result ofiterations one will obtain multiple harmoni s with k = nk0 where n is integer.The solution will onsist of a olle tion of Æ-fun tions. Correspondingly, in phy-si al spa e the solution will be a periodi fun tion of the oordinates, i.e., it willbe nonlo alized. Hen e follows the �rst sele tion rule for solitons: the di�eren e!k � k �V must be sign-de�nite, whi h is equivalent to the requirements (6) or(7). In other words, it means the absen e of Cherenkov radiation.Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 533 E.A. Kuznetsov, D.S. Agafontsev, and F. DiasIn this entire s heme, however, there is an important ex eption. Havingrepresented (5) in the form (8), we have in fa t assumed that the singularity inthe expression fk!k � k �V (9)is nonremovable. This may not be the ase � the singularity in the denominatorin (9) ould be an elled with the numerator, i.e., it ould be removable [10℄.For example, this happens for the lassi soliton of the KdV equation, for equa-tions whi h are generalizations of the KdV equation [16℄ for a ombination of theone-dimensional NLS and the MKdV equations [10, 17℄ integrated by the sameZakharov�Shabat operator [18℄ and so on. In all of these ases an ellation o ursas a result of the k dependen e of the matrix elements. However, even in these ases, the sele tion rule for solitons remains the same after the resonan e (1) isremoved � the part remaining in the denominator must be sign-de�nite.In what follows the singularities in (9) are assumed to be nonremovable in theforbidden region, and we study the behavior of the soliton solution as the solitonvelo ity approa hes the riti al value. For de�niteness, it is assumed that theplane ! = k �V is tangent to the dispersion surfa e ! = !k from below, i.e., the riterion (6) holds. Let tou hing o ur at the point k = k0. Then, instead of (8),in the allowed region k = fk!k � k �V :As the velo ity V approa hes the riti al value V r, the denominator inthis expression be omes small near the tou hing point, so that k gets a sharppeak at this point k = �12!������ + k0(V r � V )��1 fk: (10)Here !�� = �2!=�k��k� is a symmetri , positive-de�nite, tensor of the se ondderivatives, evaluated at k = k0, and � = k� k0.It is evident from (10) that as V approa hes the riti al velo ity, the width ofthe peak narrows as pV r � V , and the distribution orresponding to the mainpeak k = k0 approa hes a mono hromati wave. A ounting for nonlinearity, thespe trum ontains harmoni s whi h are multiples of k = k0. If it is assumed thatthe amplitude of the soliton vanishes gradually as V ! V r (whi h would orre-spond to a se ond-order phase transition), then the solution (x) (or, equivalently, k) an be sought as an expansion in terms of harmoni s: (x0) = 1Xh=�1 n(X)eihk0�x0 ; x0 = x�Vt: (11) 534 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 Bifur ations of Solitary WavesHere the small parameter � =p1� V=V r (12)and the �slow� oordinate X = �x0 are formally introdu ed, so that n(X) is theamplitude of the envelope of n-th harmoni . The assumption that the solitonamplitude vanishes ontinuously at V = V r means that the leading term of theseries in (11) orresponds to the �rst harmoni , and all other harmoni s are smallin the parameter �. This is the ondition under whi h the nonlinear S hr�odingerequation is derived (see, for example, [10, 19, 2℄). In the present ase, we arriveat the stationary NLS�k0V r�2 1 + 12!�� �2 1�X��X� +Bj 1j2 1 = 0 (13)at leading order in �, where B is related to the matrix ~Tk1k2k3k4 of four-waveintera tions as B = �(2�)d ~Tk0k0k0k0 : (14)In this approximation the leading term in the intera tion Hamiltonian has theformHint = ~Tk0k0k0k02 Z �k1 �k2 k3 k4Æk1+k2�k3�k4dk1dk2dk3dk4 = �B2 Z j 1j4dx;(15)and the tilde means renormalization of the vertex due to the three-wave intera -tion � in the present ase the intera tion with the zeroth and se ond harmoni s.As we have already noted, !�� in (13) is a symmetri positive-de�nite tensor.For this reason, performing a rotation to its prin ipal axes and arrying out the orresponding extensions along ea h axis, (13) an be transformed into the stan-dard form ��2 +� � �j j2 = 0; (16)where � = sign( ~T!��). Hen e it follows, in the �rst pla e, that soli-tons are possible only if � is negative (fo using nonlinearity when the produ t~T!�� is negative) and, in the se ond pla e, that the amplitude of the solitons isproportional to � =p1� V=V r;i.e., the amplitude vanishes a ording to a square-root law, the size of the solitonin reases as 1=� as the velo ity approa hes the riti al value.In order to illustrate how this me hanism works onsider the simplest example,i.e., the time-dependent one-dimensional nonlinear S hr�odinger equationi� �t + xx + 2j j2 = 0: (17)Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 535 E.A. Kuznetsov, D.S. Agafontsev, and F. DiasAs is well known, this equation, unlike the general equation (3), has one additionalsymmetry, namely the gradient symmetry ! ei�. To �nd the orrespondingsolution one should put (x; t) = ei�t (x� V t), where obeys the equationL(�i�x) � iV x + � � xx = 2j j2 : (18)For the present ase, in a ordan e with (1), the ondition for Cherenkov radiationwill be written as follows: kV = (k) or L(k) = 0 ; (19)where the dispersion law for the equation (18) takes the form (k) = � + k2: (20)Hen e one an see that for � < 0 the resonan e ondition (19) is satis�ed for anyvalue of the velo ity V ! Consequently solitons do not exist in this ase. This anbe he ked dire tly by solving equation (18): for � < 0 all solutions are periodi or quasi-periodi . Soliton solutions are possible for positive �. Their velo itieslie in the range �2p� � V � 2p�. At the points k = �p� the dispersive urve = (k) tou hes the straight line = kV r. At these points the solution mustvanish in agreement with the general onsiderations. It dire tly follows from theexa t solution of (18): = ei�t eiV x0=2�k osh(�kx0) ; x0 = x� V t; �k =p� � V 2=4: (21)Solitons exist only for � > V 2=4: The upper boundary in this inequality de�nesthe riti al velo ity V r = 2p�:It is important to note also that for � > V 2=4 the operator L in the equation (18)is positive de�nite. 3. Stability of SolitonsTo in lude the time dependen e in the averaged equations the amplitudes n in the expansion (11) must be assumed to depend not only on the �slow� oordinate X but also on the slow time T = �2t. Then a multis ale expansiongives the nonstationary analog of the NLSEi t � �2 +� � �j j2 = 0 (22)instead of the stationary NLSE (16). The soliton stability problem for this equa-tion has been well studied (see, for example, [10℄ and [21℄). We re all the basi 536 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 Bifur ations of Solitary Wavespoints in the investigation of stability. The equation (22) as an equation forenvelopes inherits the anoni al Hamiltonian form (3)i� �t = Æ ~HÆ � ; (23)where the Hamiltonian~H = �2N + Z (jr j2 � j j4)dr; (� = �1); (24)arises as a result of averaging the initial Hamiltonian. The equation (22) preserves,besides ~H, the total number N of parti les (adiabati invariant), so that solitonsare stationary points of the energy fun tional E = ~H��2N with the �xed numberof parti les Æ(E + �2N) = 0: (25)The number of parti les (or intensity) in a soliton solution as a fun tion of � hasthe form Ns = Z j j2dx = �2�d Z jg(�)j2d�; (26)where d is the dimension of the spa e, and g(�) satis�es the equation�g +�g + jgj2g = 0:In the one-dimensional ase g = p2 se h � and, orrespondingly, Ns = 4�. In thetwo-dimensional ase Ns is independent of � for the entire family of solitons, whilein the three-dimensional ase Ns de reases with in reasing �. The dependen e ofNs on �2 is ru ial from the standpoint of soliton stability. It is obvious that themost dangerous disturban es will be those having wave numbers lose to k = k0moving together with the soliton, i.e., modulation-type disturban es.A ording to (25) the envelope solitons are stationary points of the energy Efor a �xed number of waves N . Therefore su h solutions will be stable in theLyapunov sense if they realize a minimum (or a maximum) of the energy for�xed N .Consider �rst the s aling transformations leaving N un hanged, (x) = 1ad=2 s �xa� ; (27)where s is the solitoni solution. The energy E under this transformationbe omes a fun tion of the s aling parameter a:E(a) = I1sa2 � I2sad ; (28)Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 537 E.A. Kuznetsov, D.S. Agafontsev, and F. Diaswhere I1s = R jr sj2dx, I2s = R j sj4dx and � = �1. Hen e it is easy to seethat in the one-dimensional ase the energy (28) is bounded from below and hasa minimum at a = 1 orresponding to the soliton solution. In this aseEs = �2�33 and 2I1s = I2s = 4�33 :The soliton also realizes a minimum of E with respe t to another simple trans-formation, i.e., the gauge one, 0(x)! 0(x) exp[i�(x)℄, whi h also preserves N ,E = Es + Z (�x)2 20dx:Thus, for d = 1 both simple transformations yield a minimum for the energy, thusindi ating soliton stability for the one-dimensional geometry.Now we give an exa t proof of this fa t. The ru ial point of this proof is basedon integral estimations of Sobolev type. These inequalities arise as the sequen esof general imbedding theorem. The (Sobolev) theorem says that a spa e Lp anbe imbedded into the Sobolev spa e W 12 if the spa e dimensionD < 2p(p+ 4):This means that between the normskukp = �Z jujpdDx�1=p ; p > 0;kukW 12 = �Z (�2juj2 + jruj2)dDx�1=2 ; �2 > 0;there exists the following inequality (see, e.g., [23℄):kukp �MkukW 12 ; (29)where M is some positive onstant. For D = 1 and p = 4, the inequality (29) is1Z�1 j j4dx �M1 24 1Z�1 (�2j j2 + j xj2)dx352 : (30)Here it is straightforward to get a multipli ative variant of the Sobolev inequality,the so- alled Gagliardo�Nirenberg inequality (GNI) [24℄ (see also [23, 25, 26℄).Use the s aling transform x!�x in (29). Instead of (30) we have1Z�1 j j4dx �M1 24�2 1Z�1 j j2dx � �+ 1Z�1 j xj2dx � 1�352 :538 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 Bifur ations of Solitary WavesThis inequality holds for any (positive) � in luding a minimal value for the r.h.s.The omputing of its minimum yields the GNI:I2�CN3=2I1=21 ; (31)where I1 = R j xj2dx, I2 = R j j4dx, and C is a new onstant. Then this inequa-lity an be improved by �nding the best (minimal possible) onstant C.To �nd Cbest onsider all extrema of the fun tionalJf g = I2N3=2I1=21 : (32)The latter problem redu es to the solution of stationary NLSE (16):��2 + xx + 2j j2 = 0:Hen e we �nd that the best onstant Cbest is a value of Jf g on the solitonsolution: Cbest = I2sN3=2I1=21s = 2I1=21sN3=2 : (33)This inequality allows us to obtain immediately a proof of 1D soliton stability.Substituting (31) into (24) results the following estimation for the energy:E � I1 � CbestI1=21 N3=2 = Es + (I1=21 � I1=21s )2: (34)The inequality be omes pre ise on the soliton solution, thus proving its stability.Note that this provides the stability of solitons not only with respe t to smallperturbations, but also against �nite perturbations.In the three-dimensional ase, in ontrary, the fun tion E(a) in (28) hasa maximum, orresponding to the soliton solution, and is unbounded from be-low as a! 0. The gauge transformation gives a minimum of E and therefore allsoliton solutions at d = 3 represent saddle points of the energy. That indi ates apossible instability of solitons in this ase.The fa t of (linear) instability of three-dimensional solitons follows from theVakhitov-Kolokolov riterion [20℄. It is as follows: if�Ns��2 > 0; (35)then the soliton is stable and, respe tively, unstable if this derivative is negative.This riterion has a simple physi al meaning. The value ��2 for the stationarynonlinear S hr�odinger equation (16) an be interpreted as the energy of the boundstate�soliton. If we add one �parti le� to the system and the energy of this boundJournal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 539 E.A. Kuznetsov, D.S. Agafontsev, and F. Diasstate de reases, then we have a stable situation. If by adding one �parti le�the level ��2 is pushed towards the ontinuous spe trum, then su h a soliton isunstable.At d = 3 the derivative �Ns=��2 < 0 and therefore 3D solitons are unsta-ble (the modulational instability). For the two-dimensional ase the Vakhitov�Kolokolov riterion (35) gives an absen e of linear exponential instability. A moredetailed analysis in this ase yields the power type instability (for details see thesurvey [21℄ and [22℄).Thus, the solitons are stable only in the one-dimensional ase, while in thetwo-dimensional ( riti al) and three-dimensional ases they are unstable and anbe onsidered as separatrix solutions separating ollapsing solutions from the dis-persive ones [27℄.This is probably the simplest method for explaining the well-known empiri alfa t that solitons, as a rule, exist only in one-dimensional systems. For multidi-mensional systems the stable solitons are rare and an only appear as a result oftopologi al onstraints or of a me hanism that removes Cherenkov singularities(whi h is dis ussed in the present paper). The latter, as an be easily understood,is due to the existen e of a ertain lass of symmetry.4. From Super riti al to Sub riti al Bifur ationsFor sub riti al bifur ation at the riti al velo ity the soliton undergoes a jumpin its amplitude. In this ase the orresponding theory an be developed near thetransition point between sub riti al and super riti al bifur ations (in analogy withthe tri- riti al point for phase transitions). In the series of papers [11, 13, 28�30℄it was shown that in this ase the soliton behavior an be des ribed by means ofthe generalized nonlinear S hr�odinger equation (NLSE) for the envelope , whi hin the one-dimensional ase is as follows:i� �t � �2!0 + !0002 xx � �j j2 + 4i�j j2 x + bkj j2 + 3Cj j4 = 0; (36)where !0 � !(k0) and k0 are the arrying frequen y and the wave number, res-pe tively, �2 = (V r � V )=V r � 1, !000 the se ond derivative of !(k) taken atk = k0. Here the four-wave oupling oe� ient � is assumed to have additionalsmallness hara terizing the proximity to the transition from super riti al to sub- riti al bifur ations. The transition point is de�ned from the equation � = 0.For example, for the interfa ial deep-water waves propagating along the interfa ebetween two ideal �uids in the presen e of apillarity [28, 29℄� = k301 + � �A2 r �A2� ; 540 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 Bifur ations of Solitary Waveswhere � is the density ratio, A = (1 � �)=(1 + �) the Atwood number, A2 r =5=16 and � r = (21 � 8p5)=11, as it was shown in [13℄. For � < � r, the four-wave oupling oe� ient � is negative, and the orresponding nonlinearity is ofthe fo using type. In this ase, the solitary waves near the riti al velo ity V rare des ribed by the stationary (�=�t = 0) NLSE and undergo a super riti albifur ation at V = V r [13℄. For � > � r, the oupling oe� ient hanges sign and,as a result, the bifur ation be omes sub riti al. For water waves in �nite depthh the oe� ient � hanges its sign at � r = k0h � 1:363 [31℄ while !000 is alwaysnegative. Thus the nonlinearity belongs to the fo using type for �(= kh) > � r andrespe tively be omes defo using in the region � < � r [31, 32℄. In nonlinear opti s,as shown in [11℄, a de rease of � (Kerr onstant) an be provided by the intera tionof light pulses with a ousti waves (Mandelstamm�Brillouin s attering).Be ause of smallness of � we keep in (36) the following order nonlinear terms:the gradient term (� �) responsible for self-steepening of the pulse (analog ofthe Lifshitz term in phase transitions), the nonlo al term (due to presen e ofthe integral operator bk, the Fourier transform of its kernel is equal to jkj) andthe six-wave nonlinear term with oupling oe� ient C. Two additional 4-waveintera tion terms, both lo al and nonlo al, appear as a result of expansion of thefour-wave matrix element Tk1k2k3k4 in powers of the small parameters �i = ki�k0:Tk1k2k3k4 = �2� + �2� (�1 + �2 + �3 + �4) (37)� 8� (j�1 � �3j+ j�2 � �3j+ j�2 � �4j+ j�1 � �4j):The existen e of nonlo al ontribution in the expansion is onne ted with non-analyti al dependen e of the matrix element T in its arguments. For interfa ialdeep-water waves (IW) this nonanalyti ity originates from the solution of Lapla eequation for the hydrodynami potential and its redu tion to the moving inter-fa e. For instan e, for water waves (WW) with a �nite depth the nonlo al termis absent [32℄ as well as for ele tromagneti waves in nonlinear diele tri s [11℄be ause of the analyti ity of matrix elements with respe t to frequen ies, whi his a onsequen e of ausality (see, for example, Refs. [11, 33℄). In the latter asethe spatial dispersion e�e ts are relativisti ally small and an be negle ted.For both IW and WW near the transition point, !000C is positive; moreover, is also positive for IW, and therefore the orresponding nonlinearities are fo using,thus providing the existen e of lo alized solutions. Depending on the sign of �there exist two bran hes of solitons. For IW they were found numeri ally [28, 29℄using the Petviashvili s heme [34℄. Expli it solutions for both kinds of IW solitons an be obtained in the limiting ase only when V ! V r. For negative � theseare the lassi al NLS solitons with a se h shape. For the sub riti al bifur ationat V = V r the soliton amplitude remains �nite with algebrai de ay (� 1=jxj)Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 541 E.A. Kuznetsov, D.S. Agafontsev, and F. Diasat in�nity [28, 29℄. When a nonlo al nonlinearity is absent ( = 0), the solitonsolutions an be found expli itly. For both bran hes at large � the number ofwaves N = R j j2dx approa hes from below and above the same value N r whi h oin ides with the number of waves N on the solitons with � = 0. For the solitonsin �bers this property means that the energy of opti al pulse saturates tendingto the onstant value with a de rease of the pulse duration.On the other hand, all solitons of (36) are stationary points of the energy Efor �xed number of waves : Æ(E + �2N) = 0, where the energy in dimensionlessvariables is given byE = Z �j xj2 + �2 j j4 + i�( �x � x �)j j2 � 2 j j2k̂j j2 � Cj j6� dx: (38)This allows one to use the Lyapunov theorem in the analysis of their stability.Here, for the IW, � � (V r � V )1=2j�� � rj�1, C = 319=1281 and� = sign(�� � r); � = 6=p427; = 32=p427; (39)and for the WW ase ( ompare with [32℄)� = sign(� r � �); � � �0:397; C � 0:176: (40)As shown in [11, 28, 29, 32℄, for N < N r solitons orresponding to the super riti- al bran h realize the minimum values of the energy and therefore they are stablein the Lyapunov sense, i.e., stable with respe t to not only small perturbationsbut also against �nite ones. In parti ular, the boundedness of E from below anbe viewed if one onsiders the s aling transformation = (1=a)1=2 s(x=a) retain-ing the number of waves N , where = s(x) is the soliton solution. Under thistransform E be omes a fun tion of the s aling parameter a:E (a) = �1a � 12a2� �2 Z j sj4 dx: (41)It is worth noting that the dispersion term and all nonlinear terms in E, ex eptR �2 j j4dx, have the same s aling dependen e / a�2. The latter means that at� = 0 (36) an be related to the riti al NLS equation like the two-dimensional ubi NLS equation. From (41) it is also seen that for � < 0 E(a) has a min-imum orresponding to the soliton. Unlike in the super riti al ase, the s alingtransformation for another soliton bran h with � > 0 gives a maximum of E(a)on solitons and unboundedness of E as a ! 0. Under the gauge transformation = sei�, on the ontrary, the energy rea hes a minimum on soliton solutionsand onsequently the solitons with � > 0 represent saddle points. This indi atesa possible instability of solitons for the whole sub riti al bran h, at least withrespe t to �nite perturbations.542 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 Bifur ations of Solitary WavesWe onsider this question in detail and emphasize a nonlinear stage of insta-bility following to our re ent paper [30℄. This problem, indeed, is not trivial inspite of a lose similarity with the riti al NLSE. It is worth noting that (36) at� = = C = 0 represents an integrable model (the so- alled derivative NLSE)[35℄, and exponentially de aying solitons in this model are stable. It is more orless evident also that small oe� ients ; C annot break the stability of soli-tons. This means that in the spa e of parameters we may expe t the existen e ofa threshold. Above this, the threshold solitons must be unstable and the deve-lopment of this instability may lead to ollapse, i.e., the formation of a singularityin �nite time.Consider the energy (38) written in terms of amplitude r and phase ' ( =rei'): E = Z �r2x + �2 r4 � 2 r2bkr2 � 13r6 + r2 �'x + �r2�2� dx; (42)where by an appropriate hoi e of the new dimensionless variables the renormal-ized onstant ~C = C + �2 an be taken equal to 1/3. Hen e one an see that theenergy takes its minimum value when the last term in (42) vanishes, i.e. when'x + �r2 = 0: (43)The integration of this equation gives an x�dependen e for the phase that is alled the hirp in nonlinear opti s. It is interesting to note that the remainingpart of the energy does not ontain the phase at all.First, study the lo al model when = 0. Let the energy be negative insome region : E < 0. Then, following [36, 26℄, one an establish that dueto radiation of small amplitude waves E < 0 an only de rease, be oming moreand more negative, but the maximum value of j j, a ording to the mean valuetheorem, an only in rease: maxx2 j j4 � 3jE jN : (44)This pro ess is possible only for the energies whi h are unbounded from below.In a ordan e with (41) su h a situation is realized when � > 0. In this ase theradiation leads to the appearan e of in�nitely large amplitudes r. However, it isimpossible to on lude that the singularity formation develops in �nite time.For > 0 the estimations on the maximum value of j j are not as transparentas they are for the lo al ase. Instead of (44), it is possible to obtain a similarestimate maxx j j4 � 3jEjN :However, it is expressed through the total energy E and the total number ofwaves N . Besides, two inequalities must be satis�ed: E < 0 and N < 2N2 .Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 543 E.A. Kuznetsov, D.S. Agafontsev, and F. DiasFor interfa ial waves, N2 � 1:39035 > N r � 1:3521. Thus, the maximum am-plitude in this ase is bounded from below by a onservative quantity and thismaximum an never disappear during the nonlinear evolution.Now we onsider the situation where the self-steepening pro ess anbe negle ted (� = 0). In this ase (36) be omesi t + xx � �2 � �j j2 + bkj j2 + 3Cj j4 = 0:It is possible to obtain a riterion of ollapse using the virial equation (for details,see [37, 36, 38℄). This equation is written for the positive de�nite quantityR = Z x2j j2dx;whi h, up to the multiplier N , oin ides with the mean square size of the distri-bution. The se ond derivative of R with respe t to time is de�ned by the virialequation Rtt = 8�E � �4 Z j j4dx� : (45)Hen e, for � > 0 one an easily obtain the following inequality:Rtt < 8E;whi h yields, after double integration, R < 4Et2+�1t+�2. Here �1;2 are onstantswhi h are obtained from the initial onditions. Hen e, it follows that for the stateswith negative energy, E < 0, there always exists su h a moment of time t0 whenthe positive de�nite quantity R vanishes. At this moment of time the amplitudebe omes in�nite. Therefore the ondition E < 0 represents a su� ient riterion of ollapse ( ompare with [37, 36℄). However, it is ne essary to add that this riterion an be improved in the same way as in [39, 27℄ for the three-dimensional ubi NLS equation. From (45) one an see that for the stationary ase (on the solitonsolution) Es = �4 R j sj4dx, in agreement with (41). As it was shown before, for� > 0 the soliton realizes a saddle point of E for �xed N . It follows from (41) thatfor small a the energy be omes unbounded from below, but for a > 1 it de reases(this orresponds to spreading). Therefore in order to a hieve a blow-up regimethe system should pass through the energeti barrier equal to Es. Thus, for this ase the riterion E < 0 must be hanged into the sharper riterion: E < Es.This riterion an be obtained rigorously using step by step the s heme presentedin [39, 27℄ and therefore we skip its derivation. 544 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 Bifur ations of Solitary Waves −30 −20 −10 0 10 20 30 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 xFig. 1: Initial (solid line) and �nal (dashed line) at t = 1:18 distributions for j j,interfa ial waves, selfsimilar variables. The soliton amplitude was in reased by1%, � = 1, � = 1. The ratio between �nal and initial soliton amplitudes in thephysi al variables is about 11.In order to verify all the theoreti al arguments about the formation of ollapsepresented above we performed a numeri al integration of the NLSE (36) for � > 0by using the standard 4th order Runge�Kutta s heme. The initial onditions were hosen in the form of solitons but with larger amplitudes than for the stationarysolitons. In reasing in the initial amplitude varied in the interval from 0.1% up to10%. The initial phase was given by means of (43). In all runs with theses initial onditions we observed a high in rease of the soliton amplitude up to a fa tor 14with a shrinking of its width. In a peak region the pulses for both IW and WW ases behaved similarly. Near the maximum the pulse peak was almost symmetri :anisotropy was not visible. The di�eren e was observed in the asymptoti regionsfar from the pulse ore, where the pulses had di�erent asymmetries for IW andWW be ause of the opposite sign for � (see (39), (40)). For the given values of �we did not observe the simultaneous formation of two types of singularities withblowing-up amplitudes and sharp gradients as it was demonstrated in the re entnumeri al experiments for the three-dimensional ollapse of short opti al pulsesdue to self-fo using and self-steepening in the framework of the generalized NLSequation [40℄ and equations of the Kadomtsev�Petvishvili type [41℄. Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 545 E.A. Kuznetsov, D.S. Agafontsev, and F. Dias −80 −60 −40 −20 0 20 40 60 80 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 xFig. 2: Initial (solid line) and �nal (dashed line) at t = 2:7192 distributions forj j, WW solitons, selfsimilar variables. The soliton amplitude was in reased by1%, � = 1, � = 1. The ratio between �nal and initial soliton amplitudes in thephysi al variables is about 11.In our numeri al omputations we found that the amplitude and its spatial ollapsing distribution developed in a selfsimilar manner. Near the ollapse pointin the equation (with � > 0) one an negle t the term proportional to �. In thisasymptoti regime (36) admits selfsimilar solutions,r(x; t) = (t0 � t)�1=4f � x(t0 � t)1=2� ; (46)where t0 is the ollapse time.To verify that we approa hed the asymptoti behavior given by (46), at ea hmoment of time we normalized the -fun tion by the maximum (in x) of itsmodulus max j j �M and introdu ed new self-similar variables, (x; t) =M (�; �); � =M2(x� xmax); � = lnM: (47)Here xmax is the point orresponding to the maximum of j j. In omparisonwith those given by (46), new variables are more onvenient be ause they do notrequire the determination of the ollapsing time t0.Fig. 1 and Fig. 2 show typi al dependen es of j j as a fun tion of the self-similar variable � at t = 0 (solid line) and at the �nal time (dashed line) for546 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 Bifur ations of Solitary Waves 1 1.02 1.04 1.06 1.08 1.1 1.12 1.14 1.16 1.18 1.2 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 tFig. 3: Dependen e of 1=max j j4 on time. Interfa ial waves.both the IW and WW ases. In both �gures one an see a fairly good oin iden ebetween the initial soliton distribution and the �nal one at the entral ( ollapsing)part of the pulse and asymmetry of the pulse at its tails due to self-steepening.The latter demonstrates that the ollapse has a selfsimilar behavior. The form ofthe entral part of the pulse approa hes the soliton shape be ause asymptoti allythe NLS model (36) tends to the riti al NLS system. It should be mentioned thatthis has been well known for the lassi al two-dimensional NLS equation sin e thepaper by Fraiman [42℄.Fig. 3 shows how 1=max j j4 depends on time. This dependen e is almostlinear in the orresponden e with the selfsimilar law (46). If the initial amplitudeswere less than the stationary soliton values, then the distribution would spreadin time dispersively, in full orresponden e with qualitative arguments based onthe s aling transformations (41).A knowledgments. The authors thank A.I. Dya henko for valuable dis us-sions on erning the numeri al simulations. We a knowledge support from CNRSunder the framework of PICS No. 4251 and RFBR under Grant 07-01-92165.The work of DA and EK was also supported by RFBR (Grant 06-01-00665),the Program of RAS "Fundamental problems in nonlinear dynami s" and GrantNSh 7550.2006.2. Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 547 E.A. Kuznetsov, D.S. Agafontsev, and F. DiasReferen es[1℄ V.E. Zakharov, S.V. Manakov, S.P. Novikov, and L.P. Pitaevskii, Teoriya solitonov.Nauka, Mos ow, 1980. (Russian) [In Engl.: Soliton Theory, Plenum Press, New York,London, 1984.℄[2℄ A.C. Newell, Solitons in Mathemati s and Physi s. SIAM, Philadelphia, 1985.[3℄ L.D. Landau and E.M. Lifshitz, Statisti al Physi s. Part 1. Pergamon Press, NewYork, 1980. [Russian orig.: Nauka, Mos ow, 1995, 521.℄[4℄ M.S. Longuet-Higgins, Capillary-Gravity Waves of Solitary Type on Deep Water.� J. Fluid Me h. 200 (1989), 451�470.[5℄ G. Iooss and K. Kir hg�assner, Bifur ation d'Ondes Solitaires en Pr�esen e d'uneFaible Tension super�Cielle. � C.R. A ad. S i. Paris, S�er. I. 311 (1990), 265�268.[6℄ J.-M. Vanden-Broe k and F. Dias, Gravity-Capillary Solitary Waves in Water ofIn�nite Depth and Related Free-Surfa e Flows. � J. Fluid Me h. 240 (1992), 549�557.[7℄ F. Dias and G. Iooss, Gravity-Capillary Solitary Waves with Damped Os illations.� Physi a D 65 (1993), 399�423.[8℄ M.S. Longuet-Higgins, Capillary-Gravity Waves of Solitary Type and EnvelopeSolitons on Deep Water. � J. Fluid Me h. 252 (1993), 703�711.[9℄ T.R. Akylas, Envelope Solitons with Stationary Crests. � Phys. Fluids 5 (1993),789�791.[10℄ V.E. Zakharov and E.A. Kuznetsov, Opti al Solitons and Quasisolitons. � Z. �Eksp.Teor. Fiz. 113 (1998), 1892. [JETP 86 (1998), 1035�1046.℄ (Russian)[11℄ E.A. Kuznetsov, Hard Soliton Ex itation: Stability Investigation. � Z. Eksp. Teor.Fiz. 116 (1999), 299. [JETP 89 (1999), 163.℄ (Russian)[12℄ L.D. Landau, On the Problem of a Turbulen e. � Dokl. Akad. Nauk USSR 44 (1944),339; L.D. Landau and E.M. Lifshitz, Fluid me hani s. 3rd Engl. Ed. PergamonPress, New York, 1986. [Russian orig., 3rd Ed.: Nauka, Mos ow, 1986.℄[13℄ F. Dias and G. Iooss, Capillary-Gravity Interfa ial Waves in Deep Water. � Eur.J. Me h. B/Fluids 15 (1996), 367�390.[14℄ G. Iooss, Existen e d'Orbites Homo lines �a un �Equilibre Elliptique, pour un Syst�emeR�eversible. � C.R. A ad. S i. Paris 324 (1997), 933�997.[15℄ V.E. Zakharov and E.A. Kuznetsov, Hamiltonian Formalism for Nonlinear Waves. �Usp. Fiz. Nauk 167 (1997), 1137�1168. (Russian) [Phys. Usp. 40 (1997), 1087�1116.℄[16℄ J. Ny ander, Steady Vorti es in Plasmas and Geophysi al Flows. � Chaos 4 (1994),253.[17℄ E.M. Gromov and V.I. Talanov, Nonlinear Dynami s of Short Pulses in DispersiveMedia. � Z. �Eksp. Teor. Fiz. 110 (1996), 137. [JETP 83 (1996), 73℄ (Russian);Short Opt. Solitons in Fibers. � Chaos 10 (2000) 551�558.548 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 Bifur ations of Solitary Waves[18℄ V.E. Zakharov and A.B. Shabat, Exa t Theory of Two-Dimensional Self-Fo usingand One-Dimensional Self-Modulation of Waves in Nonlinear Media. � Z. Eksp.Teor. Fiz. 61 (1996), 118�134. [Sov. Phys. JETP 34 (1972), 62�69.℄ (Russian)[19℄ V.E. Zakharov, Handbook of Plasma Physi s. 2. Basi Plasma Physi s II, (A. Galeevand R. Sudan, Eds.). North�Holland, Amsterdam, 1984.[20℄ N.G. Vakhitov and A.A. Kolokolov, Stationary Solutions of the Wave Equation inMedia with Saturated Nonlinearity. � Izv. VUZ. Radio�zika 16 (1973), 1020�1028.[Radiophys. Quant. Ele tron. 16 (1973), 783.℄[21℄ E.A. Kuznetsov, A.M. Ruben hik, and V.E. Zakharov, Soliton Stability in Plasmasand Fluids. � Phys. Rep. 142 (1986), 103.[22℄ E.A. Kuznetsov and S.K. Turitsyn, Talanov Transformations for Sel�o using Prob-lems and Instability of Waveguides. � Phys. Let. 112A (1985), 273.[23℄ O.A. Ladyzhenskaya, The Mathemati al Theory of Vis ous In ompressible Flow.Fizmatgiz, Mos ow, 1961. (Russian)[24℄ L. Nirenberg, An Extended Interpolation Inequality. � Ann. S i. Norm. Sup. Pisa20 (1966), No. 4, 733�737.[25℄ M.I. Weinstein, Nonlinear S hrodinger Equations and Sharp Interpolation Esti-mates. � Comm. Math. Phys. 87 (1983), 567�576.[26℄ E.A. Kuznetsov, Wave Collapse in Plasmas and Fluids. � Chaos 6 (1996), 381�390.[27℄ E.A. Kuznetsov, J.J. Rasmussen, K. Rypdal, and S.K. Turitsyn, Sharper Criteriaof the Wave Collapse. � Physi a D 87 (1995), 273�284.[28℄ D.S. Agafontsev, F. Dias, and E.A. Kuznetsov, Bifur ations and Stability of InternalSolitary Waves. � Pis'ma v ZhETF 85 (2006), 241. (Russian) [JETP Let. 83 (2006),201.℄[29℄ D.S. Agafontsev, F. Dias, and E.A. Kuznetsov, Deep-Water Internal Solitary WavesNear Criti al Density Ratio. � Physi a D 225 (2007), 153.[30℄ D.S. Agafontsev, F. Dias, and E.A. Kuznetsov, Collapse of Solitary Wavesnear Transition from Super riti al to Sub riti al Bifur ations. � Pis'ma ZhETF87 (2008), 767�771. (Russian) [JETP Let. 87 (2008), No. 12, 667�671.℄,arXiv:0805.1620v1.[31℄ G.B. Whitham, Nonlinear Dispersive Waves. � Pro . Rod. So . A. 283 (1965), 238.[Eur. J. Me h. B/Fluids 22 (1965), 273.℄[32℄ D.S. Agafontsev, Deep-Water Internal Solitary Waves near Criti al Density Ratio.� Pis'ma ZhETF 87 (2008), 225. (Russian)[33℄ L.D. Landau and E.M. Lifshitz, Ele trodynamika Sploshnykh Sred. Nauka, Mos ow,1982. (Russian) [In Engl.: Ele trodynami s of Continuous Media, 2-nd Ed. London,Pergamon, 1984.℄[34℄ V.I. Petviashvili, On the Equation for Unusual Soliton. � Fizika Plasmy 2 (1976),469. [Sov. J. Plasma Phys. 2 (1976), 247.℄Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4 549 E.A. Kuznetsov, D.S. Agafontsev, and F. Dias[35℄ D.J. Kaup and A.C. Newell, Stationary Solutions of the Wave Equation in theMedium with Saturated Nonlinearity. � J. Math. Phys. 19 (1978), Iss.4, 798�801.[36℄ V.E. Zakharov, Collapse of Langmure Waves. � Z. Eksp. Teor. Fiz. [Sov. Phys.JETP 35 (1972), 908.℄[37℄ S.N. Vlasov, V.A. Petrish hev, and V.I. Talanov, Averaged Des ription of WaveBeams in Linear and Nonlinear Media (the Method of Moments). � Izv. VUZ. Ra-dio�zika 14 (1971), 1353. (Russian) [Radiophys. Quant. Ele tron. 14 (1974), 1062.℄[38℄ V.E. Zakharov, Collapse and Self-Fo using of Langmuir Waves. Handbook ofPlasma Physi s. 2. Basi Plasma Physi s. (A.A. Galeev and R.N. Sudan, Eds.)Elsevier, North-Holland, 1984.[39℄ S.K. Turitsyn, Nonstable Solitons and Sharp Criteria for Wave Collapse. � Phys.Rev. E47 (1993), R1316.[40℄ N.A. Zharova, A.G. Litvak, and V.A. Mironov, Self-Fo using of Wave Pa kets andEnvelope Sho k Formation in Nonlinear Dispersive Media. � Z. Eksp. Teor. Fiz.130 (2006), 21. [JETP 103 (2006), 15.℄ (Russian)[41℄ A.A. Balakin, A.G. Litvak, V.A. Mironov, and S.A. Skobelev, Stru tural Features ofthe Self-A tion Dynami s of Ultrashort Ele tromagneti Pulses. � Zh. Eksp. Teor.Fiz. 131 (2007), 408. [JETP 104 (2007), 363.℄ (Russian)[42℄ G.M. Fraiman, Asymptoti Stability of Manifold of Self-Similar Solutions in Self-Fo using. � Z. Eksp. Teor. Fiz. 88 (1985), 390. (Russian) [Sov. Phys. JETP 61(1985), 228.℄ 550 Journal of Mathemati al Physi s, Analysis, Geometry, 2008, vol. 4, No. 4