Integrability, Multi-Solitary Wave Solutions and Riemann Theta Functions Periodic Wave Solutions of the Newell Equation

This paper systematically studies the complete integrability of the Newell equation. Using generalized Bell polynomials, the corresponding bilinear equation, bilinear B äcklund transformation, Lax pair, and multi-shock wave solutions are successfully obtained. In addition, using the multidimensional Riemann theta functions, the periodic wave solutions of the Newell equation are constructed. On this basis, the asymptotic behavior of the periodic wave solution is given, which is the relationship between the periodic wave solution and the solitary wave solution.

• The Newell Equation
• Bäcklund Transformation
• Lax Pair
• Solitary Wave Solution
• Periodic Wave Solution

1. Introduction

Since the pioneering introduction of solitons into nonlinear science, many experts and scholars have studied and explored its various aspects. In a great deal of work and efforts of scientists, the main focus has been on the search for exact solutions of nonlinear evolution equations and the study of integrability [1] - [7]. The application of the bilinear method, pioneered by the Japanese mathematician Hirota is very important in nonlinear theory [8]. Nakamura proposed a synthetic approach for constructing multi-periodic wave solutions of nonlinear equations based on the Hirota bilinear method. The advantage of this method is that it only depends on the bilinear form. However, it is not easy to seek a bilinear transformation of a nonlinear evolution equation. It is necessary to make an appropriate transformation, and considerable skill is required to find this transformation. To solve this problem, Lambert et al. introduced the concept of the Bell polynomials in algebra to nonlinear differential equations, which made it simple to construct bilinear forms of nonlinear equations [9]. This method not only simplifies the computational complexity but is also practical. Furthermore, it was found that based on the bilinear form, the Bäcklund transformation, Lax pair and infinite conservation laws can be obtained [10] [11]. At present, Bell polynomial theory and the Riemann theta function are extended to superintegrable systems, discrete systems and higher-order integrable systems [12]. In addition, Wang wrote the corresponding package, so that the Bell polynomial theory method became more procedural and simplified [13].

The Newell equation is a very important equation in nonlinear mathematical physics. It is widely used in hydrodynamics, and for coastal wave simulation, nonlinear motion of charged particles in electromagnetic fields, propagation of solitary waves in media, and vibration of a one-dimensional nonlinear lattice. Its general form is

u t t − c 0 2 u x x − α u x u x x − β u x x x x = 0 , (1)

using the homogeneous balance method, the Bäcklund transformation and exact solutions of the Newell equation were presented in [14] [15]. The content of this paper is arranged as follows. In Section 2, the bilinear expression of Newell equation is constructed by using multidimensional Bell polynomials. In Section 3, multiple solitary wave solutions and Riemann theta function periodic wave solutions are obtained by using the bilinear expression of Newell equation. In Section 4, the relationship between periodic wave solution and solitary wave solution is further analyzed, and it is strictly proved that the periodic solution tends to soliton solution under the small amplitude limit. In Section 5, the Bäcklund transformation and Lax pair of Newell equation are derived.

2. Bilinear Form

First, we do the following transformation by introducing the potential function q

u = m ( t ) q x ,

where m = m ( t ) is a function to be determined, bring the above into Equation (1) you can get

m t t q x + 2 m t q x , t + m q x , 2 t − c 0 2 m q 3 x − α m t 2 q 2 x q 3 x − β m q 5 x = 0 ,

integrating it with respect to x once, and taking the integral constant to zero, can derive the following equation

m t t q + 2 m t q t + m q 2 t − c 0 2 m q 2 x − 0.5 α m t 2 q 2 x 2 − β m q 4 x = 0 , (2)

take m ( t ) = 6 β α , then Equation (2) can be written as a combination of P-polynomials,

E ( q ) = P 2 t ( q ) − c 0 2 P 2 x − β P 4 x = 0.

By using the following transformation

q = 2 ( ln f ) ↔ u = m ( t ) q 2 x = 6 β α q 2 x = 12 β α ( ln f ) x ,

the bilinear form of Equation (1) can be obtained as

( D t 2 − c 0 2 D x 2 − β D x 4 ) f ⋅ f = 0 . (3)

3. Solitary Wave Solutions and Riemann Theta Functions Periodic Wave Solutions

In this section, we use the bilinear form of the Newell equation to construct the multi-solitary wave solutions and periodic wave solutions.

3.1. Multi-Solitary Wave Solutions

Once the bilinear form of Equation (1) is given, the multi-solitary wave solution of Equation (1) can easily be obtained by using symbolic computation and Hirota method. The multi-solitary wave solutions of the Newell equation can be written as

u = 12 β α ( ln f ) x , (4)

with

f = ∑ ρ = 0 , 1 exp ( ∑ j = 1 N ρ j η j + ∑ 1 ≤ j < i ≤ N ρ i ρ j A i j )     ,

e A i j = − ( v i − v j ) 2 − c 0 2 ( μ i − μ j ) 2 − β ( μ i − μ j ) 4 ( v i + v j ) 2 − c 0 2 ( μ i + μ j ) 2 − β ( μ i + μ j ) 4 ,

η j = μ j x + v j t + c j ,     v j = ± c 0 2 μ j 2 + β μ j 4 ,     ( 1 ≤ j < i ≤ N ) , (5)

in which μ j , v j are all free constants, and ∑ ​ ρ = 0 , 1 sums all possible combinations of ρ i , ρ j = 0 , 1 .

When N = 1 , one-solitary wave solutions of Equation (1) is given by

u = 12 β α ( 1 + e η ) x ,     η = μ x ± c 0 2 μ 2 + β μ 4 t + c , (6)

where u , c are free constants (Figure 1).

3.2. Riemann Theta Function Periodic Wave Solutions

Here, we obtain one-periodic wave solution of the Newell equation and analyze its asymptotic properties. To do this, we introduce the following one-Riemann theta function:

ϑ ( ξ , τ ) = ∑ n = − ∞ + ∞ e π i n 2 τ + 2 π i n ξ , (7)

in which the phase variable ξ = k x + ω t + ε , and Im ( τ ) > 0 .

To construct the periodic wave solution of Equation (1), we consider the more general form of the bilinear equation as

( D t 2 − c 0 2 D x 2 − β D x 4 + δ ) f ⋅ f = 0 , (8)

where δ is an arbitrary nonzero constant.

Theorem 1. If ϑ ( ξ , τ ) is a one-Riemann theta function (7) and ξ = k x + ω t + ε , the Newell equation allows a one-periodic wave solution as follows:

u = 12 β α ∂ x ln ϑ ( ξ , τ ) , (9)

with

ω 2 = b 1 a 22 − b 2 a 12 a 11 a 22 − a 12 a 21 ,     δ = b 1 a 21 − b 2 a 11 a 12 a 21 − a 11 a 22 ,

and

a 11 = − ∑ n = − ∞ + ∞ 16 n 2 π 2 ℘ 2 n 2 ,     a 12 = ∑ n = − ∞ + ∞ ℘ 2 n 2 ,

a 21 = − ∑ n = − ∞ + ∞ 4 ( 2 n − 1 ) 2 π 2 ℘ 2 n 2 − 2 n + 1 ,     a 22 = ∑ n = − ∞ + ∞ ℘ 2 n 2 − 2 n + 1 ,

b 1 = ∑ n = − ∞ + ∞ ( 256 β n 4 π 4 k 4 − 16 c 0 2 n 2 π 2 k 2 ) ℘ 2 n 2 ,

b 2 = − ∑ n = − ∞ + ∞ ( 16 β ( 2 n − 1 ) 4 π 4 k 4 − 4 c 0 2 ( 2 n − 1 ) 2 π 2 k 2 ) ℘ 2 n 2 − 2 n + 1 ,     ℘ = e π i τ , (10)

in which the parameters k , ε , τ are free.

Proof. For Equation (1), we consider the Riemann theta function solutions of the form (7), according to the Theorem 1 in Ref. [10] and the bilinear equation (8), the parameters k , ω , ε satisfy the following system

∑ n = − ∞ + ∞ ( − 16 n 2 π 2 ω 2 + 16 c 0 2 n 2 π 2 k 2 − 256 β n 4 π 4 k 4 + δ ) e 2 π i n 2 τ = 0 ,

∑ n = − ∞ + ∞ ( − 4 ( 2 n − 1 ) 2 π 2 ω 2 + 4 c 0 2 ( 2 n − 1 ) 2 π 2 k 2 − 16 β ( 2 n − 1 ) 4 π 4 k 4 + δ ) e π i ( 2 n 2 − 2 n + 1 ) τ = 0.

According to the representations in (10), the upper expressions are naturally written in the following matrix form:

( a 11 a 12 a 21 a 22 ) ( ω 2 δ ) = ( b 1 b 2 ) , (11)

solving the above linear system results in one-periodic wave solution (9) of Equation (1) (Figure 2).

4. Asymptotic Analysis

In this section, we use asymptotic analysis theory to discuss the relationship between one-solitary wave solution and one-periodic wave solution. It is worth noting that if the coefficient matrix, right end vector and solution vector of the system (6.14) can be expanded by a power series with parameter ℘ , the asymptotic behavior of the periodic wave solution can be analyzed by using the method of small parameter expansion. Via the results from Ref. [10], we expand system (6.14) by parameter ℘

( a 11 a 12 a 21 a 22 ) = A 0 + A 1 ℘ + A 2 ℘ 2 + ⋯

( b 1 b 2 ) = B 0 + B 1 ℘ + B 2 ℘ 2 + ⋯

( ω 2 δ ) = X 0 + X 1 ℘ + X 2 ℘ 2 + ⋯ (12)

based on Equation (10), rewrite the functions a i j , b j ( i , j = 1 , 2 ) as series expansion with respect to ℘

a 11 = − 32 π 2 ( ℘ 2 + 4 ℘ 8 + 9 ℘ 18 + ⋯ + n 2 ℘ 2 n 2 + ⋯ ) ,

a 12 = 1 + 2 ( ℘ 2 + ℘ 8 + ℘ 18 + ⋯ + ℘ 2 n 2 + ⋯ ) ,

a 21 = − 8 π 2 ( ℘ + 9 ℘ 5 + 25 ℘ 13 + ⋯ + ( 2 n − 1 ) 2 ℘ 2 n 2 − 2 n + 1 + ⋯ ) ,

a 22 = 2 ( ℘ + ℘ 5 + ℘ 13 + ⋯ + ℘ 2 n 2 − 2 n + 1 + ⋯ ) ,

b 1 = 32 π 2 ( ( 16 β   π 2 k 4 − c 0 2 k 2 ) ℘ 2 + ( 256 β   π 2 k 4 − 4 c 0 2 k 2 ) ℘ 8 + ⋯       + ( 16 β n 4 π 2 k 4 − c 0 2 n 2 k 2 ) ℘ 2 n 2 + ⋯ ) ,

b 2 = 8 π 2 ( ( 4 β   π 2 k 4 − c 0 2 k 2 ) ℘ + ( 4 β ( 2 n − 1 ) 4 π 2 k 4               − c 0 2 ( 2 n − 1 ) 2 k 2 ) ℘ 2 n 2 − 2 n + 1 + ⋯ ) , (13)

combined with Equations (12) and (13), we get

A 1 = ( 0 0 − 8 π 2 2 ) ,     A 2 = ( − 32 π 2 2 0 0 ) ,

A 3 = A 4 = 0 , ⋯     B 1 = ( 0 8 π 2 Δ 1 ) ,

B 2 = ( 32 π 2 Δ 2 0 ) ,     B 5 = ( 0 72 π 2 Δ 3 ) ,     B 0 = B 3 = B 4 = 0 , ⋯ (14)

in which Δ i are presented as follows:

Δ 1 = 4 β   π 2 k 4 − c 0 2 k 2 ,     Δ 2 = 16 β   π 2 k 4 − c 0 2 k 2 ,     Δ 3 = 36 β   π 2 k 4 − c 0 2 k 2 .

Then, using Equation (4.15) in Ref. [10], we can calculate the following formula

X 0 = ( − Δ 1 0 ) ,     X 2 = ( − 8 Δ 1 k − 2 32 π k − 1 Δ 1 ) ,

X 4 = ( − 39 k − 2 Δ 1 − 9 k − 1 Δ 3 192 π 2 k − 1 Δ 1 ) ,     X 1 = X 3 = 0 , ⋯ (15)

combined with Equation (12) and the above expression, we can derive

ω 2 = − k − 1 Δ 1 − 8 k − 2 Δ 1 ℘ 2 + ( − 64 k − 3 Δ 1 − 25 k − 2 Δ 1 − 9 k − 1 Δ 3 ) ℘ 4 + o ( ℘ 4 ) ,

δ = 32 π k − 1 Δ 1 ℘ 2 + ( − 256 π 2 k − 2 Δ 1 + 64 π 2 k − 1 Δ 1 ) ℘ 4 + o ( ℘ 4 ) . (16)

Further, rewrite the Riemann theta function ϑ ( ξ ) to

ϑ ( ξ ) = 1 + ( e 2 π i ξ + e − 2 π i ξ ) ℘ + ( e 4 π i ξ + e − 4 π i ξ ) ℘ 4 + ⋯ (17)

with the following transformation

k = μ 2 π i ,     ε = c + π τ 2 π i , (18)

it is easy to get that Equation (17) can be rewritten to

ϑ ( ξ , τ ) = 1 + e ξ ˜ + ( e − ξ ˜ + e 2 ξ ˜ ) ℘ 2 + ( e − 2 ξ ˜ + e 3 ξ ˜ ) ℘ 6 + ⋯ → 1 + e ξ ˜ , (19)

where ξ ˜ = 2 π i ξ − π τ = μ x + 2 π i ω   t + c .

Combining Eqs. (17)-(19), we can get

ξ ˜ → μ x ± c 0 2 μ 2 + β μ 4 t + c = η ,     when     ℘ → 0 ,

ξ → η + π τ 2 π i ,     when     ℘ → 0 ,

which means

ϑ ( ξ , τ ) → 1 + e η ,     when     ℘ → 0 .

According to the above analysis, we can conclude that under the limit condition ℘ → 0 , the one-periodic wave solution (9) is exactly inclined to the one-solitary wave solution (6).

5. The Bӓcklund Transformation and Related Lax Pair

Bilinear Bäcklund transformation plays an important role in constructing solutions of nonlinear equations and characterizing the integrability of given systems. Next, we derive the Bäcklund transformation and related Lax pair of the Newell equation.

Let q = 2 ln g , q ′ = 2 ln f be two different solutions of Equation (1), then we can get

E ( q ′ ) − E ( q ) = ( q ′ − q ) 2 t − c 0 2 ( q ′ − q ) 2 x − 3 β ( q ′ − q ) 2 x ( q ′ + q ) 2 x − β ( q ′ − q ) 4 x = 0 ,

the above formula can be regarded as the bilinear Bäcklund transformation of Equation (1). Based on this, we introduce two new independent variables

v = q ′ − q 2 = ln ( f / g ) ,     w = q ′ + q 2 = ln ( f g ) ,

according to multidimensional binary Bell polynomials, namely Y -polynomials [9], the corresponding two-field condition is rewritten as

E ( q ′ ) − E ( q ) = E ( v + w ) − E ( w − v ) = v 2 t − c 0 2 v 2 x − 6 β v 2 x w 2 x − β v 4 x = − β ∂ x Y 3 x ( v , w ) + R ( v , w ) = 0 , (20)

where

R ( v , w ) = − 3 β W r o n s k i a n [ Y 2 x ( v , w ) , Y x ( v ) ] + v 2 t − c 0 2 v 2 x .

In order to represent Equation (20) as a Y -polynomial Bäcklund transformation, we need to add a constraint condition, we can take this constraint as

k Y 2 x ( v , w ) + Y t ( v , w ) = λ ,

in which k = ± − 3 β and λ are free parameters.

According to the above constraints, we can get

R ( v , w ) = ∂ x [ − k Y x , t ( v , w ) + ∂ x ( λ k − c 0 2 ) Y x ( v , w ) ] ,

by the above relations, the two-field condition (20) is decomposed into Y -polynomial type Bäcklund transformation as

k Y 2 x ( v , w ) + Y t ( v , w ) = λ , ∂ x [ − β Y 3 x ( v , w ) − k Y x , t ( v , w ) + ( λ k − c 0 2 ) Y x ( v , w ) ] = 0.

On the basis of the relation between Y -polynomials and Hirota bilinear operator [9], we can directly obtain the bilinear Bäcklund transformation of the Newell equation as

( k D x 2 + D t − λ ) f ⋅ g = 0 , ( − β D x 3 − k D x D t + ( λ k − c 0 2 ) D x + δ ) f ⋅ g = 0 , (21)

in which k = ± − 3 β and λ , δ are arbitrary parameters.

Based on related theories of Bell polynomials, and by making use of the Hopf-Cole transformation v = ln ψ , we have

Y x ( v ) = ψ x ψ ,     Y 2 x ( v , w ) = q 2 x + ψ 2 x ψ ,

Y x , y ( v , w ) = q x , y + ψ x y ψ ,     Y y ( v ) = ψ y ψ ,

Y 3 x ( v , w ) = 3 q 2 x ψ x ψ + ψ 3 x ψ ,     Y t ( v ) = ψ t ψ ,

thus, the system (21) can be linearized as a Lax pair with parameters λ , δ

( k P 1 + ∂ t ) ψ = k ( q 2 x ψ + ψ 2 x ) + ψ t = λ ψ , ( P 2 − k ∂ x t ) ψ = ( − 3 β q 2 x − c 0 2 ) ψ x − β ψ 3 x − k q x , t ψ − k ψ x , t + δ ψ = 0.

6. Conclusion

In this paper, we use the binary Bell polynomials to study the integrable properties of the Newell equation from many aspects, such as the bilinear equation, bilinear Bäcklund transformation and Lax pair. In Section 3, we use the Riemann theta functions and some of the results in Ref. [10] to obtain the periodic wave solutions of the Newell equation. In addition, by considering the small amplitude limit, we further obtain the asymptotic behaviour between the one-periodic wave solution and the one-solitary wave solution. It is noteworthy that the Hirota bilinear method, the Riemann theta function and the Bell polynomials play important roles in constructing the solutions of nonlinear evolution equations.

Acknowledgements

This research was supported by the Higher Educational Scientific Research Projects of Inner Mongolia Autonomous Region (Grant No. NJZY20248; NJZY22307).

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this article.

References

1. Lü, X. and Ma, W.X. (2016) Study of Lump Dynamics Based on a Dimensionally Reduced Hirota Bilinear Equation. Nonlinear Dynamics, 85, 1217-1222. https://doi.org/10.1007/s11071-016-2755-8

2. Mao, J.J., Tian, S.F., Zou, L. and Zhang, T.T. (2019) Bilinear Formalism, Lump Solution, Lumpoff and Instant on Rogue Wave Solution of a (3+1)-Dimensional B-Type Kadomtsev-Petviashvili Equation. Nonlinear Dynamics, 95, 3005-3017. https://doi.org/10.1007/s11071-018-04736-2

3. Ma, W.X. (2021) N-Soliton Solution and the Hirota Condition of a (2+1)-Dimensional Combined Equation. Mathematics and Computers in Simulation, 190, 270-279. https://doi.org/10.1016/j.matcom.2021.05.020

4. Rao, J.G., Chow, K.W., Mihalache, D. and He, J.S. (2021) Completely Resonant Collision of Lumps and Line Solitons in the Kadomtsev-Petviashvili I Equation. Studies in Applied Mathematics, 147, 1007-1035. https://doi.org/10.1111/sapm.12417

5. Peng, W.Q., Tian, S.F. and Zhang, T.T. (2018) Analysis on Lump, Lumpoff and Rogue Waves with Predictability to the (2+1)-Dimensional B-Type Kadomtsev-Petviashvili Equation. Physics Letters A, 382, 2701-2708. https://doi.org/10.1016/j.physleta.2018.08.002

6. Zhou, Y. (2020) Periodic Solitary Wave Solutions of the (2+1)-Dimensional Variable-Coefficient Caudrey-Dodd-Gibbon-Kotera-Sawada Equation. Open Journal of Applied Sciences, 10, 60-68. https://doi.org/10.4236/ojapps.2020.103005

7. Dong, Z.Z. and Wang, L. (2018) Solitary Wave Solutions and Rational Solutions for Modified Zakharov-Kuznetsov Equation with Initial Value Problem. Journal of Applied Mathematics and Physics, 6, 949-959. https://doi.org/10.4236/jamp.2018.65081

8. Hirota, R. (2004) Direct Methods in Soliton Theory. Springer, Berlin. https://doi.org/10.1017/CBO9780511543043

9. Lambert, F. and Springael, J. (2008) From Soliton Equations to Their Zero Curvature Formulation. Acta Applicandae Mathematicae, 102, 147-178. https://doi.org/10.1007/s10440-008-9209-3

10. Tian, S.F. and Zhang, H.Q. (2010) Riemann Theta Functions Periodic Wave Solutions and Rational Characteristics for the Nonlinear Equations. Journal of Mathematical Analysis and Applications, 371, 585-608. https://doi.org/10.1016/j.jmaa.2010.05.070

11. Ma, P.L., Tian, S.F., Tu, J.M. and Xu, M.J. (2015) On the Integrability and Quasi-Periodic Wave Solutions of the Boussinesq Equation in Shallow Water. European Physical Journal Plus, 130, 98. https://doi.org/10.1140/epjp/i2015-15098-0

12. Dong, C. and Deng, S.F. (2020) The Quasi-Periodic Solutions for the Supersymmetric Variable-Coefficient KdV Equation. Modern Physics Letters B, 34, Article ID: 2050171. https://doi.org/10.1142/S0217984920501717

13. Miao, Q., Wang, Y.H., Chen, Y. and Yang, Y.Q. (2014) PDEBellII: A Maple Package for Finding Bilinear Forms, Bilinear B&#228;cklund Transformations, Lax Pairs and Conservation Laws of the KdV-Type Equations. Computer Physics Communications, 185, 357-367. https://doi.org/10.1016/j.cpc.2013.09.005

14. Yuan, P. and Deng, W.D. (2017) B&#228;cklund Transformation and Exact Solutions of the Newell Equation. Journal of Jishou University (Natural Science Edition), 38, 1-3.

15. Shi Y.Q. ,et al. (2017)The Exact Solutions of Newell Equation Journal of Guangxi University of Science and Technology 28, 6-12.