Next Article in Journal
Generalized Čebyšev and Grüss Type Results in Weighted Lebesgue Spaces
Previous Article in Journal
Spectrum of Zariski Topology in Multiplication Krasner Hypermodules
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Breathers, Transformation Mechanisms and Their Molecular State of a (3+1)-Dimensional Generalized Yu–Toda–Sasa–Fukuyama Equation

1
School of Computer Science and Technology, China University of Mining and Technology, Xuzhou 221116, China
2
School of Mathematics, North University of China, Taiyuan 030051, China
3
School of Mathematics, China University of Mining and Technology, Xuzhou 221116, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(7), 1755; https://doi.org/10.3390/math11071755
Submission received: 13 February 2023 / Revised: 28 March 2023 / Accepted: 5 April 2023 / Published: 6 April 2023
(This article belongs to the Section Mathematical Physics)

Abstract

:
A (3+1)-dimensional generalized Yu–Toda–Sasa–Fukuyama equation is considered systematically. N-soliton solutions are obtained using Hirota’s bilinear method. The employment of the complex conjugate condition of parameters of N-soliton solutions leads to the construction of breather solutions. Then, the lump solution is obtained with the aid of the long-wave limit method. Based on the transformation mechanism of nonlinear waves, a series of nonlinear localized waves can be transformed from breathers, which include the quasi-kink soliton, M-shaped kink soliton, oscillation M-shaped kink soliton, multi-peak kink soliton, and quasi-periodic wave by analyzing the characteristic lines. Furthermore, the molecular state of the transformed two-breather is studied using velocity resonance, which is divided into three aspects, namely the modes of non-, semi-, and full transformation. The analytical method discussed in this paper can be further applied to the investigation of other complex high-dimensional nonlinear integrable systems.

1. Introduction

In the field of physics and mechanics, many scholars have become devoted to investigating the exact solutions of nonlinear integrable systems, which include solitons [1,2,3,4], breathers [5,6,7], lumps [8,9,10,11,12], quasi-periodic wave solutions [13,14,15,16], and so on. Scholars have undertaken a lot of work on exact solutions, and have summarized some effective methods, such as the Bäcklund transformation [17,18,19], inverse scattering method [20,21,22], Darboux transformation [23,24], algebraic geometry theory [25,26], Hirota’s bilinear method [27,28] and Lie symmetry analysis [29,30,31,32].
A breather can be counted as a sort of soliton that propagates periodically along the direction of the intersection with a soliton. Wang et al. investigated transformed one- and two-breathers, which are supported by a theoretical framework of a transformation mechanism of nonlinear waves [33,34,35,36]. It was found that there is a relationship of one-way transformation between breathers and nonlinear localized waves by adjusting parameters, including the quasi-anti-dark (kink) soliton, M-shaped (kink) soliton, oscillation M-shaped (kink) soliton, multi-peak (kink) soliton, and quasi-periodic wave. Furthermore, the molecular state of the transformed two-breather is discussed with the aid of velocity resonance [37]. The method of velocity resonance is widely used to study the relatively stable state of nonlinear localized waves. A lot of molecular phenomena, such as breather molecules (BMs), soliton molecules (SMs), soliton–breather molecules (SBMs), and lump-soliton molecules (LMs), have been acquired by Lou et al. in the references [38,39]. Li et al. discussed the soliton molecules of the (2+1)-dimensional B-type Kadomtsev–Petviashvili equation and the (2+1)-dimensional fifth-order Korteweg–de Vries(KdV) equation [40,41]. In addition, we studied the nonlinear superposition of the bifurcation of T-resonance Y-type solitons, lumps, breathers, and solitons with the aid of velocity resonance [42,43].
In this paper, we consider the following (3+1)-dimensional generalized Yu–Toda–Sasa– Fukuyama equation
r 1 ψ x t + r 2 ψ x x x z + r 3 ψ y y + r 4 ψ x x ψ z + r 5 ψ x ψ x z + r 6 ψ x y + r 7 ψ x z + r 8 ψ y z = 0 ,
where ψ = ψ ( x , y , z , t ) , r 1 , r 2 , r 3 , r 4 , r 5 , r 6 , r 7 and r 8 are eight arbitrary constants. When the parameters are taken as r 1 = 4 , r 2 = 1 , r 3 = 3 , r 4 = 2 , r 5 = 4 , r 6 = r 7 = r 8 = 0 , Equation (1) can be transformed into a (3+1)-dimensional potential Yu–Toda–Sasa– Fukuyama equation
4 ψ x t + ψ x x x z + 3 ψ y y + 2 ψ x x ψ z + 4 ψ x ψ x z = 0 ,
Yu et al. first proposed Equation (2) when they studied N-soliton solutions to the Bogoyavlenskii–Schiff equation [44], which can present an interfacial wave in a two-layer liquid or elastic quasi-plane wave in a lattice [45]. In recent years, many scholars in related fields have done a lot of work on the (3+1)-dimensional generalized Yu–Toda–Sasa–Fukuyama equation. Xia et al. investigated the dynamics of abundant solutions based on Hirota’s bilinear method [46]. Tian et al. obtained a bilinear form and bilinear auto-Bäcklund transformation with the aid of Hirota’s bilinear method and certain coefficient constraints. Then, breather and lump solutions were acquired [47]. Khalique et al. studied variational and non-variational approaches using Lie algebra [48].
This paper is devoted to studying the transformation mechanism of breathers and the molecular state of the transformed two-breather, which are based on parameter constraints and velocity resonance, respectively. Furthermore, it is worth noting that we need to select the appropriate parameters and carry out a lot of numerical simulations. The methods given in this paper can be extended to other high-dimensional integrable systems and further study the dynamic behaviors of ocean waves.
The organization of this paper is as follows. In Section 2, we obtain the bilinear form and N-soliton solutions using Hirota’s bilinear method. It can be concluded that Equation (3) is integrable in the sense of N-soliton solutions. In Section 3, a one-breather solution is studied, which, by taking the complex conjugate conditions to the two-soliton solution, then, with the aid of the transformation mechanism of nonlinear waves, is converted into a series of nonlinear localized waves. Furthermore, a one-lump wave is obtained by taking the long-wave limit to the one-breather solution [10,49,50,51]. In Section 4, the two-breather solution and its transformation mechanism are investigated systematically. Then, the molecular state of the transformed two-breather is discussed from the point of view of velocity resonance. Finally, some conclusions are presented in the last section.

2. Bilinear Form and the Soliton Solution

To facilitate the discussion of Equation (1), the parameters of Equation (1) are taken as r 1 = 1 , r 4 = r 5 . Then, it can be transformed into
ψ x t + r 2 ψ x x x z + r 3 ψ y y + r 4 ψ x x ψ z + r 4 ψ x ψ x z + r 6 ψ x y + r 7 ψ x z + r 8 ψ y z = 0 .
In this section, the bilinear form of Equation (3) and the N-soliton solutions are obtained. It can be noticed that the multi-soliton solutions are obtained only if a nonlinear wave equation is converted into a bilinear form through a dependent variable transformation. By means of dependent variable transformation
ψ = 6 r 2 r 4 ln ξ ( x , y , z , t ) x ,
Equation (3) can be transformed into a bilinear form [46]
D t D x + r 2 D x 3 D z + r 3 D y 2 + r 6 D x D y + r 7 D x D z + r 8 D y D z ξ · ξ = 0 ,
where the bilinear differential operators D x , D y , D z and D t are defined by
D x m D y n D z p D t q ξ ( x , y , z , t ) ζ ( x , y , z , t ) = ( x x ) m ( y y ) n ( z z ) p ( t t ) q ξ ( x , y , z , t )
× ζ ( x , y , z , t ) | x = x , y = y , z = z , t = t .
Furthermore, the N-soliton solutions of Equation (3) can be written as the following form
ξ = ξ N = μ = 0 , 1 exp j = 1 N μ j ε j + j < s N μ j μ s A j s ,
with phase variable ε j = u j ( x + v j y + h j z + w j t ) + ε j ( 0 ) , w j = ( r 2 u j 2 h j + r 3 v j 2 + r 6 v j + r 7 h j + r 8 v j h j ) ,   u j , v j , h j and ε j ( 0 ) are free constants, and
e A j s = r 2 ( u j u s ) ( ( 2 u j u s ) h j + ( u j 2 u s ) h s ) r 3 ( v j v s ) 2 r 8 ( v j v s ) ( h j h s ) r 2 ( u j + u s ) ( ( 2 u j + u s ) h j + ( u j + 2 u s ) h s ) r 3 ( v j v s ) 2 r 8 ( v j v s ) ( h j h s ) .
The symbol μ = 0 , 1 means summation over all possible combinations of μ j = 0 , 1 ( j = 1 N ) , when all μ j is zero, the corresponding term is 1; when μ j takes 0 and the rest of μ j takes 1, the corresponding term is exp ( j = 2 n ε j + 2 j < s n A j s ) , and the j < s N is the summation over all possible combinations of N elements in the specific condition j < s .

3. One-Breather Solution and Transformation Mechanism

In this section, by taking the complex conjugate conditions to the parameters of the two-soliton solution and imposing restrictions on parameters, a one-breather solution and transformation mechanism are analyzed systematically. Based on the above analysis, the two-soliton solution can be written as
ψ 2 = 6 r 2 r 4 ln ( 1 + e ε 1 + e ε 2 + e ε 1 + ε 2 + A 12 ) x .
By letting the parameters be
u 1 = a 1 + b 1 i , u 2 = a 1 b 1 i , v 1 = c 1 + d 1 i , v 2 = c 1 d 1 i , h 1 = e 1 + f 1 i ,
h 2 = e 1 f 1 i , ε 1 ( 0 ) = ln β 1 2 + γ 1 + κ 1 i , ε 2 ( 0 ) = ln β 1 2 + γ 1 κ 1 i ,
where a 1 , b 1 , c 1 , d 1 , e 1 , f 1 0 , β 1 > 0 , γ 1 and ζ 1 are arbitrary real constants, and ordering ξ 2 = ( 1 + e ε 1 + e ε 2 + e ε 1 + ε 2 + A 12 ) in (6), then, substituting (7) into ξ 2 , we obtain the one-breather solution
ξ 2 2 β 2 cosh ( θ 1 + 1 2 ln β 2 ) + β 1 cos ( Λ 1 ) ,
where
θ 1 = a 1 x + ( a 1 c 1 b 1 d 1 ) y + ( a 1 e 1 b 1 f 1 ) z + ( a 1 w 1 R b 1 w 1 I ) t + γ 1 ,
Λ 1 = b 1 x + ( a 1 d 1 + b 1 c 1 ) y + ( a 1 f 1 + b 1 e 1 ) z + ( a 1 w 1 I + b 1 w 1 R ) t + κ 1 , β 2 = β 1 2 G 1 4 ,
G 1 = 4 r 2 ( a 1 b 1 f 1 + 3 b 1 2 e 1 ) + 4 r 3 d 1 2 16 r 8 d 1 2 f 1 2 4 r 2 ( a 1 b 1 f 1 3 a 1 2 e 1 ) + 4 r 3 d 1 2 16 r 8 d 1 2 f 1 2 ,
w 1 R = r 2 ( a 1 2 e 1 b 1 2 e 1 2 a 1 b 1 f 1 ) r 3 ( c 1 2 d 1 2 ) r 6 c 1 r 7 e 1 r 8 ( c 1 e 1 d 1 f 1 ) ,
w 1 I = r 2 ( a 1 2 f 1 b 1 2 f 1 + 2 a 1 b 1 e 1 ) 2 r 3 c 1 d 1 r 6 d 1 r 7 f 1 r 8 ( c 1 f 1 + e 1 d 1 ) .
Then, substituting (8) into (6) leads to the one-breather solution of Equation (3)
ψ 2 = 6 r 2 ( 2 β 2 sinh ( θ 1 + 1 2 ln β 2 ) a 1 β 1 sin ( Λ 1 ) b 1 ) r 4 ( 2 β 2 cosh ( θ 1 + 1 2 ln β 2 ) + β 1 cos ( Λ 1 ) ) ,
which can be shown in Figure 1. The one-breather can be regarded as a one-soliton that propagates periodically along the direction of the intersection with the soliton.
Remark 1. 
Since the breathers, lump waves, and transformed localized waves to be studied have similar dynamic behaviors along the x , y , and z axes, in the following content, we take the ( x , y ) plane as an example to discuss all the subsequent problems. Similar results are obtained in the ( x , z ) plane and ( y , z ) plane.
To obtain the lump solution, we take
u 1 = a 1 + b 1 i , u 2 = a 1 b 1 i , v 1 = c 1 + d 1 i , v 2 = c 1 d 1 i , h 1 = e 1 + f 1 i ,
h 2 = e 1 f 1 i , ε 1 ( 0 ) = ln β 1 2 + π i , ε 2 ( 0 ) = ln β 1 2 π i ,
then ξ 2 can be rewritten as the form of
ξ 2 e θ 1 β 1 c o s ( Λ 1 ) + β 2 e θ 1 ,
where θ 1 , Λ 1 , β 2 are determined by (8). Then, taking a 1 0 , b 1 0 , β 1 = 2 , β 2 = 1 + ι 2 and expanding according to Taylor formula at ι = 0 , one has
ξ 2 ( θ 11 2 + Λ 11 2 + 1 ) ι 2 + O ( ι 3 ) , θ 11 = x + ( c 1 d 1 ) y + ( e 1 f 1 ) z + ( r 3 ( d 1 2 c 1 2 + 2 c 1 d 1 ) + r 6 ( d 1 c 1 ) + r 7 ( f 1 e 1 ) + r 8 ( c 1 f 1 + e 1 d 1 + d 1 f 1 c 1 e 1 ) ) t , Λ 11 = x + ( c 1 + d 1 ) y + ( e 1 + f 1 ) z + ( r 3 ( d 1 2 c 1 2 2 c 1 d 1 ) r 6 ( c 1 + d 1 ) r 7 ( e 1 + f 1 ) r 8 ( c 1 f 1 + e 1 d 1 d 1 f 1 + c 1 e 1 ) ) t .
By substituting (3) into (6), we have
ψ 2 = 6 r 2 ( 2 θ 11 + 2 Λ 11 ) r 4 ( θ 11 2 + Λ 11 2 + 1 ) .
According to the above discussion and Figure 2, we draw the following properties.
(1)
The lump solution (13) has one crest, one trough, which can be called binocular structure, and two characteristic lines that have the form of θ 11 + Λ 11 and θ 11 Λ 11 , as shown in Figure 2.
(2)
It can be noticed that the velocity of the one-lump wave (13) is available. On the ( x , y ) plane, the velocity of lump wave along the x axis is v x L = ( r 3 d 1 + r 8 f 1 2 ) ( c 1 2 + d 1 2 ) + r 7 ( c 1 f 1 e 1 d 1 ) d 1 , and the speed along the y axis is v y L = 2 r 3 c 1 d 1 + r 6 d 1 + r 7 f 1 + r 8 ( c 1 f 1 + e 1 d 1 ) d 1 .
Figure 2. (Color online) One-lump Solution (13) with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , a 1 = 0 , b 1 = 0 , c 1 = 1 , d 1 = 0.5 , e 1 = 0.3 , f 1 = 0.8 . (a) The three-dimensional stereograms when t = 0 , z = 0 . (b) the corresponding contour figure. (c) Characteristic lines figure, the red line ( y = 0 ) and the green line ( x + y = 0 ) are two characteristic lines of the one-lump wave.
Figure 2. (Color online) One-lump Solution (13) with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , a 1 = 0 , b 1 = 0 , c 1 = 1 , d 1 = 0.5 , e 1 = 0.3 , f 1 = 0.8 . (a) The three-dimensional stereograms when t = 0 , z = 0 . (b) the corresponding contour figure. (c) Characteristic lines figure, the red line ( y = 0 ) and the green line ( x + y = 0 ) are two characteristic lines of the one-lump wave.
Mathematics 11 01755 g002
Based on the analysis of the transformation mechanism of the breathers [33,34,35,36], some conclusions about the one-breather solution (9) of Equation (3) can be discussed.
(1)
It is obvious from expression (9) that the one-breather solution contains a trigonometric function ( sin , cos ) and a hyperbolic function ( cosh , sinh ) , in which the localized properties of the one-breather are controlled by the hyperbolic function, and the periodic properties are decided by the trigonometric function, so the one-breather can be considered to be the combination of the soliton wave and periodic wave.
(2)
Two characteristic lines of the one-breather have the form of θ 1 + 1 2 ln β 2 = a 1 x + ( a 1 c 1 b 1 d 1 ) y + ( a 1 e 1 b 1 f 1 ) z + ( a 1 w 1 R b 1 w 1 I ) t + γ 1 + 1 2 ln β 2 and Λ 1 = b 1 x + ( a 1 d 1 + b 1 c 1 ) y + ( a 1 f 1 + b 1 e 1 ) z + ( a 1 w 1 I + b 1 w 1 R ) t + κ 1 .
(3)
It can be noticed that the velocities of the soliton along the x axis and y axis can be written as v x s = a 1 ( a 1 w 1 R b 1 w 1 I ) a 1 2 + ( a 1 c 1 b 1 d 1 ) 2 and v y s = ( a 1 c 1 b 1 d 1 ) ( a 1 w 1 R b 1 w 1 I ) a 1 2 + ( a 1 c 1 b 1 d 1 ) 2 , respectively; the velocities of periodic wave along x axis and y axis are v x p = b 1 ( a 1 w 1 I + b 1 w 1 R ) b 1 2 + ( a 1 d 1 + b 1 c 1 ) 2 and v y p = ( a 1 d 1 + b 1 c 1 ) ( a 1 w 1 I + b 1 w 1 R ) b 1 2 + ( a 1 d 1 + b 1 c 1 ) 2 , respectively; the above results are discussed on the plane ( x , y ) , and similar conclusions can be obtained on the ( x , z ) plane and ( y , z ) plane.
(i)
If the relationship a 1 a 1 c 1 b 1 d 1 b 1 a 1 d 1 + b 1 c 1 0 is satisfied, i.e., d 1 0 , the two characteristic lines θ 1 + 1 2 ln β 2 and Λ 1 will not be parallel in the plane ( x , y ) , as shown in Figure 1.
(ii)
If the relationship a 1 a 1 c 1 b 1 d 1 b 1 a 1 d 1 + b 1 c 1 = 0 is satisfied, i.e., d 1 = 0 , the two characteristic lines θ 1 + 1 2 ln β 2 and Λ 1 will be parallel, as shown in Figure 3 and Figure 4. Under special conditions, the one-breather can be converted into a series of nonlinear waves which include quasi-kink soliton, M-shaped kink soliton, oscillation M-shaped kink soliton, multi-peak kink soliton and quasi-periodic wave. In Figure 3a–c, the one-breather will be a transformed quasi-kink soliton ( b 1 a 1 = 0.1 ) , which has one characteristic line and presents the shape of a ladder. In Figure 3d–f, the M-shaped kink soliton ( b 1 a 1 = 1 ) has two peaks and one valley, and appears in the shape of M climbing upward. With the increase of value of b 1 a 1 , the one-breather will become an oscillation M-shaped kink soliton ( b 1 a 1 = 3 ) , and the number of characteristic lines also increases. If the value of b 1 a 1 keeps growing, the periodicity will become obvious, then the one-breather will become an asymmetric multi-peak kink soliton ( b 1 a 1 = 5 ) , as shown in Figure 4a–f. When the value b 1 a 1 becomes very large, the one-breather will be transformed into a quasi-periodic wave ( b 1 a 1 = 1000 ) , as shown in Figure 4g–i. Gradually, with the values of b 1 a 1 increasing, their periodicity becomes more and more obvious and their locality almost disappears.
Remark 2. 
Figure 5 shows that quasi-kink soliton moves along the y axis at a different time. Then, the other figure of the transformation mechanism of breathers at a different time can also be obtained similarly.

4. Two-Breather Solution, Transformation Mechanism and Molecular State

4.1. Two-Breather Solution of Equation (3)

In this section, we investigate the transformation mechanism of the two-breather wave. The four-soliton solution can be written as
ξ 4 = 1 + e ε 1 + e ε 2 + e ε 3 + e ε 4 + e ε 1 + ε 2 + A 12 + e ε 1 + ε 3 + A 13 + e ε 1 + ε 4 + A 14 + e ε 2 + ε 3 + A 23 + e ε 2 + ε 4 + A 24 + e ε 3 + ε 4 + A 34 + e ε 1 + ε 2 + ε 3 + A 12 + A 23 + A 13 + e ε 1 + ε 2 + ε 4 + A 12 + A 24 + A 14 + e ε 1 + ε 3 + ε 4 + A 13 + A 34 + A 14 + e ε 2 + ε 3 + ε 4 + A 23 + A 34 + A 24 + e ε 1 + ε 2 + ε 3 + ε 4 + A 12 + A 23 + A 13 + A 23 + A 34 + A 24 .
Letting
u 3 = a 2 + b 2 i , u 4 = a 2 b 2 i , v 3 = c 2 + d 2 i , v 4 = c 2 d 2 i , h 3 = e 2 + f 2 i , h 4 = e 2 f 2 i ,
ε 3 ( 0 ) = ln β 3 2 + γ 2 + κ 2 i , ε 4 ( 0 ) = ln β 3 2 + γ 2 κ 2 i ,
where a 2 , b 2 , c 2 , d 2 , e 2 , f 2 0 , β 3 > 0 , γ 2 and κ 2 are arbitrary real constants, then substituting (10) and (15) into (14), we obtain
ξ 4 = 1 + β 1 e θ 1 cos ( Λ 1 ) + β 1 2 G 1 e 2 θ 1 4 + β 3 e θ 2 cos ( Λ 2 ) + β 3 2 G 2 e 2 θ 2 4
+ β 1 2 β 3 2 G 1 G 2 16 e 2 θ 1 + 2 θ 2 ( M 3 R 2 + M 3 I 2 ) ( M 4 R 2 + M 4 I 2 )
+ β 1 β 3 2 e θ 1 + θ 2 ( M 3 R cos ( Λ 1 + Λ 2 ) M 3 I sin ( Λ 1 + Λ 2 ) )
+ β 1 β 3 2 e θ 1 + θ 2 ( M 4 R cos ( Λ 1 Λ 2 ) M 4 I sin ( Λ 1 Λ 2 ) )
+ β 1 2 β 3 G 1 4 e 2 θ 1 + θ 2 ( ( M 3 R M 4 R + M 3 I M 4 I ) cos ( Λ 2 ) ( M 3 I M 4 R M 3 R M 4 I ) sin ( Λ 2 ) )
+ β 1 β 3 2 G 2 4 e θ 1 + 2 θ 2 ( ( M 3 R M 4 R M 3 I M 4 I ) cos ( Λ 1 ) ( M 3 I M 4 R + M 3 R M 4 I ) sin ( Λ 1 ) ) ,
where
θ 2 = a 2 x + ( a 2 c 2 b 2 d 2 ) y + ( a 2 e 2 b 2 f 2 ) z + ( a 2 w 2 R b 2 w 2 I ) t + γ 2 ,
Λ 2 = b 2 x + ( a 2 d 2 + b 2 c 2 ) y + ( a 2 f 2 + b 2 e 2 ) z + ( a 2 w 2 I + b 2 w 2 R ) t + κ 2 , β 4 = β 3 2 G 2 4 ,
G 2 = 4 r 2 ( a 2 b 2 f 2 + 3 b 2 2 e 2 ) + 4 r 3 d 2 2 16 r 8 d 2 2 f 2 2 4 r 2 ( a 2 b 2 f 2 3 a 2 2 e 2 ) + 4 r 3 d 2 2 16 r 8 d 2 2 f 2 2 ,
w 2 R = r 2 ( a 2 2 e 2 b 2 2 e 2 2 a 2 b 2 f 2 ) r 3 c 2 2 d 2 2 r 6 c 2 r 7 e 2 r 8 ( c 2 e 2 d 2 f 2 ) ,
w 2 I = r 2 ( a 2 2 f 2 b 2 2 f 2 + 2 a 2 b 2 e 2 ) 2 r 3 c 2 d 2 r 6 d 2 r 7 f 2 r 8 ( c 2 f 2 + e 2 d 2 ) .
M 3 R = Re ( e A 13 ) , M 3 I = Im ( e A 13 ) , M 4 R = Re ( e A 14 ) , M 4 I = Im ( e A 14 ) .
Unlike the one-breather, the two-breather has two sets of characteristic lines, i.e., θ 1 + 1 2 ln β 2 = a 1 x + ( a 1 c 1 b 1 d 1 ) y + ( a 1 e 1 b 1 f 1 ) z + ( a 1 w 1 R b 1 w 1 I ) t + γ 1 + 1 2 ln β 2 , Λ 1 = b 1 x + ( a 1 d 1 + b 1 c 1 ) y + ( a 1 f 1 + b 1 e 1 ) z + ( a 1 w 1 I + b 1 w 1 R ) t + κ 1 , and θ 2 + 1 2 ln β 4 = a 2 x + ( a 2 c 2 b 2 d 2 ) y + ( a 2 e 2 b 2 f 2 ) z + ( a 2 w 2 R b 2 w 2 I ) t + γ 2 + 1 2 ln β 4 , Λ 2 = b 2 x + ( a 2 d 2 + b 2 c 2 ) y + ( a 2 f 2 + b 2 e 2 ) z + ( a 2 w 2 I + b 2 w 2 R ) t + κ 2 .
Proposition 1. 
If one satisfies
a 1 a 1 c 1 b 1 d 1 b 1 a 1 d 1 + b 1 c 1 0 , a 2 a 2 c 2 b 2 d 2 b 2 a 2 d 2 + b 2 c 2 0 ,
the two sets of waves are two-breathers on the plane ( x , y ) .
(1)
If the two-breather satisfies a 1 a 1 c 1 b 1 d 1 a 2 a 2 c 2 b 2 d 2 = 0 , i.e., the two-breather is parallel, whereas they will only collide at a certain time due to different speeds, which is called the short-lived collision, as shown in Figure 6.
(2)
If the two-breather satisfies a 1 a 1 c 1 b 1 d 1 a 2 a 2 c 2 b 2 d 2 0 , i.e., the two-breather is not parallel; in other words, they will always be in a state of intersection, which is called the long-lived collision, as shown in Figure 7.
Figure 6. (Color online) The collision between parallel two-breathers with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , a 1 = 0.3 , b 1 = 0 , c 1 = 1.2 , d 1 = 0.8 , e 1 = 0.75 , f 1 = 0.75 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 , a 2 = 0.2 , b 2 = 0 , c 2 = 1.2 , d 2 = 1 , e 2 = 2 3 , f 2 = 2 3 , β 3 = 2 , γ 2 = 0 , κ 2 = 0 . (ac) are vertical view when t = 50 , t = 0 and t = 50 . (df) are contour plots.
Figure 6. (Color online) The collision between parallel two-breathers with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , a 1 = 0.3 , b 1 = 0 , c 1 = 1.2 , d 1 = 0.8 , e 1 = 0.75 , f 1 = 0.75 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 , a 2 = 0.2 , b 2 = 0 , c 2 = 1.2 , d 2 = 1 , e 2 = 2 3 , f 2 = 2 3 , β 3 = 2 , γ 2 = 0 , κ 2 = 0 . (ac) are vertical view when t = 50 , t = 0 and t = 50 . (df) are contour plots.
Mathematics 11 01755 g006
Figure 7. (Color online) The collision between non-parallel two-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , a 1 = 0.3 , b 1 = 1 , c 1 = 1.2 , d 1 = 1 , e 1 = 1 , f 1 = 1 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 , a 2 = 0.3 , b 2 = 1 , c 2 = 1.2 , d 2 = 0.5 , e 2 = 0.5 , f 2 = 1 , β 3 = 2 , γ 2 = 0 , κ 2 = 0 . (ac) are vertical view when t = 3 , t = 0 and t = 1 . (df) are contour plots.
Figure 7. (Color online) The collision between non-parallel two-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , a 1 = 0.3 , b 1 = 1 , c 1 = 1.2 , d 1 = 1 , e 1 = 1 , f 1 = 1 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 , a 2 = 0.3 , b 2 = 1 , c 2 = 1.2 , d 2 = 0.5 , e 2 = 0.5 , f 2 = 1 , β 3 = 2 , γ 2 = 0 , κ 2 = 0 . (ac) are vertical view when t = 3 , t = 0 and t = 1 . (df) are contour plots.
Mathematics 11 01755 g007

4.2. Transformation Mechanism of Two-Breather for Equation (3)

In this section, the transformation mechanism of the two-breather is studied systematically. According to the one-breather transformation mechanism, we know that the breathers can be converted into a series of nonlinear localized waves. Then, the investigation of the two-breather transformation mechanism can be divided into three aspects, including the modes of non-, semi-, and full transformation.
Regarding non-transformed modes, it is obvious that the two-breather will not be converted, as shown in Figure 6 and Figure 7. Regarding semi-transformed modes, that is one where the two-breather is transformed.
Proposition 2. 
If one satisfies
a 1 a 1 c 1 b 1 d 1 b 1 a 1 d 1 + b 1 c 1 = 0 , a 2 a 2 c 2 b 2 d 2 b 2 a 2 d 2 + b 2 c 2 0 ,
i.e., d 1 = 0 . Then, one of the two sets of waves is the breather, the other can be converted into a series of nonlinear waves including the quasi-kink soliton, M-shaped kink soliton, oscillation M-shaped kink soliton, and multi-peak kink soliton, as shown in Figure 8 and Figure 9. In Figure 8a–c, one of the two-breathers is transformed into the quasi-kink soliton. In Figure 8d–e, it is transformed into the M-shaped kink soliton. In Figure 9a–c, one of the two-breathers is transformed into the oscillation M-shaped kink soliton. In Figure 9d–e, it is transformed into the multi-peak kink soliton.
For the full-transformation modes, the two-breather will be transformed into a series of nonlinear waves.
Proposition 3. 
If one satisfies
a 1 a 1 c 1 b 1 d 1 b 1 a 1 d 1 + b 1 c 1 = 0 , a 2 a 2 c 2 b 2 d 2 b 2 a 2 d 2 + b 2 c 2 = 0 ,
i.e., d 1 = d 2 = 0 , then, the two-breather will be transformed into a series of nonlinear waves including the quasi-kink soliton, M-shaped kink soliton, oscillation M-shaped kink soliton, and multi-peak kink soliton, as shown in Figure 10. In Figure 10a–c, the breathers are all transformed into quasi-kink solitons. In Figure 10d,e, they are transformed into the quasi-kink soliton and M-shaped kink soliton. In Figure 10g–i, they are transformed into the M-shaped kink soliton and oscillation M-shaped kink soliton.
It is worth noting that the above propositions are based on a 1 a 1 c 1 b 1 d 1 a 2 a 2 c 2 b 2 d 2 = 0 , which means the characteristic lines θ 1 and θ 2 are parallel. The above two waves collide at a certain time due to their velocities being different, which is the so-called short-lived collision. If the characteristic lines θ 1 and θ 2 are not parallel, i.e., a 1 a 1 c 1 b 1 d 1 a 2 a 2 c 2 b 2 d 2 0 , which is called the long-lived collision due to two waves being in collision with time.

4.3. Molecular State of Transformed Two-Breather

Given the above discussion and analysis, we know that breathers can be transformed into a series of nonlinear localized waves that include the quasi-kink soliton, M-shaped kink soliton, oscillation M-shaped kink soliton, multi-peak kink soliton and quasi-periodic wave. Then, the different types of the molecular state of the transformed two-breather are investigated under the condition of the same speed, which can be called velocity resonance.
Then, we further consider the velocity resonance of the transformed two-breather. The propagation velocity of along the direction perpendicular to transformed two-breather is equal, i.e.,
v 1 x s 2 + v 1 y s 2 = v 2 x s 2 + v 2 y s 2 ,
then, the relative position of the transformed two-breather will not be changed by time t, where
v 1 x s = a 1 ( a 1 w 1 R b 1 w 1 I ) a 1 2 + ( a 1 c 1 b 1 d 1 ) 2 , v 2 x s = a 2 ( a 2 w 2 R b 2 w 2 I ) a 2 2 + ( a 2 c 2 b 2 d 2 ) 2 ,
v 1 y s = ( a 1 c 1 b 1 d 1 ) ( a 1 w 1 R b 1 w 1 I ) a 1 2 + ( a 1 c 1 b 1 d 1 ) 2 , v 2 y s = ( a 2 c 2 b 2 d 2 ) ( a 2 w 2 R b 2 w 2 I ) a 2 2 + ( a 2 c 2 b 2 d 2 ) 2 .
One important factor is that the molecular state is built on the condition of a 1 a 1 c 1 b 1 d 1 a 2 a 2 c 2 b 2 d 2 = 0 , i.e., the transformed two-breather is parallel. Furthermore, the distance of the transformed two-breather is not zero; otherwise, it will always be in coincidence. The distance has the form of
( γ 1 + 1 2 ln β 2 ) ( γ 2 + 1 2 ln β 4 ) a 1 2 + ( a 1 c 1 b 1 d 1 ) 2 .
Proposition 4. 
If the following conditions
a 1 a 1 c 1 b 1 d 1 a 2 a 2 c 2 b 2 d 2 = 0 , a 1 a 1 c 1 b 1 d 1 b 1 a 1 d 1 + b 1 c 1 0 , a 2 a 2 c 2 b 2 d 2 b 2 a 2 d 2 + b 2 c 2 0 ,
and
v B 1 s = v B 2 s , ( γ 1 + 1 2 ln β 2 ) ( γ 2 + 1 2 ln β 4 ) 0 ,
are satisfied, the two-breather will not be transformed, which is called the mode of non-transformation, and the distance of the two-breather will not be changed by time, as shown in Figure 11.
Proposition 5. 
If the following conditions
a 1 a 1 c 1 b 1 d 1 a 2 a 2 c 2 b 2 d 2 = 0 , a 1 a 1 c 1 b 1 d 1 b 1 a 1 d 1 + b 1 c 1 0 , a 2 a 2 c 2 b 2 d 2 b 2 a 2 d 2 + b 2 c 2 = 0 ,
and
v B 1 s = v B 2 s , ( γ 1 + 1 2 ln β 2 ) ( γ 2 + 1 2 ln β 4 ) 0 ,
are satisfied, one of the two-breathers will be transformed, the other will not be, which is called the mode of semi-transformation, and the distance of one-breather and a nonlinear localized wave will not be changed by time, as shown in Figure 12.
Proposition 6. 
If the following conditions
a 1 a 1 c 1 b 1 d 1 a 2 a 2 c 2 b 2 d 2 = 0 , a 1 a 1 c 1 b 1 d 1 b 1 a 1 d 1 + b 1 c 1 = 0 , a 2 a 2 c 2 b 2 d 2 b 2 a 2 d 2 + b 2 c 2 = 0 ,
and
v B 1 s = v B 2 s , ( γ 1 + 1 2 ln β 2 ) ( γ 2 + 1 2 ln β 4 ) 0 ,
are satisfied, the two-breather will be both transformed, which is called the mode of full transformation, and the distance of two nonlinear localized waves will not be changed by time, as shown in Figure 13.

5. Conclusions

In this paper, we focus on investigating breather solutions, transformation mechanisms, and the molecular state of the transformed two-breather. The N-soliton solutions are obtained using Hirota’s bilinear method. Then, it can be concluded that Equation (3) is integrable in the sense of N-soliton solutions. Furthermore, by using the complex conjugate conditions to the two- and four-soliton solution and imposing restrictions on the parameters, the one-breather solution, two-breather solution, and their transformation mechanism are analyzed systematically. The one-breather has two characteristic lines, θ 1 + 1 2 ln β 2 = a 1 x + ( a 1 c 1 b 1 d 1 ) y + ( a 1 e 1 b 1 f 1 ) z + ( a 1 w 1 R b 1 w 1 I ) t + γ 1 + 1 2 ln β 2 and Λ 1 = b 1 x + ( a 1 d 1 + b 1 c 1 ) y + ( a 1 f 1 + b 1 e 1 ) z + ( a 1 w 1 I + b 1 w 1 R ) t + κ 1 . If the condition a 1 a 1 c 1 b 1 d 1 b 1 a 1 d 1 + b 1 c 1 is not equal to zero, the one-breather will not be transformed, as shown in Figure 1. If the condition a 1 a 1 c 1 b 1 d 1 b 1 a 1 d 1 + b 1 c 1 is equal to zero, two characteristic lines are parallel, then the one-breather will be transformed. With the aid of the transformation mechanism of the nonlinear waves, the one-breather can be transformed into a series of nonlinear localized waves, such as the quasi-kink soliton, M-shaped kink soliton, oscillation M-shaped kink soliton, multi-peak kink soliton and quasi-periodic wave, as shown in Figure 3 and Figure 4. Then, the one-lump wave is obtained by taking the long-wave limit to the two-breather solution, as shown in Figure 2. The transformation mechanism of the two-breather is further studied similarly, as shown in Figure 8, Figure 9 and Figure 10. Furthermore, under the conditions of velocity resonance ( v 1 x s 2 + v 1 y s 2 = v 2 x s 2 + v 2 y s 2 ), the molecular state of the transformed two-breather is investigated systematically, which is shown in Figure 11, Figure 12 and Figure 13.
The phenomena presented in this paper are helpful to our further analysis of the complex dynamic behaviors in shallow-water waves, and play an important role in explaining the nonlinear phenomena existing in complex waves modeled by Equation (3). Furthermore, the dynamic behaviors of other high-dimensional integrable systems can be analyzed using the characteristic line method presented in this paper. It is worth noting that the characteristic line method cannot determine amplitude, which needs to be further studied.

Author Contributions

Methodology, J.Z., J.Y. and Y.Z.; Software, J.Y.; Validation, Y.Z.; Investigation, Z.Z.; Writing—original draft, J.Z. and J.Y.; Supervision, Z.Z.; Project administration, Z.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundations of China (No. 62206297) and the Fundamental Research Funds for the Central Universities (No. 2021QN1073).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mani Rajan, M.S.; Mahalingam, A. Nonautonomous solitons in modified inhomogeneous Hirota equation: Soliton control and soliton interaction. Nonlinear Dyn. 2015, 79, 2469–2484. [Google Scholar] [CrossRef]
  2. Ma, W.X. N-soliton solution of a combined pKP-BKP equation. J. Geom. Phys. 2021, 165, 104191. [Google Scholar] [CrossRef]
  3. Shen, Y.; Tian, B. Bilinear auto-Bäcklund transformations and soliton solutions of a (3+1)-dimensional generalized nonlinear evolution equation for the shallow water waves. Appl. Math. Lett. 2021, 122, 107301. [Google Scholar] [CrossRef]
  4. Mo, Y.F.; Ling, L.M.; Zeng, D.L. Data-driven vector soliton solutions of coupled nonlinear Schrödinger equation using a deep learning algorithm. Phys. Lett. A 2022, 421, 127739. [Google Scholar] [CrossRef]
  5. Wang, L.H.; Porsezian, K.; He, J.S. Breather and rogue wave solutions of a generalized nonlinear Schrödinger equation. Phys. Rev. E 2013, 87, 53202. [Google Scholar] [CrossRef] [Green Version]
  6. Feng, B.F.; Ma, R.Y.; Zhang, Y.J. General breather and rogue wave solutions to the complex short pulse equation. Phys. D 2022, 439, 133360. [Google Scholar] [CrossRef]
  7. Yusuf, A.; Sulaiman, T.A.; Alshomrani, A.S.; Baleanu, D. Breather and lump-periodic wave solutions to a system of nonlinear wave model arising in fluid mechanics. Nonlinear Dyn 2022, 110, 3655–3669. [Google Scholar] [CrossRef]
  8. Zhao, Z.L.; Chen, Y.; Han, B. Lump soliton, mixed lump stripe and periodic lump solutions of a (2+1)-dimensional asymmetrical Nizhnik-Novikov-Veselov equation. Modern Phys. Lett. B 2017, 31, 1750157. [Google Scholar] [CrossRef]
  9. Ma, W.X.; Zhou, Y. Lump solutions to nonlinear partial differential equations via Hirota bilinear forms. J. Differ. Equ. 2018, 264, 2633–2659. [Google Scholar] [CrossRef] [Green Version]
  10. Zhao, Z.L.; Yue, J.; He, L.C. New type of multiple lump and rogue wave solutions of the (2+1)-dimensional Bogoyavlenskii-Kadomtsev-Petviashvili equation. Appl. Math. Lett. 2022, 133, 108294. [Google Scholar] [CrossRef]
  11. Zhao, Z.L.; He, L.C.; Wazwaz, A.M. Dynamics of lump chains for the BKP equation describing propagation of nonlinear waves. Chin. Phys. B 2023, 32, 040501. [Google Scholar] [CrossRef]
  12. Zhao, Z.L.; He, L.C. Multiple lump molecules and interaction solutions of the Kadomtsev-Petviashvili I equation. Commun. Theor. Phys. 2022, 74, 105004. [Google Scholar] [CrossRef]
  13. Fan, E.G.; Hon, Y.C. Quasiperiodic waves and asymptotic behavior for Bogoyavlenskii’s breaking soliton equation in (2+1) dimensions. Phys. Rev. E 2008, 78, 36607. [Google Scholar] [CrossRef]
  14. Luo, L.; Fan, E.G. Bilinear approach to the quasi-periodic wave solutions of Modified Nizhnik-Novikov-Vesselov equation in (2+1) dimensions. Phys. Lett. A 2010, 374, 3001–3006. [Google Scholar] [CrossRef]
  15. Yue, J.; Zhao, Z.L. Solitons, breath-wave transitions, quasi-periodic waves and asymptotic behaviors for a (2+1)-dimensional Boussinesq-type equation. Eur. Phys. J. Plus 2022, 137, 914. [Google Scholar] [CrossRef]
  16. Li, J.; Sun, G.Q.; Jin, Z. Interactions of time delay and spatial diffusion induce the periodic oscillation of the vegetation system. Discrete Contin. Dyn. Syst. B 2022, 27, 2147–2172. [Google Scholar] [CrossRef]
  17. Zhou, T.Y.; Tian, B.; Chen, Y.Q.; Shen, Y. Painlevé analysis, auto-Bäcklund transformation and analytic solutions of a (2+1)-dimensional generalized Burgers system with the variable coefficients in a fluid. Nonlinear Dyn. 2022, 108, 2417–2428. [Google Scholar] [CrossRef]
  18. Dong, S.; Lan, Z.Z.; Gao, B.; Shen, Y.J. Bäcklund transformation and multi-soliton solutions for the discrete Korteweg-de Vries equation. Appl. Math. Lett. 2022, 125, 107747. [Google Scholar] [CrossRef]
  19. Yin, Y.H.; Lü, X.; Ma, W.X. Bäcklund transformation, exact solutions and diverse interaction phenomena to a (3+1)-dimensional nonlinear evolution equation. Nonlinear Dyn. 2022, 108, 4181–4194. [Google Scholar] [CrossRef]
  20. Mikhailov, A.V. The reduction problem and the inverse scattering method. Phys. D 1981, 3, 73–117. [Google Scholar] [CrossRef]
  21. Ning, T.K.; Chen, D.Y.; Zhang, D.J. The exact solutions for the nonisospectral AKNS hierarchy through the inverse scattering transform. Phys. A 2004, 339, 248–266. [Google Scholar] [CrossRef]
  22. Chen, M.S.; Fan, E.G.; He, J.S. l2-Sobolev space bijectivity of the scattering-inverse scattering transforms related to defocusing Ablowitz-Ladik systems. Phys. D 2023, 443, 133565. [Google Scholar] [CrossRef]
  23. Guo, B.L.; Ling, L.M.; Liu, Q.P. Nonlinear Schrödinger equation: Generalized Darboux transformation and rogue wave solutions. Phys. Rev. E 2012, 85, 26607. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Song, J.Y.; Xiao, Y.; Zhang, C.P. Darboux transformation, exact solutions and conservation laws for the reverse space-time Fokas-Lenells equation. Nonlinear Dyn. 2022, 107, 3805–3818. [Google Scholar] [CrossRef]
  25. Krichever, I.M. Methods of algebraic geometry in the theory of nonlinear equations. Russ. Math. Surv. 1977, 32, 185–213. [Google Scholar] [CrossRef]
  26. Cao, C.W.; Xiao, Y. Algebraic-geometric solution to (2+1)-dimensional Sawada-Kotera equation. Commun. Theor. Phys. 2008, 49, 34–36. [Google Scholar]
  27. Hirota, R. Exact solution of the Korteweg-de Vries equation for multiple collisions of solitons. Phys. Rev. Lett. 1971, 27, 1192–1194. [Google Scholar] [CrossRef]
  28. Hirota, R. The Direct Method in Soliton Theory; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  29. Bluman, G.W.; Yang, Z.Z. A symmetry-based method for constructing nonlocally related partial differential equation systems. J. Math. Phys. 2013, 54, 93504. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, H.Z.; Li, J.B.; Zhang, Q.X. Lie symmetry analysis and exact explicit solutions for general Burgers’ equation. J. Comput. Appl. Math. 2009, 228, 1–9. [Google Scholar] [CrossRef] [Green Version]
  31. Zhao, Z.L.; Han, B. Lie symmetry analysis, Bäcklund transformations, and exact solutions of a (2+1)-dimensional Boiti-Leon-Pempinelli system. J. Math. Phys. 2017, 58, 101514. [Google Scholar] [CrossRef]
  32. Zhao, Z.L.; He, L.C. Lie symmetry, nonlocal symmetry analysis, and interaction of solutions of a (2+1)-dimensional KdV-mKdV equation. Theor. Math. Phys. 2021, 206, 142–162. [Google Scholar] [CrossRef]
  33. Zhang, X.; Wang, L.; Liu, C.; Li, M.; Zhao, Y.C. High-dimensional nonlinear wave transitions and their mechanisms. Chaos 2020, 30, 113107. [Google Scholar] [CrossRef]
  34. Zhang, X.; Wang, L.; Chen, W.Q.; Yao, X.M.; Wang, X.; Zhao, Y.C. Dynamics of transformed nonlinear waves in the (3+1)-dimensional B-type Kadomtsev-Petviashvili equation I: Transitions mechanisms. Commun. Nonlinear Sci. Numer. Simul. 2022, 105, 106070. [Google Scholar] [CrossRef]
  35. Yin, Z.Y.; Tian, S.F. Nonlinear wave transitions and their mechanisms of (2+1)-dimensional Sawada-Kotera equation. Phys. D 2021, 427, 133002. [Google Scholar] [CrossRef]
  36. Zhang, D.D.; Wang, L.; Liu, L.; Liu, T.X.; Sun, W.R. Shape-changed propagations and interactions for the (3+1)-dimensional generalized Kadomtsev-Petviashvili equation in fluids. Commun. Theor. Phys. 2021, 73, 95001. [Google Scholar] [CrossRef]
  37. Yao, X.M.; Wang, L.; Zhang, X.; Zhang, Y.B. Dynamics of transformed nonlinear waves in the (3+1)-dimensional B-type Kadomtsev-Petviashvili equation II: Interactions and molecular waves. Nonlinear Dyn. 2023, 111, 4613–4629. [Google Scholar] [CrossRef]
  38. Jia, M.; Lin, J.; Lou, S.Y. Soliton and breather molecules in few-cycle-pulse optical model. Nonlinear Dyn. 2020, 100, 3745–3757. [Google Scholar] [CrossRef]
  39. Yan, Z.W.; Lou, S.Y. Special types of solitons and breather molecules for a (2+1)-dimensional fifth-order KdV equation. Commun. Nonlinear Sci. Numer. Simul. 2020, 91, 105425. [Google Scholar] [CrossRef]
  40. Wang, M.M.; Qi, Z.Q.; Chen, J.C.; Li, B. Resonance Y-shaped soliton and interaction solutions in the (2+1)-dimensional B-type Kadomtsev-Petviashvili equation. Internat. J. Modern Phys. B 2021, 35, 2150222. [Google Scholar] [CrossRef]
  41. Zhang, Z.; Yang, S.X.; Li, B. Soliton molecules, asymmetric solitons and hybrid solutions for (2+1)-dimensional fifth-order KdV equation. Chinese Phys. Lett. 2019, 36, 120501. [Google Scholar] [CrossRef]
  42. Zhao, Z.L.; He, L.C. Nonlinear superposition between lump waves and other waves of the (2+1)-dimensional asymmetrical Nizhnik-Novikov-Veselov equation. Nonlinear Dyn. 2022, 108, 555–568. [Google Scholar] [CrossRef]
  43. Yue, J.; Zhao, Z.L. Interaction solutions and molecule state between resonance Y-type solitons and lump waves, and transformed 2-breather molecular waves of a (2+1)-dimensional asymmetrical Nizhnik-Novikov-Veselov equation. Nonlinear Dyn. 2023, 111, 7565–7589. [Google Scholar] [CrossRef]
  44. Yu, S.J.; Toda, K.; Sasa, N.; Fukuyama, T. N soliton solutions to the Bogoyavlenskii-Schiff equation and a quest for the soliton solution in (3+1) dimensions. J. Phys. A: Math. Gen. 1998, 31, 3337–3347. [Google Scholar] [CrossRef]
  45. Yin, H.M.; Tian, B.; Chai, J.; Wu, X.Y.; Sun, W.R. Solitons and bilinear Bäcklund transformations for a (3+1)-dimensional Yu-Toda-Sasa-Fukuyama equation in a liquid or lattice. Appl. Math. Lett. 2016, 58, 178–183. [Google Scholar] [CrossRef]
  46. Guo, H.D.; Xia, T.C.; Hu, B.B. Dynamics of abundant solutions to the (3+1)-dimensional generalized Yu-Toda-Sasa-Fukuyama equation. Appl. Math. Lett. 2020, 105, 106301. [Google Scholar] [CrossRef]
  47. Shen, Y.; Tian, B.; Zhao, X.; Shan, W.R.; Jiang, Y. Bilinear form, bilinear auto-Bäcklund transformation, breather and lump solutions for a (3+1)-dimensional generalised Yu-Toda-Sasa-Fukuyama equation in a two-layer liquid or a lattice. Pramana-J. Phys. 2021, 95, 137. [Google Scholar] [CrossRef]
  48. Adeyemo, O.D.; Khalique, C.M.; Gasimov, Y.S. Variational and non-variational approaches with Lie algebra of a generalized (3+1)-dimensional nonlinear potential Yu-Toda-Sasa-Fukuyama equation in engineering and physics. Alex. Eng. J. 2023, 63, 17–43. [Google Scholar] [CrossRef]
  49. Ablowitz, M.J.; Satsuma, J. Solitons and rational solutions of nonlinear evolution equations. J. Math. Phys. 1978, 19, 2180–2186. [Google Scholar] [CrossRef]
  50. Satsuma, J.; Ablowitz, M.J. Two-dimensional lumps in nonlinear dispersive systems. J. Math. Phys. 1979, 20, 1496–1503. [Google Scholar] [CrossRef]
  51. Zhang, Z.; Guo, Q.; Li, B.; Chen, J.C. A new class of nonlinear superposition between lump waves and other waves for Kadomtsev-Petviashvili I equation. Commun. Nonlinear Sci. Numer. Simul. 2021, 101, 105866. [Google Scholar] [CrossRef]
Figure 1. (Color online) One-breather solution (9) with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , a 1 = 0.3 , b 1 = 1 , c 1 = 1.2 , d 1 = 1 , e 1 = 1 , f 1 = 1 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 . (a) The three-dimensional stereograms when t = 0 , z = 0 . (b) the corresponding contour figure. (c) Characteristic lines figure, the red line ( 0.3 x + 1.36 y + 0.3659935070 = 0 ) and the green line ( x + 0.9 y = 0 ) are two characteristic lines of one-breather.
Figure 1. (Color online) One-breather solution (9) with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , a 1 = 0.3 , b 1 = 1 , c 1 = 1.2 , d 1 = 1 , e 1 = 1 , f 1 = 1 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 . (a) The three-dimensional stereograms when t = 0 , z = 0 . (b) the corresponding contour figure. (c) Characteristic lines figure, the red line ( 0.3 x + 1.36 y + 0.3659935070 = 0 ) and the green line ( x + 0.9 y = 0 ) are two characteristic lines of one-breather.
Mathematics 11 01755 g001
Figure 3. (Color online) The transformation of one-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 . In (a), a 1 = 1 , b 1 = 0.1 , c 1 = 1 , d 1 = 0 , e 1 = 1.2 , f 1 = 2.5 . In (d), a 1 = 0.2 , b 1 = 0.2 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 . These two figures are three-dimensional stereograms when t = 0 , z = 0 . (b,e) are the corresponding contour figure. (c,f) show the wave moves along the y axis when x = 0 , t = 0 , z = 0 . .
Figure 3. (Color online) The transformation of one-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 . In (a), a 1 = 1 , b 1 = 0.1 , c 1 = 1 , d 1 = 0 , e 1 = 1.2 , f 1 = 2.5 . In (d), a 1 = 0.2 , b 1 = 0.2 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 . These two figures are three-dimensional stereograms when t = 0 , z = 0 . (b,e) are the corresponding contour figure. (c,f) show the wave moves along the y axis when x = 0 , t = 0 , z = 0 . .
Mathematics 11 01755 g003
Figure 4. (Color online) The transformation of one-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 . In (a), a 1 = 0.1 , b 1 = 0.3 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.1 , f 1 = 3.5 . In (d), a 1 = 0.04 , b 1 = 0.2 , c 1 = 0.5 , d 1 = 0 , e 1 = 0.01 , f 1 = 1 . In (g), a 1 = 0.004 , b 1 = 4 , c 1 = 1 , d 1 = 0 , e 1 = 0.01 , f 1 = 5 . These three figures are three-dimensional stereograms when t = 0 , z = 0 . (b,e,h) are the corresponding contour figure. (c,f,i) show the wave moves along the y axis when x = 0 , t = 0 , z = 0 . .
Figure 4. (Color online) The transformation of one-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 . In (a), a 1 = 0.1 , b 1 = 0.3 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.1 , f 1 = 3.5 . In (d), a 1 = 0.04 , b 1 = 0.2 , c 1 = 0.5 , d 1 = 0 , e 1 = 0.01 , f 1 = 1 . In (g), a 1 = 0.004 , b 1 = 4 , c 1 = 1 , d 1 = 0 , e 1 = 0.01 , f 1 = 5 . These three figures are three-dimensional stereograms when t = 0 , z = 0 . (b,e,h) are the corresponding contour figure. (c,f,i) show the wave moves along the y axis when x = 0 , t = 0 , z = 0 . .
Mathematics 11 01755 g004
Figure 5. (Color online) Quasi-kink soliton with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 , a 1 = 1 , b 1 = 0.1 , c 1 = 1 , d 1 = 0 , e 1 = 1.2 , f 1 = 2.5 . (ac) show quasi-kink soliton moves along the y axis when t = 1 , t = 0 and t = 1 , respectively.
Figure 5. (Color online) Quasi-kink soliton with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = 2 , γ 1 = 0 , κ 1 = 0 , a 1 = 1 , b 1 = 0.1 , c 1 = 1 , d 1 = 0 , e 1 = 1.2 , f 1 = 2.5 . (ac) show quasi-kink soliton moves along the y axis when t = 1 , t = 0 and t = 1 , respectively.
Mathematics 11 01755 g005
Figure 8. (Color online) The transformation of two-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , a 2 = 0.2 , b 2 = 0 , c 2 = 1.2 , d 2 = 1 , e 2 = 2 3 , f 2 = 2 3 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 , γ 1 = 20 , γ 2 = 10 . In (a), a 1 = 1 , b 1 = 0.1 , c 1 = 1.2 , d 1 = 0 , e 1 = 1.5 , f 1 = 2.5 . In (d), a 1 = 0.1 , b 1 = 0.1 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 . These two figures are three-dimensional stereograms of two-breather transformation when t = 0 , z = 0 . (b,e) are the corresponding contour figure. (c,f) show the wave moves along the x axis when y = 50 , t = 0 , z = 0 . .
Figure 8. (Color online) The transformation of two-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , a 2 = 0.2 , b 2 = 0 , c 2 = 1.2 , d 2 = 1 , e 2 = 2 3 , f 2 = 2 3 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 , γ 1 = 20 , γ 2 = 10 . In (a), a 1 = 1 , b 1 = 0.1 , c 1 = 1.2 , d 1 = 0 , e 1 = 1.5 , f 1 = 2.5 . In (d), a 1 = 0.1 , b 1 = 0.1 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 . These two figures are three-dimensional stereograms of two-breather transformation when t = 0 , z = 0 . (b,e) are the corresponding contour figure. (c,f) show the wave moves along the x axis when y = 50 , t = 0 , z = 0 . .
Mathematics 11 01755 g008
Figure 9. (Color online) The transformation of two-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 . In (a), a 1 = 0.3 , b 1 = 1 , c 1 = 1.2 , d 1 = 0 , e 1 = 1 , f 1 = 1 , γ 1 = 10 , a 2 = 0.2 , b 2 = 0 , c 2 = 1.2 , d 2 = 1 , e 2 = 2 3 , f 2 = 2 3 , γ 2 = 10 . In (d), a 1 = 0.04 , b 1 = 0.2 , c 1 = 0.5 , d 1 = 0 , e 1 = 0.01 , f 1 = 1 , γ 1 = 5 , a 2 = 0.2 , b 2 = 0 , c 2 = 0.5 , d 2 = 1 , e 2 = 0.005 , f 2 = 1 , γ 2 = 0 . These two figures are three-dimensional stereograms of two-breather transformation when t = 0 , z = 0 . (b,e) are the corresponding contour figure. (c,f) show the wave moves along the x axis when y = 50 , t = 0 , z = 0 . .
Figure 9. (Color online) The transformation of two-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 . In (a), a 1 = 0.3 , b 1 = 1 , c 1 = 1.2 , d 1 = 0 , e 1 = 1 , f 1 = 1 , γ 1 = 10 , a 2 = 0.2 , b 2 = 0 , c 2 = 1.2 , d 2 = 1 , e 2 = 2 3 , f 2 = 2 3 , γ 2 = 10 . In (d), a 1 = 0.04 , b 1 = 0.2 , c 1 = 0.5 , d 1 = 0 , e 1 = 0.01 , f 1 = 1 , γ 1 = 5 , a 2 = 0.2 , b 2 = 0 , c 2 = 0.5 , d 2 = 1 , e 2 = 0.005 , f 2 = 1 , γ 2 = 0 . These two figures are three-dimensional stereograms of two-breather transformation when t = 0 , z = 0 . (b,e) are the corresponding contour figure. (c,f) show the wave moves along the x axis when y = 50 , t = 0 , z = 0 . .
Mathematics 11 01755 g009
Figure 10. (Color online) The transformation of two-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 . In (a), a 1 = 0.5 , b 1 = 0.005 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 , γ 1 = 30 , a 2 = 0.3 , b 2 = 0.003 , c 2 = 1.2 , d 2 = 0 , e 2 = 0.5 , f 2 = 3.5 , γ 2 = 0 . In (d), a 1 = 0.5 , b 1 = 0.005 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 , γ 1 = 30 , a 2 = 0.3 , b 2 = 0.3 , c 2 = 1.2 , d 2 = 0 , e 2 = 0.5 , f 2 = 3.5 , γ 2 = 0 . In (g), a 1 = 0.5 , b 1 = 0.05 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 , γ 1 = 20 , a 2 = 0.3 , b 2 = 0.5 , c 2 = 1.2 , d 2 = 0 , e 2 = 0.5 , f 2 = 3.5 , γ 2 = 0 . These three figures are three-dimensional stereograms of the two-breather transformation when t = 0 , z = 0 . (b,e,h) are the corresponding contour figure. (c,f,i) show the wave moves along the x axis when y = 0 , t = 0 , z = 0 . .
Figure 10. (Color online) The transformation of two-breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 . In (a), a 1 = 0.5 , b 1 = 0.005 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 , γ 1 = 30 , a 2 = 0.3 , b 2 = 0.003 , c 2 = 1.2 , d 2 = 0 , e 2 = 0.5 , f 2 = 3.5 , γ 2 = 0 . In (d), a 1 = 0.5 , b 1 = 0.005 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 , γ 1 = 30 , a 2 = 0.3 , b 2 = 0.3 , c 2 = 1.2 , d 2 = 0 , e 2 = 0.5 , f 2 = 3.5 , γ 2 = 0 . In (g), a 1 = 0.5 , b 1 = 0.05 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 , γ 1 = 20 , a 2 = 0.3 , b 2 = 0.5 , c 2 = 1.2 , d 2 = 0 , e 2 = 0.5 , f 2 = 3.5 , γ 2 = 0 . These three figures are three-dimensional stereograms of the two-breather transformation when t = 0 , z = 0 . (b,e,h) are the corresponding contour figure. (c,f,i) show the wave moves along the x axis when y = 0 , t = 0 , z = 0 . .
Mathematics 11 01755 g010
Figure 11. (Color online) The molecular state between breather and breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 , a 1 = 0.3 , b 1 = 0 , c 1 = 1.2 , d 1 = 0.8 , e 1 = 1 , f 1 = 2.279 , γ 1 = 10 , a 2 = 0.2 , b 2 = 0 , c 2 = 1.2 , d 2 = 1 , e 2 = 2 3 , f 2 = 2 3 , γ 2 = 10 . (ac) are vertical view when t = 30 , t = 0 and t = 30 . (df) are contour plots.
Figure 11. (Color online) The molecular state between breather and breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 , a 1 = 0.3 , b 1 = 0 , c 1 = 1.2 , d 1 = 0.8 , e 1 = 1 , f 1 = 2.279 , γ 1 = 10 , a 2 = 0.2 , b 2 = 0 , c 2 = 1.2 , d 2 = 1 , e 2 = 2 3 , f 2 = 2 3 , γ 2 = 10 . (ac) are vertical view when t = 30 , t = 0 and t = 30 . (df) are contour plots.
Mathematics 11 01755 g011
Figure 12. (Color online) The molecular state between quasi-kink soliton and breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 , a 1 = 1 , b 1 = 0.1 , c 1 = 1.2 , d 1 = 0 , e 1 = 1.5 , f 1 = 2.5 , γ 1 = 10 , a 2 = 0.2 , b 2 = 0 , c 2 = 1.2 , d 2 = 1 , e 2 = 2 3 , f 2 = 5.559 , γ 2 = 10 . (ac) are vertical view when t = 8 , t = 0 and t = 10 . (df) are contour plots of (ac), respectively.
Figure 12. (Color online) The molecular state between quasi-kink soliton and breather with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 , a 1 = 1 , b 1 = 0.1 , c 1 = 1.2 , d 1 = 0 , e 1 = 1.5 , f 1 = 2.5 , γ 1 = 10 , a 2 = 0.2 , b 2 = 0 , c 2 = 1.2 , d 2 = 1 , e 2 = 2 3 , f 2 = 5.559 , γ 2 = 10 . (ac) are vertical view when t = 8 , t = 0 and t = 10 . (df) are contour plots of (ac), respectively.
Mathematics 11 01755 g012
Figure 13. (Color online) The molecular state between the quasi-kink soliton and quasi-kink soliton with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 , a 1 = 0.5 , b 1 = 0.005 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 , γ 1 = 10 , a 2 = 0.3 , b 2 = 0.003 , c 2 = 1.2 , d 2 = 0 , e 2 = 0.542 , f 2 = 3.5 , γ 2 = 10 . (ac) are vertical view when t = 30 , t = 0 and t = 30 . (df) are contour plots of (ac), respectively.
Figure 13. (Color online) The molecular state between the quasi-kink soliton and quasi-kink soliton with r 2 = r 3 = r 4 = r 6 = r 7 = r 8 = 1 , β 1 = β 3 = 2 , κ 1 = κ 2 = 0 , a 1 = 0.5 , b 1 = 0.005 , c 1 = 1.2 , d 1 = 0 , e 1 = 0.5 , f 1 = 3.5 , γ 1 = 10 , a 2 = 0.3 , b 2 = 0.003 , c 2 = 1.2 , d 2 = 0 , e 2 = 0.542 , f 2 = 3.5 , γ 2 = 10 . (ac) are vertical view when t = 30 , t = 0 and t = 30 . (df) are contour plots of (ac), respectively.
Mathematics 11 01755 g013
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, J.; Yue, J.; Zhao, Z.; Zhang, Y. Breathers, Transformation Mechanisms and Their Molecular State of a (3+1)-Dimensional Generalized Yu–Toda–Sasa–Fukuyama Equation. Mathematics 2023, 11, 1755. https://doi.org/10.3390/math11071755

AMA Style

Zhang J, Yue J, Zhao Z, Zhang Y. Breathers, Transformation Mechanisms and Their Molecular State of a (3+1)-Dimensional Generalized Yu–Toda–Sasa–Fukuyama Equation. Mathematics. 2023; 11(7):1755. https://doi.org/10.3390/math11071755

Chicago/Turabian Style

Zhang, Jian, Juan Yue, Zhonglong Zhao, and Yufeng Zhang. 2023. "Breathers, Transformation Mechanisms and Their Molecular State of a (3+1)-Dimensional Generalized Yu–Toda–Sasa–Fukuyama Equation" Mathematics 11, no. 7: 1755. https://doi.org/10.3390/math11071755

APA Style

Zhang, J., Yue, J., Zhao, Z., & Zhang, Y. (2023). Breathers, Transformation Mechanisms and Their Molecular State of a (3+1)-Dimensional Generalized Yu–Toda–Sasa–Fukuyama Equation. Mathematics, 11(7), 1755. https://doi.org/10.3390/math11071755

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop