## Mach-Type Soliton in the Novikov-Veselov Equation

Jen-Hsu CHANG

Department of Computer Science and Information Engineering, National Defense University, Tauyuan County 33551, Taiwan E-mail: jhchang@ndu.edu.tw

Received September 18, 2014, in final form December 10, 2014; Published online December 18, 2014 http://dx.doi.org/10.3842/SIGMA.2014.111

Abstract. Using the reality condition of the solutions, one constructs the Mach-type soliton of the Novikov-Veselov equation by the minor-summation formula of the Pfaffian. We study the evolution of the Mach-type soliton and find that the amplitude of the Mach stem wave is less than two times of the one of the incident wave. It is shown that the length of the Mach stem wave is linear with time. One discusses the relations with V -shape initial value wave for different critical values of Miles parameter.

Key words: Pfaffian; Mach-type soliton; Mach stem wave; V -shape wave

2010 Mathematics Subject Classification: 35C08; 35A22

## 1 Introduction

Recently, the resonance theory of line solitons of KP-(II) equation (shallow water wave equation)

$$\partial _ { x } ( - 4 u _ { t } + u _ { x x x } + 6 u u _ { x } ) + 3 u _ { y y } = 0$$

has attracted much attractions using the totally non-negative Grassmannian [1, 3, 5, 6, 16, 22, 23], that is, those points of the real Grassmannian whose Plucker coordinates are all nonnegative. For the KP-(II) equation case, the τ -function is described by the Wronskian form with respect to x . The Mach reflection problem describes the resonant interaction of solitary waves appearing in the reflection of an obliquely incident waves onto a vertical wall. John Miles discussed an oblique interaction of solitary waves and found a resonant interaction to describe the Mach reflection phenomena [26]. In this end, he predicts an extraordinary fourfold application of the wave at the wall. The Miles theory in terms of the KP equation and the Mach-type solution in KP observed experimentally are investigated in [16, 17, 18, 20, 31] (and references therein). The point is that irregular reflection can be described by the (3142)-type soliton and the stem in the middle part should be a Mach stem wave. Inspired by their works, one can consider the Novikov-Veselov equation similarly.

One considers the Novikov-Veselov (NV) equation [4, 11, 15, 28, 30] with real solution U :

$$U _ { t } = \text{Re} \left [ \partial _ { z } ^ { 3 } U + 3 \partial _ { z } ( Q U ) - 3 \epsilon \partial _ { z } Q \right ], \\ \bar { \partial } _ { z } Q = \partial _ { z } U, \quad \ t \in R,$$

where glyph[epsilon1] is a real constant. The NV equation (1.1) is one of the natural generation of the famous KdV equation and can have the Manakov's triad representation [24]

$$L _ { t } = [ A, L ] + B L,$$

where L is the two-dimension Schr¨dinger operator o

$$L = \partial _ { z } \bar { \partial } _ { z } + U - \epsilon$$

and

$$A = \partial _ { z } ^ { 3 } + Q \partial _ { z } + \bar { \partial } _ { z } ^ { 3 } + \bar { Q } \bar { \partial } _ { z }, \quad B = Q _ { z } + \bar { Q } _ { \bar { z } }.$$

It is equivalent to the linear representation

$$L \phi = 0, \quad \partial _ { t } \phi = A \phi.$$

We remark that when glyph[epsilon1] →±∞ , the Veselov-Novikov equation reduces to the KP-I ( glyph[epsilon1] →-∞ ) and KP-(II) ( glyph[epsilon1] →∞ ) equation respectively [10]. To make a comparison with KP-(II) equation, we only consider glyph[epsilon1] &gt; 0.

Let φ 1 , φ 2 be any two independent solutions of (1.2). Then one can construct the extended Moutard transformation using the skew product [2, 13, 25]

$$W ( \phi _ { 1 }, \phi _ { 2 } ) & = \int ( \phi _ { 1 } \partial \phi _ { 2 } - \phi _ { 2 } \partial \phi _ { 1 } ) d z - ( \phi _ { 1 } \bar { \partial } \phi _ { 2 } - \phi _ { 2 } \bar { \partial } \phi _ { 1 } ) d \bar { z } + [ \phi _ { 1 } \partial ^ { 3 } \phi _ { 2 } - \phi _ { 2 } \partial ^ { 3 } \phi _ { 1 } \\ & \quad + \phi _ { 2 } \bar { \partial } ^ { 3 } \phi _ { 1 } - \phi _ { 1 } \bar { \partial } ^ { 3 } \phi _ { 2 } + 2 ( \partial ^ { 2 } \phi _ { 1 } \partial \phi _ { 2 } - \partial \phi _ { 1 } \partial ^ { 2 } \phi _ { 2 } ) - 2 ( \bar { \partial } ^ { 2 } \phi _ { 1 } \bar { \partial } \phi _ { 2 } - \bar { \partial } \phi _ { 1 } \bar { \partial } ^ { 2 } \phi _ { 2 } ) \\ & \quad + 3 Q ( \phi _ { 1 } \partial \phi _ { 2 } - \phi _ { 2 } \partial \phi _ { 1 } ) - 3 \bar { Q } ( \phi _ { 1 } \bar { \partial } \phi _ { 2 } - \phi _ { 2 } \bar { \partial } \phi _ { 1 } ) ] d t,$$

such that

$$\hat { U } ( t, z, \bar { z } ) = U ( t, z, \bar { z } ) + 2 \partial \bar { \partial } \ln W ( \phi _ { 1 }, \phi _ { 2 } ), \quad \hat { Q } ( t, z, \bar { z } ) = Q ( t, z, \bar { z } ) + 2 \partial \partial \ln W ( \phi _ { 1 }, \phi _ { 2 } ),$$

is also a solution of the NV equation (1.1).

For fixed potential U 0 ( z, z, t ¯ ) and Q 0 ( z, z, t ¯ ) of the NV equation (1.1), we can take any 2 N wave functions φ , φ 1 2 , φ 3 , . . . , φ 2 N (or their linear combinations) of (1.2). Then the 2 N -step successive extended Moutard transformation can be expressed as the Pfaffian form [2, 27] (also see [12, 29])

$$U = U _ { 0 } + 2 \partial \bar { \partial } [ \ln \mathbf P f ( \phi _ { 1 }, \phi _ { 2 }, \phi _ { 3 }, \dots, \phi _ { 2 N } ) ], \quad Q = Q _ { 0 } + 2 \partial \partial [ \ln \mathbf P f ( \phi _ { 1 }, \phi _ { 2 }, \phi _ { 3 }, \dots, \phi _ { 2 N } ) ],$$

where Pf φ ,φ , φ , . . . , φ ( 1 2 3 2 N ) is the Pfaffian defined by

$$\text{Pf} ( \phi _ { 1 }, \phi _ { 2 }, \phi _ { 3 }, \dots, \phi _ { 2 N } ) = \sum _ { \sigma } \epsilon ( \sigma ) W _ { \sigma _ { 1 } \sigma _ { 2 } } W _ { \sigma _ { 3 } \sigma _ { 4 } } \cdots W _ { \sigma _ { 2 N - 1 } \sigma _ { 2 N } },$$

and W σ σ i j = W φ ( σ i ( ) , φ σ j ( ) ) is the extended Moutard transformation (1.3), σ being some permutations.

To construct the N -solitons solutions, we take V = U = 0 in (1.1) and then (1.2) becomes

$$\partial \bar { \partial } \phi = \epsilon \phi, \quad \phi _ { t } = \phi _ { z z z } + \phi _ { z z z z },$$

where glyph[epsilon1] is non-zero real constant. The general solution of (1.4) can be expressed as

$$\phi ( z, \bar { z }, t ) = \int _ { \Gamma } e ^ { ( i \lambda ) z + ( i \lambda ) ^ { 3 } t + \frac { c } { i \lambda } \bar { z } + \frac { c } { ( i \lambda ) ^ { 3 } } t } \nu ( \lambda ) d \lambda,$$

where ν λ ( ) is an arbitrary distribution and Γ is an arbitrary path of integration such that the r.h.s. of (1.5) is well defined. One takes ν m ( λ ) = δ λ ( -p m ), where p m is a complex number. Define

$$\phi _ { m } = \frac { \phi ( p _ { m } ) } { \sqrt { 3 } } = \frac { 1 } { \sqrt { 3 } } e ^ { F ( p _ { m } ) },$$

where

$$F ( \lambda ) = ( i \lambda ) z + ( i \lambda ) ^ { 3 } t + \frac { \epsilon } { i \lambda } \bar { z } + \frac { \epsilon ^ { 3 } } { ( i \lambda ) ^ { 3 } } t.$$

Plugging ( φ m , φ n ) into the extended Moutard transformation (1.3), we obtain

$$W ( \phi _ { m }, \phi _ { n } ) = i \frac { p _ { n } - p _ { m } } { p _ { n } + p _ { m } } e ^ { F ( p _ { m } ) + F ( p _ { n } ) }.$$

$$W ( \phi _ { m }, \phi _ { n } ) = i \frac { p _ { n } - p _ { m } } { p _ { n } + p _ { m } } e ^ { F ( \mathcal { I } }$$

To study resonance, we introduce the real Grassmannian (or the 2 N × M matrix) to construct N solitons. To this end, one considers linear combination of φ n . Let

$$\vec { \Psi } = ( \phi _ { 1 }, \phi _ { 2 }, \phi _ { 3 }, \dots, \phi _ { M } ) ^ { T }$$

and H be an 2 N × M (2 N ≤ M ) of real constant matrix (or Grassmannian). Suppose that

$$H \vec { \Psi } = \vec { \Psi } ^ { * } = ( \Psi _ { 1 } ^ { * }, \Psi _ { 2 } ^ { * }, \Psi _ { 3 } ^ { * }, \dots, \Psi _ { 2 N } ^ { * } ) ^ { T },$$

that is,

$$\Psi _ { n } ^ { * } = h _ { n 1 } \phi _ { 1 } + h _ { n 2 } \phi _ { 2 } + h _ { n 3 } \phi _ { 3 } + \dots + h _ { n M } \phi _ { M }, \quad 1 \leq n \leq 2 N.$$

Then one has by the minor-summation formula [14, 19]

$$\tau _ { N } = \Pr f ( \Psi _ { 1 } ^ { * }, \Psi _ { 2 } ^ { * }, \Psi _ { 3 } ^ { * }, \dots, \Psi _ { 2 N } ^ { * } ) = \Pr f \left ( H W _ { M } H ^ { T } \right ) = \sum _ { I \subset [ M ], \, \dots I = 2 N } \Pr f \left ( H _ { I } ^ { I } \right ) \det ( H _ { I } ), \quad ( 1. 7 )$$

where the M × M matrix W M is defined by the element (1.6) and H I I denote the 2 N × M submatrix of H obtained by picking up the rows and columns indexed by the same index set I . By this formula, the resonance of real solitons of the Novikov-Veselov equation can be investigated just like the resonance theory of KP-(II) equation [16, 21, 23]. Finally, the N -solitons solutions are defined by [7, 9]

$$U ( z, \bar { z }, t ) = 2 \partial \bar { \partial } \ln \tau _ { N } ( z, \bar { z }, t ), \quad V ( z, \bar { z }, t ) = 2 \partial \partial \ln \tau _ { N } ( z, \bar { z }, t ).$$

To obtain the real potential U , the following reality conditions [9] for resonance have to be considered

$$| p _ { k } | ^ { 2 } = | q _ { k } | ^ { 2 } = \epsilon > 0, \quad k = 1, 2, 3, \dots, m,$$

given m pairs of complex numbers ( p , q 1 1 ) , ( p , q 2 2 ) , . . . , ( p m , q m ). Letting p m = √ glyph[epsilon1]e iα m and removing i factor from (1.6) afterwards , one has

$$W ( \phi _ { m }, \phi _ { n } ) = - \tan \frac { \alpha _ { n } - \alpha _ { m } } { 2 } e ^ { \phi _ { m n } },$$

where

$$\phi _ { m n } & = F ( p _ { m } ) + F ( p _ { n } ) = - 2 \sqrt { \epsilon } [ x ( \sin \alpha _ { m } + \sin \alpha _ { n } ) + y ( \cos \alpha _ { m } + \cos \alpha _ { n } ) ] \\ & \quad + 2 t \epsilon \sqrt { \epsilon } ( \sin 3 \alpha _ { m } + \sin 3 \alpha _ { n } ).$$

Therefore, given a 2 N × M matrix H , the associated τ H -function can be written as by (1.7) [8]

$$\tau _ { H } = \sum _ { I \subset [ M ], \, \natural I = 2 N } \Gamma _ { I } \Lambda _ { I } ( x, y, t ),$$

where

$$\Lambda _ { I } ( x, y, t ) = \text{Pf} ( W _ { 2 N } ) = ( - 1 ) ^ { N } \left ( \prod _ { i = 2, \, i > j } ^ { 2 N } \tan \frac { \alpha _ { i } - \alpha _ { j } } { 2 } \right ) e ^ { \sum _ { m = 1 } ^ { 2 N } F ( p _ { m } ) },$$

Γ I being the 2 N × 2 N minor for the columns with the index set I = { i 1 , i 2 , i 3 , . . . , i 2 N } . Also, to keep τ H totally positive (or totally negative), we assume that the matrix H belongs to the totally non-negative Grassmannian [22, 23] and the angle α n satisfies the following condition:

$$- \frac { \pi } { 2 } \leq \alpha _ { 1 } < \alpha _ { 2 } < \alpha _ { 3 } < \dots < \alpha _ { M - 1 } < \alpha _ { M } \leq \frac { \pi } { 2 }.$$

For one-soliton solution, we have, -π 2 ≤ α &lt; α i j &lt; α k ≤ π 2

$$\text{or one-solution, we have, } - \frac { \frac { } { 2 } \geq \alpha _ { i } < \alpha _ { j } < \alpha _ { k } \geq \frac { } { 2 } } { \tau _ { 1 } } & = \tan \frac { \alpha _ { i } - \alpha _ { j } } { 2 } e ^ { \phi _ { i j } } + a \tan \frac { \alpha _ { i } - \alpha _ { k } } { 2 } e ^ { \phi _ { i k } } \\ & = a e ^ { \phi _ { i k } } \tan \frac { \alpha _ { i } - \alpha _ { k } } { 2 } \left [ 1 + \frac { 1 } { a } \frac { \tan \frac { \alpha _ { i } - \alpha _ { j } } { 2 } } e ^ { F ( p _ { j } ) - F ( p _ { k } ) } \right ] \\ & = a e ^ { \phi _ { i k } } \tan \frac { \alpha _ { i } - \alpha _ { k } } { 2 } \left [ 1 + e ^ { F ( p _ { j } ) - F ( p _ { k } ) + \theta _ { j k } } \right ], \\ \text{where } a \text{ is a constant and the phase shift}$$

where a is a constant and the phase shift

$$\theta _ { j k } = \ln \frac { 1 } { a } \frac { \tan \frac { \alpha _ { i } - \alpha _ { j } } { 2 } } { \tan \frac { \alpha _ { i } - \alpha _ { k } } { 2 } } = \ln \frac { \tan \frac { \alpha _ { i } - \alpha _ { j } } { 2 } } { \tan \frac { \alpha _ { i } - \alpha _ { k } } { 2 } } - \ln a.$$

Hence the real one-soliton solution is [8]

$$\text{Hence the real one-soliton solution is } [ 8 ] \\ U & = 2 \partial _ { z } \partial _ { \bar { z } } \ln a e ^ { \phi _ { k } } \tan \frac { \alpha _ { i } - \alpha _ { k } } { 2 } \left [ 1 + e ^ { F ( p _ { j } ) - F ( p _ { k } ) + \theta _ { j k } } \right ] = 2 \partial _ { z } \partial _ { \bar { z } } \left [ 1 + e ^ { F ( p _ { j } ) - F ( p _ { k } ) + \theta _ { j k } } \right ] \\ & = \frac { 1 } { 2 } | p _ { k } - p _ { j } | ^ { 2 } \, \text{sech} ^ { 2 } \left [ \frac { F ( p _ { j } ) - F ( p _ { k } ) + \theta _ { j k } } { 2 } \right ] \\ & = 2 \epsilon \sin ^ { 2 } \left ( \frac { \alpha _ { k } - \alpha _ { j } } { 2 } \right ) \text{sech} ^ { 2 } \left [ \frac { F ( p _ { j } ) - F ( p _ { k } ) + \theta _ { j k } } { 2 } \right ] \\ & = A _ { \left [ j, k \right ] } \, \text{sech} ^ { 2 } \frac { 1 } { 2 } ( \vec { K } _ { \left [ j, k \right ] } \cdot \vec { x } - \Omega _ { \left [ j, k \right ] } t + \theta _ { j k } ). \\ \text{From } ( 1. 9 ) \, \text{the amplitude } A _ { \left [ i, k \right ] }, \, \text{the wave vector } \vec { K } _ { \left [ i, k \right ] } \, \text{and the frequency } \Omega _ { \left [ i, k \right ] } \, \text{are defined by}$$

From (1.9) the amplitude A [ j,k ] , the wave vector glyph[vector] K [ j,k ] and the frequency Ω [ j,k ] are defined by

$$\cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cd. \\ A _ { [ j, k ] } & = 2 \epsilon \sin ^ { 2 } \left ( \frac { \alpha _ { k } - \alpha _ { j } } { 2 } \right ), \\ \vec { \mathbf K } _ { [ j, k ] } & = 2 \sqrt { \epsilon } ( - \sin \alpha _ { j } + \sin \alpha _ { k }, - \cos \alpha _ { j } + \cos \alpha _ { k } ), \\ \Omega _ { [ j, k ] } & = 2 \epsilon \sqrt { \epsilon } [ - \sin 3 \alpha _ { j } + \sin 3 \alpha _ { k } ], \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{array}$$

The direction of the wave vector glyph[vector] K [ j,k ] = ( K x [ j,k ] , K y [ j,k ] ) is measured in the clockwise sense from the y -axis and it is given by

$$\frac { K _ { [ j, k ] } ^ { y } } { K _ { [ j, k ] } ^ { x } } = \frac { - \cos \alpha _ { j } + \cos \alpha _ { k } } { - \sin \alpha _ { j } + \sin \alpha _ { k } } = - \tan \frac { \alpha _ { j } + \alpha _ { k } } { 2 },$$

that is, α j + α k 2 gives the angle between the line soliton and the y -axis in the clockwise sense. In addition, the soliton velocity V [ j k , ] is [8]

$$\mathbf V _ { [ j, k ] } = \frac { \epsilon } { 4 } \frac { \sin 3 \alpha _ { k } - \sin 3 \alpha _ { j } } { \sin ^ { 2 } \frac { \alpha _ { j } - \alpha _ { k } } { 2 } } \left ( \sin \alpha _ { k } - \sin \alpha _ { j }, \cos \alpha _ { k } - \cos \alpha _ { j } \right ).$$

The paper is organized as follows. In Section 2, one investigates Mach-type or (3142)-type soliton for the Novikov-Veselov equation. One shows the evolution of the Mach-type soliton and obtains the relation of the amplitude of the Mach stem wave ([1,4]-soliton) with the one of the incident wave ([1,3]-soliton). Furthermore, the length of the Mach stem wave is linear with time. In Section 3, we discuss the relations with V -shape initial value wave for different critical value of Miles parameter κ . It is shown that the amplitude of the Mach stem wave is less than two times of the one of the incident wave. In Section 4, we conclude the paper with several remarks.

## 2 Mach type soliton

In this section, we investigate the Mach-type or (3142)-type soliton. The corresponding totally non-negative Grassmannian is the the matrix [16]

$$H _ { M } = \begin{bmatrix} 1 & a & 0 & - c \\ 0 & 0 & 1 & b \end{bmatrix},$$

where a , b , c are positive numbers. When c = 0, one has the O -type soliton for the NovikovVeselov equation. For V -shape initial value wave, one can introduce parameter κ to determine the evolution into Mach-type or O -type soliton (see next section). We remark that the Y -shape, O -type, and P -type solitons for the Novikov-Veselov equation are investigated in [8]. Now,

$$H _ { M } \frac { 1 } { \sqrt { 3 } } [ \phi ( p _ { 1 } ), \phi ( p _ { 2 } ), \phi ( p _ { 3 } ), \phi ( p _ { 4 } ) ] ^ { T } = \frac { 1 } { \sqrt { 3 } } \begin{bmatrix} \phi ( p _ { 1 } ) + a \phi ( p _ { 2 } ) \\ \phi ( p _ { 3 } ) + b \phi ( p _ { 4 } ) \end{bmatrix} = \begin{bmatrix} \Psi _ { 1 } ^ { * } \\ \Psi _ { 2 } ^ { * } \end{bmatrix}.$$

A direct calculation yields by (1.7), (1.8) and (1.9)

$$A \text{ direct calculation yields by } ( 1. 7 ), ( 1. 8 ) \text{ and } ( 1. 9 ) \\ \tau _ { M } = W ( \Psi _ { 1 } ^ { * }, \Psi _ { 2 } ^ { * } ) = W ( \phi _ { 1 }, \phi _ { 3 } ) + b W ( \phi _ { 1 }, \phi _ { 4 } ) + a W ( \phi _ { 2 }, \phi _ { 3 } ) + a b W ( \phi _ { 2 }, \phi _ { 4 } ) + c W ( \phi _ { 3 }, \phi _ { 4 } ) \\ = \tan \frac { \alpha _ { 1 } - \alpha _ { 3 } } { 2 } e ^ { - 2 \sqrt { \left [ z ( \sin \alpha _ { 1 } + s \sin \alpha _ { 3 } ) + y ( \cos \alpha _ { 1 } + \cos \alpha _ { 3 } ) ] + 2 t e \sqrt { \left ( s \sin 3 \alpha _ { 1 } + s \sin 3 \alpha _ { 3 } } ) } \\ + b \tan \frac { \alpha _ { 1 } - \alpha _ { 4 } } { 2 } e ^ { - 2 \sqrt { \left [ z ( \sin \alpha _ { 1 } + s \sin \alpha _ { 4 } ) + y ( \cos \alpha _ { 1 } + \cos \alpha _ { 4 } ) ] + 2 t e \sqrt { \left ( s \sin 3 \alpha _ { 1 } + s \sin 3 \alpha _ { 4 } } ) } \\ + a \tan \frac { \alpha _ { 2 } - \alpha _ { 3 } } { 2 } e ^ { - 2 \sqrt { \left [ z ( \sin \alpha _ { 2 } + s \sin \alpha _ { 3 } ) + y ( \cos \alpha _ { 2 } + \cos \alpha _ { 3 } ) ] + 2 t e \sqrt { \left ( s \sin 3 \alpha _ { 2 } + s \sin 3 \alpha _ { 3 } } ) } \\ + a b \tan \frac { \alpha _ { 2 } - \alpha _ { 4 } } { 2 } e ^ { - 2 \sqrt { \left [ z ( \sin \alpha _ { 2 } + s \sin \alpha _ { 4 } ) + y ( \cos \alpha _ { 2 } + \cos \alpha _ { 4 } ) ] + 2 t e \sqrt { \left ( s \sin 3 \alpha _ { 2 } + s \sin 3 \alpha _ { 4 } } ) } \\ + c \tan \frac { \alpha _ { 3 } - \alpha _ { 4 } } { 2 } e ^ { - 2 \sqrt { \left [ z ( \sin \alpha _ { 3 } + s \sin \alpha _ { 4 } ) + y ( \cos \alpha _ { 3 } + \cos \alpha _ { 4 } ) ] + 2 t e \sqrt { \left ( s \sin 3 \alpha _ { 3 } + s \sin 3 \alpha _ { 4 } } ) },$$

where

$$- \frac { \pi } { 2 } \leq \alpha _ { 1 } < \alpha _ { 2 } < \alpha _ { 3 } < \alpha _ { 4 } \leq \frac { \pi } { 2 }.$$

To investigate the asymptotic behavior for | y | → ∞ , we use the notation [16], considering the line x = -cy ,

$$\eta _ { m } ( c ) = - c \sin \alpha _ { m } + \cos \alpha _ { m }.$$

When η m ( c ) = η n ( c ), one gets

$$c = \frac { \cos \alpha _ { m } - \cos \alpha _ { n } } { \sin \alpha _ { m } - \sin \alpha _ { n } } = - \tan \frac { \alpha _ { m } + \alpha _ { n } } { 2 }.$$

Since

$$[ \eta _ { m } ( c ) - \eta _ { i } ( c ) ] \Big | _ { c = - \tan \frac { \alpha _ { i } + \alpha _ { j } } { 2 } } & = \cos \alpha _ { m } - \cos \alpha _ { i } + \tan \frac { \alpha _ { i } + \alpha _ { j } } { 2 } ( \sin \alpha _ { m } - \sin \alpha _ { i } ) \\ & = ( \sin \alpha _ { m } - \sin \alpha _ { i } ) \left [ \tan \frac { \alpha _ { i } + \alpha _ { j } } { 2 } - \tan \frac { \alpha _ { i } + \alpha _ { m } } { 2 } \right ],$$

we have the following order relations among the other η m ( c ) ′ s along c = -tan α i + α j 2

$$\begin{cases} \eta _ { i } = \eta _ { j } < \eta _ { m } & \text{if} \quad i < m < j, \\ \eta _ { i } = \eta _ { j } > \eta _ { m } & \text{if} \quad m < i \quad \text{or} \quad m > j. \end{cases}$$

Then by a similar argument in [16], one knows that by (1.10):

(a) For y glyph[greatermuch] 0, there are two unbounded line solitons, whose types from left to right are

$$[ 1, 3 ], \quad [ 3, 4 ].$$

(b) For y glyph[lessmuch] 0, there are two unbounded line solitons, whose types from left to right are

$$[ 4, 2 ], \quad [ 2, 1 ].$$

It can be verified by the Maple software.

Now, we can discuss the relations between the parameters a , b , c and phase shifts of these line solitons. Let us first consider the line solitons in x glyph[greatermuch] 0. There are two solitons which are [3,4]soliton and [2,1]-soliton. The [3,4]-soliton is obtained by the balance between the exponential terms W φ ,φ ( 1 3 ) and bW φ ,φ ( 1 4 ), and the [2,1]-soliton is obtained by the balance between the exponential terms W φ ,φ ( 1 3 ) and aW φ ,φ ( 2 3 ). Therefore, the phase shifts of [3,4]-soliton and [2,1]-soliton for x glyph[greatermuch] 0 are given by

$$\theta _ { [ 3, 4 ] } ^ { + } = \ln \frac { \tan \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 } } { \tan \frac { \alpha _ { 4 } - \alpha _ { 1 } } { 2 } } - \ln b, \quad \theta _ { [ 2, 1 ] } ^ { + } = \ln \frac { \tan \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 } } { \tan \frac { \alpha _ { 3 } - \alpha _ { 2 } } { 2 } } - \ln a.$$

For the line solitons in x glyph[lessmuch] 0, there are two solitons, which are [1,3]-soliton and [4,2]-soliton. The [1,3]-soliton is obtained by the balance between the exponential terms cW φ ,φ ( 3 4 ) and bW φ ,φ ( 1 4 ), and the [4,2]-soliton is obtained by the balance between the exponential terms cW φ ,φ ( 3 4 ) and aW φ ,φ ( 2 3 ). Therefore, the phase shifts of [1,3]-soliton and [4,2]-soliton for x glyph[lessmuch] 0 are given by

$$\theta _ { [ 1, 3 ] } ^ { - } = \ln \frac { \tan \frac { \alpha _ { 4 } - \alpha _ { 1 } } { 2 } } { \tan \frac { \alpha _ { 4 } - \alpha _ { 3 } } { 2 } } + \ln \frac { b } { c }, \quad \theta _ { [ 4, 2 ] } ^ { - } = \ln \frac { \tan \frac { \alpha _ { 3 } - \alpha _ { 2 } } { 2 } } { \tan \frac { \alpha _ { 4 } - \alpha _ { 3 } } { 2 } } + \ln \frac { a } { c }.$$

So one can see that

$$\theta _ { [ 1, 3 ] } ^ { - } + \theta _ { [ 3, 4 ] } ^ { + } = \theta _ { [ 4, 2 ] } ^ { - } + \theta _ { [ 2, 1 ] } ^ { + } = \text{total phase shift} = \ln \frac { \tan \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 } } { \tan \frac { \alpha _ { 4 } - \alpha _ { 3 } } { 2 } } - \ln c.$$

We define the parameter s (representing the total phase shift) by s = e -θ -[4 2] , -θ + [2 1] , , which leads to

$$a = \frac { \tan \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 } } { \tan \frac { \alpha _ { 3 } - \alpha _ { 2 } } { 2 } } s e ^ { \theta _ { [ 4, 2 ] } ^ { - } }, \quad b = \frac { \tan \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 } } { \tan \frac { \alpha _ { 4 } - \alpha _ { 1 } } { 2 } } s e ^ { \theta _ { [ 1, 3 ] } ^ { - } }, \quad c = \frac { \tan \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 } } { \tan \frac { \alpha _ { 4 } - \alpha _ { 3 } } { 2 } } s.$$

Hence we know that the three parameters a , b , c can be used to determine the locations of three asymptotic line solitons, that is, two in x glyph[lessmuch] 0 and one in x &gt; 0. The s -parameter represents the

Figure 1. The middle portion, having maximum amplitude, is the [1,4]-soliton (stem wave). The y -axis is slightly enlarged to make the middle portion longer.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0588v3-figure-102edd81-ec08-4285-a04b-91310f13ca5a.jpg)

Figure 2. α 1 = -23 50 π , α 2 = -1 5 π , α 3 = 1 5 π , α 4 = 23 50 π , a = b = c = 1, glyph[epsilon1] = 5.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0588v3-figure-8ef6310e-237b-4e90-9d63-289158afd6ca.jpg)

relative locations of the intersection point of the [1,3]-soliton and [3,4]-soliton with the x -axis. Especially, when s = 1, θ -[4 2] , = 0, θ -[1 3] , = 0, all of the four solitons will intersect at (0 , 0) when t = 0. One remarks that the bounded line soliton [1,4] (Mach stem wave), obtained by the balance between the exponential terms W φ ,φ ( 1 3 ) and cW φ ,φ ( 3 4 ), has the maximal amplitude among all the solitons by (1.11) (Fig. 1) and the velocity is obtained by (1.12). Furthermore, when t &lt; 0, there is a bounded line [2,3]-soliton (Fig. 2, the left side of the triangle), obtained by the balance between the exponential terms abW φ ,φ ( 2 4 ) and cW φ ,φ ( 3 4 ).

Figure 3. Initial wave.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0588v3-figure-9527c74a-5d71-4a0f-a2cd-c88a23d9b366.jpg)

Now, we consider the case α 3 = -α 2 ≥ 0, α 4 = -α 1 ≥ 0, and the amplitude

$$A = A _ { [ 1, 3 ] } = A _ { [ 4, 2 ] } \leq 2 \epsilon$$

is fixed. Then one can see that [1 3]-soliton and [4 2]-soliton is symmetric to the , , x -axis and similarly for [3 , 4]-soliton and [2 , 1]-soliton. By (1.11), one knows

$$\frac { \alpha _ { 3 } + \alpha _ { 1 } } { 2 } \leq \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 } = \frac { \alpha _ { 3 } + \alpha _ { 4 } } { 2 } = \arcsin \sqrt { \frac { A } { 2 \epsilon } }.$$

Therefore the angle between the [1,3]-soliton and the y -axis (counter-clockwise) is less than the the angle between the [3,4]-soliton and the y -axis (clockwise). We see that given A and 2 glyph[epsilon1] there is a critical angle ϕ C = arcsin √ A 2 glyph[epsilon1] for the angle between the [1,3]-soliton and the y -axis (counter-clockwise). Then one can introduce the following Miles-parameter [5, 8, 26] to describe the interaction for the Mach-type solution, noticing α 3 + α 1 2 ≤ 0,

$$\kappa = \frac { | \tan \frac { \alpha _ { 3 } + \alpha _ { 1 } } { 2 } | } { \tan \frac { \alpha _ { 3 } + \alpha _ { 4 } } { 2 } } = \frac { | \tan \frac { \alpha _ { 3 } + \alpha _ { 1 } } { 2 } | } { \tan \varphi _ { C } } = \frac { | \tan \frac { \alpha _ { 3 } + \alpha _ { 1 } } { 2 } | } { \sqrt { \frac { A } { 2 \epsilon - A } } } \leq 1.$$

From (1.11), we have thus using κ

$$& \text{from (1.11)}, \, \text{we have thus using } \kappa \\ & A = A _ { [ 1, 3 ] } = A _ { [ 4, 2 ] } = \frac { 2 \epsilon ( \tan \varphi _ { C } ) ^ { 2 } } { 1 + ( \tan \varphi _ { C } ) ^ { 2 } }, \quad A _ { [ 3, 4 ] } = A _ { [ 2, 1 ] } = \frac { 2 \epsilon ( \tan \varphi _ { C } ) ^ { 2 } } { \frac { 1 } { \kappa ^ { 2 } } + ( \tan \varphi _ { C } ) ^ { 2 } } \leq A, \\ & A _ { [ 1, 4 ] } = 2 \epsilon \sin ^ { 2 } \frac { \alpha _ { 4 } - \alpha _ { 1 } } { 2 } = 2 \epsilon \left [ \sin \left ( \frac { \alpha _ { 4 } - \alpha _ { 3 } } { 2 } + \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 } \right ) \right ] ^ { 2 } \\ & \quad = \frac { 2 \epsilon ( \tan \varphi _ { C } ) ^ { 2 } ( \kappa + 1 ) ^ { 2 } } { [ 1 + ( \tan \varphi _ { C } ) ^ { 2 } ] [ 1 + \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } ] } = A \frac { ( \kappa + 1 ) ^ { 2 } } { [ 1 + \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } ] } < 4 A. \\ & \text{Remark} \,. \, \text{To make a common/on with } \text{KP-(II)}. \, \text{we see that}.$$

Remark. To make a comparison with KP-(II), we see that

$$2 \epsilon - A = 2 \epsilon \left ( 1 - \sin ^ { 2 } \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 } \right ) = 2 \epsilon \cos ^ { 2 } \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 }.$$

Figure 4. The stem wave moves to the left.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0588v3-figure-43528e91-7eb0-401e-b822-1dd090088349.jpg)

When glyph[epsilon1] →∞ , α 1 →-π 2 and α 3 → π 2 such that

$$\epsilon \cos ^ { 2 } \frac { \alpha _ { 3 } - \alpha _ { 1 } } { 2 } = \frac { 1 } { 4 }.$$

Then

$$\kappa \rightarrow \frac { | \tan \frac { \alpha _ { 3 } + \alpha _ { 1 } } { 2 } | } { \sqrt { 2 A } },$$

which is the Miles parameter of KP-(II) to describe the interactions of water wave solitons [16, 17, 18].

Since the [1,4]-soliton (Mach stem wave) is increasing its length with time but its end points will lie in a line (see Figs. 4 and 5), we can obtain them as follows.

We choose s = 1, θ -[4 2] , = 0, θ -[1 3] , = 0 such that [1,3]-soliton and [1,4]-soliton will intersect at (0 , 0) when t = 0 (see Fig. 3). From (1.10)), the ridges of [1,3]-soliton and [1,4]-soliton are given by F p ( 1 ) -F p ( 3 ) = 0, F p ( 1 ) -F p ( 4 ) = 0, which lead to

$$x ( - \sin \alpha _ { 1 } + \sin \alpha _ { 3 } ) + y ( - \cos \alpha _ { 1 } + \cos \alpha _ { 3 } ) + t \epsilon ( \sin 3 \alpha _ { 1 } - \sin 3 \alpha _ { 3 } ) = 0,$$

$$x ( - \sin \alpha _ { 1 } + \sin \alpha _ { 4 } ) + y ( - \cos \alpha _ { 1 } + \cos \alpha _ { 4 } ) + t \epsilon ( \sin 3 \alpha _ { 1 } - \sin 3 \alpha _ { 4 } ) = 0.$$

$$\begin{smallmatrix} 2 & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & & \\ & \\ & & \\ & & \\ & \\ & & \\ & \\ & & \\ & \\ & & \\ & \\ & \\ & & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ & \\ \\ & \\ & \\ & \\ & \\ \\ & \\ & \\ & \\
\end{smallmatrix}$$

Noticing that α 3 = -α 2 ≥ 0, α 4 = -α 1 ≥ 0, one gets

$$x = \frac { t \epsilon \sin 3 \alpha _ { 4 } } { \sin \alpha _ { 2 } } = t \epsilon ( 4 \cos ^ { 2 } \alpha _ { 4 } - 1 ),$$

$$x & = \frac { t \epsilon \sin 3 \alpha _ { 4 } } { \sin \alpha _ { 4 } } = t \epsilon ( 4 \cos ^ { 2 } \alpha _ { 4 } - 1 ), & ( 2. 5 ) \\ y & = \frac { x ( \sin \alpha _ { 1 } - \sin \alpha _ { 3 } ) + t \epsilon ( - \sin 3 \alpha _ { 1 } + \sin 3 \alpha _ { 3 } ) } { - \cos \alpha _ { 1 } + \cos \alpha _ { 3 } } \\ & = t \epsilon \frac { ( 4 \cos ^ { 2 } \alpha _ { 4 } - 1 ) ( \sin \alpha _ { 1 } - \sin \alpha _ { 3 } ) + ( - \sin 3 \alpha _ { 1 } + \sin 3 \alpha _ { 3 } ) } { - \cos \alpha _ { 1 } + \cos \alpha _ { 3 } } \\ & = 4 t \epsilon \frac { \sin \alpha _ { 3 } ( \sin \alpha _ { 3 } + \sin \alpha _ { 4 } ) ( \sin \alpha _ { 4 } - \sin \alpha _ { 3 } ) } { \cos \alpha _ { 3 } - \cos \alpha _ { 4 } } = 4 t \epsilon \sin \alpha _ { 3 } ( \sin \alpha _ { 3 } + \sin \alpha _ { 4 } ) \cot \frac { \alpha _ { 3 } + \alpha _ { 4 } } { 2 }.$$

Figure 5. The length of the stem wave is increasing.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0588v3-figure-fe7a1e38-848f-48be-8005-2f3ddf228d58.jpg)

Using (2.3), one has

$$\text{Using (2.3), one has} \\ \sin \alpha _ { 3 } = \sin \left ( \frac { \alpha _ { 3 } + \alpha _ { 4 } } { 2 } - \frac { \alpha _ { 4 } - \alpha _ { 3 } } { 2 } \right ) = \frac { ( 1 - \kappa ) \tan \varphi _ { C } } { \sqrt { [ 1 + ( \tan \varphi _ { C } ) ^ { 2 } ] [ 1 + \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } ] } }, \\ \cos \alpha _ { 4 } = \cos \left ( \frac { \alpha _ { 3 } + \alpha _ { 4 } } { 2 } + \frac { \alpha _ { 4 } - \alpha _ { 3 } } { 2 } \right ) = \frac { 1 - \kappa ( \tan \varphi _ { C } ) ^ { 2 } } { \sqrt { [ 1 + ( \tan \varphi _ { C } ) ^ { 2 } ] [ 1 + \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } ] } }, \\ \sin \alpha _ { 3 } + \sin \alpha _ { 4 } = 2 \sin \frac { \alpha _ { 3 } + \alpha _ { 4 } } { 2 } \cos \frac { \alpha _ { 3 } - \alpha _ { 4 } } { 2 } = \frac { 2 \tan \varphi _ { C } } { \sqrt { [ 1 + ( \tan \varphi _ { C } ) ^ { 2 } ] [ 1 + \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } ] } }, \\ 4 \cos ^ { 2 } \alpha _ { 4 } - 1 = \frac { 4 [ 1 - \kappa ( \tan \varphi _ { C } ) ^ { 2 } ] ^ { 2 } } { [ 1 + ( \tan \varphi _ { C } ) ^ { 2 } ] [ 1 + \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } ] } - 1 \\ = \frac { 3 + ( \tan \varphi _ { C } ) ^ { 2 } [ 3 \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } - \kappa ^ { 2 } - 8 \kappa - 1 ] } { [ 1 + ( \tan \varphi _ { C } ) ^ { 2 } ] [ 1 + \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } ] }. \\ \text{A simple calculation yields using (2.4)}$$

A simple calculation yields using (2.4)

$$& y = 8 t \epsilon \tan \varphi _ { C } \frac { 1 - \kappa } { [ 1 + ( \tan \varphi _ { C } ) ^ { 2 } ] [ 1 + \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } ] } = 4 t A _ { [ 1, 4 ] } \frac { 1 - \kappa } { ( 1 + \kappa ^ { 2 } ) ( \tan \varphi _ { C } ) }, \\ & \tan \chi = \frac { y } { x } = \frac { 8 ( 1 - \kappa ) \tan \varphi _ { C } } { 3 + ( \tan \varphi _ { C } ) ^ { 2 } [ 3 \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } - \kappa ^ { 2 } - 8 \kappa - 1 ] }.$$

Hence one knows that the length of [1,4]-soliton is linear with time and its end points will lie in a line having slope ± tan χ (see Figs. 4 and 5). Furthermore, from (2.5), one gets that the [1,4]-soliton moves to the right if α 4 &lt; π 3 , and moves to the left if α 4 &gt; π 3 . In particular, if α 4 = π 3 or by (2.7)

$$3 + ( \tan \varphi _ { C } ) ^ { 2 } [ 3 \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } - \kappa ^ { 2 } - 8 \kappa - 1 ] = 0.$$

then [1,4]-soliton's length is increasing along the y -axis. When κ = 1 (or α 3 = 0), one has A = glyph[epsilon1] 2 by (2.8) and α 4 = π 3 . In this special case, the soliton is fixed. It is different from the KP-(II) case [18, 20].

## 3 Relations with V -shape initial value waves

In this section, we investigate some relations with the V -shape initial value wave for the NovikovVeselov equation (1.1), glyph[epsilon1] being fixed, as compared with the KP-(II) case [16, 17, 18, 20]. The main purpose is to study the interactions between line solitons, especially for the meaning of the critical angle ϕ C .

Recalling the one-soliton solution (1.10) and (1.11), one considers the initial data given in the shape of V with amplitude A and the oblique angle ϕ I &lt; 0 (measured in the clockwise sense from the y -axis):

$$A \, \text{sech} ^ { 2 } \left [ \sqrt { 2 A } \cos \varphi _ { I } ( x - | y | \tan \varphi _ { I } ) \right ].$$

For simplicity, one considers A ≤ 2 . glyph[epsilon1] We notice here the V -shape initial wave is in the negative x -region. The main idea is that we can think the initial value wave as a part of Mach-type soliton (2.1) or O -type soliton [8], that is, c = 0 in (2.1). In order to identify those soliton solutions from the V -shape (3.1), we denote them as [ i + , j + ]-soliton for y glyph[greatermuch] 0 and [ i -, j -]-soliton for y glyph[lessmuch] 0. Solitons for y →±∞ have by (1.10)

$$A & = 2 \epsilon \sin ^ { 2 } \frac { \alpha _ { j + } - \alpha _ { i + } } { 2 } = 2 \epsilon \sin ^ { 2 } \frac { \alpha _ { i - } - \alpha _ { j - } } { 2 }, \\ \varphi _ { I } & = \frac { \alpha _ { j + } + \alpha _ { i + } } { 2 } = - \left ( \frac { \alpha _ { i - } + \alpha _ { j - } } { 2 } \right ).$$

Assume that i + &lt; j + and i -&gt; j -. Then symmetry gives

$$\alpha _ { i ^ { + } } = - \alpha _ { i _ { = } }, \quad \alpha _ { j ^ { + } } = - \alpha _ { j _ { = } }.$$

Using the parameter (2.3) [5, 8, 26]

$$\kappa = \frac { | \tan \varphi _ { I } | } { \sqrt { \frac { A } { 2 \epsilon - A } } } = \frac { | \tan \varphi _ { I } | } { \tan \varphi _ { C } },$$

one can yield, noticing that ϕ C = α j + -α i + 2 = α i --α j -2 = arctan √ A 2 glyph[epsilon1] -A from (3.2),

- · κ ≥ 1 ⇒| ϕ I | ≥ ϕ C ⇒- ≤ π 2 α i + &lt; α j + &lt; α j -&lt; α i -≤ π 2 ( O -type),
- · 0 &lt; κ &lt; 1 ⇒| ϕ I | &lt; ϕ C ⇒- ≤ π 2 α i + &lt; α j -&lt; α j + &lt; α i -≤ π 2 (Mach-type).

We remark here that if κ = 1 (or α 3 = 0), then it is of O -type by (2.4) and (2.6). One can see that if the angle ϕ I is small, then an intermediate wave called the Mach stem ([1,4]-soliton) appears. The Mach stem, the incident wave ([1,3]-soliton) and the reflected wave ([3,4]-soliton) interact resonantly, and those three waves form a resonant triple. It is similar to the KP-(II) case [16].

Let's compute the maximal amplitude of the Mach stem ([1,4]-soliton) for fixed amplitude A and glyph[epsilon1] . By (2.4), a simple calculation shows that

$$\frac { d A _ { [ 1, 4 ] } } { d \kappa } = A \frac { 2 ( \kappa + 1 ) [ 1 - \kappa ( \tan \varphi _ { C } ) ^ { 2 } ] } { [ 1 + \kappa ^ { 2 } ( \tan \varphi _ { C } ) ^ { 2 } ] ^ { 2 } }.$$

Hence one can see that when κ = 1 (tan / ϕ C ) 2 , that is,

$$( \tan \varphi _ { C } ) ( | \tan \varphi _ { I } | ) = 1,$$

the Mach stem has the maximal amplitude. Consequently, if

$$\varphi _ { C } + \varphi _ { I } = \frac { \pi } { 2 },$$

then one obtains by (3.4), recalling that 0 &lt; κ &lt; 1 (or tan ϕ C &gt; 1, i.e., A &gt; glyph[epsilon1] ),

$$A _ { [ 1, 4 ] } ^ { \max } = A \left ( 1 + \frac { 1 } { ( \tan \varphi _ { C } ) ^ { 2 } } \right ) = 2 \epsilon < 2 A.$$

Therefore one sees that from (3.4), A and glyph[epsilon1] being fixed,

- · 0 &lt; κ &lt; 1 (tan ϕ C ) 2 , the amplitude A [1 4] , (stem wave) is increasing;
- · κ = 1 (tan ϕ C ) 2 (or (3.5)), the amplitude A [1 4] , has the maximal value 2 ; glyph[epsilon1]
- · 1 (tan ϕ C ) 2 &lt; κ &lt; 1, the amplitude A [1 4] , is decreasing.

It is noteworthy that the maximal amplitude is independent of A . Also, we know that the maximal amplitude of Mach stem for NV equation is less than twice of the incident wave's one; however, for the KP equation (shallow water waves), the Mach stem's amplitude can be four times of the incident wave's one [18]. This is the different point from the case of the KP equation.

On the other hand, one can see that for κ &gt; 1 ( O -type) we have 0 ≤ α j + -α i + 2 ≤ π 4 by (3.3), that is, A ≤ glyph[epsilon1] . Thus, if we choose A such that

$$\epsilon < A \leq 2 \epsilon,$$

we get π 2 &lt; α i -≤ π ; therefore, under the condition (3.7), the initial value wave (3.1) would develop into a singular O -type soliton by (2.1) ( c = 0) when glyph[epsilon1] is fixed. On the other hand, when | ϕ I | ≤ π 2 , A and κ are fixed, one can choose

$$\epsilon = \frac { A } { 2 } \left [ 1 + \left ( \frac { \kappa } { \tan \varphi _ { I } } \right ) ^ { 2 } \right ] \geq \frac { A } { 2 }.$$

Then we can obtain regular soliton solutions.

Finally, from (2.5) one remarks that the [1,4]-soliton (stem wave) moves to the right if α 4 &lt; π 3 , and moves to the left if α 4 &gt; π 3 . The former case is different from the KP equation (shallow water waves); i.e., the stem wave moves with the same side of incident wave for the KP equation. On the other hand, if we replace the condition (2.2) by A = A [3 4] , = A [2 1] , ≤ 2 , then by (1.11) the glyph[epsilon1] [3,4]-soliton (the incident wave) has smaller amplitude than the [1,3]-soliton's one (the reflected wave). But this is not physically interesting.

## 4 Concluding remarks

One investigates the Mach-type (or (3142)-type) soliton of the Novikov-Veselov equation. The Mach stem ([1,4]-soliton), the incident wave ([1,3]-soliton) and the reflected wave ([3,4]-soliton) form a resonant triple. From (3.6), we see that the amplitude of Mach stem is less than two times of the one of the incident wave, which is different from the KP equation [18]; moreover, the length of the Mach stem is computed and show it is linear with time (2.6). On the other hand, one uses the parameter κ (2.3) to describe the critical behavior for the O -type and Mach-type solitons and notices that it depends on the the fixed parameter glyph[epsilon1] . We see that the amplitude A of the incident wave is small than 2 ; glyph[epsilon1] furthermore, if glyph[epsilon1] &lt; A &lt; 2 , glyph[epsilon1] then the soliton will be singular. Now, a natural question is: what happens if A &gt; 2 glyph[epsilon1] when glyph[epsilon1] is fixed in (1.1)? Another

question is the minimal completion [20]. It means the resulting chord diagram has the smallest total length of the chords. This minimal completion can help us study the asymptotic solutions and estimate the maximum amplitude generated by the interaction of those initial waves. A numerical investigation of these issues will be published elsewhere.

## Acknowledgements

The author thanks the referees for their valuable suggestions. This work is supported in part by the National Science Council of Taiwan under Grant No. NSC 102-2115-M-606-001.

## References

- [1] Ablowitz M.J., Baldwin D.E., Nonlinear shallow ocean-wave soliton interactions on flat beaches, Phys. Rev. E 86 (2012), 036305, 5 pages, arXiv:1208.2904.
- [2] Athorne C., Nimmo J.J.C., On the Moutard transformation for integrable partial differential equations, Inverse Problems 7 (1991), 809-826.
- [3] Biondini G., Chakravarty S., Soliton solutions of the Kadomtsev-Petviashvili II equation, J. Math. Phys. 47 (2006), 033514, 26 pages, nlin.SI/0511068.
- [4] Bogdanov L.V., Veselov-Novikov equation as a natural two-dimensional generalization of the Korteweg-de Vries equation, Theoret. and Math. Phys. 70 (1987), 219-223.
- [5] Chakravarty S., Kodama Y., Soliton solutions of the KP equation and application to shallow water waves, Stud. Appl. Math. 123 (2009), 83-151, arXiv:0902.4423.
- [6] Chakravarty S., Lewkow T., Maruno K.I., On the construction of the KP line-solitons and their interactions, Appl. Anal. 89 (2010), 529-545, arXiv:0911.2290.
- [7] Chang J.H., On the N -solitons solutions in the Novikov-Veselov equation, SIGMA 9 (2013), 006, 13 pages, arXiv:1206.3751.
- [8] Chang J.H., The interactions of solitons in the Novikov-Veselov equation, arXiv:1310.4027.
- [9] Dubrovsky V.G., Topovsky A.V., Basalaev M.Y., New exact multi line soliton and periodic solutions with constant asymptotic values at infinity of the NVN integrable nonlinear evolution equation via dibar-dressing method, arXiv:0912.2155.
- [10] Grinevich P.G., The scattering transform for the two-dimensional Schr¨dinger operator with a potential o that decreases at infinity at fixed nonzero energy, Russ. Math. Surv. 55 (2000), 1015-1083.
- [11] Grinevich P.G., Manakov S.V., Inverse problem of scattering theory for the two-dimensional Schr¨dinger o operator, the ¯ -method and nonlinear equations, ∂ Funct. Anal. Appl. 20 (1986), 94-103.
- [12] Hu H.-C., Lou S.-Y., Construction of the Darboux transformaiton and solutions to the modified NizhnikNovikov-Veselov equation, Chinese Phys. Lett. 21 (2004), 2073-2076.
- [13] Hu H.-C., Lou S.-Y., Liu Q.-P., Darboux transformation and variable separation approach: the NizhnikNovikov-Veselov equation, Chinese Phys. Lett. 20 (2003), 1413-1415, nlin.SI/0210012.
- [14] Ishikawa M., Wakayama M., Applications of minor-summation formula. II. Pfaffians and Schur polynomials, J. Combin. Theory Ser. A 88 (1999), 136-157.
- [15] Kazeykina A.V., Novikov R.G., Large time asymptotics for the Grinevich-Zakharov potentials, Bull. Sci. Math. 135 (2011), 374-382, arXiv:1011.4038.
- [16] Kodama Y., KP solitons in shallow water, J. Phys. A: Math. Theor. 43 (2010), 434004, 54 pages, arXiv:1004.4607.
- [17] Kodama Y., KP solitons and Mach reflection in shallow water, arXiv:1210.0281.
- [18] Kodama Y., Lectures delivered at the NSF/CBMS Regional Conference in the Mathematical Sciences 'Solitons in Two-Dimensional Water Waves and Applications to Tsunami' (UTPA, May 20-24, 2013), available at http://faculty.utpa.edu/kmaruno/nsfcbms-tsunami.html .
- [19] Kodama Y., Maruno K.-I., N -soliton solutions to the DKP equation and Weyl group actions, J. Phys. A: Math. Gen. 39 (2006), 4063-4086, nlin.SI/0602031.
- [20] Kodama Y., Oikawa M., Tsuji H., Soliton solutions of the KP equation with V -shape initial waves, J. Phys. A: Math. Theor. 42 (2009), 312001, 9 pages, arXiv:0904.2620.

- [21] Kodama Y., Williams L., KP solitons, total positivity, and cluster algebras, Proc. Natl. Acad. Sci. USA 108 (2011), 8984-8989, arXiv:1105.4170.
- [22] Kodama Y., Williams L., The Deodhar decomposition of the Grassmannian and the regularity of KP solitons, Adv. Math. 244 (2013), 979-1032, arXiv:1204.6446.
- [23] Kodama Y., Williams L., KP solitons and total positivity for the Grassmannian, Invent. Math. 198 (2014), 637-699, arXiv:1106.0023.
- [24] Manakov S.V., The method of the inverse scattering problem, and two-dimensional evolution equations, Russian Math. Surveys 31 (1976), no. 5, 245-246.
- [25] Matveev V.B., Salle M.A., Darboux transformations and solitons, Springer Series in Nonlinear Dynamics , Springer-Verlag, Berlin, 1991.
- [26] Miles J.W., Resonantly interacting solitary waves, J. Fluid Mech. 79 (1977), 171-179.
- [27] Nimmo J.J.C., Darboux transformations in (2 + 1)-dimensions, in Applications of Analytic and Geometric Methods to Nonlinear Differential Equations (Exeter, 1992), NATO Adv. Sci. Inst. Ser. C Math. Phys. Sci. , Vol. 413, Kluwer Acad. Publ., Dordrecht, 1993, 183-192.
- [28] Novikov S.P., Veselov A.P., Two-dimensional Schr¨dinger operator: inverse scattering transform and evoluo tional equations, Phys. D 18 (1986), 267-273.
- [29] Ohta Y., Pfaffian solutions for the Veselov-Novikov equation, J. Phys. Soc. Japan 61 (1992), 3928-3933.
- [30] Veselov A.P., Novikov S.P., Finite-gap two-dimensional potential Schr¨dinger operators. Explicit formulas o and evolution equations, Sov. Math. Dokl. 30 (1984), 588-591.
- [31] Yeh H., Li W., Kodama Y., Mach reflection and KP solitons in shallow water, Eur. Phys. J. ST 185 (2010), 97-111, arXiv:1004.0370.
Jen-Hsu Chang
Mach-Type Soliton in the Novikov-Veselov Equation
Using the reality condition of the solutions, one constructs the Mach-type
soliton of the Novikov-Veselov equation by the minor-summation formula of the
Pfaffian. We study the evolution of the Mach-type soliton and find that the
amplitude of the Mach stem wave is less than two times of the one of the
incident wave. It is shown that the length of the Mach stem wave is linear with
time. One discusses the relations with V -shape initial value wave for
different critical values of Miles parameter.