threephase freak waves  arxiv · threephase freak waves a. o. smirnov, s. g. matveenko, s. k....
Embed Size (px)
TRANSCRIPT

THREEPHASE FREAK WAVES
A. O. SMIRNOV, S. G. MATVEENKO, S. K. SEMENOV, AND E. G. SEMENOVA
Abstract. Threephase finitegap solutions of the focusing nonlinear Schrödinger, KadomtsevPetviashvili and Hirota equations are constructed. These solutions have a behavior of almost
periodic “freak waves”.
Introduction
This study was motivated by the intention to demonstrate the behavior of threephase extremewaves. In the last time it was realized that the simplest and most universal model for such wavesis the focusing nonlinear Schrödinger equation (NLS)
(1) ipt + pxx + 2 p2 p = 0, i2 = −1,Equation (1) is used when describing distribution on the surface of the ocean weakly nonlinearquasimonochromatic wave packets with a relatively steep fronts since 1968 [49]. An application ofthis equation to the problems of nonlinear optics was known earlier [7]. Since the equation (1) is amodel of first approximation, it will appear in simulations of many weakly nonlinear phenomena.Fields of application of this equation are from plasma physics [32] to financial markets [48].
One of equation (1) property is a modulation instability, leading to the appearance of socalled“freak waves” (in hydrodynamics known as “rogue waves”) [3]. These waves are localized in spaceand time amplitude’s peaks. In the last 20 years, first in hydrodynamics, and then in nonlinearoptics, these waves have been the object of numerous theoretical and experimental studies [2].Such attention to the problem of “freak waves” is due, in particular, losses from destruction “roguewaves” oil platforms, tankers, container ships and other large vessels.
There are many more precise and more complex models, which give a more fine description of“freak waves” [2]. These models can be divided into two classes. For some models, like to equation(1) can be applied analytical methods. Other models are nonintegrable and can be solved bynumerical methods only. Analytical methods include:
• inverse scattering transform method;• finitegap integration method;• Bäcklund transform method;• Darboux transform method;• Hirota method.
In the present work we use a finitegap integration method. This method was created in theworks of Dubrovin, Novikov, Marchenko, Lax, McKean, van Moerbeke, Matveev, Its, Krichever[11–15, 24, 25, 30, 33, 35, 37, 39] (see also review article [36]). It should be mentioned that another
2010 Mathematics Subject Classification. 35Q55, 37C55.Key words and phrases. nonlinear Schrödinger equation; Hirota equation; freak waves; theta function; reduction;
covering; spectral curve.
1
arX
iv:1
412.
1562
v1 [
mat
hph
] 4
Dec
201
4

2 A. O. SMIRNOV, S. G. MATVEENKO, S. K. SEMENOV, AND E. G. SEMENOVA
method of constructing finitegap solutions of integrable nonlinear equations exists [27, 28, 38, 40].Let us remark that first method is based on BakerAkhiezer function but second method is basedon some Fays identities [16]. In our paper we use first method and Its’ and Kotlyarov’s classicformulas [23,26] (see also [6]).
Our aim here is to show a behavior of threephase algebrogeometric solutions of NLS, KPI andHirota equations. The first section of this paper contains the basic notations and classic formulasfor algebrogeometric solutions of considered integrable nonlinear equations. The second section ofthe paper is devoted to the periodicity of threephase solutions of NLS, KPI and Hirota equations.In the third section we consider an example of threephase algebrogeometric solution of KPI andHirota equations for different values of parameters.
1. Finitegap multiphase solutions of the NLS equation
The nonlinear differential equations that are integrated by methods of the algebraic geometry,can be obtained as a compatibility condition of system of the ordinary linear differential equationswith spectral parameter [6,19,20]. In particular, let us consider the following equations [19,22,42]
Yx = UY,(2a)
Yz = VY,(2b)
Yt = WY,(2c)
where
U = −λ(i 00 −i
)+
(0 iψ−iφ 0
),
V = 2λU + V0, W = 4λ2U + 2λW0 + W1,
λ is a spectral parameter. Using these equations and additional relations
(Yx)z = (Yz)x, (Yx)t = (Yt)x
one can easy obtain socalled equations of zero curvature
(3) Uz −Vx + UV−VU = 0 and Ut −Wx + UW−WU = 0,which should be valid for all values of spectral parameter λ. Respectively, it follows from eqs. (3)that matrixes V0,W0,W1 have to look like
(4) W0 = V0 =
(−iψφ −ψx−φx iψφ
), W1 =
(ψxφ− ψφx 2iψ2φ− iψxx−2iψφ2 + iφxx ψφx − ψxφ
),
Also that, accordingly, W = 2λV+W1. Conditions (3) lead also to additional system of equations(parities). First system is the coupled nonlinear Schrödinger equation
(5)
{iψz + ψxx − 2ψ2φ = 0,iφz − φxx + 2ψφ2 = 0,
and the second system is the coupled modified Kortewegde Vries equation
(6)
{ψt + ψxxx − 6ψφψx = 0,φt + φxxx − 6ψφφx = 0.

THREEPHASE FREAK WAVES 3
These two systems of the nonlinear differential equations are closely connected with two otherones. Namely, by differentiating eqs. (5) on x and substituting them in (6), one obtains the coupledmodified twodimensional nonlinear Schrödinger equation in a cone coordinates [31]
(7)
{iψt + ψxz + 2i(ψφx − φψx)ψ = 0,iφt − φxz + 2i(φψx − ψφx)φ = 0,
Also the functions ψ(x, t,−αt) and φ(x, t,−αt) are solutions of the coupled integrable Hirota equation (α ∈ R)
(8)
{iψt + ψxx − 2ψ2φ− iα(ψxxx − 6ψφψx) = 0,iφt − φxx + 2ψφ2 − iα(φxxx − 6ψφφx) = 0,
if ψ(x, z, t) and φ(x, z, t) are solutions of (5) and (6).Systems of the nonlinear differential equations (5), (6) are the first two integrable systems from
the AKNS hierarchy [19]. One of features of finitegap multiphase solutions of the integrablenonlinear equations is that fact that in some sense they are the solutions of all hierarchy. Particulary,our solutions can be used for constructing solutions of generalized nonlinear Schrödinger equation[47]. Substitutions φ = ±ψ into eq. (5) give us a standard form of the nonlinear Schrödingerequation. Particularly, for φ = −ψ equations (5) transform to (1) [13, 22, 26] and equations (8)transform to the integrable Hirota equation [4, 8, 21,34]
(9) iψt + ψxx + 2 ψ2 ψ − iα(ψxxx + 6 ψ2 ψx) = 0.It is also easy to check that for any ψ and φ, that satisfy both (5) and (6) simultaneously, the
function u(x, z, t) = −2ψφ is a solution of the KadomtsevPetviashviliI equation (KPI)(10) 3uzz = (4ut + uxxx + 6uux)x.
In the case φ = ±ψ this solution is a real function.Finitegap solutions of systems (5), (6) are parameterized by the hyperelliptic curve Γ = {(χ, λ)}
of the genus g [19, 42]:
Γ : χ2 =
2g+2∏j=1
(λ− λj),
The branch points (λ = λj , j = 1, . . . , 2g+ 2) of this curve are the endpoints of the spectral arcs ofcontinuous spectrum of Dirac operator (2a). Infinitely far point of the spectrum corresponds twodifferent points P±∞ on the curve Γ. In the case φ = −ψ the curve Γ has the form
(11) Γ : χ2 =
g+1∏j=1
(λ− λj)(λ− λj) = λ2g+2 +2g+2∑j=1
χjλ2g+2−j , =χj = 0, =(λj) 6= 0.
Following a standard procedure of constructing a finitegap solutions [6, 13, 42], let us to chooseon Γ a canonical basis of cycles γt = (a1, . . . , ag, b1, . . . , bg) with matrix of intersection indices
C0 =
(0 I−I 0
).
From the condition φ = −ψ it following that this basis of cycles satisfies the transformation relations[6, 13]
(12) τ̂1a = −a, τ̂1b = b +Ka.

4 A. O. SMIRNOV, S. G. MATVEENKO, S. K. SEMENOV, AND E. G. SEMENOVA
where τ1 is antiholomorphic involution, τ1 : (χ, λ)→ (χ, λ).Let us also consider normalized holomorphic differentials dUj :
(13)
∮ak
dUj = δkj , k, j = 1, . . . , g,
and a matrix of periods B of the curve Γ:
(14) Bkj =
∮bk
dUj , k, j = 1, . . . , g.
It is well known (see, for example [5, 13]) that the matrix B is a symmetric matrix with positivelydefined imaginary part.
Let us introduce in consideration gdimensional Riemann theta function with characteristicsη, ζ ∈ Rg [5, 13,16]:
(15)Θ[ηt; ζt](pB) =
∑m∈Zg
exp{πi(m + η)tB(m + η) + 2πi(m + η)t(p + ζ)},
Θ[0t; 0t](pB) ≡ Θ(pB),where B is a matrix of periods, p ∈ Cg and summation passes over an integer gdimensional lattice.
Let us also define normalized Abelian integrals on the Γ: of the second kind, Ωj(P), j = 1, 2, 3and the third kind, ω0(P), with the following asymptotic at infinitely far points P±∞:∮
ak
dΩ1 =
∮ak
dΩ2 =
∮ak
dΩ3 =
∮ak
dω0 = 0, k = 1, . . . , g,
Ω1(P) = ∓i(λ−K1 +O
(λ−1
)), P → P±∞,
Ω2(P) = ∓i(2λ2 −K2 +O
(λ−1
)), P → P±∞,
Ω3(P) = ∓i(4λ3 −K3 +O
(λ−1
)), P → P±∞,
ω0(P) = ∓(lnλ− lnK0 +O
(λ−1
)), P → P±∞,
χ = ±(λg+1 +O (λg )
), P → P±∞.
Let us denote by 2πiU, 2πiV, 2πiW the vectors of bperiods of Abelian integrals of the secondkind Ω1(P), Ω2(P), Ω3(P), respectively.
Theorem 1 ( [6, 42]). Function
Y (P, x, z, t) =(y1(P, x, z, t) y1(τ0P, x, z, t)y2(P, x, z, t) y2(τ0P, x, z, t)
),
where τ0 is hyperelliptic involution, τ0 : (χ, λ)→ (−χ, λ),
(16) y1(P, x, z, t) =Θ(U(P) + Ux+ Vz + Wt−X)Θ(Z)Θ(U(P)−X)Θ(Ux+ Vz + Wt+ Z) ×
× exp{Ω1(P)x+ Ω2(P)z + Ω3(P)t+ iΦ(x, z, t)},
y2(P, x, z, t) = ρΘ(U(P) + Ux+ Vz + Wt+ ∆−X)Θ(Z−∆)
Θ(U(P)−X)Θ(Ux+ Vz + Wt+ Z) ×
× exp{Ω1(P)x+ Ω2(P)z + Ω3(P)t− iΦ(x, z, t) + ω0(P)},

THREEPHASE FREAK WAVES 5
is the eigenfunction of the Dirac operator (2a) with functions
(17)
ψ(x, z, t) =2K0ρ
Θ(Z)Θ(Ux+ Vz + Wt+ Z−∆)Θ(Z−∆)Θ(Ux+ Vz + Wt+ Z) exp{2iΦ(x, z, t)},
φ(x, z, t) = 2ρK0Θ(Z−∆)Θ(Ux+ Vz + Wt+ Z + ∆)
Θ(Z)Θ(Ux+ Vz + Wt+ Z)exp{−2iΦ(x, z, t)},
for any z, t and ρ 6= 0. The functions (17) satisfy to the equations (5) and (6). Here ∆ is vectorof holomorphic Abelian integrals, calculated along a path, connecting P−∞ and P+∞, without crossingany of basic’s cycles,
∆ = U(P+∞)− U(P−∞), Φ(x, z, t) = K1x+K2z +K3t,
X = K +g∑j=1
U(Pj), Z = U(P+∞)−X,
K is a vector of Riemann constants [5, 13, 16, 29]; Pj, j = 1, . . . , g is a nonspecial divisor. If thespectral curve Γ satisfies the condition (11) then the following equalities hold
(18) ψ2 = −4K20Θ(Ux+ Vz + Wt+ Z−∆)Θ(Ux+ Vz + Wt+ Z + ∆)
Θ2(Ux+ Vz + Wt+ Z),
=U = =V = =W = =Z = 0, K20 < 0.It is easy to see that corresponding solution of KPI equation (10) has the form
(19) u(x, z, t) = −8K20Θ(Ux+ Vz + Wt+ Z−∆)Θ(Ux+ Vz + Wt+ Z + ∆)
Θ2(Ux+ Vz + Wt+ Z),
and that the square of amplitude of solution of Hirota equation (9) equals
(20) ψH 2 (x, t) = −4K20Θ(Ux+ (V − αW)t+ Z−∆)Θ(Ux+ (V − αW)t+ Z + ∆)
Θ2(Ux+ (V − αW)t+ Z) .
2. Features of threephase solutions
In a case g = 3 basis of normalized holomorphic differentials is defined by the formula [6, 13]:
(21) dUk = (ck1λ2 + ck2λ+ ck3)dλ
χ,
where
C = (At)−1, Ajm =
∮aj
λ3−mdλ
χ.
It follows from equation (` is an arbitrary path on Γ)∫τ̂ `
dω =
∫`
τ∗dω,
that
Ajm =
∮aj
λ3−mdλ
χ=
∮aj
τ∗1
(λ3−m
dλ
χ
)=
=
∮τ̂1aj
λ3−mdλ
χ= −
∮aj
λ3−mdλ
χ= −Ajm.

6 A. O. SMIRNOV, S. G. MATVEENKO, S. K. SEMENOV, AND E. G. SEMENOVA
Therefore A = −A and C = −C. Dealing similarly with integrals on bcycles, we obtain
(22) B = −B −K or

THREEPHASE FREAK WAVES 7
in the difference of a case of the twophase solution [43,45,46], where change of initial phase Z ledto trivial shift of the solution on plane XOZ, the amplitude of the threephase solution (17)) ofequations (1) depends on a choice of initial phase Wt0 + Z a little bit more complicated.
3. An example of threephase solution
Let us consider a spectral curve Γ3 = {χ, λ} of genus g = 3:(24) Γ3 : χ
2 = ((λ− λ0)4 − 2a2(λ− λ0)2 cos 2ϕ+ a4)((λ− λ0)4 − 2b2(λ− λ0)2 cos 2ϕ+ b4),where 0 < a < b, π/4 < ϕ < π/2.
Let us choose the basis of cycles on Γ3 as it is shown on fig. 1.
ℜ(λ− λ0)
ℑ(λ− λ0)
aeiϕ
be−iϕ
a3a1 a2
b1 b2 b3
Figure 1. Canonical basis of cycles on Γ3
It is easy to check that the antiholomorphic involution τ1 transforms the canonical basis of cyclesaccordingly to relations (12) with the matrix
K =
0 1 11 0 11 1 0
.Also there are three holomorphic involutions on Γ3:
τ0 : (χ, λ)→ (−χ, λ),τ2 : (χ, λ)→ (χ, 2λ0 − λ),τ3; (χ, λ)→ (a2b2(λ− λ0)−4χ, λ0 + ab(λ− λ0)−1).
As a corollary the curve Γ3 covers two following curves:
(1) Γ1 = Γ3/τ2 of genus g = 1
Γ1 : χ2+ = (t
2 − 2a2t cos 2ϕ+ a4)(t2 − 2b2t cos 2ϕ+ b4),(2) Γ2 = Γ3/(τ0τ2) of genus g = 2
Γ2 : χ2− = t(t
2 − 2a2t cos 2ϕ+ a4)(t2 − 2b2t cos 2ϕ+ b4),

8 A. O. SMIRNOV, S. G. MATVEENKO, S. K. SEMENOV, AND E. G. SEMENOVA
where t = (λ− λ0)2, χ+ = χ, χ− = (λ− λ0)χ, anddt
χ+=
2(λ− λ0)dλχ
,tdt
χ−=
2(λ− λ0)2dλχ
,dt
χ−=
2dλ
χ.
The curves Γ1 and Γ2 are shown on fig. 2 and 3, where t1 = b2e2iϕ, t2 = a
2e2iϕ.
t2
t∗2
t1
t∗1
ℑt
ℜt
a1
b1
Figure 2. The curve Γ1
t2
t∗2
t1
t∗1
ℑt
ℜta21
b21
b22
a22
Figure 3. The curve Γ2
The coverings generate the following transformations of cycles:a1a2a3
→ Sa1a21a22
+ Pb1b21b22
,b1b2b3
→ Qa1a21a22
+Rb1b21b22
,where
S =
−1 1 01 0 −11 0 1
, P =0 −2 00 0 2
0 0 −2
, Q =−1 1 00 0 −1
0 0 1
, R =0 0 01 1 1
1 1 −1
.Let us remind that these matrices should satisfy to relations
StQ = QtS, RtP = P tR, StR−QtP = nI,where I is identity matrix, n = 2 is number of covering sheets.
Due to involution τ3, the curve Γ2 covers two elliptic curves Γ± (fig. 4 and 5):
Γ± : ν2± = (s± 2ab)(s2 − 2(a2 + b2)s cos 2ϕ+ a4 + b4 + 2a2b2 cos 4ϕ),
where
s = t+a2b2
t, ν± =
t± abt2
χ−,ds
ν±=
(t∓ ab)dtχ−
.
The coverings of Γ2 on Γ± generate the following mappings of cycles(a21a22
)→(
1 1−1 1
)(a+a−
),
(b21b22
)→(
1 1−1 1
)(b+b−
),

THREEPHASE FREAK WAVES 9
b
b
ℜs
ℑs
b2 − a2
b2 + a2−2ab
s1
s∗1
b+
a+
Figure 4. The curve Γ+
b
b
ℜs
ℑs
b2 − a2
b2 + a22ab
s1
s∗1
b−
a−
Figure 5. The curve Γ−
As a result we havea1a2a3
→−1 1 11 1 −1
1 −1 1
a1a+a−
+0 −2 −20 −2 2
0 2 −2
b1b+b−
,(25)b1b2b3
→−1 1 10 1 −1
0 −1 1
a1a+a−
+0 0 01 0 2
1 2 0
b1b+b−
.(26)From (25), (26) and from the relations
dλ
χ=
1
4ab
ds
ν−− 1
4ab
ds
ν+,
λdλ
χ=
1
2
dt
χ++
λ04ab
ds
ν−− λ0
4ab
ds
ν+,
λ2dλ
χ= λ0
dt
χ++λ20 + ab
4ab
ds
ν−− λ
20 − ab4ab
ds
ν+
it follows that the matrices C and B equal
C =
c1 + c3 −2λ0(c1 + c3) (λ20 − ab)c1 + (λ20 + ab)c3c1 c2 − 2λ0c1 (λ20 − ab)c1 − λ0c2c3 c2 − 2λ0c3 (λ20 + ab)c3 − λ0c2
,B =
i(b1 + b3) ib1 − 1/2 ib3 − 1/2ib1 − 1/2 i(b1 + b2) ib2 − 1/2ib3 − 1/2 ib2 − 1/2 i(b2 + b3)
,

10 A. O. SMIRNOV, S. G. MATVEENKO, S. K. SEMENOV, AND E. G. SEMENOVA
where
c1 =1
2(α1 − 2β1), c2 =
1
2α2, c3 =
1
2(α3 − 2β3),
ib1 =α1
2(α1 − 2β1), ib2 =
β22α2
, ib3 =α3
2(α3 − 2β3),
α1 =1
2
∮a+
ds
ν+, α2 =
1
2
∮a1
dt
χ+, α3 =
1
2
∮a−
ds
ν−,
β1 =1
2
∮b+
ds
ν+, β2 =
1
2
∮b1
dt
χ+, β3 =
1
2
∮b−
ds
ν−.
From the structure of Bmatrix and from matrix version of Appel’s theorem [41] it follows thatthe Riemann theta function of curve Γ3 equals:
(27)
Θ(pB) = f(p̃1, p̃2, p̃3) == ϑ3(p̃1h1)ϑ3(p̃2h2)ϑ3(p̃3h3) + ϑ4(p̃1h1)ϑ1(p̃2h2)ϑ1(p̃3h3)++ ϑ1(p̃1h1)ϑ4(p̃2h2)ϑ1(p̃3h3) + ϑ1(p̃1h1)ϑ1(p̃2h2)ϑ4(p̃3h3),
where p̃j = pj + pj+1 − pj+2, pj+3 ≡ pj , hj = exp(−4bj).The functions ϑj(ph) are Jacoby elliptic theta functions [1]:
ϑ1(ph) = 2∞∑m=1
(−1)m−1h(m−1/2)2 sin[(2m− 1)πp],
ϑ2(ph) = 2∞∑m=1
h(m−1/2)2
cos[(2m− 1)πp],
ϑ3(ph) = 1 + 2∞∑m=1
hm2
cos(2mπp),
ϑ4(ph) = 1 + 2∞∑m=1
(−1)mhm2 cos(2mπp).
Using the reduced form of thetafunction (27) and values for vectors of periods, one obtains thefollowing formula for the square of absolute value of the threephase solution (18) of the focusingNLS equation (1)
(28) ψ2 = −4K20f(k1x+ κ1t+ δ1, k2z + δ2, k3x+ κ3t+ δ3)×× f(k1x+ κ1t− δ1, k2z − δ2, k3x+ κ3t− δ3)× {f(k1x+ κ1t, k2z, k3x+ κ3t)}−2,
where the function f(p̃1, p̃2, p̃3) is defined by equation (27), and
k1 = −4ic1, k2 = −8ic2, k3 = −4ic3,κ1 = 4k1(3λ
20 − ab+ (a2 + b2) cos(2ϕ)),
κ3 = 4k3(3λ20 + ab+ (a
2 + b2) cos(2ϕ)).
It follows from (27), (28) that for λ0 = 0 the amplitude of the constructed solution of NLS equation
(1) is a periodic function in z, and for λ0 = 0, ϕ =1
2arccos
( ±aba2 + b2
)it is a periodic function in
z and in t.

THREEPHASE FREAK WAVES 11
Let us recall that the threephase solution u(x, z, t) of the KPI equation (10) and the square of
amplitude ψH(x, t)2 of threephase solution of Hirota equation (9) can be constructed from (28)by relations u(x, z, t) = 2 ψ(x, z, t)2 and ψH(x, t)2 = ψ(x, t,−αt)2.
The threephase solution of KPI equation for ab = 1,√b/a = 1.3, ϕ = 0.3π at the different
moment of time t and for λ0 = 0 is presented on the fig. 69. The same solution for λ0 = k2/(4k1)one can see on the fig. 1013. It is easy to see all three phase of solution on figures 613. Twophase are shortwave and third phase is a longwave envelope. One can see also on fig. 69 that thesolution for λ0 = 0 is periodic in z, and that the longwave envelope move to the right side.
Figure 6. Threephase solution of KPI equation forλ0 = 0, t = 0
Figure 7. Threephase solution of KPI equation forλ0 = 0, t = 0.1
Figure 8. Threephase solution of KPI equation forλ0 = 0, t = 0.2
Figure 9. Threephase solution of KPI equation forλ0 = 0, t = 0.3

12 A. O. SMIRNOV, S. G. MATVEENKO, S. K. SEMENOV, AND E. G. SEMENOVA
Figure 10. Threephase solution of KPI equation forλ0 = k2/(4k1), t = 0
Figure 11. Threephase solution of KPI equation forλ0 = k2/(4k1), t = 0.1
Figure 12. Threephase solution of KPI equation forλ0 = k2/(4k1), t = 0.2
Figure 13. Threephase solution of KPI equation forλ0 = k2/(4k1), t = 0.3
The threephase solution of Hirota equation for ab = 1,√b/a = 1.3, ϕ = 0.3π, α = 0.1 and for
different values of λ0 is presented on fig. 1417. It is easy to see all three phase of solution only onfig. 15.
Concluding remarks
Let us remark that the considered solution can not be used for a degeneration to quasirationalone because odd number of phases. I.e. very popular solutions can not be obtained as a limit ofthreephase solution. Analyzing a family of quasirational solutions from [9,10,17,18] we make thefollowing conjecture.

THREEPHASE FREAK WAVES 13
Figure 14. Amplitude ofthreephase solution of Hirota equation for λ0 = 0
Figure 15. Amplitude ofthreephase solution of Hirota equation for λ0 = 4
Figure 16. Amplitude ofthreephase solution of Hirota equation for λ0 =k2/(4k1)
Figure 17. Amplitude ofthreephase solution of Hirota equation for λ0 =k2/(4k3)
(1) Quasirational rang n solution of NLS equation can be obtained as a limit of 2nphaseelliptic solution.
(2) Spectral curve of corresponding elliptic solution should be a 2sheet unramified coveringover spectral curve of nphase elliptic solution of KdV equation.
The rang 1 quasirational solution (Peregrine soliton) was obtained as limit of twophase ellipticone in [44], and for obtaining rang 2 quasirational solution we should use a spectral curve of genusg = 4.
Authors acknowledge Prof. V.B. Matveev for his support and discussion over this paper andquasirational solutions of NLS equation. This work performed within the framework of the state

14 A. O. SMIRNOV, S. G. MATVEENKO, S. K. SEMENOV, AND E. G. SEMENOVA
order of the Ministry of education of Russia, and partially supported by RFBR, research project140100589 a.
References
[1] N. I. Akhiezer, Elements of the theory of elliptic functions, American Mathematical Society, Providence, RI,1990, Translated from the second Russian edition by H. H. McFaden.
[2] N. Akhmediev and E. Pelinovsky (eds.), Discussion & debate: Rogue waves  towards a unifying concept?, Eur.
Phys. J. Special Topics, vol. 185, Springer, 2010.[3] N. N. Akhmediev and A. Ankiewicz, Solitons, nonlinear pulses and beams, CHAPMAN & HALL, 1997.
[4] A. Ankiewicz, J. M. SotoCrespo, and N. Akhmediev, Rogue waves and rational solutions of the Hirota equation,
Phys. Rev. E 81 (2010).[5] H. F. Baker, Abel’s theorem and the allied theory including the theory of theta functions, Cambridge, 1897.
[6] E. D. Belokolos, A. I. Bobenko, V. Z. Enol’skii, A. R. Its, and V. B. Matveev, Algebrogeometrical approach to
nonlinear evolution equations, Springer Ser. Nonlinear Dynamics, Springer, 1994.[7] R. Y. Chiao, E. Garmire, and C. H. Townes, Selftrapping of optical beams, Phys. Rev. Lett. 13 (1964), no. 15,
479–482.
[8] C. Q Dai and J. F. Zhang, New solitons for the Hirota equation and generalized higherorder nonlinearSchrödinger equation with variable coefficients, J. Phys. A 39 (2006), 723–737.
[9] P. Dubard and V. B. Matveev, Multirogue waves solutions: from the NLS equation to the KPI equation,
Nonlinearity 26 (2013), no. 12, R93–R125.[10] Ph. Dubard, Multirogue solutions to the focusing nls equation, Ph.D. Thesis, 2010.
[11] B. A. Dubrovin, Inverse problem for periodic finitezoned potentials in the theory of scattering, Funct. Anal.Appl. 9 (1975), no. 1, 61–62.
[12] , Periodic problems for the Kortewegde Vries equation in the class of finite band potentials, Funct. Anal.
Appl. 9 (1975), no. 3, 215–223.[13] , Theta functions and nonlinear equations, Russ. Math. Surv. 36 (1981), no. 2, 11–92.
[14] B. A. Dubrovin, V. B. Matveev, and S. P. Novikov, Nonlinear equations of Kortewegde Vries type, finiteband
linear operators and Abelian varietes, Russ. Math. Surv. 31 (1976), no. 1, 59–146.[15] B. A. Dubrovin and S. P. Novikov, A periodicity problem for the Kortewegde Vries and SturmLiouville equa
tions. Their connection with algebraic geometry, Sov. Math., Dokl. 15 (1974), 1597–1601.[16] J. D. Fay, Thetafunctions on Riemann surfaces, Lect. Notes in Math., vol. 352, Springer, 1973.
[17] P. Gaillard, Deformations of thirdorder Peregrine breather solutions of the nonlinear Schrödinger equation
with four parameters, Phys. Rev. E 88 (2013).[18] , Tenparameter deformations of the sixthorder Peregrine breather solutions of the NLS equation, Phys
ica Scripta 89 (2014).
[19] F. Gesztesy and H. Holden, Soliton equation and their algebrogeometric solutions: Vol. 1, (1+1)dimensionalcontinuous models, Cambridge Stud. in Adv. Math., vol. 79, Cambridge University Press, 2003.
[20] F. Gesztesy, H. Holden, J. Michor, and G. Teschl, Soliton equation and their algebrogeometric solutions: Vol.
2, (1+1)dimensional discrete models, Cambridge Stud. in Adv. Math., vol. 114, Cambridge University Press,2008.
[21] B. Guo, L. Ling, and Q. P. Liu, Nonlinear Schrödinger equation: Generalized Darboux transformation and rogue
wave solutions, Phys. Rev. E 85 (2012).[22] A. R. Its, “isomonodromy” solutions of equations of zero curvature, Mathematics of the USSRIzvestiya 26
(1986), no. 3, 497–529.[23] A. R. Its and V. P. Kotlyarov, On a class of solutions of the nonlinear Schrödinger equation, Dokl. Akad. Nauk
Ukrain. SSR, Ser. A 11 (1976), 965–968, (Russian).
[24] A. R. Its and V. B. Matveev, Hill’s operator with finitely many gaps, Funct. Anal. Appl. 9 (1975), no. 1, 65–66.[25] , Schrödinger operators with finitegap spectrum and nsoliton solutions of the KdV equation, Theor.
Math. Phys. 23 (1975), no. 1, 343–355.
[26] A.R. Its, Inversion of hyperelliptic integrals and integration of nonlinear differential equations, Vestn. Leningr.Univ. (Mat. Mekh. Astron.) 7 (1976), no. 2, 39–46, (Russian).
[27] C. Kalla and C. Klein, New construction of algebrogeometric solutions to the CamassaHolm equation andtheir numerical evaluation, Proc. R. Soc. A 468 (2012), no. 2141, 1371–1390.

THREEPHASE FREAK WAVES 15
[28] , On the numerical evaluation of algebrogeometric solutions to integrable equations, Nonlinearity 25(2012), no. 3, 569–596.
[29] A. Krazer, Lehrbuch der thetafunktionen, Teubner, Leipzig, 1903.
[30] I. M. Krichever, Methods of algebraic geometry in the theory of nonlinear equations, Russ. Math. Surv. 32(1977), no. 6, 185–213.
[31] A. Kundu, A. Mukherjee, and T. Naskar, Modelling rogue waves through exact dynamical lump soliton controlledby ocean currents, Preprint, arXiv:1204.0916, 2012, 5p.
[32] E. A. Kuznetsov, Solitons in a parametricaly unstable plasma, Sov. Phys., Dokl. 22 (1977), 507–508.
[33] P. D. Lax, Periodic solutions of the KdV equations, Lect. in Appl. Math. 15 (1974), 85–96.[34] Ch. Zh. Li and J. S. He, Darboux transformation and positons of the inhomogeneous Hirota and the Maxwell
Bloch equation, Preprint, arXiv:1210.2501, 2012, 5p.
[35] V. A. Marchenko, The periodic Kortewegde Vries problem, Math. USSR Sb. 24 (1974), no. 3, 319–344.[36] V. B. Matveev, 30 years of finitegap integration theory, Phil. Trans. R. Soc. A 366 (2008), 837–875.
[37] H. P. McKean and P. van Moerbeke, The spectrum of Hill’s equation, Invent. Math. 30 (1975), 217–274.
[38] D. Mumford, Tata lectures on theta. II, Progress in Math., vol. 43, Birkhäuser Boston Inc., Boston, MA, 1984.[39] S. P. Novikov, The periodic problem for the Kortewegde Vries equation. i, Funct. Anal. Appl. 8 (1974), no. 3,
236–246.
[40] E. Previato, Hyperelliptic quasiperiodic and soliton solutions of the nonlinear Schrödinger equation, DukeMath. J. 52 (1985), no. 2, 281–545.
[41] A. O. Smirnov, A matrix analogue of the Appell theorem and reduction of multidimensional Riemann theta
functions, Math. USSR Sb. 61 (1988), no. 2, 379–388.[42] , Elliptic solutions of the nonlinear Schrödinger equation and the modified Kortewegde Vries equation,
Russ. Acad. Sci. Sb. Math. 82 (1995), no. 2, 461–470.[43] , Solution of a nonlinear Schrödinger equation in the form of twophase freak waves, Theor. Math. Phys.
173 (2012), no. 1, 1403–1416.
[44] , The elliptic breather for the nonlinear Schrödinger equation, J. Math. Sci. 192 (2013), no. 1, 117–125.[45] , Periodic twophase “rogue waves”, Math. Notes 94 (2013), no. 6, 897–907.
[46] A. O. Smirnov, E. G. Semenova, V. Zinger, and N. Zinger, On a periodic solution of the focusing nonlinear
Schrödinger equation, Preprint, arXiv:1407.7974, 2014, 24p.[47] L. H. Wang, K. Porsezian, and J. S. He, Breather and rogue wave solutions of a generalized nonlinear Schrödinger
equation, Preprint, arXiv:1304.8085, 2013, 21p.[48] Z. Yan, Financial rogue waves appearing in the coupled nonlinear volatility and option pricing model, Preprint,
arXiv:1101.3107, 2011.
[49] V. E. Zakharov, Stability of periodic waves of finite amplitude on the surface of a deep fluid, J. Appl. Mech.Tech. Phys. 9 (1968), no. 2, 190–194.
Email address: [email protected]
Email address: [email protected]
Email address: [email protected]
Email address: [email protected]
St.Petersburg State University of Aerospace Instrumentation (SUAI), 67 Bolshaya Morskaya str.,
St.Petersburg, 190000,Russia
http://arxiv.org/abs/1204.0916http://arxiv.org/abs/1210.2501http://arxiv.org/abs/1407.7974http://arxiv.org/abs/1304.8085http://arxiv.org/abs/1101.3107mailto:[email protected]:[email protected]:[email protected]:[email protected]
Introduction1. Finitegap multiphase solutions of the NLS equation2. Features of threephase solutions3. An example of threephase solutionConcluding remarksReferences