arxiv: v1 [math.ap] 23 Apr 2008

Transcription

1 ON UNIVERSALITY OF CRITICAL BEHAVIOUR IN HAMILTONIAN PDES arxiv: v1 [math.ap] 3 Apr 008 BORIS DUBROVIN Dedicated to Sergei Petrovich Novikov on the occasion o his 70th birthday. Abstract. Our main goal is the comparative study o singularities o solutions to the systems o irst order quasilinear PDEs and their perturbations containing higher derivatives. The study is ocused on the subclass o Hamiltonian PDEs with one spatial dimension. For the systems o order one or two we describe the local structure o singularities o a generic solution to the unperturbed system near the point o gradient catastrophe in terms o standard objects o the classical singularity theory; we argue that their perturbed companions must be given by certain special solutions o Painlevé equations and their generalizations. Contents 1. Introduction. Nonlinear wave equation: Hamiltonian structure, irst integrals and solutions 6.1. Main examples o nonlinear wave equations 6.. Hamiltonian ormulation o nonlinear wave equation 7.3. Domain o hyperbolicity and Riemann invariants 9.4. Solutions o nonlinear wave equations Conserved quantities and ininitesimal symmetries o nonlinear wave equation Commutative Lie algebra o ininitesimal symmetries o nonlinear wave equation 1 3. Singularities o solutions to nonlinear wave equation Critical points o Type I and Whitney singularities Critical points o Type II and elliptic umbilic catastrophe Generalization to the multicomponent case 4. Perturbations o nonlinear wave equations Deormation problem-i: main deinitions and the programme o research Deormation problem-ii: extension o ininitesimal symmetries and D-operator 31 1

2 B. DUBROVIN 4.3. Deormation problem-iii: an integrability test Deormation problem-iv: on classiication o integrable perturbations Quasitriviality and solutions to the perturbed nonlinear wave equations Special Painlevé-type transcendents and Universality Conjectures 4 5. Main examples Boussinesq equation Toda hierarchy and its extension Nonlinear Schrödinger equation Ablowitz - Ladik equation 53 Reerences Introduction In [3] we developed a perturbative approach to classiication o systems o integrable Hamiltonian partial dierential equations PDEs). The systems under consideration are written in the orm o long wave expansion also called small dispersion expansion) 1.1) u t = Au) u x + B u; u x, u xx ) + B 3 u; u x, u xx, u xxx ) Here u = u 1 x, t),..., u n x, t) ) T is the unknown vector unction, the entries o the coeicient matrix Au) are smooth on some domain D R n ; the characteristic roots λ 1 u),..., λ n u) will be assumed to be pairiwise distinct 1.) det Au) λ 1) = 0, λ i u) λ j u) or i j, u D. The terms B, B 3,..., B k,... o the expansion 1 in the rhs o 1.1) are polynomials in the jet coordinates u x = u 1 x,..., u n x), uxx = u 1 xx,..., u n xx),... graded homogeneous o the degrees, 3,..., k,..., 1.3) deg B k u; ux,..., u k)) = k deg m u i x m = m, m > 0, deg ui = 0, i = 1,..., n. The coeicients o these polynomials are smooth unctions on the same domain D. 1 In order to uniormize the notations the leading term could be redenoted as Au) u x =: B 1 u; u x ). However we will not ollow this system o notations. We emphasize that the terms o degree zero thus, depending only on u but not on the jets) must not appear in the rhs o 1.1).

3 UNIVERSALITY IN HAMILTONIAN PDES 3 The system 1.1) can be considered as a perturbation o the irst order quasilinear system 1.4) u t = Au) u x when considering slow varying solutions, i.e., working with vector unctions ux) changing by 1 on the spatial scale o order L 1. The magnitude ux) will not be assumed to be small; the natural small parameter 1.5) ɛ = 1 L appears in the estimates or the derivatives 1.6) u x ɛ, u xx ɛ,..., u k) ɛ k,.... It is convenient to introduce slow variables by rescaling 1.7) x ɛ x, t ɛ t. The system 1.1) rewrites in the orm 1.8) u t = Au) u x + ɛ B u; u x, u xx ) + ɛ B 3 u; u x, u xx, u xxx ) In these notations the leading term 1.4) results rom the limit ɛ 0. For suiciently small ɛ one expects to see no major dierences in the behaviour o solutions to the perturbed and unperturbed equations 1.8) and 1.4) within the regions where the x-derivatives are bounded. However the dierences become quite serious near the critical point also called the point o gradient catastrophe) where derivatives o the solution to the unperturbed equation tend to ininity. Although the case o small viscosity perturbations has been well studied and understood see [5] and reerences therein), the critical behaviour o solutions to general conservative perturbations 1.8) to our best knowledge has not been investigated in a systematic way beore the author s paper [17] see the papers [9, 3, 36, 38, 39] or the study o various particular cases). The main goal o this paper is to extend the technique o [17] to Hamiltonian perturbations o systems o PDEs. A serious dierence rom the scalar case treated in [17] can be seen already in the case o second order systems considered in the present paper. For the scalar case any Hamiltonian perturbation preserves integrability up to a high order o the small parameter ɛ. For systems this is not the case: or generic perturbation o a second order quasilinear system integrability is destroyed already at the order ɛ. Due to this dierence we begin our study o critical behaviour in perturbed systems o Hamiltonian PDEs with considering integrable perturbations only. Our idea is that the local structure o the critical behaviour should depend only on the order ɛ 4 perturbation near a weak singularity o a hyperbolic system on the order ɛ perturbation near a weak singularity o an elliptic system. We postpone to a separate publication the study o singularities on the boundary between domains o ellipticity and hyperbolicity.) Because o this it suices to have integrability only at low orders Recall that any second order quasilinear PDE 1.4) can be integrated by the method o characteristics.

4 4 B. DUBROVIN o the perturbative expansion. We give a precise meaning to such an approximate integrability, explain the basic tools involved extension o ininitesimal symmetries, D-operator, quasitriviality and string equation), and ormulate the Universality Conjectures about the behaviour o a generic solution to the perturbed Hamiltonian system near the point o gradient catastrophe o the unperturbed solution. We argue that, up to shits, Galilean transormations and rescalings this behaviour essentially does not depend on the choice o solution provided certain genericity assumptions hold true). Moreover, this behaviour near the critical point is given by the ollowing ormulae: near a critical point within the domain o hyperbolicity the solution, ater an aine transormation o the independent variables x, t) x +, x ) and a nonlinear change o the dependent variables u, v) r +, r ) the solution has the ollowing asymptotic representation r + r c x + + α + ɛ 4/7 U a ɛ 6/7 x ; b ɛ 4/7 x + ) + O ɛ 6/7 ) 1.9) r r 0 + α ɛ /7 U a ɛ 6/7 x ; b ɛ 4/7 x + ) + O ɛ 4/7 ) where U = UX; T ) is a particular real smooth or all X R solution to the ourth order analogue o the Painlevé-I equation [ ) X = T U 6 U U + U U ) + 1 ] 40 U IV, U = du dx etc. depending on the real parameter T the so-called PI equation). Here a, b, c, α ± are some constants that depend on the choice o the equation and the solution. For the original variables u, v the above asymptotic behaviour implies u u c + λ ɛ /7 U ) a ɛ 6/7 x ; b ɛ 4/7 x ) + + O ɛ 4/7 1.11) with some constants λ, µ. v v c + µ ɛ /7 U a ɛ 6/7 x ; b ɛ 4/7 x + ) + O ɛ 4/7 ) The solution UX, T ) was completely characterized by T.Claeys and M.Vanlessen in [9] via the isomonodromy deormations technique see [33]) where a long standing problem o existence o a pole-ree real solution to PI has been settled3. In the very recent paper [8] by T.Claeys and T.Grava the Universality Conjecture ormulated in [17] was rigorously established, with the help o Deit - Zhou steepest descent method [13], or the case o critical behaviour o solutions to the KdV equation with analytic initial data. For the systems the somewhat weaker asymptotics 1.11) was proved in [10] or a particular class o solution to the hierarchy o Toda equations with analytic initial data. These are the solutions appearing in the study o Hermitean matrix integrals. They depend on one arbitrary analytic unction that determines the probability 3 The equation 1.10) appeared in [6] or the particular value o the parameter T = 0) in the study o the double scaling limit or the matrix model with the multicritical index m = 3. It was observed that the generic solution to 1.10) blows up at some point o real line; the conjecture about existence o a unique smooth solution has been ormulated. Our original conjecture ormulated in [17] says that the solution to 1.10) smooth on the entire real axis exists and is unique or any real T. The uniqueness part o this conjecture remains open.

5 UNIVERSALITY IN HAMILTONIAN PDES 5 distribution; the general solution to Toda hierarchy equations should depend on two arbitrary unctions. The critical behaviour 1.11) appears, according to [10], in the large N asymptotic behaviour o orthogonal polynomials near the so-called singular edge points. The small dispersion parameter in this case is ɛ = 1/N. near a critical point within the domain o ellipticity a suitable complex combination w = wu, v) o the solution ux, t), vx, t) as unction o an aine complex combination z o the independent variables behaves as ollows w w 0 + α ɛ /5 W 0 ɛ 4/5 z ) + O ɛ 4/5) 1.1) z = a x + b t + z 0, a, b C. Here W 0 = W 0 Z) is the so-called tritronquée solution to the Painlevé-I equation P I ) W = 6 W Z. This solution discovered by Boutroux in [4] is uniquely characterized by the ollowing asymptotic behaviour Z W Z) 6, Z, arg Z < 4π 5. In particular it has no poles or suiciently large Z in the above sector. In [18] it was conjectured that the tritronquée solution has no poles inside the sector 1.13) arg Z < 4π 5. Although a rigorous proo o this conjecture remains an open problem, numerical evidences ound in [18] support the conjectur. Importance o absence o poles in the sector 1.13) has the ollowing explanation. It turns out that the complex parameters a and b in 1.1) are such that or suiciently small t and ɛ > 0 the argument Z = a x + b t + z 0 ɛ 4/5 o the tritronquée solution belongs to the sector 1.13) or all real x. So the Universality Conjecture 1.1) makes sense only i no poles occur within the sector 1.13). The Universality Conjecture 1.1) appeared irst in the recent paper [18] or the case o ocusing nonlinear Schrödinger NLS) equation. In this paper we will explain that the same critical behaviour arises also or other Hamiltonian perturbations o second order quasilinear PDEs o elliptic type. The paper is organized as ollows. In Section we recall some basic acts about second order quasilinear PDEs, their Hamiltonian structures, solutions, conserved quantities and symmetries. In Section 3 we give a very simple local classiication o weak singularities o solutions, or both hyperbolic and elliptic cases. An attempt to extend this local classiication o singularities to quasilinear hyperbolic systems o order higher than is also given; only the case o integrable hyperbolic PDEs the so-called semihamiltonian systems in the terminology o [48]), or systèmes rich according to [45] is

6 6 B. DUBROVIN considered. In Section 4 the main tools or the study o Hamiltonian perturbations are developed. Applications o these techniques to the study o concrete examples is given in Section 5. Acknowledgments. This work is partially supported by European Science Foundation Programme Methods o Integrable Systems, Geometry, Applied Mathematics MISGAM), Marie Curie RTN European Network in Geometry, Mathematical Physics and Applications ENIGMA), and by Italian Ministry o Universities and Researches MUR) research grant PRIN 006 Geometric methods in the theory o nonlinear waves and their applications.. Nonlinear wave equation: Hamiltonian structure, irst integrals and solutions We now proceed to considering the two-component systems o Hamiltonian PDEs. In this section we start with the irst order quasilinear systems 1.4). Moreover, we restrict ourselves with considering an important subclass o such systems associated with nonlinear wave equation.1) u tt xp u) = 0 or a given smooth unction P u). The equation.1) is linear or a quadratic unction P u); we assume thereore that P u) Main examples o nonlinear wave equations. The equation.1) o course is o interest per se; it also arises in the study o dispersionless limits o various PDEs o higher order. In particular, or P u) = 1 6 u3 one obtains the dispersionless limit.) u tt + u u x ) x = 0 o Boussinesq equation.3) u tt + u u x ) x + u xxxx = 0. For P u) = e u.1) reduces to the long-wave limit.4) u tt = xe u o Toda equations.5) ü n = e u n+1 u n e un u n 1 More generally, one-dimensional system o particles with neighboring interaction.6) H = 1 p n + P q n q n 1 ) with the potential P u) generalized Fermi - Pasta - Ulam system) yields, ater interpolation.7) q n t) q n 1 t) = wɛ n, ɛ t) p n t) = vɛ n, ɛ t)

7 UNIVERSALITY IN HAMILTONIAN PDES 7 and substitution.8) u = the ollowing system ɛ x w 1 e ɛ x.9) u t = v x [ e v t = ɛ 1 P ɛ x 1 ɛ x ) u 1 e P ɛ x )] u ɛ x = x P u) + ɛ [ ] P u) u xxx + 4 P u) u x u xx + P IV u) u 3 x + Oɛ 4 ). 4 The above ormulae are understood as ormal power series in ɛ: w = eɛ x 1 ɛ x u = 1 ɛ x+ɛ x us) ds = u + k 1 ɛ k k + 1)! uk),.10) u = w + 1 ɛ w + k>1 B k k! ɛk w k), B k are the Bernoulli numbers. In the limit ɛ 0.9) reduces to the nonlinear wave equation written as a system.11) u t = v x v t = x P u). Further examples will be considered in Section 5... Hamiltonian ormulation o nonlinear wave equation. The system.11) can be written in the Hamiltonian orm δh u t = x δvx).1) δh v t = x δux) with the Hamiltonian.13) H = The associated Poisson bracket reads [ ] 1 v + P u) dx..14) {ux), vy)} = δ x y) The Poisson bracket o two local unctionals F = u, v; u x, v x,... ) dx, G = gu, v; u x, v x,... ) dx

8 8 B. DUBROVIN is again a local unctional given by [ δf.15) {F, G} = δux) δg x δvx) + δf ] δvx) δg x dx δux) where δf δux) = E u, δf δvx) = E v, δg δux) = E ug δg δvx) = E vg are the variational derivatives,.16) E u := u x + x... u x u xx E v := v x + x... v x v xx are the Euler - Lagrange operators. Recall that the densities o the local unctionals are considered modulo total x- derivatives. So, the Poisson bracket.15) vanishes i the integrand in the rhs is a total x-derivative. The ollowing important criterion see, e.g., [11]) gives a necessary and suicient condition or a dierential polynomial to be a total x-derivative o another dierential polynomial. Theorem.1. Let h = M A i,j u, v) i 1 x u... i k x u j 1 x v... j l x v k,l 0 i,j be a dierential polynomial with coeicients A i,j u, v) smooth on a disk D R. Here i = i 1,..., i k ), j = j 1,..., j l ), 1 i 1,..., i k, j 1,..., j l N are multiindices o the summation, M, N are two nonnegative integers. Then there exists another dierential polynomial h o the same orm such that h = x h, x = ) u m+1) + vm+1) u m) v m) m 0 i δh δux) = 0, δh δvx) = 0, H = h dx.

9 UNIVERSALITY IN HAMILTONIAN PDES 9 Corollary.. Two local unctionals F = dx and G = g dx commute with respect to the Poisson bracket.14) i δf E u δux) δg x δvx) + δf ) δvx) δg x = 0 δux).17) δf E v δux) δg x δvx) + δf ) δvx) δg x = 0. δux) Alternatively.1) can be recast into the orm o Euler - Lagrange equations [ ] 1 δs = 0, S = φ t V φ x ) dx dt with by introducing a potential The canonical Poisson bracket yields.14). P u) = V u) φ x = u. {φx), φ t y)} = δx y).3. Domain o hyperbolicity and Riemann invariants. From.11) it readily ollows that the nonlinear wave equation is hyperbolic on the domain o convexity o P u),.18) u, v) R such that P u) > 0 and elliptic when P u) becomes concave. On the domain o hyperbolicity.11) can be reduced to the diagonal orm.19) t r ± = ± P u) x r ± introducing Riemann invariants.0) r ± = v ± Qu), where Q u) = P u). We also remind the ollowing property o solutions written in the Riemann invariants: introducing the characteristic variables.1) x ± = x ± P u)t the unction r + /r does not depend on x /x +,.) r + x = r x + = 0.

10 10 B. DUBROVIN.4. Solutions o nonlinear wave equations. Let us irst describe the local structure o general solution to the system.11). Lemma.3. Let = u, v) be an arbitrary solution to the linear PDE.3) uu = P u) vv. Let u 0, v 0 ) be a point such that.4) [ uv u 0, v 0 )] P u 0 ) [ vv u 0, v 0 )] 0. Then the system o equations.5) x = u u, v) t = v u, v) deines a unique pair o unctions ux, t), vx, t)) near the point such that x 0, t 0 ) := u u 0, v 0 ), v u 0, v 0 )) ux 0, t 0 ) = u 0, vx 0, t 0 ) = v 0. The unctions ux, t), vx, t)) solve.11). Their x-derivatives at the point x 0, t 0 ) do not vanish simultaneously. Conversely, any solution to.11) satisying the condition [.6) v x P u)ux] 0 x 0,t 0 ) can be locally obtained by the above procedure starting rom a suitable solution u, v) to the linear PDE.3). This is a reormulation o the classical method o characteristics also called hodograph transormation). Remark.4. The Riemann invariants.0) are also the characteristic variables or the linear PDE.3), i.e., in these variables the PDE or the unction = r +, r ) reads.7) r + r 1 8 P ) = 0. P ) 3/ r + r Remark.5. The exceptional solutions to.11) that cannot be written in the orm.5) can be represented in the orm v = Qu) + c or v = Qu) + c or a constant c where u = ux, t) satisies the scalar quasilinear PDE u t = P u) u x or u t = P u) u x. On these solutions one o the Riemann invariants takes constant values.

11 UNIVERSALITY IN HAMILTONIAN PDES Conserved quantities and ininitesimal symmetries o nonlinear wave equation. An alternative interpretation o solutions to the linear PDE.3) is given by Lemma.6. Given a unction = u, v), consider the Hamiltonian.8) H = u, v) dx. This Hamiltonian commutes with the Hamiltonian.13) o nonlinear wave equation i the unction satisies.3). Proo. The Poisson bracket.14) o H and H reads [ δh {H, H} = δux) δh x δvx) + δh ] δvx) δh x dx = δux) The integrand is a total x-derivative i u u ) = v v P u)). This gives.3). [ u v x + v P u)u x ] dx. Recall that the dierential polynomial = u, v; u x, v x,... ) is called conserved quantity or density o a irst integral) o the PDE.11) i the ollowing identity holds.9) t = x g, t = [ ] u t ) m) u + v t) m) m) v m) m 0 or some dierential polynomial g. It is understood that in.9) one has to substitute rom.11) u t v x, v t P u)u x. Corollary.7. The unction = u, v) is a conserved quantity or the nonlinear wave equation.11) i it satisies the linear PDE.3). Proo. For a solution to.3) one has t = x g where g = gu, v) is determined by a quadrature The converse statement is straightorward. Corollary.8. The Hamiltonian PDE.30) commutes with the wave equation dg = P u) v du + u dv. u s = x v u, v) v s = x u u, v) u t ) s = u s ) t, v t ) s = v s ) t i the unction u, v) satisies the linear PDE.3).

12 1 B. DUBROVIN Proo. Indeed, the low.30) is generated by the Hamiltonian H : u s = x δh δvx) v s = x δh δux). As it is well known rom Hamiltonian mechanics, commutativity o Hamiltonians implies commutativity o the lows..6. Commutative Lie algebra o ininitesimal symmetries o nonlinear wave equation. We conclude that the solutions to the linear PDE.3) are associated with ininitesimal symmetries o the nonlinear wave equation.1). A somewhat stronger statement says that all Hamiltonians o the orm H commute pairwise, i.e., the Lie algebra o symmetries is commutative. Lemma.9. Given two arbitrary solutions, g to the linear PDE.3), the Hamiltonians H = u, v) dx and H g = gu, v) dx commute: {H, H g } = 0. We leave the proo o this lemma as an exercise or the reader. For a given P u) one obtains thereore an ininite amily o commuting Hamiltonians parametrized by solutions to the linear PDE.3). The associated Hamiltonian systems commute pairwise. Remark.10. All equations.30) o the commutative amily have the same Riemann invariants.0):.31) s r ± = h uv ± P u)h vv ) x r ±. Example.11. Given a real number κ 0, ±1, choose.3) P u) = The unction.33) = 1 u v + u κ κκ 1). u κ κκ + 1) satisies the linear equation.3). The associated Hamiltonian low reads.34) u t = x u v) ) v v t = x + uκ. κ

13 UNIVERSALITY IN HAMILTONIAN PDES 13 These equations coincide ater the change o sign o the time variable) with equations o motion o one-dimensional isentropic gas with the equation o state p ρ κ+1. Here u = ρ the mass density), v the velocity. Thereore the equations.34) o gas dynamic can be considered as an ininitesimal symmetry o the nonlinear wave equation u tt = uκ 1 κ 1 u xx + u κ u x. Example.1. For the exceptional value κ = 1 take.35) P u) = u log u 1). The choice.36) = 1 u u v ) o a solution to the linear PDE.3) produces the dispersionless limit.37) u t + u v) x = 0 v t + v v x u x = 0 o the ocusing/deocusing nonlinear Schrödinger equation.38) i ψ t + 1 ψ xx ± ψ ψ = 0 written in the coordinates i.e., u = ψ, v = 1 i u t + u v) x = 0 v t + v v x u x = 1 4 ψx ψ ψ ), x ψ uxx u 1 ) u x. u x So, the dispersionless NLS is an ininitesimal symmetry o the ollowing nonlinear wave equation ) 1 u tt ± x = 0. u Remark.13. In [14] a hierarchy o the irst order commuting quasilinear PDEs has been associated with an arbitrary solution F = F v 1,..., v n ) o equations o associativity.39) 3 F 3 F v i v j v p ηpq = v q v k v l 3 F v 1 v i v j = η ij. 3 F 3 F v l v j v p ηpq v q v k v, i i, j, k, l = 1,..., n

14 14 B. DUBROVIN Here η ij are entries o a constant symmetric nondegenerate matrix, η ij are entries o the inverse matrix. The quasihomogeneity axiom has not been used or the construction o the hierarchy. For n = the solution F = F v, u) must have the ollowing orm or some unction Φu). For F = 1 v u + Φu) Φ u) = P u) this amily o commuting PDEs coincides with the one associated with solutions to.3). The Hamiltonians o the hierarchy o [14] are solutions to.3) polynomial in v. Beore proceeding to the discussion o the critical phenomena let us extend the method o characteristics described in Lemma.3 to solving any equation o the commutative amily.40) u s = x h v with the Hamiltonian H h = v s = x h u hu, v) dx, h uu = P u)h vv. Lemma.14. Let = u, v) be any solution to.3) satisying the condition ) 0 v h 0 uvv uvh 0 0 vv P u 0 ) ) v 0 h 0 vvv vvh 0 0 vv 0 or some u 0, v 0 ). Then the system.41) x + s h uv u, v) = u u, v) s h vv u, v) = v u, v) ux 0, s 0 ) = u 0, vx 0, s 0 ) = v 0 deines a unique pair o unctions ux, s), vx, s)) near x 0, s 0 ) determined rom x 0 + s 0 h uv u 0, v 0 ) = u u 0, v 0 ) s 0 h vv u 0, v 0 ) = v u 0, v 0 ). These unctions solve the system.40). Conversely, any solution to.40) satisying certain nondegeneracy conditions can be obtained by this procedure.

15 UNIVERSALITY IN HAMILTONIAN PDES Singularities o solutions to nonlinear wave equation In this section we develop a local classiication o singularities o solutions to the nonlinear wave equation. Actually, the technique based on the implicit unction theorem works equally well also or any PDE o the orm.40) rom the commuting amily o ininitesimal symmetries o the wave equation. So, our study o singularities will be developed or the solutions to the general system.40). The unctions ux, s), vx, s)) determined by the system.41) are smooth provided the conditions o the implicit unction theorem hold true. At the points x c, s c ) where the implicit unction theorem ails to be applicable the solution to.40) has a weak singularity, i.e., there exist the limits lim ux, s) = u c, x x c, s s c lim x x c, s s c vx, s) = v c but the derivatives u x, v x blow up. The main goal o this Section is to describe the local behaviour o a solution near the generic weak singularity. The precise deinition o genericity will be given below. Let us irst subdivide weak singularities into three classes. Type I. The limiting value u c, v c ) belongs to the domain o hyperbolicity o the system.40), i.e. 3.1) P u c ) > 0. Type II. The limiting value u c, v c ) belongs to the domain o ellipticity o the system.40), i.e. 3.) P u c ) < 0. Type III. The limiting value u c, v c ) belongs to the boundary between the domains o hyperbolicity and ellipticity, 3.3) P u c ) = 0. Observe that the linear PDE.3) belongs to the hyperbolic/elliptic type near the critical point u c, v c ) o Type I/II, while the type o this PDE changes near a generic i.e., P u c ) 0) critical point o Type III. In this Section we will study the local behaviour o the solution near the weak singularity ollowing the above classiication. Only Type I and Type II singularities will be considered; the study o Type III singularities is postponed to a subsequent publication. Let us do irst some general simpliication o notations. First, since the partial derivative h v satisies the linear PDE.3) i the unction h does so, it suices to study singularities o implicit unctions ux, s) and vx, s) determined by the system 3.4) x + s h u u, v) = u u, v) s h v u, v) = v u, v)

16 16 B. DUBROVIN where both unctions h and satisy.3). Next, the shits can be absorbed by a redeinition x x x c, s s s c + x c u + s c h o the unction. So, without loss o generality we may assume that x c = s c = 0. We will also adopt the ollowing system o notations: the values o the unctions u, v), hu, v) and o their derivatives u u, v), v u, v) etc. at the critical point u c, v c ) will be denoted 0, h 0, u, 0 v 0 etc. We will also denote etc. P 0 := P u c ), P 0 := P u c ) According to the above conventions we will have 3.5) 0 u = 0 v = 0. Violation o the implicit unction theorem conditions in our notations reads ) uu s h 3.6) det uu uv s h uv = ) uv s h uv vv s h uv 0 ) P 0 0 vv = 0. vv x=0, s= Critical points o Type I and Whitney singularities. As it ollows rom the representation.5) o solutions, the inverse map to x, t) ux, t), vx, t)) is well deined also on a neighborhood o the critical point. Moreover, it is given by the gradient 3.7) u, v) u u, v), v u, v)) = x, t) o a unction satisying the linear PDE.3) that, in the case o Type I singularity at u c, v c ) belongs to the hyperbolic type. The points o blow-up o the derivatives o the solution ux, t), vx, t)) are exactly the points where the Jacobian o the gradient map 3.7) vanishes. Classiication o singularities o smooth maps o plane to plane obtained in the celebrated paper o H.Whitney [49] was one o the irst successes o the singularity theory. Recall that Whitney obtained the complete list o normal orms o smooth maps R, 0) R, 0) u, v) x, t) stable with respect to small perturbations. He proved that, up to dieomorphisms x, t) x, t) u, v) ũ, ṽ) these maps belong to one o the ollowing three normal orms:

17 UNIVERSALITY IN HAMILTONIAN PDES 17 Type W 1 3.8) Type W 3.9) x = ṽ t = ũ x = ṽ t = ũ Type W ) x = ũ ṽ 1 6ṽ3 t = ũ For the W 1 case the inverse unctions u = ux, t), v = vx, t) are well deined and smooth. We will show below that also in the general situation the solutions to the quasilinear system.40) near a generic Type I singularity locally has the universal shape 3.9). Another natural question is about the structure o solutions near the irst point o singularity. We say that x c, t c ) is the irst singularity o the solution to the nonlinear wave equation i the unctions ux, t) and vx, t) are smooth or suiciently small x x c and or t < t c. We will see that the local behaviour o a generic solution near a irst point singularity, within the hyperbolicity domain, is described by the W 3 singularity on the Whitney list. This sounds like being almost obvious; an additional point is that only aine transormation o independent variables x, t will be allowed near the critical point. Let proceed now to the study o singularities o solutions to the general system.40). Due to the assumption P 0 > 0 the jacobian 3.6) actorizes ) ) 0 uv P 0. 0 vv ) = 0 uv + P 0 0 vv ) 0 uv P 0 0 vv The irst o the genericity assumption says that, at the critical point only one o the actors vanishes. Without loss o generality we may assume that 3.11) 0 uv P 0 0 vv = 0. Deinition 3.1. We say that the critical point 3.11) o Type 1 is generic i: 3.1) 0 uv + P 0 0 vv 0 and 3.13) uvv 0 P 0 vvv 0 + vv 0 P P 0

18 18 B. DUBROVIN or but 3.14) uvv 0 P 0 vvv 0 + vv 0 P 0 4 P 0 A 0 := h 0 uv P 0 h 0 vv + h 0 P 0 v 4 P 0 0 B 0 := uvvv 0 P 0 vvvv 0 + vvv 0 P 0 4 P and 0 vv = 0 [ 5P P 0 ) 5/ 0 ) 4P IV 0 P 0 ] 0. The ormula 3.15) x r ± = 1 0 uv ± P 0 0 vv s h 0 uv ± ) P 0 h 0 vv or the x-derivatives o Riemann invariants clariies the meaning o the irst genericity assumption 3.1): at the critical point x c, s c ) = 0, 0) only one Riemann invariant blows up under our assumptions r does so) while another one remains smooth. We now show that the case 3.13) is incompatible with the assumption that x c, s c ) = 0, 0) is locally the irst critical point. Lemma 3.. Let x c, s c ) = 0, 0) be a critical point o Type I satisying 3.11), 3.1) such that the unctions ux, s), vx, s) determined by 3.4) along with normalization ux c, s c ) = u c, vx c, s c ) = v c are smooth or s < s c or suiciently small x x c, s s c. Then 3.16) uvv 0 P 0 vvv 0 + vv 0 P 0 4 P 0 The meaning o the constraint 3.16) is the ollowing: the graph o the Riemann invariant r x, s c ) must have an inlection point at x = x c. = 0. Proo. For simplicity o calculations we will assume that h 0 v 0. Under this assumption one can solve the system 3.4) in the orm s = su, v) = v h v, x = xu, v) = uh v v h u h v. Expanding these unctions in Taylor series in r +, r at u = u c, v = v c we obtain the ollowing expressions or the linear combinations ) x ± = x + h 0 u ± ) P 0 h 0 v s 4 Clearly these correspond to the characteristics o the hyperbolic PDE.

19 UNIVERSALITY IN HAMILTONIAN PDES 19 o independent variables 3.18) x + = vv 0 P [ h 0 v + 1 h 0 v x = 0 vv + 1 P 0 8P 0 0 r + 0 uvv + P 0 0 vvv r + 0 vv h 0 v 0 uvv P 0 0 vvv + 0 vv where we introduced the notations ) vv 0 h 0 uv P 0 h 0 vv + h 0 v ) P 0 4 P ) r ± = r ± r 0 ± h 0 uv + )] P 0 h 0 vv r + + O w 3 ) ) P 0 r + r 4P 0 r + O w 3 ) or the shited Riemann invariants. Assuming 3.13) do a rescaling x ± k x ±, r + k r +, r k 1/ r with a small positive k. Ater simple calculations one obtains 3.0) x + = r + + Ok 1/ ) x = r + Ok 1/ ) where x + = x + 0 vv P 0 x = x 0 uvv P 0 0 vvv + 0 vv P 0 4 P 0 observe that 0 vv 0 due to the genericity assumption 3.1)). We obtain Whitney normal orm o type W. Thereore the solution to the implicit unction system 3.4) near the critical point is deined on a hal plane with the boundary x = 0 passing through the critical point. Clearly such a behaviour is incompatible with the assumptions o the lemma. Theorem 3.3. Under the genericity assumption 3.14) the local structure o the solution to the implicit unction equations is described by the Whitney normal orm o type W 3. Proo We have to include the degree 3 terms in the Taylor expansions o xu, v), su, v). The rescaling 3.1) x + k /3 x +, x k x, r + k /3 r +, r k 1/3 r

20 0 B. DUBROVIN yields 3.) x + = 0 vv P 0 r + + Ok 1/3 ) x = vv 0 A h 0 0 r + r + 1 v 6 B 0 r 3 + Ok 1/3 ) where A 0 and B 0 are deined in 3.14). The theorem is proved. 3.. Critical points o Type II and elliptic umbilic catastrophe. Recall [47] that elliptic umbilic catastrophe can be represented as the singularity at X = Y = 0 o the unction 3.3) F X, Y ) = 1 3 X3 X Y. It can also be identiied see []) as one o the real orms o the codimension singularity o type D 4. The ollowing interpretation will be useul in sequel. Let ΦX, Y ) be a harmonic unction o two variables, Consider singularities o the gradient map Φ XX + Φ Y Y = ) X, Y ) Φ X X, Y ), Φ Y X, Y )) = U, V ) i.e. the points where the Jacobian o 3.4) vanishes, det U/ X U Y = [ ] Φ XX + Φ XY = 0. V/ X V/ Y As it ollows rom the above ormula or the Jacobian, at this point all the second derivatives o Φ vanish. So, the 3-jet o the unction must be a cubic harmonic polynomial. It must coincide with 3.3) up to a linear change o variables. The gradient map o 3.3) is given by the ollowing system U = X Y 3.5) V = X Y. In the complex coordinates this gives the quadratic map W = U i V, 3.6) W = Z. Z = X + i Y We will now show that the Type II singularity o a generic solution to the quasilinear system.40) is locally described by the elliptic umbilic catastrophe. As above we assume that the singular point is x C = s C = 0; it will also be assumed, as above, that 0 u = 0 v = 0.

21 UNIVERSALITY IN HAMILTONIAN PDES 1 Denote ū = u u C, v = v v C. Consider a solution to.40) determined by the system.41) near the critical point. Lemma 3.4. At the critical point all second derivatives o vanish: 3.7) 0 uu = 0 vv = 0 uv = 0. Proo ollows rom 3.6) due to negativity o P 0. Deinition 3.5. The Type II critical points is called generic i the 3-jet o u, v) does not vanish at this point. Lemma 3.6. Ater the rescaling x k x s k s ū k 1/ ū v k 1/ v at the limit k 0+ the solution to the system 3.6) near a generic Type II critical point x = x C = 0, s = s C = 0, u = u C, v = v C ) tends to the solution o the equation 3.8) z = 1 a 0w z := x + s h 0 u + i c 0 s h 0 v, w := v + i c 0 ū, a 0 := 0 uvv + i c 0 0 vvv where c 0 = P 0. Observe that, due to.3) and 3.7) at the critical point the third derivatives o the unction u, v) satisy 0 uuu = P 0 0 uvv, 0 uuv = P 0 0 vvv. So, at the generic Type II critical point we have necessarily a 0 0. Proo o Lemma is similar to that o Lemma 3.. Expanding the let hand sides o eqs. 3.4) in Taylor series and using equation.3) valid or both unctions and h one obtains x + s h 0 u + s c 0h 0 vvū + h 0 uv v) + s [ h 0 uvv c 0ū + v ) c 0h 0 vvvū v ] = 1 [ 0 uvv c 0ū + v ) c 0 0 vvvū v ] + O ū + v ) 3/) s h 0 v + s h 0 uvū + h 0 vv v) + s [ h 0 vvv c 0ū + v ) + h 0 uvvū v ] = 1 [ 0 vvv c 0ū + v ) + uvvū 0 v ] + O ū + v ) 3/)

22 B. DUBROVIN Ater the rescaling and division by k one arrives at the system x + s h 0 u = 1 [ 0 uvv c 0ū + v ) c 0 0 vvvū v] + O k 1/) s h 0 v = 1 [ 0 vvv c 0ū + v ) + 0 uvvū v] + O k 1/). This converges to the quadratic equation 3.8). The Lemma is proved. We see that, ater a linear change o independent variables the local structure o the solution near a generic Type II singularity is described by the biurcation diagram 3.6) o the elliptic umbilic catastrophe Generalization to the multicomponent case. Let us consider an n-component system o PDEs written in the diagonal orm 3.9) t u i = a i u) x u i, i = 1,..., n, n 3. Thus the dependent variables u 1,..., u n are the Riemann invariants or the system 3.9); the characteristic velocities a 1 u),..., a n u) will be assumed to be smooth unctions taking pairwise distinct values, Recall the ollowing deinition [48]. a i u) a j u), i j or all u D. Deinition 3.7. The system 3.9) is called semihamiltonian i the coeicients satisy the ollowing dierential equations ) ) ai,j ai,k 3.30) k = j a i a j a i a k or all triples o pairwise distinct indices i, j, k. Here i :=, a i,j := j a i ; u i we will also use similar notations or higher derivatives, i.e. etc. a i,jk := j k a i, A i,jkl := j k l A i The meaning o this deinition is clariied by the ollowing theorem [48]. Theorem 3.8. i) Let A 1 u),..., A n u) be an arbitrary solution to the overdetermined system o PDEs 3.31) A i,j = a i,j a i a j A i A j ), i j. Then the system o PDEs 3.3) s u i = A i u) x u i, i = 1,..., n commutes with 3.9), 3.33) s t u i t s u i.

23 UNIVERSALITY IN HAMILTONIAN PDES 3 ii) For a semihamiltonian system 3.9) the general solution to 3.3) depends on n arbitrary unctions o one variable. iii) Given a solution to the system 3.31) such that n A 1 u 0 ) = = A n u 0 ) = 0, A i,i u 0 ) 0 or some u 0 D, then the the system o equations 3.34) x = a i u) t + A i u), i = 1,..., n deines a unique smooth vector unction ux, t) or suiciently small x, t such that u0, 0) = u 0. This unction satisies the system 3.9). Conversely, any solution ux, t) to the semihamiltonian system 3.9) satisying i=1 x u 1 0, 0) 0,..., x u n 0, 0) 0 can be locally obtained by the above procedure. Remark 3.9. For a semihamiltonian system 3.9) there also exists an ininite amily o conserved quantities hu) such that the ollowing identity holds true on the solutions to 3.9) 3.35) t hu) = x u) or some unction u). The conserved quantities are determined rom the ollowing overdetermined system o PDEs 3.36) i j h = a i,j i h + a j,i j h, i j. a j a i a i a j The solutions to this system also depend on n arbitrary unctions o one variable [ibid]. The conserved quantities hu) were called entropies in [45]; the systems 3.9) possessing a rich amilies o entropies i.e., amilies depending on n arbitrary unctions o one variable) were called systèmes rich in [45]. This class coincides with the class o semihamiltonian systems. Let us describe the local structure o singularities o solutions to the semihamiltonian systems near the point o gradient catastrophe. Without loss o generality we may assume that the gradient catastrophe takes place at the point x 0 = 0, t 0 = 0. At the point o catastrophe the assumptions o the implicit unction theorem applied to the system o equations 3.34) ail to hold true. Explicitly this means that n 3.37) det a i,j u) t + A i,j u)) x=x0, t=t 0, u=u 0 = A i,i u 0 ) = 0 at the point o catastrophe. Deinition The point x = x 0 = 0, t = t 0 = 0, u = u 0 o gradient catastrophe is called generic i at this point 3.38) i=1 A i,i u 0 ) = 0, A j,j u 0 ) 0 or j i, A i,ii u 0 ) = 0, A i,iii u 0 ) 0.

24 4 B. DUBROVIN The irst part o this deinition says that only i-th Riemann invariant breaks down at the point o catastrophe. The assumption A i,ii u 0 ) = 0 is not a restriction; it says that the graph o the i-th Riemann invariant u i x, 0) has an inlection with vertical tangent at the point o catastrophe. As we will see the assumption A i,iii u 0 ) 0 guarantees that the inlection point does not degenerate. Let us assume that x = x 0 = 0, t = t 0 = 0, u = u 0 = 0 is a generic point o gradient catastrophe with i = n. We will denote a 0 i := a i u 0 ), A 0 i,i := A i,i u 0 ) etc. Observe that A 0 i,j = 0 or all i j due to the system 3.31). Perorm a linear change in the space o independent variables putting 3.39) z i := x a 0 i t, i = 1,..., n. Theorem Ater the rescaling 3.40) z m k /3 z m, z n k z n t k /3 t u m k /3 u m, u n k 1/3 u n m n m n the system 3.31) converges to 3.41) z m = A 0 m,mu m + O k 1/3), m n z n = a 0 n,nt u n A0 n,nnnu 3 n + O k 1/3). Proo. Spelling the equations 3.30) out one obtains 3.4) a i,jk = a i a j a k a i a j )a i a k ) a i,ja i,k + a i a j ) a i,k a k,j a i a k ) a i,j a j,k a i a j )a i a k )a j a k ) or all triples i, j, k o pairwise distinct indices. Using these equations along with the system 3.31) yields

25 UNIVERSALITY IN HAMILTONIAN PDES 5 Lemma 3.1. At the point o catastrophe one has A 0 m,nn = 0, A 0 m,ni = 0 or m n, i m, n 3.43) A 0 m,nnn = 0, A 0 m,nni = 0, A 0 m,nij = 0, A 0 m,ijk = 0 or m n, i, j, k m, n, i, j, k distinct A 0 n,ni = 0, A 0 n,ij = 0, or i, j n, i j A 0 n,nni = 0, A 0 n,nij = 0, A 0 n,ijk = 0 or i, j, k n, i, j, k distinct Corollary The Taylor expansions o the equations 3.34) near the point o catastrophe read [ x a 0 mt = t a 0 m,i u i + a 0 m,mu m + a 0 m,nu n + O u )] + A 0 m,mu m 1 i a 0 m,i a 0 m a 0 i i A 0 i,iu i + u m i a 0 m,i a 0 m a 0 i A 0 m,mu i 3.44) A 0 m,m + 1 a 0 A0 m,mmu m,n m + u m u n a 0 m a 0 n + 1 A 0 m,iij u i u j + 1 A 0 6 m,iii u 3 i i j + 1 u m A 0 m,mij u i u j + 1 u m A 0 m,mii u i i j + 1 u m Am,mmi u i A0 m,mmmu 3 m + O u 4). i i i 3.45) [ x a 0 nt = t a 0 n,i u i + a 0 n,nu n + O u )] 1 a 0 n,i a 0 n a 0 i + 1 i i A 0 i,iu i A 0 n,iij u i u j + 1 A 0 6 n,iii u 3 i + 1 u n A 0 n,nii u i A0 n,nnnu 3 n + O u 4). i j i It is understood that the summation indices i, j etc. in the sums denoted are all distinct rom n and in 3.44)) rom m. Applying the rescaling 3.40) to the equations 3.44), 3.45) one completes the proo o the Theorem. i

26 6 B. DUBROVIN 4. Perturbations o nonlinear wave equations 4.1. Deormation problem-i: main deinitions and the programme o research. Given a system o the irst order quasilinear PDEs 4.1) u t = Au) u x, u = u 1 x, t),..., u n x, t) ) admitting a Hamiltonian description 4.) u t = {ux), H 0 } 0, H 0 = h 0 u) dx with respect to a Poisson bracket o hydrodynamic type written in the lat coordinates in the orm 4.3) {u i x), u j y)} 0 = η ij δ x y), η ji = η ij = const, det η ij) 0, we say that a system o the orm 1.8) is a Hamiltonian deormation o 4.1) i it can be represented in the orm 4.4) u t = {ux), H} with a perturbed Hamiltonian 4.5) H = H 0 + ɛ H 1 + ɛ H +... H k = h k u; u x,..., u k) ) dx, k 1 and a perturbed Poisson bracket 4.6) deg h k u; u x,..., u k) ) = k {u i x), u j y)} = {u i x), u j y)} 0 + ɛ {u i x), u j y)} 1 + ɛ {u i x), u j y)} +... k+1 {u i x), u j y)} k = A ij k s ux); ux x),..., u s) x) ) δ k s+1) x y), k 1 s=0 deg A ij k s u; ux,..., u s)) = s. In the above ormulae the notations like u; u x,..., u k) ) are used or dierential polynomials u; u x,..., u k) ) C M)[u x,..., u k) ], u M, i.e., or polynomial unctions on the jet bundle J k M). The degree o a dierential polynomial is deined by the rule deg u i x = 1, deg u i xx =,..., i = 1,..., n. The sums in 4.5) and 4.6) are allowed to be ininite. However every term o the ɛ- expansions o 4.4) is a well-deined dierential polynomial. Antisymmetry and Jacobi identity or the Poisson bracket are understood as identities or ormal power series

27 UNIVERSALITY IN HAMILTONIAN PDES 7 in ɛ. The delta-unction symbol can be spelled out in the ollowing way. Introduce a matrix o linear dierential operators depending on the parameter ɛ 4.7) One clearly has Π ij = Π ij 0 + ɛ Π ij 1 + ɛ Π ij +... Π ij 0 = η ij x, k+1 Π ij k = A ij k s u; ux,..., u s)) x k s+1, k 1. s=0 4.8) {u i x), u j y)} = Π ij δx y). Then the perturbed Hamiltonian system reads 4.9) u i t = Π ij δh δu j x) = m 0 ɛ m k+l=m Π ij k δh l δu j x). The assumptions about the structure o the perturbative expansions in 4.5) and 4.6) imply that the terms o the expansion 1.8) B i m u; ux,... u m+1)) = k+l=m Π ij k are dierential polynomials o the degree m + 1. δh l, m 0, i = 1,..., n δu j x) The class o Hamiltonian deormations is invariant with respect to special) Miuratype transormations o the dependent variables 4.10) u ũ = u + k 1 ɛ k F k u; ux,..., u k)), deg F k u; ux,..., u k)) = k. The transormation o the Hamiltonian is deined by the direct substitution; the Poisson bracket is transormed by the rule 4.11) {ũ i x), ũ j y)} = Π ij δx y) Π ij = L i p Π pq L j q where L and L are respectively the operator o linearization o the transormation 4.10) and its ormal adjoint: 4.1) L i k = s ũ i u k,s s x, L i k = s x ) s ũi u k,s. The transormations o the orm 4.10) orm a group. Two Hamiltonian deormations o the quasilinear system 4.1) related by a transormation 4.10) are called equivalent. In particular the Hamiltonian deormation is called trivial i it is equivalent to the unperturbed system 4.1).

28 8 B. DUBROVIN Example 4.1. The ocusing/deocusing nonlinear Schrödinger equation 4.13) i ɛ ψ t + 1 ɛ ψ xx ± ψ ψ = 0 written in the coordinates u = ψ, v = ɛ i ψx ψ ψ ), x ψ takes the orm o a Hamiltonian perturbation o the system.37): 4.14) u t + u v) x = 0 v t + v v x u x = 1 4 ɛ uxx u 1 The Poisson bracket gets no perturbation: ) u x. u x {ux), vy)} = δ x y), {ux), uy)} = {vx), vy)} = 0. The perturbed Hamiltonian reads 4.15) H = H 0 + ɛ H, H 0 = 1 u u v ) dx, H = u x 8 u dx. Example 4.. The generalized FPU system with the Hamiltonian.6) in the translation invariant variables 4.16) and interpolation w n = q n q n 1 v n = p n w n = wɛ n, ɛ t) v n = vɛ n, ɛ t) yields the ollowing equations or the unctions w = wx, t), v = vx, t) 4.17) w t x) = 1 ɛ [vx) vx ɛ)] = v x 1 ɛ v xx ɛ v xxx +... v t x) = 1 ɛ [P wx + ɛ)) P wx))] = x P w) + 1 ɛ xp w) ɛ 3 xp w) This system is a Hamiltonian perturbation o the nonlinear wave equation.11). The Poisson brackets o the variables wx), vx) have a nontrivial dependence on ɛ. Indeed, the brackets o the discrete variables read {w n, v m } = δ n,m δ n,m+1, other brackets vanish. Ater interpolation and division by ɛ one obtains the perturbed bracket 4.18) {wx), vy)} = 1 ɛ [δx y) δx y ɛ)] = δ x y) 1 ɛ δ x y) ɛ δ x y) +....

29 UNIVERSALITY IN HAMILTONIAN PDES 9 One can reduce the Poisson bracket 4.18) to the standard orm {ux), vy)} = δ x y) by the Miura-type transormation.8). Indeed, introducing the shit operator Λ = e ɛ x rewrite the bracket 4.18) as ollows {wx), vy)} = 1 ) 1 Λ 1 δx y). ɛ The substitution.8) represented in the orm u = ɛ x 1 Λ w 1 implies {ux), vy)} = ɛ x 1 Λ {wx), vy)} = xδx y). 1 In the new variables the deormed system is described by the perturbed Hamiltonian: δh u t = x δvx) H = h dx = δh v t = x δux) [ 1 v x) + P wx) wx ɛ)) ] dx 4.19) h = 1 v + P u) ɛ 4 P u) u x + ɛ4 [ ] 8 P u) u xx P IV u) u 4 x + Oɛ 6 ) 5760 modulo inessential total derivatives). The series in ɛ in this case are all ininite. In act the ollowing general statement about deormations o the Poisson bracket 4.3) holds true. Theorem ). Any deormation 4.5) o the Poisson bracket 4.3) is trivial. ). Any special Miura-type transormation 4.10) canonical with respect to the bracket {, } = {, } 0 {ũ i x), ũ j y)} = {u i x), u j y)} is the time-ɛ shit generated by a Hamiltonian Φ: 4.0) ũ i = u i + ɛ {u i x), Φ} + ɛ {{ui x), Φ}, Φ} +... Φ = ϕ dx ϕ = k 0 ɛ k ϕ k u; u x,..., u k) ), deg ϕ k u; u x,..., u k) ) = k, k 0

30 30 B. DUBROVIN with some dierential polynomials ϕ 0 u), ϕ 1 u; u x ),.... Proo ollows rom triviality in positive degrees o Poisson cohomology o the bracket 4.3) established in [5], [1] see details in [3]). Thanks to the Theorem one reduces the study o Hamiltonian deormations o the Hamiltonian quasilinear equations to the study o systems written in the standard orm 4.1) u i t = η ij x δh δu j x), i = 1,..., n just perturbing the Hamiltonian, 4.) H = H [0] + ɛ H [1] + ɛ H [] +... H [k] = h [k] u; u x,..., u k) ) dx, deg h [k] u; u x,..., u k) ) = k, k 0 with some dierential polynomials h [0] u), h [1] u; u x ), h [] u; u x, u xx ),.... Here H [0] = h [0] u) dx is the Hamiltonian o the unperturbed equation. The equivalence o perturbed Hamiltonians is established by the canonical transormations 4.0): 4.3) H H = H + ɛ {H, Φ} + ɛ {{H, Φ}, Φ} In the subsequent sections we will discuss the main problems o the deormation theory o Hamiltonian PDEs, namely, Problem 1: describe irst integrals/ininitesimal symmetries o a given deormed Hamiltonian system by studying the procedure o extension o integrals/symmetries o the unperturbed system to integrals/symmetries o the perturbed system. Problem : classiy integrable Hamiltonian perturbations or which all irst integrals/ininitesimal symmetries o the unperturbed system can be extended to irst integrals/ininitesimal symmetries o the perturbed one. Problem 3: comparative study o the properties o solutions o the perturbed versus unperturbed systems. In a more speciic way, characterize the perturbations or which the solutions to the perturbed Cauchy problem, in a suitable class o initial data, exists or the times beyond the moment o gradient catastrophe or the solution to the unperturbed Cauchy problem. For such a class o perturbed Hamiltonian systems establish the universal properties o solutions near the point o gradient catastrophe o the unperturbed system, that is, those properties independent rom the choice o the initial data. A more ambitious project is in establishing universality independently on the choice o Hamiltonian perturbation.

31 UNIVERSALITY IN HAMILTONIAN PDES Deormation problem-ii: extension o ininitesimal symmetries and D- operator. Let 4.4) H pert = H [0] + ɛ H [1] + ɛ H [] +... H [k] = h [k] u; u x,..., u k) ) dx, deg h [k] u; u x,..., u k) ) = k, k 0 h [0] = 1 v + P u) be a perturbation o the Hamiltonian o the nonlinear wave equation. We know that, given a solution = u, v) to the linear PDE.3) the Hamiltonians H [0] and 4.5) H [0] = dx commute: {H [0], H [0] } = 0. The goal is to construct a deormation 4.6) H = H [0] + ɛ H [1] + ɛ H [] +... H [k] = h [k] u; u x,..., u k) ) dx, h [0] = u, v) or an arbitrary solution o.3) such that {H pert, H } = 0. Deinition 4.4. We say that the linear dierential operator D = D [0] + ɛ D [1] + ɛ D [] +... deg h [k] u; u x,..., u k) ) = k, k 0 4.7) D [0] = id, mk) D [k] = b [k] i 1...i mk) u; u x,..., u k) ) u i 1... u i mk) deg b [k] i 1...i mk) u; u x,..., u k) ) = k, k 1 is a D-operator or the perturbation 4.4) i or any solution to.3) the Hamiltonian 4.8) H = D dx commutes with H pert. In the ormula 4.7) mk) is some positive integer depending on k. It is easy to see that [ ] 3k 4.9) mk) =.

32 3 B. DUBROVIN The summation is taken over all indices i 1,..., i mk) rom 1 to. As usual the coeicients b [k] i 1...i mk) u; u x,..., u k) ) are dierential polynomials. As the D-operator makes sense only acting on the kernel o 4.30) u P u) v, the coeicients are not determined uniquely. One can choose a suitable normal orm or this operator, e.g., keeping only derivatives o the orm m v and u m 1 v. We will show that the D-operator is determined uniquely just rom the commutativity with H pert, as it ollows rom the ollowing Theorem 4.5. The D-operator, i exists, is determined uniquely up to a total x- derivative, i.e., i D 1 and D are two D-operators then D D 1 Im x Ker u P u) v). Proo ollows rom the ollowing result saying that the nonlinear wave equation has no irst integrals depending nontrivially on the derivatives. Lemma 4.6. Let P u) be a potential such that P u) 0. Let g = gu; u x,..., u k) ) be a dierential polynomial o the degree k > 0 such that the unctional G = gu; u x,..., u k) ) dx commutes with H [0], 4.31) {G, H [0] } = 0 Then g is a total x-derivative: gu; u x,..., u k) ) = x gu; u x,..., u k 1) ). The proo o this lemma is straightorward and will appear elsewhere. Corollary 4.7. Let the perturbed Hamiltonian 4.4) admit a D-operator. Then or any pair o solutions, g to.3) the Hamiltonians H and H g commute: {H, H g } = 0, H = D dx, H g = Dg dx 4.3) uu = P u) vv, g uu = P u)g vv. Proo From Jacobi identity it ollows that {{H, H g }, H pert } = 0. Assume the Poisson bracket {H, H g } to have the orm {H, H g } = ɛ k F [k] + O ɛ k+1)

33 UNIVERSALITY IN HAMILTONIAN PDES 33 or some k 1. The unctional F [k] must commute with H [0]. Due to the Lemma it must vanish.. Remark 4.8. One can consider a more general class o the irst order quasilinear Hamiltonian systems written in the orm 4.33) u t = x h v v t = x h u or a given smooth unction h = hu, v). The conservation laws = u, v) or this system are determined rom the linear PDE 4.34) h vv uu = h uu vv. The associated Hamiltonians H [0] = commute pairwise. We say that the system 4.33) is regular 5 i it has no nontrivial local conservation laws depending on higher derivatives c. Lemma 4.6). It would be interesting to characterize all regular quasilinear Hamiltonian systems. We end this Section with considering a simple example o Hamiltonian perturbations o the scalar transport equation u t + u u x = 0 classiied in [17]. The Hamiltonian structure o the transport equation reads dx u t + {ux), H [0] δh [0] } u t + x δux) = 0 H [0] = 1 u 3 dx, {ux), uy)} = δ x y). 6 All perturbations o the Hamiltonian H [0] up to the order ɛ 4 are parameterized by two arbitrary unctions o one variable c = cu) and p = pu) as ollows: [ ] u 3 H pert = 6 ɛ 4 cu) u x + ɛ 4 pu) u xx dx. This perturbation turns out to be integrable with the order ɛ 4 approximation {H, H g } = O ɛ 6), H = D c,p dx, H g = D c,p g dx, = u), g = gu). The associated D-operator has the orm ) D c,p = ɛ 4 c u x + ɛ [p 4 + c 4) u xx ) c c 4) c c 5) c 6) p 4) + p ) ] 5) u 4 x + O ɛ 6), Similarly to the deinition o regular elements in Lie algebras.

34 34 B. DUBROVIN the unction in this case is arbitrary Deormation problem-iii: an integrability test. Existence o a D-operator or a given perturbed Hamiltonian gives rise to a simple integrability test, as we will now explain or the example o the generalized FPU system. The perturbed Hamiltonian reads [ ] ) H pert = v + P u) ɛ 4 P u) u x dx + O ɛ 4) see 4.19) above). Under what condition or the potential P u) a D-operator exists within the ɛ approximation? We will show that this is possible only under the constraint 4.37) P P IV = P ). To derive 4.37) let us look or the perturbed Hamiltonian density with indetermined coeicients H = h dx, {H, H pert } = O ɛ 3) h = + ɛ h [1] u; u x) + ɛ h [] u; u x, u xx ). The irst correction must be linear in u x, v x. Adding i necessary a total x-derivative one can reduce to the study o the irst perturbation o the orm h [1] = pu, v) v x or some unction p = pu, v). Computation o the bracket gives [ {H, H pert } = ɛ p u v x P u)ux] ) dx + O ɛ. It is easy to see that the integrand is never a total x-derivative unless p u = 0, i.e. p = pv). Hence h [1] must be itsel a total x-derivative. Let us now consider the order terms. Up to a total x-derivative they can be written in the orm h [] = 1 ) au, v)u x + bu, v)u x v x + cu, v)vx. Computation o the order ɛ terms in the Poisson bracket gives the ollowing unctional: {H, H pert } = = ɛ { 1 1 vp u xxx + [ ) ] 1 6 vp b P u x a v x u xx c P u x + b v x ) v xx + 1 ) v P IV + 1P a v 4P b u u 3 4 x 1 a u + P c u ) u xv x 1 a v + P c v ) u x vx + 1 } c u b v ) vx 3 dx + O ɛ 4). Denote I the integrand. One must veriy the equations E u I = 0, E v I = 0.

35 UNIVERSALITY IN HAMILTONIAN PDES 35 Vanishing o the coeicient o u xxx in E v I gives a = c 1 ) 1 vv P. Next, vanishing o the coeicient o v xx v x in E v I implies b v = c u. Hence it exists a potential λ = λu, v) such that b = λ u, c = λ v. Ater the substitution we obtain, collecting the coeicients o u xx u x, u xx v x, v xx u x, v xx v x, the ollowing system o equations It ollows that Equating the mixed derivatives 6P λ v + P uvv = 0 P λ u P λ uu + P λ vv = 0 6λ uu 6P λ vv + P vvv = 0. λ u = 1 P 6 P vvv P λ v = 1 6 P uvv. implies This proves 4.37). λ u ) v = λ v ) u P ) P P IV 6 P ) uvv = 0. The equation 4.37) or the potential can be easily solved: dividing by P P yields Integration gives Hence or some constants a, b, c, k. P P = P IV P. P = c P. P u) = k e c u + a u + b We have proved that the only generalized FPU system passing the integrability test must have exponential potential. This gives Toda lattice. Remark 4.9. It would be interesting to compare the above arguments with the theory o small amplitude normal orms or inite FPU chains developed recently by A.Henrici and T.Kappeler [30, 31].

36 36 B. DUBROVIN 4.4. Deormation problem-iv: on classiication o integrable perturbations. Let = u, v) and g = gu, v) be two irst integrals o the nonlinear wave equation.11), i.e., two solutions to the linear PDE.3). The Hamiltonians H [0] = dx, H [0] g = commute: {H [0], H[0] g } = 0. We obtain a commutative Lie algebra o Hamiltonians parameterized by solutions to.3). The goal is to classiy deormations o this algebra o the orm 4.4) preserving the commutativity. So ar the classiication is available only within the ɛ approximation. Let us describe this preliminary result and outline the classiication procedure. Theorem Any -integrable perturbation o the nonlinear wave equation.11) is parameterized by two arbitrary unctions o one variable ρ ± = ρ ± r ± ). The associated D-operator reads {[P ρ u vvv + ρ v uvv ) + 1 ] P ρ v vv g dx D ρ,r = + ɛ + P ρ v vvv + ρ u uvv + P ) 4P ρ u vv u x v x + ρ u vvv + ρ v uvv ) v x} + Oɛ 3 ), where ρ = ρu, v) = ρ + r + ) + ρ r ), the Riemann invariants r ± = r ± u, v) are deined in.0). u x Proo Let = u, v), g = gu, v) be two arbitrary solutions to the linear PDE.3). Any Oɛ)-perturbation o the commuting Hamiltonians H [0] and H g [0] can be written in the orm H = H [0] + ɛ pu, v)v x dx + O ɛ ), H g = H g [0] + ɛ qu, v)v x dx + O ɛ ) or some unctions p and q. From the commutativity {H, H g } = O ɛ ) it ollows, calculating the coeicient o u xx in u p 4.38) vv Consider now a canonical transormation generated by the Hamiltonian = uq g vv. δ {H δux), H g }, that u u + ɛ{ux), K} + O ɛ ) K = ku, v) dx.

37 UNIVERSALITY IN HAMILTONIAN PDES 37 The Hamiltonians H and H g will transorm according to the ollowing rule: H H + ɛ{h [0], K} + O ɛ ), H g H g + ɛ{h [0] g, K} + O ɛ ). Choosing the unction k = ku, v) satisying k uu P u)k vv = up vv = uq g vv one eliminates the terms linear in ɛ rom both perturbed Hamiltonians H and H g. At the next step one has to consider the perturbation o order ɛ that can be chosen in the ollowing orm: H = H [0] + ɛ a ) u x + b u x v x + c v ) x dx + O ɛ 3 H g = H [0] g + ɛ ag ) u x + b g u x v x + c g v ) x dx + O ɛ 3 with some unctions a i = a i u, v), b i = b i u, v), c i = c i u, v), i =, g. From the commutativity {H, H g } = O ɛ 3) it ollows, calculating the coeicient o v xxx in Hence or some unction ku, v). At the next step we substitute δ {H δux), H g }, that a c P u)) g vv = a g c g P u)) vv. b = r 1 vvv + r uvv + r 3 vv + r 4 uv, b g = r 1 g vvv + r g uvv + r 3 g vv + r 4 g uv, a = c P u) + ku, v) vv a g = c g P u) + ku, v)g vv c = s 1 vvv + s uvv + s 3 vv + s 4 uv c g = s 1 g vvv + s g uvv + s 3 g vv + s 4 g uv with some coeicients r i = r i u, v), s i = s i u, v), i = 1,..., 4 the coeicients o the D-operator). The calculations similar to above yield the ollowing expressions or the coeicients: r 1 = P ρ v, r = ρ u, r 3 = P s 1 = ρ u, s = ρ v, s 4 = 0 k = 1 P ρ v s 3 P 4P ρ u, r 4 = 0 where the unction s 3 = s 3 u, v) is arbitrary and the unction ρ = ρu, v) satisies the linear PDE ρ uu P ρ vv = P P ρ u.

38 38 B. DUBROVIN It is readily seen that this equation is equivalent to the D Alembert equation in Riemann invariants: ρ uu P ρ vv P P ρ u = 4P ρ. r + r Finaly, we obtain the ɛ -correction o the Hamiltonian densities in the ollowing orm h = + {[P ɛ ρ u vvv + ρ v uvv ) + 1 ] P ρ v vv u x + P ρ v vvv + ρ u uvv + P ) 4P ρ u vv u x v x + ρ u vvv + ρ v uvv ) v x + s 3 v x P u x) vv } + O ɛ 3 ), a similar expression holds or h g replacing g. It is easy to see that the term containing the arbitrary unction s 3 can be eliminated by a canonical transormation generated by the Hamiltonian Φ = ɛ φu, v)v x dx with φ = φu, v) determined rom φ u = s Quasitriviality and solutions to the perturbed nonlinear wave equations. Recall see Section 4.1 above) that the ɛ-perturbation o a irst order quasilinear system is called trivial i it can be eliminated by a special) Miura-type transormation 4.10). All terms in the ɛ-expansion o this transormation are dierential polynomials. Deinition We say that the perturbation is quasitrivial i it can be eliminated by a transormation o the orm 4.39) u ũ = u + k 1 ɛ k F k u; ux,..., u mk))), deg F k u; ux,..., u mk))) = k where the terms F k or k 1 are rational unctions in the jet variables u x, u xx,... u mk)) or some integer mk). As it was proved in [19], all bihamiltonian perturbations are quasitrivial. Moreover, the order mk) o the highest derivative in F k can be estimated as [ ] 3k 4.40) mk) =. It would be interesting to extend the Quasitriviality Theorem o [19] to all Hamiltonian perturbations 6. Existence o a quasitriviality transormation gives a universal ormula or the perturbative solutions to the deormed systems with very general classes o initial data. 6 The quasitriviality theorem was proved by S.-Q.Liu and Y.Zhang in [41] or a more general class o perturbations o scalar PDEs. In the recent paper [40] the quasitriviality technique was applied to the analysis o existence o a local Hamiltonian structure or scalar evolutionary PDEs.

39 UNIVERSALITY IN HAMILTONIAN PDES 39 Namely, substituting a solution u 0 = u 0 x, t) o the unperturbed system into 4.39) one obtains a ormal series in ɛ satisying the deormed PDE in all orders. In the examples o Section 5 below all the rational unctions in the quasitriviality transormation have the denominators that are powers o 4.41) = v x P u)u x. In the Hamiltonian case the quasitriviality transormation can be represented as a time-ɛ canonical transormation see [17]). That means that there exists a Hamiltonian K = m 1 ɛ m K m 4.4) K m = k m u; u x,..., u ) ) dx, m 1 such that the canonical transormation 4.43) u ũ = u + ɛ{ux), K} + ɛ {{ux), K}, K} +... transorms all unperturbed commuting Hamiltonians into the perturbed ones, 4.44) H = ũ) dx or any satisying uu = P u) vv. The generating Hamiltonian K has to be ound rom the equation 4.45) H [0] + ɛ{h [0], K} + ɛ {{H[0], K}, K} + = H H [0] + ɛ H [1] + ɛ H [] The Hamiltonians constructed in this way will automatically commute; the highly nontrivial set o cancellations ensures that the densities o these Hamiltonians up to total x-derivatives are polynomials in the jet coordinates u x, u xx etc. The above perturbative expansions can be used or asymptotic integration o the deormed Hamiltonian systems or suiciently small times although rigorous analytic results justiying validity o these asymptotic expansions are still to be established). They clearly lose their validity when approaching the time o gradient catastrophe since the denominators become all equal to zero. In order to study the behaviour o the perturbed solutions near the critical point we will now describe an alternative procedure o constructing solutions to the perturbed Hamiltonian systems. Let us irst prove the ollowing Lemma 4.1. Let the perturbed Hamiltonian o the nonlinear wave equation have the orm [ ) H pert = v + P u) + ] ɛ k h [k] u; u x, v x,..., u k), v k) ) dx k 1 i.e., no explicit dependence on v in the higher corrections). The the D-operator, i exists, has v-independent coeicients.

40 40 B. DUBROVIN Proo Indeed, under the above assumptions the coeicients o the commutativity equation {H, H pert } = 0 or H = D dx [ δh δ E u v x + ɛ H ) pert x + δh P δ u)u x + ɛ H )] pert x = 0 δux) δvx) δvx) δux) [ δh δ E v v x + ɛ H ) pert x + δh P δ u)u x + ɛ H )] pert x = 0 δux) δvx) δvx) δux) do not depend on v see.17)). Here H pert = ɛ k 1 h [k] u; u x, v x,..., u k), v k) ) dx. Moreover, the shit k 1 v v + c maps to itsel the space o solutions to the linear PDE.3). Thereore the same shit applied to the D-operator produces another D-operator. Due to the uniqueness theorem 4.5 the shited operator coincides with D ater adding, i necessary, a total x-derivative. Remark The class o Hamiltonan perturbations considered in the Lemma corresponds to the ollowing subclass o perturbed nonlinear wave equations: 4.47) u tt = P u)u xx + ɛ x au; u x, u t, u xx, u xt,... ; ɛ). Given a solution h = hu, v) to the PDE.3), consider a Hamiltonian perturbation o the quasilinear system u t = x h v v t = x h u we have changed notation s t or the time variable). Solutions to this system are parameterized in the orm.40) by unctions = u, v) satisying the same PDE.3) according to the method o characteristics. Let us assume that the perturbed Hamiltonian equation admits a v-independent D-operator. We will show that the solutions to the perturbed equation close to.40) can be obtained rom a system o ininite order ODEs written in a variational orm. This system o ODEs will be called string equation. Theorem The solutions to the string equation 4.48) x + t δh h δux) = δh δux) t δh h δvx) = δh δvx)

41 UNIVERSALITY IN HAMILTONIAN PDES 41 also solve the Hamiltonian equations 4.49) u t = x δh h δvx) v t = x δh h δux) where = u, v) is another solution to uu = P u) vv, h := h v. Proo Since the coeicients o D-operator are v-independent, one has an identity 4.50) v Dh = Dh. Due to obvious commutativity [ ] [ ] v, E u = v, E v = 0 in the assumptions o the Theorem rom 4.50) one derives 4.51) δh h v δux) = δh h δux), δh h v δvx) = δh h δvx). Let us use the above identities in order to prove that the low 4.5) u s = t x δh h δvx) v s = 1 + t x δh h δux) is an ininitesimal symmetry o the system 4.49). Indeed, δh h u t ) s = t { x δvx), H h } + δh h x v δvx) = t {{ux), H δh h h}, H h } + x δvx) v t ) s = t {{vx), H h }, H h } + x δh h δux). Computation o the derivatives in the opposite order yields δh h u s ) t = x δvx) + t { δh h x δvx), H δh h h} = x δvx) + {{ux), H h }, H h} v s ) t = x δh h δux) + {{vx), H h }, H h}. The Jacobi identity and commutativity o the Hamiltonians {H h, H h } = 0 imply {{ux), H h }, H h } = {{ux), H h }, H h } {{vx), H h }, H h } = {{vx), H h }, H h }.

42 4 B. DUBROVIN So u t ) s = u s ) t, v t ) s = v s ) t. The commutativity remain valid adding to the symmetry 4.5) an arbitrary low o the hierarchy, e.g., replacing 4.5) with 4.53) δh h u s = t x δvx) δh x δvx) δh h v s = 1 + t x δux) δh x δux) or some solution to the linear PDE.3). It is wellknown see, e.g., [44, 37]) that the set o stationary points o an ininitesmal symmetry is invariant or the low 4.49). Integrating one obtains the string equation 4.48). In the dispersion limit the string equation coincides with the equation.41) o the method o characteristics: x + t h uv = u t h vv = v Special Painlevé-type transcendents and Universality Conjectures. We have seen in Section 3 above that locally near the point o catastrophe all generic solutions to a nonlinear wave equation and, more generally, to any Hamiltonian PDE commuting with the wave equation, have a very standard behaviour. Namely, i the catastrophe takes place in the domain o hyperbolicity then the local behaviour is described by the Whitney singularity equations. For the catastrophe in the domain o ellipticity the local behaviour is given in terms o elliptic umbilic singularity. In both cases the shape o the solution is locally universal, i.e., independent rom the choice o the solution, up to a simple aine transormation o the independent variables and a nonlinear transormation o the dependent variables. The Universality Conjecture or the perturbed systems says that, similarly, the leading term in the asymptotic expansion o the solution at ɛ 0 near a critical point does not depend on the choice o a generic solution, up to simple transormations o the independent/dependent variables. The universal shapes that occur in such a local description are given in terms o certain very particular solutions to Painlevé equations and their generalizations. Let us irst describe these solutions. The description will be given separately or catastrophes taking place in the domain o hyperbolicity or domain o ellipticity. Type I singularity. Consider the ollowing 4th order ODE or the unction U = UX) depending on the parameter T [ ) X = T U 6 U U + U U ) + 1 ] 40 U IV, U = du dx etc. The equation 1.10) can be considered as a higher order analogue o the Painlevé-I equation. It is oten identiied as the PI equation. As it was proved by T.Claeys and M.Vanlessen [9], this equation has a particular solution without singularities on the real axis X R or all real values o the parameter T. Denote this solution UX, T ).

43 UNIVERSALITY IN HAMILTONIAN PDES 43 We reer to the paper [9] or the precise characterization o the particular solution in terms o the Stokes multipliers o the associated with 1.10) linear dierential operator with polynomial coeicients. For large X the solution grows like UX, T ) 6 X) 1/3. It is monotone decreasing or T << 0; or positive T it develops a zone o oscillations see Fig.1). Fig.1: The solution to 1.10) or T > 0 We are ready to ormulate the Universality Conjecture describing the local behaviour o solutions to the perturbed Hamiltonian system 4.49) near the Type I critical point x c, t c, u c, v c ). Denote r ± the Riemann invariants.0) o the unperturbed system. Let x ± = x x c ) + c ± t t c ), c ± = h u u c, v c ) ± P u c )h v u c, v c ) be the characteristic variables at the point o catastrophe. Let us assume that at the point o catastrophe o the unperturbed system the second Riemann invariant r breaks down, while the invariant r remains smooth see 3.) above). Then one should expect the ollowing asymptotic behaviour o solutions to 4.49) near the critical point: r + = r + u c, v c ) + c + x + + a + ɛ 4/7 U x b ɛ, b x ) + 6/7 + + O ɛ 6/7) ɛ 4/7 4.55) r = r u c, v c ) + a ɛ /7 U x b ɛ, b 6/7 + x ) + ɛ 4/7 + O ɛ 4/7) or some constants a ±, b ±, c + depending on the choice o solution and on the perturbed Hamiltonian. Here U = UX, T ) is the particular solution to the 4th order ODE 1.10) described above. In the next Section we will explain how to compute the constants or the particular examples o integrable systems.

44 44 B. DUBROVIN Type II singularity. Here we need to speciy a particular solution W = W 0 Z) to the classical Painlevé-I P I ) equation 4.56) W ZZ = 6 W Z. The particular solution in question was discovered by Boutroux in the beginning o the last century [4]. It is characterized by the asymptotic behaviour Z W 6, arg Z < 4π 5, Z the principal branch o the square root is chosen). In particular the solution W 0 Z) is ree o poles within the sector arg Z < 4π 5 or suiciently large Z. Boutroux called it the tritronquée solution. A.Kapaev [34] gave a complete characterization o the tritronquée solution in terms o particular values o the Stokes multipliers o the linear dierential operator with polynomial coeicients associated with 4.56). Fig.: The real let) and imaginary parts o the tritronquée solution to Painlevé-I equation in the sector arg Z < 4π 5. In [18] it was conjectured that all poles o the tritronquée solution to P I are conined within the sector arg Z 4π. The Fig. borrowed rom [18] shows the graphs o the 5 real and imaginary parts o the tritronquée solution suggesting that, indeed there are no poles in the sector arg Z < 4π. 5 The Universality Conjecture or the Type II critical behaviour suggests the ollowing structure o the leading term o the asymptotic behaviour o complex) Riemann