Next Article in Journal
New Modification of Polar Nonlinear Optical Iodate Fluoride PbF(IO3), the Family MX(IO3), M = Bi, Ba, Pb, X = O, F, (OH) Related to Aurivillius Phases and Similar Iodates
Previous Article in Journal
Symmetry in Σ Hyperon Decay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Construction of Novel Bright-Dark Solitons and Breather Waves of Unstable Nonlinear Schrödinger Equations with Applications

1
School of Energy and Power Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
2
Department of Mathematics and Statistics, Sub-Campus Depalpur, University of Agriculture, Faisalabad 38000, Pakistan
3
Faculty of Science, Jiangsu University, Zhenjiang 212013, China
4
Institute of Mathematics, Khawaja Fareed University of Engineering and Information Technology, Rahim Yar Khan 64200, Pakistan
5
Department of Mathematics Near East Boulevard, Mathematics Research Center, Near East University, Mersin 99138, Turkey
6
Department of Computer Science and Mathematics, Lebanese American University, Beirut 1102 2801, Lebanon
7
Department of Mathematics, Art and Science Faculty, Siirt University, Siirt 56100, Turkey
8
Department of Computer Engineering, Biruni University, Topkapı, Istanbul 34010, Turkey
9
Department of Mathematics, Azarbaijan Shahid Madani University, Tabriz 53751-71379, Iran
10
Department of Medical Research, China Medical University Hospital China Medical University, Taichung 40402, Taiwan
11
Department of Electricity and Electronics, Institute of Research and Development of Processes, Faculty of Science and Technology, University of the Basque Country, 48940 Leioa, Spain
*
Authors to whom correspondence should be addressed.
Submission received: 30 November 2022 / Revised: 20 December 2022 / Accepted: 26 December 2022 / Published: 30 December 2022
(This article belongs to the Section Mathematics)

Abstract

:
The unstable nonlinear Schrödinger equations (UNLSEs) are universal equations of the class of nonlinear integrable systems, which reveal the temporal changing of disruption in slightly stable and unstable media. In current paper, an improved auxiliary equation technique is proposed to obtain the wave results of UNLSE and modified UNLSE. Numerous varieties of results are generated in the mode of some special Jacobi elliptic functions and trigonometric and hyperbolic functions, many of which are distinctive and have significant applications such as pulse propagation in optical fibers. The exact soliton solutions also give information on the soliton interaction in unstable media. Furthermore, with the assistance of the suitable parameter values, various kinds of structures such as bright-dark, multi-wave structures, breather and kink-type solitons, and several periodic solitary waves are depicted that aid in the understanding of the physical interpretation of unstable nonlinear models. The various constructed solutions demonstrate the effectiveness of the suggested approach, which proves that the current technique may be applied to other nonlinear physical problems encountered in mathematical physics.

1. Introduction

The nonlinear Schrödinger equations (NLSEs) are fundamental nonlinear canonical evolution equations that are extensively used in several disciplines of both theoretical and applied sciences. The nonlinear equations are often employed in the modeling of many physical events, including the transmission line across optical fibers [1,2,3,4,5,6,7], attraction-repulsion chemotaxis models [8,9], the traveling of Langmuir waves in hot plasma [10], the propagation of rogue waves [11,12,13], depth water waves in oceans [14], ecology, and numerous others [15]. Partial differential equations (PDEs) are also derived for a population whose individuals randomly experience births, deaths, and size changes. PDEs are among the best mathematical tools for analyzing spatiotemporal processes. These nonlinear equations deal with non-linearity and dispersive effects at the same time and are one of the simplest forms of evolution equations. Due to their integrable characteristics, these nonlinear equations are widely used for the description of dynamics of localized stationary along with pulsating waves envelops. These nonlinear equations are known as universal equations because they explain how waves move through many physical systems. These equations are crucial tools for comprehending the physical analogies and variations in the nonlinear behaviour of dispersive waves. Rogue waves in nonlinear dispersive media are very thought provoking, for example, [16,17,18,19,20]. These rogue waves were initially examined in the background of water waves, and after that, they became an interesting topic for in-depth studies of nonlinear optics. The NLSEs are a basic part of the explanation of the construction of rogue waves that have been found in fiber optics [21,22], as well as in deep water [23,24,25,26].
With the achievements in the field of nonlinear science, many systematic and efficient methods have been formulated and improved for the finding of exact and approximate solutions of nonlinear evolution equations (NLEEs). These equations include the modified extended direct algebraic method [27], the modified simple equation method [28], the Jacobi elliptic function method [29,30], the modified extended mapping method [31], the semi-inverse variation method [32], the generalized Kudryashov method [33,34], the Bilinear residual network method [35], the modified extended tanh method [36], the sine-Gordon equation method [37], the tri-prong scheme [38], the Darboux transformation [39,40], the explicit spectral collection method [41], the generalized extended rational expansion method [42], the binary bell polynomials [43], the exponential fitting method [44], the auxiliary equation method [45], the improved auxiliary equation method [46], the modified F-expansion method [47], and many more [48,49].
The UNLSEs are NLSEs with space and time shuffled. Their role arises in plasma and elaborates the two-layer baroclinic instabilities [50,51]. In [50], the author showed that in the absence of dissipation of the equilibrated finite amplitude state causes oscillation, both of the mean flow and the baroclinic wave. The authors of [51] presented a nonlinear theory of the propagation of acoustic waves in piezoelectric semiconductors, based on an asymptotic expansion in multiple time and space scales. The UNLSE [51] has the following form:
i U t + U x x + 2 ν | U | 2 U 2 γ U = 0 ,
and the modified unstable nonlinear Schrödinger equation (mUNLSE) illustrates certain instabilities of modulated wave-trains [52], which have the following form:
i U t + U x x + 2 ν | U | 2 U 2 γ U x t = 0 ,
where real valued arbitrary constants are γ and ν . These UNLSEs depict a temporal evolution of disruption in slightly stable or unstable media. Different methods are used to find exact wave solutions of unstable NLSE such as the modified extended mapping method [53], the ( G / G ) -expansion method [54], and many others. Researchers in [28] used the extended simple equation technique to construct the exact solutions of unstable NLSEs and to modify unstable NLSEs. The modified exponential rational function method was used by the researchers in [55] to obtain the traveling wave results of unstable NLSEs.
In this article, bright-dark solitons, breather-type waves, and other wave solutions of UNLSE and mUNLSE are achieved with the help of the suggested improved auxiliary equation technique. Several novel wave solutions of these models have been constructed by using the proposed technique in the form of the Jacobi elliptic function and trigonometric and hyperbolic functions, and many of them have important applications in mathematical physics. Additionally, plotted visuals are used to further explain the various innovative structures of the solutions that were found, including bright, dark, and kink-type solitons, and numerous periodic solitary waves. The main advantage of this method is that it provides solutions for different types of waves as well as the Jacobi elliptic function. Any nonlinear PDEs, including fractional-order nonlinear PDEs, can be solved using this method. The many successful findings demonstrate the effectiveness of the suggested method, demonstrating its applicability to various nonlinear dynamical models.
The rest of the paper is arranged as follows: In Section 2, the description of proposed method is given. The application of the suggested approach on UNLSE and mUNLSE is described in Section 3 and Section 4. The findings, discussion, and physical interpretation are covered thoroughly in Section 5. Finally, Section 6 presents the conclusion.

2. Proposed Improved Auxiliary Equation Method

In this section, a brief explanation of the improved auxiliary equation technique [46] for the construction of wave solutions of nonlinear evolution equations is stated. Consider a general nonlinear evolution equation in the following form as
G ( U , U x , U t , U x t , U x x , U t t , ) = 0 ,
where polynomial G contains U ( x , t ) along with its partial derivatives U x , U t , U x t , U x x , U t t and so on.
Step 1: 
Consider the traveling wave transformation that converts U ( x , t ) into Φ ( η ) as
U ( x , t ) = Φ ( η ) , η = κ x ± λ t ,
where Equation (4) has arbitrary constants κ and λ . Using Equation (4) on Equation (3), it changes into an ordinary differential equation (ODE) as
H Φ , Φ , Φ , = 0 ,
where ( ) symbolizes the derivatives of Φ ( η ) with respect to η .
Step 2: 
Assume that Equation (5) has the following solution as
Φ ( η ) = i = 0 2 n c i Z i ( η ) ,
where the arbitrary constants n have a positive integer value; this value can be determined with the help of the homogeneous balancing principle on Equation (5), and Z ( η ) satisfies the new ansatz equation
Z ( η ) = a 0 + a 2 Z ( η ) 2 + a 4 Z ( η ) 4 + a 6 Z ( η ) 6 ,
where the constants a 0 , a 2 , a 4 , a 6 are arbitrary. Assume that Equation (7) has solutions as
Z ( η ) = 1 2 a 4 a 6 ( 1 ± z ( η ) ) ,
where z ( η ) can be written in terms of Jacobi elliptic functions s n ( η , m ) , c n ( η , m ) , and d n ( η , m ) . These Jacobi elliptic functions have m as its modulus; its value is 0 < m < 1 . As m achieves the value 0 or 1, the functions transform into trigonometric and hyperbolic functions.
Step 3: 
Use Equation (6) along with Equation (7) into Equation (5), collecting the coefficient of Z i ( η ) and setting equal to zero, it will give algebraic equations that have the parameters a 0 , a 2 , a 4 , a 6 , c i , κ , and λ . These algebraic equations will be solved by software Mathematica, and the values of the parameters are obtained.
Step 4: 
Put the parameters along with Equation (8) into Equation (6), and the solution of Equation (3) is achieved.

3. Implementation of the Improved Auxiliary Equation Method on UNLSE

The UNLSE is a sort of NLSE with space and time shuffled. The role of this kind arises in plasma and elaborates the baroclinic instability of two layers and their lossless symmetrices [50,51]. As Equation (1) is complex, suppose the transformation in the following form as
U ( x , t ) = Φ ( η ) e i ζ , η = κ x + λ t a n d ζ = μ x + σ t ,
where κ , λ , μ , and σ are the arbitrary constants. Applying the complex wave transformation Equation (9) on Equation (1) and separating into real and imaginary parts yields
κ 2 Φ ( η ) 2 γ + μ 2 + σ Φ ( η ) + 2 ν Φ 3 ( η ) = 0 , and λ = 2 k μ .
Using the homogeneous balancing principal on Equation (10) gives the positive integer of value n = 1 . Assume the solution in the following form:
Φ ( η ) = c 0 + c 1 Z ( η ) + c 2 Z ( η ) 2 .
Substituting Equation (11) along Equation (7) into Equation (10) yields a system of algebraic equations in parameters a 0 , a 2 , a 4 , a 6 , c 0 , c 1 , c 2 , κ , λ , μ , σ , γ and ν . By solving the system of algebraic equations with the help of Mathematica software, one obtains the values of parameters as
ν = a 2 2 2 a 0 a 4 a 2 κ 2 2 c 0 2 , μ = a 2 κ 2 2 γ σ + 3 κ 2 a 2 2 2 a 0 a 4 , a 6 = a 4 a 2 a 2 2 2 a 0 a 4 4 a 0 , c 1 = 0 , c 2 = a 2 c 0 + a 2 2 c 0 2 2 a 0 a 4 c 0 2 a 0 .
Using the parameters values of Equation (12) in Equation (8) generates
Z ( η ) = c 0 ( 1 ± z ( η ) ) a 2 + a 2 2 2 a 0 a 4 c 0 a 2 a 2 2 2 a 0 a 4 .
By substituting the values z ( η ) in the mode of the Jacobi elliptic functions, and by using Equation (11), the following cases are generated as
Case 1: 
If a 0 = a 4 3 m 2 1 32 a 6 2 m 2 , a 2 = a 4 2 5 m 2 1 16 a 6 m 2 , a 6 > 0 , then
U 11 ( η ) = ( c 0 a 2 + a 2 2 2 a 0 a 4 c 0 ( 1 ± sn ( τ η ) ) a 2 a 2 2 2 a 0 a 4 ) e i ζ ,
or
U 12 ( η ) = ( c 0 a 2 + a 2 2 2 a 0 a 4 c 0 1 ± 1 m sn ( τ η ) a 2 a 2 2 2 a 0 a 4 ) e i ζ .
Substituting m 1 , in Equation (15)
U 13 ( η ) = ( c 0 a 2 + a 2 2 2 a 0 a 4 c 0 ( 1 ± coth ( τ η ) ) a 2 a 2 2 2 a 0 a 4 ) e i ζ ,
where η = κ ( x 2 μ t ) and ζ = ( a 2 κ 2 2 γ σ + 3 κ 2 a 2 2 2 a 0 a 4 ) x + σ t .
Case 2: 
If a 0 = a 4 3 1 m 2 32 a 6 2 , a 2 = a 4 2 5 m 2 16 a 6 and a 6 > 0 , then
U 21 ( η ) = ( c 0 a 2 + a 2 2 2 a 0 a 4 c 0 ( 1 ± m sn ( τ η ) ) a 2 a 2 2 2 a 0 a 4 ) e i ζ ,
or
U 22 ( η ) = ( c 0 a 2 + a 2 2 2 a 0 a 4 c 0 1 ± 1 sn ( τ η ) a 2 a 2 2 2 a 0 a 4 ) e i ζ .
Substituting m 1 , in Equation (17)
U 23 ( η ) = ( c 0 a 2 + a 2 2 2 a 0 a 4 c 0 ( 1 ± tanh ( τ η ) ) a 2 a 2 2 2 a 0 a 4 ) e i ζ ,
where η = κ ( x 2 μ t ) and ζ = ( a 2 κ 2 2 γ σ + 3 κ 2 a 2 2 2 a 0 a 4 ) x + σ t .
Case 3: 
If a 0 = a 4 3 32 a 6 2 1 m 2 , a 2 = a 4 2 4 m 2 5 16 a 6 m 2 1 , a 6 > 0 , then
U 31 ( η ) = ( c 0 a 2 + a 2 2 2 a 0 a 4 c 0 1 ± 1 cn ( τ η ) a 2 a 2 2 2 a 0 a 4 ) e i ζ .
U 32 ( η ) = ( c 0 a 2 + a 2 2 2 a 0 a 4 c 0 1 ± dn ( τ η ) 1 m 2 sn ( τ η ) a 2 a 2 2 2 a 0 a 4 ) e i ζ .
Substituting m 0 in Equation (21)
U 33 ( η ) = ( c 0 a 2 + a 2 2 2 a 0 a 4 c 0 ( 1 ± csc ( τ η ) ) a 2 a 2 2 2 a 0 a 4 ) e i ζ ,
where η = κ ( x 2 μ t ) and ζ = ( a 2 κ 2 2 γ σ + 3 κ 2 a 2 2 2 a 0 a 4 ) x + σ t .
The auxiliary Equation (7) has some special solutions, which are
Case 4: 
If a 0 = a 2 2 m 2 1 m 2 a 4 2 m 2 1 , a 6 = 0 , a 4 < 0 , a 2 2 m 2 1 > 0 , then the parameters values are as follows:
a 2 = 2 γ + μ 2 + σ κ 2 , a 4 = c 1 2 ν κ 2 , c 0 = 0 , c 2 = 0 .
Using the parameters values in Equation (23) generates
U 41 ( η ) = m 2 2 γ + μ 2 + σ 2 m 2 1 ν cn η 2 γ + μ 2 + σ κ 2 2 m 2 1 .
Case 5: 
If a 0 = 0 , a 6 = 0 and a 2 > 0 , a 4 < 0 , then the parameter values are as follows
σ = a 2 κ 2 2 γ μ 2 , a 4 = c 1 2 ν κ 2 , c 0 = c 2 = 0 .
Using the parameter values in Equation (25) generates
U 51 ( η ) = 2 γ + μ 2 + σ ν sec h η 2 γ + μ 2 + σ κ 2 e i ζ .
Case 6: 
If a 0 = a 2 2 4 a 4 , a 6 = 0 , and a 2 > 0 , a 4 > 0 , then the following parameters values are solved:
a 2 = 2 γ + μ 2 + σ κ 2 , c 0 = 0 , c 1 = ± i a 4 κ ν , c 2 = 0 .
Using parameters values Equation (27) generates
U 61 ( η ) = ± i 2 γ + μ 2 + σ 2 ν tan η 2 γ + μ 2 + σ 2 κ 2 e i ζ .

4. Implementation of the Improved Auxiliary Equation Method on mUNLSE

The mUNLSE describes certain instabilities of modulated wave-trains and depicts a temporal evolution of disruption in slightly stable or unstable media [52]. As Equation (2) is also complex, suppose that the traveling wave transformation has the following form:
U ( x , t ) = Ψ ( η ) e i ζ , η = κ x + λ t and ζ = ρ x + ω t ,
where arbitrary constants in the above equation are κ , ρ , λ , and ω real. Applying the complex wave transformation Equation (29) on Equation (2) and separating it into real and imaginary parts yields
κ ( κ γ λ ) Ψ ( η ) + γ ρ ω ρ 2 ω Ψ ( η ) + 2 ν Ψ 3 ( η ) = 0 , λ = κ ( γ ω 2 ρ ) 1 ρ γ .
Using the homogeneous balancing principal on Equation (30) gives a positive integer of the value n = 1 . Assume the solution in the following form:
Ψ ( η ) = c 0 + c 1 Z ( η ) + c 2 Z ( η ) 2 .
Substituting Equation (31) along Equation (7) into Equation (30) yields a system of algebraic in parameters a 0 , a 2 , a 4 , a 6 , c 0 , c 1 , c 2 , κ , λ , ρ , ω , γ , and ν . On solving the system of algebraic equations with the aid of Mathematica software, the values of parameters are obtained as
ν = 2 a 4 2 κ ( κ γ λ ) a 2 2 2 a 0 a 4 a 2 c 2 2 , ρ = 4 a 2 κ ( κ γ λ ) + ω γ 2 ω 4 + 12 a 2 2 2 a 0 a 4 κ ( κ γ λ ) + γ ω 2 , a 6 = a 4 2 2 a 2 a 2 2 2 a 0 a 4 , c 0 = a 2 + a 2 2 2 a 0 a 4 c 2 2 a 4 , c 1 = 0 .
Using the parameter values of Equation (32) in Equation (8) generates
Z ( η ) = a 2 a 2 2 2 a 0 a 4 ( 1 ± z ( η ) ) a 4 2 .
By substituting the values z ( η ) in the mode of Jacobi elliptic functions, and by using Equation (31), the following cases are generated as
Case 1: 
If a 0 = a 4 3 m 2 1 32 a 6 2 m 2 , a 2 = a 4 2 5 m 2 1 16 a 6 m 2 , a 6 > 0 , then
U 11 ( η ) = ± sn ( τ η ) a 2 2 2 a 0 a 4 a 2 c 2 2 a 4 e i ζ ,
or
U 12 ( η ) = ± 1 m sn ( τ η ) a 2 2 2 a 0 a 4 a 2 c 2 2 a 4 e i ζ .
Substituting m 1 in Equation (35), we obtain
U 13 ( η ) = ± coth ( τ η ) a 2 2 2 a 0 a 4 a 2 c 2 2 a 4 e i ζ ,
where η = κ x + ( γ ω 2 ρ ) 1 ρ γ t and ζ = 4 a 2 κ ( κ γ λ ) + ω γ 2 ω 4 + 12 a 2 2 2 a 0 a 4 κ ( κ γ λ ) + γ ω 2 x + ω t .
Case 2: 
If a 0 = a 4 3 1 m 2 32 a 6 2 , a 2 = a 4 2 5 m 2 16 a 6 and a 6 > 0 , then
U 21 ( η ) = ± ( m sn ( τ η ) ) a 2 2 2 a 0 a 4 a 2 c 2 2 a 4 e i ζ ,
or
U 22 ( η ) = ± 1 sn ( τ η ) a 2 2 2 a 0 a 4 a 2 c 2 2 a 4 e i ζ .
Substituting m 1 in Equation (37), we obtain
U 23 ( η ) = ± tanh ( τ η ) a 2 2 2 a 0 a 4 a 2 c 2 2 a 4 e i ζ ,
where η = κ x + ( γ ω 2 ρ ) 1 ρ γ t and ζ = 4 a 2 κ ( κ γ λ ) + ω γ 2 ω 4 + 12 a 2 2 2 a 0 a 4 κ ( κ γ λ ) + γ ω 2 x + ω t .
Case 3: 
If a 0 = a 4 3 32 a 6 2 1 m 2 , a 2 = a 4 2 4 m 2 5 16 a 6 m 2 1 , a 6 > 0 , then
U 31 ( η ) = ± 1 cn ( τ η ) a 2 2 2 a 0 a 4 a 2 c 2 2 a 4 e i ζ ,
or
U 32 ( η ) = ± dn ( τ η ) 1 m 2 sn ( τ η ) a 2 2 2 a 0 a 4 a 2 c 2 2 a 4 e i ζ .
Substituting m 0 in Equation (41), we obtain
U 33 ( η ) = ± csc ( τ η ) a 2 2 2 a 0 a 4 a 2 c 2 2 a 4 e i ζ ,
where η = κ x + ( γ ω 2 ρ ) 1 ρ γ t and ζ = 4 a 2 κ ( κ γ λ ) + ω γ 2 ω 4 + 12 a 2 2 2 a 0 a 4 κ ( κ γ λ ) + γ ω 2 x + ω t .
The auxiliary Equation (7) has some special solutions, which are as
Case 4: 
If a 0 = a 2 2 m 2 1 m 2 a 4 2 m 2 1 , a 6 = 0 , a 4 < 0 , a 2 2 m 2 1 > 0 then the parameters values are as follows:
ρ = 1 γ , a 2 = 1 γ 2 κ ( κ γ λ ) , a 4 = c 1 2 ν κ ( κ γ λ ) , c 0 = 0 , c 2 = 0 .
Using the parameter values in Equation (43) generates
U 41 ( η ) = m 2 γ 2 2 m 2 1 ν cn η 1 γ 2 κ 2 m 2 1 ( κ γ λ ) .
Case 5: 
If a 0 = 0 , a 6 = 0 , and a 2 > 0 , a 4 < 0 , then the parameter values are as follows:
ω = a 2 γ κ λ a 2 κ 2 + ρ 2 γ ρ 1 , a 4 = c 1 2 ν κ ( κ γ λ ) , c 0 = 0 , c 2 = 0 .
Using the parameter values in Equation (45) generates
U 51 ( η ) = a 2 κ ( κ γ λ ) ν sec h a 2 η e i ζ .
Case 6: 
If a 0 = a 2 2 4 a 4 , a 6 = 0 and a 2 > 0 , a 4 > 0 , then the parameters values are as follows
ρ = 1 γ , a 2 = 1 γ 2 κ ( κ γ λ ) , c 0 = 0 , c 1 = ± i a 4 κ κ γ λ ν , c 2 = 0 .
Using parameters values in Equation (47) generates
U 61 ( η ) = ± i 1 ν γ 2 tan η 1 2 γ 2 κ ( κ γ λ ) 2 e i ζ .

5. Physical Interpretation

It has been seen that the obtained solutions of UNLSE and mUNLSE using the currently suggested method, which is the improved auxiliary equation method, are different from those obtained by using alternative methods. Our findings are innovative and straightforward. The noticeable features of obtaining the results of the proposed method are unique a structure of solutions and Equation (7), which produce distinct solutions. Some works on unstable NLSEs and their exact wave solutions are as follows. The UNLSE was solved by utilizing the modified extended mapping method [53]. The ( G / G ) -expansion scheme was used in [54] to obtain wave results for the UNLSE. Researchers in [28] used the extended simple equation technique to construct the exact solutions of unstable and modified unstable NLSEs. The modified exponential rational function technique was used by the researchers in [55] to obtain traveling wave results of UNLSE and mUNLSE. The extended mapping method was used by author [56] for solutions of mUNLSE. Therefore, we have accomplished many ingenious results that have not been stated before.
The suggested approach has produced a number of wave solutions in the mode of special Jacobi elliptic, trigonometric, and hyperbolic functions; all of them have important applications in applied sciences. Furthermore, with the assistance of the suitable parameter values, various kinds of bright-dark, multi-wave structures, breather and kink-type solitons, and several periodic solitary waves structures are obtained that aid in the comprehension understanding of nonlinear models. Figure 1, Figure 2, Figure 3 and Figure 4 illustrate the solitary waves solutions in different shapes of UNLSE. In Figure 1, the waves solution (19) in different architectures is drawn with the help of relevant parameters: (a) the kink-type soliton with its 2D in (d), combine dark-bright solitary waves in (b) and (c), and their 2D in (e) and (f). In Figure 2, the waves solution (22) in different architectures is drawn with the help of relevant parameters: (a) is multi peaks mixed waves with its 2D in (d), (b) is Breather-type waves of a strange structure with its 2D in (e), and (c) is the periodic soliton with its 2D in (f). In Figure 3, the waves solution (26) in different architectures is drawn with the help of relevant parameters: (a) is Bright solitons and its 2D in (d), (b) and (c) are multi-peakon solitons that have different amplitude with their 2D in (e) and (f). In Figure 4, the waves solution (48) in different architectures is drawn with the help of relevant parameters: (a) the periodic bright-dark solitary wave with its 2D in (d), (b), and (c) are Breather-type waves of strange shapes and their 2D in (e) and (f).
Figure 5, Figure 6, Figure 7 and Figure 8 illustrate the solitary wave solutions in various shapes of (mULSE). In Figure 5, the waves solution (36) in different architectures is drawn with the help of relevant parameters: (a) A Peakon soliton of amplitude 30 with its 2D in (d), (b) is a dark-bright solitary wave with its 2D in (e), and (c) is a dark solitary wave with its 2D in (f). Similarly, in Figure 6, the waves solution (39) in different architectures is drawn with the help of relevant parameters: (a) the multi-peak soliton with its 2D in (d), (b) is the Kink type solitary wave with its 2D in (e), and (c) is the dark soliton with its 2D in (f). In Figure 7, the waves solution (46) in different architectures is drawn with the help of relevant parameters: (a) is the Bright soliton with its 2D in (d), and (b) and (c) are multi peakon solitons with their 2D in (e) and (f). In Figure 8, the waves solution (48) in different architectures is drawn with the help of relevant parameters: (a), (b) and (c) are Kinky-breathers-type wave structures with their 2D in (d), (e), and (f), respectively.

6. Conclusions

The suggested improved auxiliary equation approach has been effectively used in this study to obtain innovative wave solutions of UNLSEs. These unstable nonlinear Schrödinger equations are universal equations of the nonlinear integrable system class that depict the temporal evolution of disruption in slightly stable or unstable mediums. Numerous varieties of results are generated in the form of the Jacobi elliptic function, and trigonometric and hyperbolic functions, many of which are distinctive and have significant applications in the fields of applied sciences and mathematical physics. Moreover, with the assistance of the suitable parameter values, various kinds of bright-dark, multi-wave structures, breather and kink-type solitons, and several periodic solitary waves are plotted, which help to describe the physical interpretation of unstable nonlinear equations models. The various constructed solutions demonstrate the utility of the suggested approach, which proves that the proposed technique can be utilized in other non-linear physical models occurring in mathematical physics.

Author Contributions

All authors have equal contributions. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to the Basque Government for its support through Grants IT1555-22 and KK-2022/00090, and to MCIN/AEI 269.10.13039/501100011033 for Grant PID2021-1235430B-C21/C22.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mitschke, F.; Mahnke, C.; Hause, A. Soliton Content of Fiber-Optic Light Pulses. Appl. Sci. 2017, 7, 635. [Google Scholar] [CrossRef]
  2. Chabchoub, A.; Kibler, B.; Finot, C.; Millot, G.; Onorato, M.; Dudley, J.M.; Babanin, A.V. The nonlinear Schrödinger equation and the propagation of weakly nonlinear waves in optical fibers and on the water surface. Ann. Phys. 2015, 361, 490–500. [Google Scholar] [CrossRef]
  3. Kumar, S.; Kumar, A.; Wazwaz, A.-M. New exact solitary wave solutions of the strain wave equation in microstructured solids via the generalized exponential rational function method. Eur. Phys. J. Plus 2020, 135, 870. [Google Scholar] [CrossRef]
  4. Arshad, M.; Seadawy, A.R.; Lu, D. Modulation stability and dispersive optical soliton solutions of higher order nonlinear Schrödinger equation and its applications in mono-mode optical fibers. Superlattices Microstruct. 2018, 113, 419–429. [Google Scholar] [CrossRef]
  5. Akbar, M.A.; Ali, N.; Wazwaz, A.M. Closed form traveling wave solutions of non-linear fractional evolution equations through the modified simple equation method. Therm. Sci. 2018, 22, 341352. [Google Scholar]
  6. Yin, H.M.; Pan, Q.; Chow, K.W. Doubly periodic solutions and breathers of the Hirota equation: Recurrence, cascading mechanism and spectral analysis. Nonlinear Dyn. 2022, 110, 3751–3768. [Google Scholar] [CrossRef]
  7. Yin, H.M.; Chow, K.W. Breathers, cascading instabilities and Fermi-Pasta-Ulam-Tsingou recurrence of the derivative nonlinear Schrödinger equation: Effects of ‘self-steepening’ nonlinearity. Phys. D 2021, 428, 133033. [Google Scholar] [CrossRef]
  8. Frassu, S.; Galván, R.R.; Viglialoro, G. Uniform in time L-estimates for an attraction-repulsion chemotaxis model with double saturation. Discrete Contin. Dyn. Syst. Ser. B 2023, 3, 1886–1904. [Google Scholar] [CrossRef]
  9. Frassu, S.; Li, T.; Viglialoro, G. Improvements and generalizations of results concerning attraction-repulsion chemotaxis models. Math. Methods Appl. Sci. 2022, 45, 11067–11078. [Google Scholar] [CrossRef]
  10. Dvornikov, M. Stable Langmuir solitons in plasma with diatomic ions. Nonlin. Process. Geophys. 2013, 20, 581–588. [Google Scholar] [CrossRef] [Green Version]
  11. He, J.S.; Charalampidis, E.G.; Kevrekidis, P.G.; Frantzeskakis, D.J. Rogue waves in nonlinear Schrödinger models with variable coefficients: Application to Bose–Einstein condensates. Phys. Lett. A 2014, 378, 577–583. [Google Scholar] [CrossRef]
  12. Bona, J.L.; Saut, J.-C. Dispersive Blow-Up II. Schrödinger-Type Equations, Optical and Oceanic Rogue Waves. Chin. Ann. Math. Ser. B 2010, 31, 793–818. [Google Scholar] [CrossRef]
  13. Grinevich, P.G.; Santini, P.M. The exact rogue wave recurrence in the NLS periodic setting via matched asymptotic expansions, for 1 and 2 unstable modes. Phys. Lett. A 2018, 382, 973–979. [Google Scholar] [CrossRef] [Green Version]
  14. Chabchoub, A.; Grimshaw, R.H.J. The Hydrodynamic Nonlinear Schrödinger Equation: Space and Time. Fluids 2016, 1, 23. [Google Scholar] [CrossRef] [Green Version]
  15. Kumar, S.; Kumar, D.; Kharbanda, H. Lie symmetry analysis, abundant exact solutions and dynamics of multisolitons to the (2 + 1)-dimensional KP-BBM equation. Pramana 2021, 95, 33. [Google Scholar] [CrossRef]
  16. Yin, H.M.; Pan, Q.; Chow, K.W. Four-wave mixing and coherently coupled Schrödinger equations: Cascading processes and Fermi–Pasta–Ulam–Tsingou recurrence. Chaos 2021, 31, 083117. [Google Scholar] [CrossRef]
  17. Solli, D.R.; Ropers, C.; Koonath, P.; Jalali, B. Optical rogue waves. Nature 2007, 450, 1054–1057. [Google Scholar] [CrossRef]
  18. Viotti, C.; Dutykh, D.; Dudley, J.M.; Dias, F. Emergence of coherent wave groups in deep-water random sea. Phys. Rev. E 2013, 87, 063001. [Google Scholar] [CrossRef] [Green Version]
  19. Kumar, S.; Kumar, D.; Kumar, A. Lie symmetry analysis for obtaining the abundant exact solutions, optimal system and dynamics of solitons for a higher-dimensional Fokas equation. Chaos Solitons Fractals 2021, 142, 110507. [Google Scholar] [CrossRef]
  20. Khatun, M.A.; Arefin, M.A.; Uddin, M.H.; Baleanu, D.; Akbar, M.A.; Inc, M. Explicit wave phenomena to the couple type fractional order nonlinear evolution equations. Results Phys. 2021, 28, 104597. [Google Scholar] [CrossRef]
  21. Kibler, B.; Fatome, J.; Finot, C.; Millot, G.; Dias, F.; Genty, G.; Akhmediev, N.; Dudley, J.M. The Peregrine soliton in nonlinear fibre optics. Nat. Phys. 2010, 6, 790–795. [Google Scholar] [CrossRef]
  22. Frisquet, B.; Kibler, B.; Millot, G. Collision of Akhmediev Breathers in Nonlinear Fiber Optics. Phys. Rev. X 2013, 3, 041032. [Google Scholar] [CrossRef] [Green Version]
  23. Chabchoub, A.; Hoffmann, N.; Onorato, M.; Akhmediev, N. Super Rogue Waves: Observation of a Higher-Order Breather in Water Waves. Phys. Rev. X 2012, 2, 011015. [Google Scholar] [CrossRef] [Green Version]
  24. Zakharov, V.E.; Ostrovsky, L.A. Modulation instability: The beginning. Phys. D 2009, 238, 540–548. [Google Scholar] [CrossRef]
  25. Hosseini, K.; Mirzazadeh, M.; Baleanu, D.; Salahshour, S.; Akinyemi, L. Optical solitons of a high-order nonlinear Schrödinger equation involving nonlinear dispersions and Kerr effect. Opt. Quantum Electron. 2022, 54, 177. [Google Scholar] [CrossRef]
  26. Habib, M.A.; Ali, H.M.S.; Miah, M.M.; Akbar, M.A. The generalized Kudryashov method for new closed form traveling wave solutions to some NLEEs. AIMS Math. 2019, 4, 896–909. [Google Scholar] [CrossRef]
  27. Arshad, M.; Seadawy, A.R.; Lu, D. Optical soliton solutions of the generalized higher-order nonlinear Schrödinger equations and their applications. Opt. Quantum Electron. 2017, 49, 421. [Google Scholar] [CrossRef]
  28. Lu, D.; Seadawy, A.R.; Arshad, M. Applications of extended simple equation method on unstable nonlinear Schrödinger equations. Optik 2017, 140, 136–144. [Google Scholar] [CrossRef]
  29. Elboree, M.K. The Jacobi elliptic function method and its application for two component BKP hieracy equations. Comput. Math. Appl. 2011, 62, 4402–4414. [Google Scholar] [CrossRef] [Green Version]
  30. Khan, K.; Akbar, M.A.; Islam, M. Exact solutions for (1 + 1)-dimensional nonlinear dispersive modified Benjamin-Bona-Mahony equation and coupled Klein-Gordon equations. SpringerPlus 2014, 3, 724. [Google Scholar] [CrossRef] [Green Version]
  31. Arshad, M.; Seadawy, A.R.; Lu, D. Modulation stability and optical soliton solutions of nonlinear Schrödinger equation with higher order dispersion and nonlinear terms and its applications. Superlattices Microstruct. 2017, 112, 422–434. [Google Scholar] [CrossRef]
  32. Biswa, A.; Arshed, S. Application of semi-inverse variational principle to cubic-quartic optical solitons with kerr and power law nonlinearity. Optik 2018, 172, 847–850. [Google Scholar] [CrossRef]
  33. Gaber, A.A.; Aljohani, A.F.; Ebaid, A.; Machado, J.T. The generalized Kudryashov method for nonlinear space–time fractional partial differential equations of Burgers type. Nonlinear Dyn. 2019, 95, 361–368. [Google Scholar] [CrossRef]
  34. Miah, M.M.; Ali, H.M.S.; Akbar, M.A.; Wazwaz, A.M. Some applications of the (G/G,1/G)-expansion method to find new exact solutions of NLEEs. Eur. Phys. J. Plus 2017, 132, 252. [Google Scholar] [CrossRef]
  35. Zhang, R.-F.; Li, M.-C. Bilinear residual network method for solving the exactly explicit solutions of nonlinear evolution equations. Nonlinear Dyn. 2022, 108, 521–531. [Google Scholar] [CrossRef]
  36. Seadawy, A.R.; Lu, D.; Khater, M.M.A. Bifurcations of traveling wave solutions for Dodd-Bullough-Mikhailov equation and coupled Higgs equation and their applications. Chin. J. Phys. 2017, 55, 1310–1318. [Google Scholar] [CrossRef]
  37. Yasar, E.; Yildirim, Y.; Yasar, E. New optical solitons of space-time conformable fractional perturbed Gerdjikov-Ivanov equation by sine-Gordon equation method. Results Phys. 2018, 9, 1666–1672. [Google Scholar] [CrossRef]
  38. Hosseini, M.M.; Ghaneai, H.; Mohyud-din, S.T.; Usman, M. Tri-prong scheme for regularized long wave equation. J. Assoc. Arab Univ. Basic Appl. Sci. 2016, 20, 68–77. [Google Scholar] [CrossRef]
  39. Yin, H.M.; Pan, Q.; Chow, K.W. The Fermi–Pasta–Ulam–Tsingou recurrence for discrete systems: Cascading mechanism and machine learning for the Ablowitz-Ladik equation. Commun. Nonlinear Sci. Numer. Simul. 2022, 114, 106664. [Google Scholar] [CrossRef]
  40. Yang, Y.; Suzuki, T.; Cheng, X. Darboux transformations and exact solutions for the integrable nonlocal Lakshmanan-Porsezian-Daniel equation. Appl. Math. Lett. 2020, 99, 105998. [Google Scholar] [CrossRef]
  41. Fang, J.; Wu, B.; Liu, W. An explicit spectral collocation method for the linearized Korteweg–de Vries equation on unbounded domain. Appl. Numer. Math. 2018, 126, 34–52. [Google Scholar] [CrossRef]
  42. Zeng, X.; Wang, D.S. A generalized extended rational expansion method and its application to (1 + 1)-dimensional dispersive long wave equation. Appl. Math. Comput. 2009, 212, 296–304. [Google Scholar] [CrossRef]
  43. Wang, Y.-H.; Chen, Y. Binary Bell Polynomials, Bilinear Approach to Exact Periodic Wave Solutions of (2 + 1)-Dimensional Nonlinear Evolution Equations. Commun. Theor. Phys. 2011, 56, 672–678. [Google Scholar] [CrossRef] [Green Version]
  44. Rizea, M. Exponential fitting method for the time-dependent Schrödinger equation. J. Math. Chem. 2010, 48, 55–65. [Google Scholar] [CrossRef]
  45. Zayed, E.M.E.; Alurrfi, K.A.E. New extended auxiliary equation method and its applications to nonlinear Schrödinger-type equations. Optik 2016, 127, 9131–9151. [Google Scholar] [CrossRef]
  46. Sarwar, A.; Gang, T.; Arshad, M.; Ahmed, I. Construction of brightdark solitary waves and elliptic function solutions of space-time fractional partial differential equations and their applications. Phys. Scr. 2020, 95, 045227. [Google Scholar] [CrossRef]
  47. Nasreen, N.; Lu, D.; Arshad, M. Optical soliton solutions of nonlinear Schrödinger equation with second order spatiotemporal dispersion and its modulation instability. Optik 2018, 161, 221–229. [Google Scholar] [CrossRef]
  48. Kumar, S.; Ma, W.-X.; Kumar, A. Lie symmetries, optimal system and group-invariant solutions of the (3 + 1)-dimensional generalized KP equation. Chin. J. Phys. 2021, 69, 1–23. [Google Scholar] [CrossRef]
  49. Hafiz Uddin, M.; Zaman, U.H.M.; Arefin, M.A.; Akbar, M.A. Nonlinear dispersive wave propagation pattern in optical fiber system. Chaos Solitons Fractals 2022, 164, 112596. [Google Scholar] [CrossRef]
  50. Pedlosky, J. Finite-amplitude baroclinic waves. J. Atmos. Sci. 1970, 27, 15–30. [Google Scholar] [CrossRef]
  51. Pawlik, M.; Rowlands, G. The propagation of solitary waves in piezoelectric semiconductors. J. Phys. C Solid State Phys. 1975, 8, 1189–1204. [Google Scholar] [CrossRef]
  52. Wadati, M.M.; Segur, H.; Ablowitz, M.J. A new Hamiltonian amplitude equation governing modulated wave instabilities. J. Phys. Soc. Jpn. 1992, 61, 1187–1193. [Google Scholar] [CrossRef]
  53. Arshad, M.; Seadawy, A.R.; Lu, D.; Jun, W. Optical soliton solutions of unstable nonlinear Schrödinger dynamical equation and stability analysis with applications. Optik 2018, 157, 597–605. [Google Scholar] [CrossRef]
  54. Arbabi, S.; Najafi, M. Exact solitary wave solutions of the complex nonlinear Schrödinger equations. Optik 2016, 127, 4682–4688. [Google Scholar] [CrossRef]
  55. Yue, L.; Lu, D.; Arshad, M.; Xu, X. New exact traveling wave solutions of the unstable nonlinear Schrödinger equations and their applications. Optik 2021, 226, 165386. [Google Scholar]
  56. Arshad, M.; Seadawy, A.R.; Lu, D.; Jun, W. Modulation instability analysis of modify unstable nonlinear Schrödinger dynamical equation and its optical soliton solutions. Results Phys. 2017, 7, 4153–4161. [Google Scholar] [CrossRef]
Figure 1. The solutions (19) are drawn in different shapes at c 0 = a 0 = 0.1 , a 2 = 1.5 , a 4 = 0.1 , κ = λ = μ = σ = τ = 1 : (a) Kink-type soliton with its 2D in (d), combining dark-bright solitary waves in (b) and (c), and their 2D in (e) and (f), respectively.
Figure 1. The solutions (19) are drawn in different shapes at c 0 = a 0 = 0.1 , a 2 = 1.5 , a 4 = 0.1 , κ = λ = μ = σ = τ = 1 : (a) Kink-type soliton with its 2D in (d), combining dark-bright solitary waves in (b) and (c), and their 2D in (e) and (f), respectively.
Symmetry 15 00099 g001
Figure 2. The solutions (22) are drawn in different shapes at: (a) is multi peaks mixed waves with its 2D in (d), (b) is Breather-type waves of strange structure with its 2D in (e), and (c) is periodic soliton with its 2D in (f), respectively, at c 0 = a 0 = 0.1 , a 2 = 1.25 , a 4 = 0.1 , κ = λ = μ = σ = τ = 1 .
Figure 2. The solutions (22) are drawn in different shapes at: (a) is multi peaks mixed waves with its 2D in (d), (b) is Breather-type waves of strange structure with its 2D in (e), and (c) is periodic soliton with its 2D in (f), respectively, at c 0 = a 0 = 0.1 , a 2 = 1.25 , a 4 = 0.1 , κ = λ = μ = σ = τ = 1 .
Symmetry 15 00099 g002
Figure 3. The solution (26) in different shapes are drawn at γ = μ = 0.1 , σ = 1.5 , κ = λ = ν = 1 = τ = 1 : (a) is Bright solitons, and its 2D in (bd) are multi-peakon solitons that have different amplitudes with their 2D in (e) and (f), respectively.
Figure 3. The solution (26) in different shapes are drawn at γ = μ = 0.1 , σ = 1.5 , κ = λ = ν = 1 = τ = 1 : (a) is Bright solitons, and its 2D in (bd) are multi-peakon solitons that have different amplitudes with their 2D in (e) and (f), respectively.
Symmetry 15 00099 g003
Figure 4. The solution (28) in different shapes is drawn at γ = μ = 0.1 , σ = 0.1 , κ = λ = ν = 1 = τ = 1 : (a) Periodic bright-dark solitary wave with its 2D in (bd) are Breather-type waves of strange shapes and their 2D in (e) and (f), respectively.
Figure 4. The solution (28) in different shapes is drawn at γ = μ = 0.1 , σ = 0.1 , κ = λ = ν = 1 = τ = 1 : (a) Periodic bright-dark solitary wave with its 2D in (bd) are Breather-type waves of strange shapes and their 2D in (e) and (f), respectively.
Symmetry 15 00099 g004
Figure 5. The solution (36) in different shapes is drawn at c 2 = 0.5 , a 0 = a 4 = 0.1 , a 2 = 0.35 , ρ = ω = 0.1 , τ = κ = λ = 1 : (a) Peakon soliton with its 2D in (d), (b) is dark-bright solitary wave with its 2D in (e), and (c) is dark solitary wave with its 2D in (f), respectively.
Figure 5. The solution (36) in different shapes is drawn at c 2 = 0.5 , a 0 = a 4 = 0.1 , a 2 = 0.35 , ρ = ω = 0.1 , τ = κ = λ = 1 : (a) Peakon soliton with its 2D in (d), (b) is dark-bright solitary wave with its 2D in (e), and (c) is dark solitary wave with its 2D in (f), respectively.
Symmetry 15 00099 g005
Figure 6. The solution (39) in different shapes are drawn at c 2 = 1.75 , a 0 = a 4 = 0.1 , a 2 = τ = 1.5 , ρ = ω = 0.1 , κ = λ = 1 : (a) Multi-peak dark soliton with its 2D in (d), (b) Kink type solitary wave with its 2D in (e), and (c) dark soliton with its 2D in (f), respectively.
Figure 6. The solution (39) in different shapes are drawn at c 2 = 1.75 , a 0 = a 4 = 0.1 , a 2 = τ = 1.5 , ρ = ω = 0.1 , κ = λ = 1 : (a) Multi-peak dark soliton with its 2D in (d), (b) Kink type solitary wave with its 2D in (e), and (c) dark soliton with its 2D in (f), respectively.
Symmetry 15 00099 g006
Figure 7. The solution (46) in different shapes is drawn at c 1 = 1.5 , a 2 = 1.25 , ρ = 2.5 , γ = ω = ν = τ = 1 , κ = 2.5 , λ = 1 : (a) is Bright soliton with its 2D in (bd) are multi-peakon solitons with their 2D in (e) and (f), respectively.
Figure 7. The solution (46) in different shapes is drawn at c 1 = 1.5 , a 2 = 1.25 , ρ = 2.5 , γ = ω = ν = τ = 1 , κ = 2.5 , λ = 1 : (a) is Bright soliton with its 2D in (bd) are multi-peakon solitons with their 2D in (e) and (f), respectively.
Symmetry 15 00099 g007
Figure 8. The solution (48) in different shapes is drawn at ρ = 1.5 , γ = ω = ν = τ = 1 , κ = 2 , λ = 1 : (ac) are Kinky-breathers-type wave structures with their 2D in (df), respectively.
Figure 8. The solution (48) in different shapes is drawn at ρ = 1.5 , γ = ω = ν = τ = 1 , κ = 2 , λ = 1 : (ac) are Kinky-breathers-type wave structures with their 2D in (df), respectively.
Symmetry 15 00099 g008
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

Sarwar, A.; Arshad, M.; Farman, M.; Akgül, A.; Ahmed, I.; Bayram, M.; Rezapour, S.; De la Sen, M. Construction of Novel Bright-Dark Solitons and Breather Waves of Unstable Nonlinear Schrödinger Equations with Applications. Symmetry 2023, 15, 99. https://doi.org/10.3390/sym15010099

AMA Style

Sarwar A, Arshad M, Farman M, Akgül A, Ahmed I, Bayram M, Rezapour S, De la Sen M. Construction of Novel Bright-Dark Solitons and Breather Waves of Unstable Nonlinear Schrödinger Equations with Applications. Symmetry. 2023; 15(1):99. https://doi.org/10.3390/sym15010099

Chicago/Turabian Style

Sarwar, Ambreen, Muhammad Arshad, Muhammad Farman, Ali Akgül, Iftikhar Ahmed, Mustafa Bayram, Shahram Rezapour, and Manuel De la Sen. 2023. "Construction of Novel Bright-Dark Solitons and Breather Waves of Unstable Nonlinear Schrödinger Equations with Applications" Symmetry 15, no. 1: 99. https://doi.org/10.3390/sym15010099

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