Abstract

In the present paper, we have obtained an exact biperiodic, one-phase solution of the Kawahara evolution equation. Two classes of real periodic waves generated by the biperiodic solution have been analyzed. A modification of the bilinear-transformation method has been applied allowing to provide a single solution of the residual equation derived from the bidifferential reduction of the considered nonintegrable equation. It is shown that the spatial displacements are individual for each separate harmonic of the real periodic solutions.

1. Introduction

The nonlinear evolution equation 𝑒𝑑+𝛼𝑒𝑒π‘₯+𝛽𝑒π‘₯π‘₯π‘₯=𝛾𝑒π‘₯π‘₯π‘₯π‘₯π‘₯(1.1) was introduced by Kawahara [1] as a model describing one-dimensional waves for which the angle between the front and the gradient of the external field tends to the critical angle: πœ‘π‘βˆš=arctg(π‘š1/π‘š0βˆ’βˆšπ‘š0/π‘š1). Here, 𝛼, 𝛽, and𝛾 are real parameters (𝛼≠0), 𝑒=𝑒(π‘₯,𝑑) is the elevation of dispersion medium, and π‘š0, π‘š1are the ion and electron masses, respectively. Within the zone of the critical angle πœ‘π‘, the coefficient in front of 𝑒π‘₯π‘₯π‘₯ in (1.1) decreases causing disbalance between the nonlinear and dispersion effects, which is balanced by the member with the highest-order derivative. Therefore, three varieties of the dispersion medium are distinguished: medium with positive dispersion (𝛽<0, 𝛾<0); medium with negative dispersion (𝛽>0, 𝛾>0) and medium with mixed dispersion (𝛽𝛾<0). Yamamoto [2] obtained the same evolution equation (1.1) in describing the dynamics of a nonviscous fluid in the vicinity of the critical depth. In this case, dispersion of the lowest order caused by the gravitation is balanced by the one caused by the surface tension.

The Kawahara evolution equation (1.1) is nonintegrable, which can be easily established by applying the modified criterion of Ablowtz and Segur [3]. Although it does not possess any conservation law and 𝑁-soliton solutions (𝑁β‰₯2), (1.1) has been a subject of research by a number of authors employing both analytic and numerical methods.

Kano and Nakayama [4] have found elliptic solution of (1.1), from which by means of appropriate phase modulations they obtained periodic cnoidal solutions, as well as M-type and W-type solitary waves. Yamamoto and Takizawa [5] have found a stable progressive-pulse solution 𝑒(π‘₯,𝑑)=105𝛽2169𝛾𝛼secβ„Ž412ξƒŽπ›½ξ‚΅13𝛾π‘₯βˆ’36𝛽2𝑑,169𝛾(1.2) which actually represents one-soliton impulse that can be generated only in a medium with positive or negative dispersion (i.e., for sign𝛽=sign𝛾). In the particular case 𝛽=0, (1.1) is called FKdV (Five Korteweg-de Vries) and has been thoroughly studied numerically and analytically by Boyd [6]. Summarizing the numerical analysis accomplished by a number of authors [1, 2, 7], we can say that the Kawahara equation has two types of localized solitary waves: compressed and fictitiously scattered corresponding to a medium with negative and positive dispersion, respectively.

Studies on the existence of periodic solutions of (1.1) are considerably fewer. In 1997 Berloff and Howard [8], by applying the singular manifold method adapted for partial differential equations [9], obtained solitary and approximate periodic solutions of the Kawahara equation.

In the present paper, two families of real periodic solutions of the aforementioned equation are obtained. These exact periodic solutions were proved to be dynamically equivalent. A β€œspatial” version [10] of the bilinear-transformation method of Hirota [11] and Matsuno [12] has been applied to derive them. In its classic form, this analytic model is inapplicable to nonintegrable partial differential equations.

2. Preliminaries

We could reduce the Kawahara equation in the so-called normal form if instead of the parameter 𝛽 we put πœ€π›½, that is, π›½β†’πœ€π›½, where πœ€=Β±1, 𝛽>0, and assume that 𝛾>0.

As will be seen below, the last condition does not lead to loss of generality since after the following rescaling: ξ‚™π‘₯→𝛾𝛽𝛾π‘₯,𝑑→𝛽2𝛾𝛽𝛽𝑑,𝑒→2𝛼𝛾𝑒(2.1) (1.1) is reduced to its normal form 𝑒𝑑+𝑒𝑒π‘₯+πœ€π‘’π‘₯π‘₯π‘₯=𝑒π‘₯π‘₯π‘₯π‘₯π‘₯,πœ€=Β±1.(2.2)

Hence, it is easy to realize that for πœ€=+1, we have a medium with negative dispersion, for πœ€=+1 along with replacing the variables π‘’β†’βˆ’π‘’, π‘₯β†’βˆ’π‘₯, π‘‘β†’βˆ’π‘‘, the medium possesses positive dispersion, while for πœ€=βˆ’1, we have a medium with mixed dispersion, that is, we have all available variations of the dispersion medium.

The Kawahara equation (2.2) is invariant with respect to the Galilean transformation π‘₯ξ…ž=π‘₯+πœ†π‘’0, π‘‘ξ…ž=𝑑, and π‘ˆ=π‘’βˆ’π‘’0, where πœ† = const, 𝑒0 = const, that is, if π‘ˆ(π‘₯βˆ’π‘£0𝑑) is a solution of (2.2) with constant phase velocity 𝑣0, then the function 𝑒(π‘₯,𝑑)=𝑒0+π‘ˆ[π‘₯βˆ’(𝑣0+πœ†π‘’0)𝑑] is also a solution of the same equation but with increased phase velocity 𝑣0→𝑣0+πœ†π‘’0.

3. Biperiodic Solution

We will represent the solution of (2.2) by means of the Hirota-Satsuma [13] transformation 𝑒(π‘₯,𝑑)=π‘Ž+2πœ‡(ln𝑓(π‘₯,𝑑))π‘₯π‘₯,(3.1) where π‘Ž, πœ‡ are unknown parameters so far, and 𝑓(π‘₯,𝑑) is unknown function, but on the assumption periodic and continuously differentiable to seventh order (with respect to π‘₯) in the semi-infinite domain Ξ©={(π‘₯,𝑑)βˆˆπ‘…2,βˆ’βˆž<π‘₯<∞,𝑑>0}. By substituting 𝑒(π‘₯,𝑑) from (3.1) into (2.2) and employing the bidifferential identities for the Hirota operator (see Appendix A) π·π‘›π‘‘π·π‘šπ‘₯ξ€·π‘₯πœ‘(π‘₯,𝑑)β‹…πœ“ξ…ž,π‘‘ξ…žξ€Έ=ξ‚€πœ•βˆ’πœ•πœ•π‘‘πœ•π‘‘ξ…žξ‚π‘›ξ‚€πœ•βˆ’πœ•πœ•π‘₯πœ•π‘₯ξ…žξ‚π‘šπœ‘ξ€·π‘₯(π‘₯,𝑑)πœ“ξ…ž,π‘‘ξ…žξ€Έ||||𝑑=𝑑′π‘₯=π‘₯β€²,(3.2) we will obtain the following bidifferential form of this equation: 𝑓4𝐷𝑑𝐷π‘₯+πœ€π·4π‘₯βˆ’π·6π‘₯ξ€»π‘“βˆ’8𝐢𝑓⋅𝑓+22𝐷(πœ‡βˆ’6πœ€)2π‘₯𝑓⋅𝑓2+𝐷2π‘₯𝑓⋅𝑓15𝑓2𝐷4π‘₯𝑓⋅𝑓+𝛼𝑓4ξ€·π·βˆ’302π‘₯𝑓⋅𝑓2ξ‚„=0,(3.3) where 𝐢 is an integration constant assumed to be nonzero. This constant plays a major role in constructing the periodic solutions though it does not have any dynamic features. If in (3.3) we put πœ‡=6πœ€, then a sufficient condition for the function 𝑓(π‘₯,𝑑) to be its solution is to satisfy the following two equations accordingly: 𝐷𝑑𝐷π‘₯+πœ€π·4π‘₯βˆ’π·6π‘₯ξ€Έβˆ’4𝐢𝑓⋅𝑓=0,(3.4)π‘Žπ‘“4+15𝑓2𝐷4π‘₯𝐷𝑓⋅𝑓=302π‘₯𝑓⋅𝑓2.(3.5)

The first of these equations is called bidifferential (since its structure is such), and the second one will be called residual. We will search for the solution of (3.4) in the form 𝑓(π‘₯,𝑑)=πœƒ3(πœ‰,π‘ž)=βˆžξ“π‘›=βˆ’βˆžπ‘žπ‘›2𝑒2π‘–π‘›πœ‰,(3.6) where πœƒ3 is the Jacobi biperiodic function [14] with phase variable πœ‰=π‘˜π‘₯+πœ”π‘‘+𝛿. It is possible that the parameters π‘˜, πœ”, 𝛿 are complex and π‘ž=π‘’π‘–πœ‹πœ, and Im𝜏>0 (i.e., 0<|π‘ž|<1) is the perturbation parameter. The function πœƒ3(πœ‰,π‘ž) has a real period πœ‹/π‘˜ and an imaginary one 𝜏/π‘˜. We substitute 𝑓(π‘₯,𝑑) from (3.6) into the bidifferential (3.4). We obtain an infinite system of algebraic equations βˆžξ“π‘š=βˆ’βˆžπΉ(π‘š)𝑒2π‘–π‘šπœ‰=0,(3.7) where πΉπ‘š(π‘š)=βˆžξ“π‘›=βˆ’βˆžξ€Ίβˆ’4π‘˜πœ”(2π‘›βˆ’π‘š)2+16πœ€π‘˜4(2π‘›βˆ’π‘š)4βˆ’64π‘˜6(2π‘›βˆ’π‘š)6ξ€»π‘ž+8𝐢𝑛22+(π‘›βˆ’π‘š)(3.8) and having π‘š=0,Β±1,Β±2,….

The bilinear structure of (3.7) makes it possible to apply the index parity principle therein, which means that if in (3.8) we substitute 𝑛→𝑛+1, we will have the following relations: 𝐹(π‘š)=𝐹(π‘šβˆ’2)π‘ž2(π‘šβˆ’1)=𝐹(π‘šβˆ’4)π‘ž2(2π‘šβˆ’4)=β‹―=𝐹(0)π‘žπ‘š2/2ifπ‘šisanevennumber𝐹(π‘š)=𝐹(π‘šβˆ’2)π‘ž2(π‘šβˆ’1)=𝐹(π‘šβˆ’4)π‘ž2(2π‘šβˆ’4)=β‹―=𝐹(1)π‘ž(π‘š2βˆ’1)/2ifπ‘šisanoddnumber.(3.9) Summing up separately the even and odd addends in (3.7), we will obtain the following more compact form of this equation: 𝐹(0)πœƒ3ξ€·2πœ‰,π‘ž2ξ€Έ+π‘žβˆ’1/2𝐹(1)πœƒ2ξ€·2πœ‰,π‘ž2ξ€Έ=0,(3.10) where πœƒ2(𝑧,π‘ž) is the second Jacobi biperiodic function [14] defined by the equality πœƒ2βˆ‘(𝑧,π‘ž)=βˆžπ‘›=βˆ’βˆžπ‘ž(𝑛+1/2)2𝑒𝑖(2π‘›βˆ’1)𝑧. Accounting for the linear independence of the functions πœƒ2(2πœ‰,π‘ž2) and πœƒ3(2πœ‰,π‘ž2) in (3.10), it is reduced to two equations: 𝐹(0)=0 and 𝐹(1)=0. In the context of equality (3.8) and the functional identities for the πœƒ-functions, given in Appendix B, the last two equations are reduced to the following algebraic linear system: ξ€·π‘˜πœƒξ…ž3ξ€Έξ‚΅πœƒπœ”+3π‘žξ‚ΆπΆ=8πœ€π‘˜4ξ€·πœƒξ…ž3+π‘žπœƒ3ξ…žξ…žξ€Έβˆ’64π‘˜6ξ€·πœƒξ…ž3+3π‘žπœƒ3ξ…žξ…ž+π‘ž2πœƒ3ξ…žξ…žξ…žξ€Έ,ξ€·π‘˜πœƒξ…ž2ξ€Έξ‚΅πœƒπœ”+2π‘žξ‚ΆπΆ=8πœ€π‘˜4ξ€·πœƒξ…ž2+π‘žπœƒ2ξ…žξ…žξ€Έβˆ’64π‘˜6ξ€·πœƒξ…ž2+3π‘žπœƒ2ξ…žξ…ž+π‘ž2πœƒ2ξ…žξ…žξ…žξ€Έ,(3.11) with respect to πœ” and 𝐢. This system is compatible and definite since Ξ”=(π‘˜/π‘ž)(πœƒ2πœƒξ…ž3βˆ’πœƒ3πœƒξ…ž2)=(π‘˜/π‘ž)π‘Š(πœƒ2,πœƒ3)β‰ 0, where π‘Š is the Wronskian determinant, and the wave number is different from zero. We have denoted for convenience πœƒπ‘—=πœƒπ‘—(0,π‘ž2), 𝑗=2,3, so that all derivatives of these functions in the system (3.11) are with respect to the parameter π‘ž. The solution of system (3.11) is as follows: πœ”=8πœ€π‘˜3ξƒ¬π‘Š1+π‘žξ…žξ€·πœƒ2,πœƒ3ξ€Έπ‘Šξ€·πœƒ2,πœƒ3ξ€Έξƒ­βˆ’64π‘˜5π‘Šξ€·πœƒ2,πœƒ3ξ€Έξ€½π‘Šξ€·πœƒ2,πœƒ3ξ€Έ+3π‘žπ‘Šξ…žξ€·πœƒ2,πœƒ3ξ€Έ+π‘ž2ξ€Ίπ‘Šξ…žξ…žξ€·πœƒ2,πœƒ3ξ€Έξ€·πœƒβˆ’π‘Šξ…ž2,πœƒξ…ž3,ξ€Έξ€»ξ€Ύ(3.12)𝐢=βˆ’8πœ€π‘˜4π‘ž2π‘Šξ€·πœƒξ…ž2,πœƒξ…ž3ξ€Έπ‘Šξ€·πœƒ2,πœƒ3ξ€Έ+64π‘˜6π‘ž2π‘Šξ€·πœƒ2,πœƒ3ξ€Έξ€Ίπ‘Šξ€·πœƒξ…ž2,πœƒξ…ž3ξ€Έ+3π‘žπ‘Šξ…žξ€·πœƒξ…ž2,πœƒξ…ž3ξ€Έξ€».(3.13)

The residual equation (3.5) does not possess a bilinear structure, which is a typical feature for the nonintegrable partial differential equations. This means that we cannot apply the index parity principle to (3.5), that is, we cannot reduce the infinite system generated by this equation (after substituting 𝑓(π‘₯,𝑑) with its equal function πœƒ3(πœ‰,𝑑)) to a system of two equations. Consequently, we will represent the parameter π‘Ž in the formal series π‘Ž=π‘˜4βˆžξ“π‘š=βˆ’βˆžπ‘Žπ‘š,(3.14) by the terms π‘Žπ‘š, which are real parameters unknown at this stage. Substituting (3.6) and (3.14) into the residual equation (3.5), we will obtain the infinite system π‘Žπ‘šβˆžξ“π‘›=βˆ’βˆžπ‘ž2(π‘›βˆ’π‘š)2+(2π‘›βˆ’3π‘š)2=240βˆžξ“π‘›=βˆ’βˆžξ€Ί2𝑛2(3π‘›βˆ’π‘š)2βˆ’(2π‘›βˆ’π‘š)4ξ€»π‘ž2𝑛2+(2π‘›βˆ’π‘š)2,π‘š=0,Β±1,Β±2,….(3.15)

The infinite series on both sides of the system (3.15) are absolutely convergent for 0<|π‘ž|<1. It is easy to deduce that in summation term by term of the equalities from (3.15) for each whole number π‘š, the series (3.14) is an absolutely and uniformly convergent series under the same condition for the perturbation parameter. Thus, we obtain that the terms of the functional series (3.14) are determined for every integer π‘š by the formula π‘Žπ‘šβˆ‘(π‘ž)=240βˆžπ‘›=βˆ’βˆžξ€Ί2𝑛2(3π‘›βˆ’π‘š)2βˆ’(2π‘›βˆ’π‘š)4ξ€»π‘ž2𝑛2+(2π‘›βˆ’π‘š)2βˆ‘βˆžπ‘›=βˆ’βˆžπ‘ž2(π‘›βˆ’π‘š)2+(2π‘›βˆ’3π‘š)2,π‘š=0,Β±1,Β±2,….(3.16)

The terms π‘Žπ‘š(π‘ž) of the uniformly convergent series (3.14) will be called spatial displacements. This will be explained further.

And finally in this section, we can make the conclusion that the continuously differentiable function 𝑒(π‘₯,𝑑)=π‘˜4βˆžξ“π‘š=βˆ’βˆžπ‘Žπ‘š(π‘ž)+12πœ€π‘˜2π‘‘ξƒ¬πœƒπ‘‘πœ‰ξ…ž3(πœ‰,π‘ž)πœƒ3ξƒ­(πœ‰,π‘ž)(3.17) is a family of localized biperiodic solutions of the Kawahara equation (2.2). For convenience, in (3.17), we have denoted πœƒξ…ž3𝑑(πœ‰,π‘ž)=πœƒπ‘‘πœ‰3(πœ‰,π‘ž).(3.18)

Figure 1 shows the discrete values of the spatial displacements π‘Žπ‘š(π‘ž) for a mean value of the perturbation parameter π‘ž=π‘’βˆ’πœ€πœ‹, πœ€=1, π‘˜=1, and π‘š=0,1,2,…. In addition to the foregoing, the function 𝑒(π‘₯,𝑑), expressed by (3.17), is a solution of the considered equation if the phase velocity πœ”(π‘˜,π‘ž) satisfies the dispersion relation (3.12), and the integration constant 𝐢 is as in (3.13). The wave number π‘˜β‰ 0 is complex in the general case, but we have to exclude those values for which the constant 𝐢 is zero, that is, π‘˜2β‰ πœ€8ξƒ¬π‘Š1+3π‘žξ…žξ€·πœƒξ…ž2ξ€·0,π‘ž2ξ€Έ,πœƒξ…ž3ξ€·0,π‘ž2ξ€Έξ€Έπ‘Šξ€·πœƒξ…ž2ξ€·0,π‘ž2ξ€Έ,πœƒξ…ž3ξ€·0,π‘ž2ξƒ­ξ€Έξ€Έβˆ’1.(3.19)

4. Analicity Condition and Real Biperiodic Solutions

The obtained exact biperiodic solution (3.17) of the Kawahara evolution equation (2.2) is generally speaking a rational complex function of the phase variable 𝝃=π‘˜π‘₯+πœ”π‘‘+𝛿. This function has twofold poles in the lattice of complex points πœ‰π‘šπ‘›=πœ‹(π‘š+1/2)+π‘–πœπœ‹(𝑛+1/2),π‘š,π‘›βˆˆπ‘, which are zeros of the function πœƒ3(πœ‰,π‘ž). As we are interested in the physical relevance of the obtained biperiodic solution, we can suitably choose the free parameters π‘˜, π‘ž, 𝛿, so that we can achieve real biperiodic analogues of the solution and at the same time avoid the twofold poles. For this purpose, let us assume that 𝜏=𝑖𝑠,𝑠>0,thatis,π‘ž=π‘’βˆ’πœ‹π‘ <1.(4.1)

Under these conditions, to avoid the poles of the solution 𝑒(π‘₯,𝑑) at the points πœ‰π‘šπ‘›, it is sufficient to limit the phase variable within the horizontal strip βˆ’πœ‹π‘ <Im(πœ‰)<πœ‹π‘ ,(4.2) which is actually the analicity condition for the solution. To provide real biperiodic solutions from (3.17), we will consider two possible options for the wave number π‘˜, under the condition of (3.19).(i)The wave number π‘˜ is real. Under this assumption and provided the condition of (4.1) holds, the phase variable πœ‰=π‘˜π‘₯+πœ”π‘‘+𝛿 is real for real phase shift 𝛿. The logarithmic derivative from (3.17) in this case can be expressed in a Fourier series [14] by the equality πœƒξ…ž3(πœ‰,π‘ž)πœƒ3(πœ‰,π‘ž)=4βˆžξ“π‘š=βˆ’βˆž(βˆ’1)π‘šπ‘žπ‘š1βˆ’π‘ž2π‘šsin(2π‘šπœ‰).(4.3) Thus, from (3.17), we obtain a well-defined, real function which is periodical on the real straight line 𝑒(π‘₯,𝑑)=π‘˜2βˆžξ“π‘š=βˆ’βˆžξ€Ίπ‘˜2π‘Žπ‘š(π‘ž)+48πœ€π‘š(βˆ’1)π‘šξ€».coscβ„Ž(π‘šπœ‹π‘ )β‹…cos(2π‘šπœ‰)(4.4) Figure 2 illustrates the forms of the real periodic solution for a real value of the wave number π‘˜, (π‘˜=1). It can be seen from the last formula that the spatial displacements π‘Žπ‘š(π‘ž) contribute to each separate harmonic of the biperiodic solution.(ii)The wave number π‘˜ is imaginary. Let π‘˜β†’π‘–π‘˜, and without loss of generality, we can assume that π‘˜>0. Under this hypothesis, it is possible to generate real periodic solutions if the phase velocity is an imaginary number and the phase shift is properly chosen. It is obvious from the dispersion equality (3.12) that the phase velocity πœ” is an imaginary number (since πœ”(π‘˜) contains only odd powers of π‘˜). For the phase shift 𝛿, we choose 𝛿→𝛿+πœ‹πœ, and hence π‘–πœ‰β†’π‘–πœ‰+πœ‹πœ. With this transition, the dispersion relation (3.12) remains the same. Applying the quasiperiodic property for πœƒ3, πœƒ2(π‘–πœ‰,π‘ž)=π‘ž1/4π‘’π‘–πœ‰πœƒ3ξ‚€π‘–πœ‰+πœ‹πœ2,,π‘ž(4.5) and from formula (3.17), we have for the logarithmic derivative (see [2]) 𝑑2π‘‘πœ‰2lnπœƒ2(π‘–πœ‰,π‘ž)=βˆ’βˆžξ“π‘š=βˆ’βˆžsecβ„Ž2[].𝑖(π‘–πœ‰βˆ’π‘–π‘šπœ‹π‘ )(4.6) In this case, we also obtain a real periodic function which is well defined in the strip (4.2): 𝑒(π‘₯,𝑑)=βˆžξ“π‘š=βˆ’βˆžξ€Ίπ‘˜4π‘Žπ‘š(π‘ž)+12πœ€π‘˜2secβ„Ž2(ξ€».πœ‰βˆ’π‘šπœ‹π‘ )(4.7)

Figure 3 visualizes the periodic solitary-wave profiles of the secβ„Ž2 type for imaginary values of the wave number π‘˜βˆΆπ‘–,𝑖/2,𝑖/4, and π‘ž=π‘’βˆ’πœ‹. Despite their structural differences, the two types of real biperiodic solutions of the Kawahara equation obtained in (4.4) and (4.7) are dynamically equivalent. This comes to confirm the nonlinear superposition principle established by Toda [15] for the KdV evolution equation. In conclusion, let us note that the biperiodic real solutions (4.4) and (4.7) actually describe the dynamics of dispersing waves resulting from the dispersion relation (3.12). Indeed, πœ”ξ…žξ…ž(π‘˜)β‰ 0, which means that these waves are propagating in two directions, as periodically repeating groups traveling with group velocity 𝑉=πœ”ξ…ž(π‘˜).

5. Conclusions

Regardless of the purely technical difficulties arising in the application of the bilinear transformation method, it turns out to be a versatile theoretical instrument in analyzing nonintegrable, nonlinear partial differential equations, such as the Kawahara equation (1.1). The major intricacy arises from the residual equation after the bilinear reduction of the initial equation. Harmonization of the bilinear and the residual equation practically means to satisfy an algebraic system with infinite number of equations, but with a finite number of unknowns (usually one or two). The solution of a similar infinite system is feasible if a suitable unknown quantity is presented in the form of a convergent numerical or functional series with unknown termsβ€”a procedure often applied in the Fourier method for linear partial differential equations. What is most surprising here is that the terms of these numerical or functional series have unambiguous physical interpretation. They describe spatial displacements individually for each harmonic. For this reason, we called this modification of the bilinear-transformation method, as well in [10], β€œspatial.” Note that for some nonintegrable nonlinear evolution equations, such as the Kuramoto-Sivashinsky equation 𝑒𝑑+𝑒𝑒π‘₯+𝛼1𝑒π‘₯π‘₯+𝛼2𝑒π‘₯π‘₯π‘₯+𝛼3𝑒π‘₯π‘₯π‘₯π‘₯=0,(5.1) spatial displacements are absent. In this case, their role is taken by the β€œwave” displacements, that is, the wave number is presented as a convergent infinite series.

Appendices

A. Logarithmic Derivatives Expressed by the Hirota’s Bilinear Differential Operators 𝐷𝑑, 𝐷π‘₯

Consider the following: (ln𝜁)π‘₯π‘₯=𝐷2π‘₯πœβ‹…πœ2𝜁2;(ln𝜁)𝑑π‘₯=𝐷𝑑𝐷π‘₯πœβ‹…πœ2𝜁2,(ln𝜁)π‘₯π‘₯π‘₯π‘₯=𝐷4π‘₯πœβ‹…πœ2𝜁2ξ‚΅π·βˆ’62π‘₯πœβ‹…πœ2𝜁2ξ‚Ά2,(ln𝜁)π‘₯π‘₯π‘₯π‘₯π‘₯π‘₯=𝐷6π‘₯πœβ‹…πœ2𝜁2ξ‚΅π·βˆ’302π‘₯πœβ‹…πœ2𝜁2𝐷4π‘₯πœβ‹…πœ2𝜁2𝐷+1202π‘₯πœβ‹…πœ2𝜁2ξ‚Ά3.(A.1)

B. Identities for the Jacobi πœƒ-Functions

Consider the following: βˆžξ“π‘›=βˆ’βˆžπ‘ž2𝑛2=πœƒ3=πœƒ3ξ€·0,π‘ž2ξ€Έ,βˆžξ“π‘›=βˆ’βˆžπ‘žπ‘›2+(π‘›βˆ’1)2=π‘ž1/2πœƒ2=π‘ž1/2πœƒ2ξ€·0,π‘ž2ξ€Έ;βˆžξ“π‘›=βˆ’βˆžπ‘›2π‘ž2𝑛2=π‘žπœƒξ…ž32,βˆžξ“π‘›=βˆ’βˆž(2π‘›βˆ’1)2π‘žπ‘›2+(π‘›βˆ’1)2=2π‘ž3/2πœƒξ…ž2;βˆžξ“π‘›=βˆ’βˆžπ‘›4π‘ž2𝑛2=π‘žξ€·πœƒξ…ž3+π‘žπœƒ3ξ…žξ…žξ€Έ4,βˆžξ“π‘›=βˆ’βˆž(2π‘›βˆ’1)4π‘žπ‘›2+(π‘›βˆ’1)2=4π‘ž3/2ξ€·πœƒξ…ž2+π‘žπœƒ2ξ…žξ…žξ€Έ;βˆžξ“π‘›=βˆ’βˆžπ‘›6π‘ž2𝑛2=π‘žξ€·πœƒξ…ž3+3π‘žπœƒ3ξ…žξ…ž+π‘ž2πœƒ3ξ…žξ…žξ…žξ€Έ8,βˆžξ“π‘›=βˆ’βˆž(2π‘›βˆ’1)6π‘žπ‘›2+(π‘›βˆ’1)2=8π‘ž3/2ξ€·πœƒξ…ž2+3π‘žπœƒ2ξ…žξ…ž+π‘ž2πœƒ2ξ…žξ…žξ…žξ€Έ.(B.1)

Acknowledgment

The authors would like to thank the Academic Editors for the useful remarks that helped in improving the presentation of the paper.