7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Correspondence
Corrigendum
Editorial
Full Length Article
Invited review
Letter to the Editor
Original Article
Retraction notice
REVIEW
Review Article
SHORT COMMUNICATION
Short review
7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Correspondence
Corrigendum
Editorial
Full Length Article
Invited review
Letter to the Editor
Original Article
Retraction notice
REVIEW
Review Article
SHORT COMMUNICATION
Short review
View/Download PDF

Translate this page into:

Original article
01 2020
:33;
101180
doi:
10.1016/j.jksus.2020.09.007

New exact solitary wave solutions, bifurcation analysis and first order conserved quantities of resonance nonlinear Schrödinger’s equation with Kerr law nonlinearity

Department of Mathematics, Namal Institute, 30KM Talagang Road, Mianwali 42250, Pakistan
Department of Mathematics and Science Education, Faculty of Education, Harran University, Sanliurfa, Turkey
Faculty of Educational Sciences, Final International University, Kyrenia, Mersin 10, Turkey
School of Information Science and Technology, Yunnan Normal University, Yunnan, China

⁎Corresponding author. gaowei@ynnu.edu.cn (Wei Gao)

Disclaimer:
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.

Peer review under responsibility of King Saud University.

Abstract

This paper anatomizes the exact solutions of the resonant non-linear Schrödinger’s equation (R-NLSE) with the Kerr law non-linearity with the assistance of the new extended direct algebraic technique. The secured soliton erections are newfangled and unreservedly invigorating for investigators. The graphically comprehensive report of some specific solutions is embellished with the well-judged values of parameters to illustrate their propagation. Then a planer dynamical system is introduced and the bifurcation analysis has been executed to figure out the bifurcation structures of the non-linear and super non-linear traveling wave solutions of the heeded model. All possible phase portraits are exhibited with specific values of parameters. Furthermore, a precise class of non-trivial and first-order conserved quantities is enumerated by the intervention of the multiplier approach.

Keywords

Schrödinger’s equation
Bifurcation theory
Conservation laws
1

1 Introduction

The non-linear Schrödinger’s equation (NLSE) has an imperative influence and pivotal significance in different branches of science. Quantum mechanics, non-linear optics, fluid dynamics, and plasma physics are some of the fields where NLSE has appeared. The R-NLSE has availed oneself of the study of Madelung fluids and solitonic dynamics in various non-linear processes (Baleanua et al., 2017 and references therein):

(1)
i Φ t + α Φ xx + β Ω ( | Φ | 2 ) Φ + γ | Φ xx | Φ | Φ = 0 , i = - 1 .

The independent variables x and t depicts the non-dimensional distance in Eq. (1), along the fiber and temporal variables, respectively. The complex valued dependent variable Φ ( t , x ) represents wave profile, while, the constants parameters α , β and γ are the coefficients of group-velocity dispersion, non-Kerr non-linearity and resonant non-linearity respectively. In the Eq. (1), real valued and n-times differentiable function Ω is defined as: Ω ( | Φ | 2 ) Φ : C C . where the two-dimensional linear space of R 2 represents by Argand plane C . Here we are assuming Ω ( y ) = y , which arises in water waves and the non-linear fiber optics when the refractive index is directly proportional to intensity (Biswas and Konat, 2006).

The R-NLSE equation has been discussed by many means in literature. Over the last few years, scientists have been quantified the exact solutions of Eq. (1) by utilizing numerous distinct approaches. The Jacobi elliptic tool has been exploited to quantifying the exact structures of Eq. (1) and also simplest equation approach took under consideration (Eslami et al., 2013). To extract the exact solitary wave solutions of Eq. (1), D. Baleanu has been prosecuted the Ricatti-Bernoulli sub-ODE technique (Baleanua et al., 2017). One of the aspirations of this exploration is to deliberate a futuristic class of exact solutions. On that account, a new direct extended algebraic approach is maneuvered to come across the exact solutions of Eq. (1). Some of the classical contributions on numerical aspects of the PDEs are reported in Cattani et al. (2013), Cattani et al. (2012), Heydari et al. (2015), Kumar et al. (2018), Kumar et al. (2020), Rushchitsky et al. (2004), Singh and Srivastava (2020), Singh (xxxxx), Singh (2020), Singh et al. (2019), Singh et al. (2020) and Yadav et al. (2019). In recent years, the study of differential equations employing bifurcation analysis is a hot topic of research. Dubinov et. al. in Dubinov et al. (2012) characterized a new class of nonlinear and super-nonlinear waves. In Zhang et al. (2013), a class of solutions is obtained for Klein–Gordon–Zakharov equations by using bifurcation theory. More recently, Sharma-Tasso-Olver equation is dealt with utilizing bifurcation theory (Ali et al., 2018) and the complete classification of waves is presented. By the best of our knowledge, any study related to the behavior of nonlinear and super-nonlinear traveling waves for R-NLSE equation is not done before. Consequently, a profundity anatomy of Eq. (1) along with these trajectories is fascinating and is propounded here.

Conservation laws have miraculous contributions to unfold the partial differential equations and in various applications. The suggestion about conservation laws ejaculated by the perception of physical laws like mass, energy, and momentum. Initially, The ideology and exploitation of the Noether’s theorem was authenticated by German mathematician Emmy Noether and also certified as a well-ordered approach for uncovering the conservation laws Bessel-Hagen, 1921. Noether theorem states that each Euler - Lagrange equation’s Noether symmetry corresponds to a difference equation’s conservation law. Noether’s theorem works only for a variational differential equation, yet there are differential equations which have no Lagrangian equations which can be dealt with different approaches available in the literature, some of them are (Anco and Bluman, 2002; Ibragimov, 2007; Kara and Mahomed, 2006) while computer package for construction of conserved quantities by using (Anco and Bluman, 2002) is also developed (Cheviakov, 2007) and utilized in this research. Here we search out the first order nontrivial conserved quantities of Eq. (1) by using a multiplier approach (Anco and Bluman, 2002). Some latest works related to exact solutions and conservation laws are given in Baskonus et al. (2019), Eskitaolu et al. (2019), Khalique et al. (2018), Khalique and Mhlanga (2018) and Moleleki et al. (2018).

This report is designated as segment 2 is customized for the exact solutions of Eq. (1) by using the new direct extended algebraic technique. Bifurcation analysis of Eq. (1) is presented in Section 3. Section 4 is devoted to conserved quantities while concluding remarks are asserted at the termination.

2

2 Travelling wave solutions

2.1

2.1 Specification of the lodged approach

Suppose that the non-linear partial differential equation:

(2)
Q ( Φ , Φ t , Φ x , Φ y , Φ tt , Φ xx , ) = 0 , remoulded into non-linear ordinary differential equation:
(3)
G ( u , u , u , ) = 0 ,
operating the complex transformation
(4)
Φ ( t , x ) = u ( ξ ) e iv ( t , x ) ,
where, ξ = K ( x + st ) and v ( t , x ) = - kx + ω t + θ and here prime in Eq. (3) demonstrates the notation of differentiation concerning to ξ .

Assume Eq. (3) has solution of the such pattern:

(5)
u ( ξ ) = a 0 + i = 1 m a i W ( ξ ) , where,
(6)
W ( ξ ) = ln ( ρ ) μ + ν W ( ξ ) + ζ W 2 ( ξ ) , ρ 0 , 1 ,
where μ , ν and ζ are the real constants.

Eq. (6) has general solutions in respect of parameters μ , ν and ζ are prescribed as Rezazadeh, 2018:

  1. When ν 2 - 4 μ ζ < 0 with ζ 0 ,

    (7)
    W 1 ( ξ ) = - ν 2 ζ + - ( ν 2 - 4 μ ζ ) 2 ζ tan ρ - ( ν 2 - 4 μ ζ ) 2 ξ ,
    (8)
    W 2 ( ξ ) = - ν 2 ζ - - ( ν 2 - 4 μ ζ ) 2 ζ cot ρ - ( ν 2 - 4 μ ζ ) 2 ξ ,
    (9)
    W 3 ( ξ ) = - ν 2 ζ + - ( ν 2 - 4 μ ζ ) 2 ζ tan ρ - ν 2 - 4 μ ζ ξ ± mn sec ρ - ν 2 - 4 μ ζ ξ ,
    (10)
    W 4 ( ξ ) = - ν 2 ζ + - ( ν 2 - 4 μ ζ ) 2 ζ cot ρ - ν 2 - 4 μ ζ ξ ± mn csc ρ - ν 2 - 4 μ ζ ξ ,
    (11)
    W 5 ( ξ ) = - ν 2 ζ + - ( ν 2 - 4 μ ζ ) 4 ζ tan ρ - ( ν 2 - 4 μ ζ ) 4 ξ - cot ρ - ( ν 2 - 4 μ ζ ) 4 ξ .

  2. When ν 2 - 4 μ ζ > 0 with ζ 0 ,

    (12)
    W 6 ( ξ ) = - ν 2 ζ - ν 2 - 4 μ ζ 2 ζ tanh ρ ν 2 - 4 μ ζ 2 ξ ,
    (13)
    W 7 ( ξ ) = - ν 2 ζ - ν 2 - 4 μ ζ 2 ζ coth ρ ν 2 - 4 μ ζ 2 ξ ,
    (14)
    W 8 ( ξ ) = - ν 2 ζ + ν 2 - 4 μ ζ 2 ζ - tanh ρ ν 2 - 4 μ ζ ξ ± i mn sech ρ ν 2 - 4 μ ζ ξ ,
    (15)
    W 9 ( ξ ) = - ν 2 ζ + ν 2 - 4 μ ζ 2 ζ - coth ρ ν 2 - 4 μ ζ ξ ± mn csch ρ ν 2 - 4 μ ζ ξ ,
    (16)
    W 10 ( ξ ) = - ν 2 ζ - ν 2 - 4 μ ζ 4 ζ tanh ρ ν 2 - 4 μ ζ 4 ξ + coth ρ ν 2 - 4 μ ζ 4 ξ .

  3. When μ ζ > 0 with ν = 0 ,

    (17)
    W 11 ( ξ ) = μ ζ tan ρ μ ζ ξ ,
    (18)
    W 12 ( ξ ) = - μ ζ cot ρ μ ζ ξ ,
    (19)
    W 13 ( ξ ) = μ ζ tan ρ 2 μ ζ ξ ± mn sec ρ 2 μ ζ ξ ,
    (20)
    W 14 ( ξ ) = μ ζ - cot ρ 2 μ ζ ξ ± mn csc ρ 2 μ ζ ξ ,
    (21)
    W 15 ( ξ ) = 1 2 μ ζ tan ρ μ ζ 2 ξ - cot ρ μ ζ 2 ξ .

  4. When μ ζ < 0 with ν = 0 ,

    (22)
    W 16 ( ξ ) = - - μ ζ tanh ρ - μ ζ ξ ,
    (23)
    W 17 ( ξ ) = - - μ ζ coth ρ - μ ζ ξ ,
    (24)
    W 18 ( ξ ) = - μ ζ - tanh ρ 2 - μ ζ ξ ± i mn sech ρ 2 - μ ζ ξ ,
    (25)
    W 19 ( ξ ) = - μ ζ - coth ρ 2 - μ ζ ξ ± mn csch ρ 2 - μ ζ ξ ,
    (26)
    W 20 ( ξ ) = - 1 2 - μ ζ tanh ρ - μ ζ 2 ξ + coth ρ - μ ζ 2 ξ .

  5. When ν = 0 with μ = ζ ,

    (27)
    W 21 ( ξ ) = tan ρ μ ξ ,
    (28)
    W 22 ( ξ ) = - cot ρ μ ξ ,
    (29)
    W 23 ( ξ ) = tan ρ 2 μ ξ ± mn sec ρ 2 μ ξ ,
    (30)
    W 24 ( ξ ) = - cot ρ 2 μ ξ ± mn csc ρ 2 μ ξ ,
    (31)
    W 25 ( ξ ) = 1 2 tan ρ μ 2 ξ - cot ρ μ 2 ξ .

  6. When ν = 0 with ζ = - μ ,

    (32)
    W 26 ( ξ ) = - tanh ρ μ ξ ,
    (33)
    W 27 ( ξ ) = - coth ρ μ ξ ,
    (34)
    W 28 ( ξ ) = - tanh ρ 2 μ ξ ± i mn sech ρ 2 μ ξ ,
    (35)
    W 29 ( ξ ) = - cot ρ 2 μ ξ ± mn csch ρ 2 μ ξ ,
    (36)
    W 30 ( ξ ) = - 1 2 tanh ρ μ 2 ξ + cot ρ μ 2 ξ .

  7. When ν 2 = 4 μ ζ ,

    (37)
    W 31 ( ξ ) = - 2 μ ( ν ξ ln ρ + 2 ) ν 2 ξ ln ρ .

  8. When ν = p and μ = pq , ( q 0 ) but ζ = 0 ,

    (38)
    W 32 ( ξ ) = ρ p ξ - q .

  9. When ν = ζ = 0 ,

    (39)
    W 33 ( ξ ) = μ ξ ln ρ .

  10. When ν = 0 and μ = 0 ,

    (40)
    W 34 ( ξ ) = - 1 ζ ξ ln ρ .

  11. When μ = 0 but ν 0 ,

    (41)
    W 35 ( ξ ) - m ν ζ cosh ρ ν ξ - sinh ρ ν ξ + m ,
    (42)
    W 36 ( ξ ) = - ν sinh ρ ν ξ + cosh ρ ν ξ ζ sinh ρ ν ξ + cosh ρ ν ξ + n .

  12. When ν = p and ζ = pq , ( q 0 but μ = 0 ) ,

    (43)
    W 37 ( ξ ) = - m ρ p ξ m - qn ρ p ξ .

sinh ρ ( ξ ) = m ρ ξ - n ρ - ξ 2 , cosh ρ ( ξ ) = m ρ ξ + n ρ - ξ 2 , tanh ρ ( ξ ) = m ρ ξ - n ρ - ξ m ρ ξ + n ρ - ξ , coth ρ ( ξ ) = m ρ ξ + n ρ - ξ m ρ ξ - n ρ - ξ , sech ρ ( ξ ) = 2 m ρ ξ + n ρ - ξ , csch ρ ( ξ ) = 2 m ρ ξ - n ρ - ξ , sin ρ ( ξ ) = m ρ i ξ - n ρ - i ξ 2 i , cos ρ ( ξ ) = m ρ i ξ + n ρ - i ξ 2 , tan ρ ( ξ ) = - i m ρ i ξ - n ρ - i ξ m ρ i ξ + n ρ - i ξ , cot ρ ( ξ ) = i m ρ i ξ + n ρ - i ξ m ρ i ξ - n ρ - i ξ , sec ρ ( ξ ) = 2 m ρ ξ + n ρ - ξ , csc ρ ( ξ ) = 2 i m ρ ξ - n ρ - ξ , where, m and n are deformation parameters, which are arbitrary positive constants.

2.2

2.2 Applications to Eq. (1)

The R-NLSE with the Kerr law nonlinearity is considered for estimation of exact solutions via method purposed in Rezazadeh (2018). For this let us substitute a complex envelope (4) into Eq. (1) and splitting into the real and imaginary portion, respectively. We will attain:

(44)
K 2 ( γ + α ) u - ( k 2 α + ω ) u + β u 3 = 0 , s = 2 k α .

After stabilizing the highest order derivative expressions with contact to the highest power of non-linear expressions in Eq. (44), one can take the solution of such type:

(45)
u = a 0 + a 1 W ( ξ ) , where W ( ξ ) satisfies Eq. (6). On exchange Eq. (45) into Eq. (44) and equating the coefficients of non-identical powers of W ( ξ ) , leads to a system of an algebraic equations. W ( ξ ) 0 : - k 2 α a 0 - ω a 0 + β a 0 3 + K 2 μ ν α a 1 log ( ρ ) 2 + K 2 μ ν γ a 1 log ( ρ ) 2 = 0 , W ( ξ ) 1 : - k 2 α a 1 - ω a 1 + 3 β a 0 2 a 1 + K 2 ν 2 α a 1 log ( ρ ) 2 + 2 K 2 μ ζ α a 1 log ( ρ ) 2 + K 2 ν 2 γ a 1 log ( ρ ) 2 + 2 K 2 μ ζ γ a 1 log ( ρ ) 2 = 0 , W ( ξ ) 2 : 3 β a 0 a 1 2 + 3 K 2 ν ζ α a 1 log ( ρ ) 2 + 3 K 2 ν ζ γ a 1 log ( ρ ) 2 = 0 , W ( ξ ) 3 : β a 1 3 + 2 K 2 ζ 2 α a 1 log ( ρ ) 2 + 2 K 2 ζ 2 γ a 1 log ( ρ ) 2 = 0 .

To extricate the solution of an algebraic equations by means of Mathematica, following set of solution is achieved:

(46)
a 0 = ± ν Θ β Π , a 1 = ± 2 ζ Θ β Π , K = ± - 2 Θ Π ( α + γ ) log ( ρ ) 2 , where Θ = k 2 α + ω , Π = ν 2 - 4 μ ζ

Case 1. If Π < 0 with ζ 0 , then.

By mentioning the values of a 0 and a 1 from (46) into Eq. (45) we get: u 1 ± ( ξ ) = - Θ β tan ρ - Π 2 ξ , where, along with complex transformation (4), u 1 ± ( ξ ) capitulates as solutions of Eq. (1): Φ 1 ± ( t , x ) = ± - Θ β tan ρ - Π 2 ξ e i ( - kx + ω t + θ ) .

Thus working on the same line following solutions are obtained. Φ 2 ± ( t , x ) = ± - Θ β cot ρ - Π 2 ξ e i ( - kx + ω t + θ ) . Φ 3 ± ( t , x ) = ± - Θ β tan ρ - Π ξ ± mn sec ρ - Π ξ e i ( - kx + ω t + θ ) . Φ 4 ± ( t , x ) = ± - Θ β cot ρ - Π ξ ± mn csc ρ - Π ξ e i ( - kx + ω t + θ ) . Φ 5 ± ( t , x ) = ± - Θ β tan ρ - Π 4 ξ - cot ρ - Π 4 ξ e i ( - kx + ω t + θ ) .

Case 2. If Π > 0 with ζ 0 , then Φ 6 ± ( t , x ) = Θ β tanh ρ Π 2 ξ e i ( - kx + ω t + θ ) . Φ 7 ± ( t , x ) = Θ β coth ρ Π 2 ξ e i ( - kx + ω t + θ ) . Φ 8 ± ( t , x ) = Θ β - tanh ρ Π ξ ± i mn sech ρ Π ξ e i ( - kx + ω t + θ ) . Φ 9 ± ( t , x ) = Θ β - coth ρ Π ξ ± mn csch ρ Π ξ e i ( - kx + ω t + θ ) . Φ 10 ± ( t , x ) = Θ β tanh ρ Π 4 ξ + coth ρ Π 4 ξ e i ( - kx + ω t + θ ) .

Case 3. If μ ζ > 0 with ν = 0 , then Φ 11 ( t , x ) = ± Θ β Π ν + 2 μ ζ tan ρ μ ζ ξ e i ( - kx + ω t + θ ) . Φ 12 ( t , x ) = ± Θ β Π ν - 2 μ ζ cot ρ μ ζ ξ e i ( - kx + ω t + θ ) . Φ 13 ( t , x ) = ± Θ β Π ν + 2 μ ζ tan ρ 2 μ ζ ξ ± mn sec ρ 2 μ ζ ξ e i ( - kx + ω t + θ ) . Φ 14 ( t , x ) = ± Θ β Π ν + 2 μ ζ - cot ρ 2 μ ζ ξ ± mn csc ρ 2 μ ζ ξ e i ( - kx + ω t + θ ) . Φ 15 ( t , x ) = ± Θ 2 β Π ν + 2 μ ζ tan ρ μ ζ 2 ξ - cot ρ μ ζ 2 ξ e i ( - kx + ω t + θ ) .

Case 4. If μ ζ < 0 with ν = 0 , then Φ 16 ( t , x ) = ± Θ β Π ν - 2 - μ ζ tanh ρ - μ ζ ξ e i ( - kx + ω t + θ ) . Φ 17 ( t , x ) = ± Θ β Π ν - 2 - μ ζ coth ρ - μ ζ ξ e i ( - kx + ω t + θ ) . Φ 18 ( t , x ) = ± Θ β Π ν + 2 μ ζ - tanh ρ 2 - μ ζ ξ ± i mn sech ρ 2 - μ ζ ξ e i ( - kx + ω t + θ ) . Φ 19 ( t , x ) = ± Θ β Π ν + 2 μ ζ - coth ρ 2 - μ ζ ξ ± mn csch ρ 2 - μ ζ ξ e i ( - kx + ω t + θ ) . Φ 20 ( t , x ) = ± Θ 2 β Π ν - 2 - μ ζ tanh ρ - μ ζ 2 ξ + coth ρ - μ ζ 2 ξ e i ( - kx + ω t + θ ) .

Case 5. If When ν = 0 with μ = ζ , then Φ 21 ( t , x ) = ± Θ β Π ν + 2 ζ tan ρ μ ξ e i ( - kx + ω t + θ ) . Φ 22 ( t , x ) = ± Θ β Π ν - 2 ζ cot ρ μ ξ e i ( - kx + ω t + θ ) . Φ 23 ( t , x ) = ± Θ β Π ν + 2 ζ tan ρ 2 μ ξ ± mn sec ρ 2 μ ξ e i ( - kx + ω t + θ ) . Φ 24 ( t , x ) = ± Θ β Π ν + 2 ζ - cot ρ 2 μ ξ ± mn csc ρ 2 μ ξ e i ( - kx + ω t + θ ) . Φ 25 ( t , x ) = ± Θ 2 β Π ν + 2 ζ tan ρ μ 2 ξ - cot ρ μ 2 ξ e i ( - kx + ω t + θ ) .

Case 6. If ν = 0 with ζ = - μ , then Φ 26 ( t , x ) = ± Θ 2 β Π ν - 2 ζ tanh ρ μ ξ e i ( - kx + ω t + θ ) . Φ 27 ( t , x ) = ± Θ 2 β Π ν - 2 ζ coth ρ μ ξ , Φ 28 ( t , x ) = ± Θ 2 β Π ν + 2 ζ - tanh ρ 2 μ ξ ± i mn sech ρ 2 μ ξ e i ( - kx + ω t + θ ) . Φ 29 ( t , x ) = ± Θ 2 β Π ν + 2 ζ - coth ρ 2 μ ξ ± mn csch ρ 2 μ ξ e i ( - kx + ω t + θ ) . Φ 30 ( t , x ) = Θ 4 β Π ν + 2 ζ tanh ρ μ 2 ξ + coth ρ μ 2 ξ e i ( - kx + ω t + θ ) .

Case 7. If ν 2 = 4 μ ζ , then Φ 31 ( t , x ) = ± Θ 2 β Π ν + 2 ζ - 2 μ ( ν ξ ln ρ + 2 ) ν 2 ξ ln ρ e i ( - kx + ω t + θ ) .

Case 8. If ν = p and μ = pq , ( q 0 ) but ζ = 0 , then Φ 32 ( t , x ) = ± Θ 2 β Π ν + 2 ζ ρ p ξ - q e i ( - kx + ω t + θ ) .

Case 9. If ν = ζ = 0 , then Φ 33 ( t , x ) = ± Θ 2 β Π ν + 2 ζ μ ξ ln ρ e i ( - kx + ω t + θ ) .

Case 10. If ν = 0 and also μ = 0 , then Φ 34 ( t , x ) = ± Θ 2 β Π ν + 2 ζ - 1 ζ ξ ln ρ e i ( - kx + ω t + θ ) .

Case 11. If μ = 0 but ν 0 then Φ 35 ( t , x ) = ± Θ 2 β Π ν + 2 ζ - m ν ζ cosh ρ ν ξ - sinh ρ ν ξ + m e i ( - kx + ω t + θ ) . Φ 36 ( t , x ) = ± Θ 2 β Π ν + 2 ζ - ν sinh ρ ν ξ + cosh ρ ν ξ ζ sinh ρ ν ξ + cosh ρ ν ξ + n e i ( - kx + ω t + θ ) .

Case 12. If ν = p and ζ = pq , ( q 0 but μ = 0 ) then Φ 37 ( t , x ) = ± Θ 2 β Π ν + 2 ζ - m ρ p ξ m - qn ρ p ξ e i ( - kx + ω t + θ ) .

3D-graphics, 2D-graphics and contour plots of different solutions | Φ i | of Eq. (1) with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 ands = 0.2 are presented to describe their behaviour in Fig. 1–6.

Different graphical representations of ∣ Φ 1 ∣ with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Fig. 1
Different graphical representations of Φ 1 with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Different graphical representations of ∣ Φ 8 ∣ with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Fig. 2
Different graphical representations of Φ 8 with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Different graphical representations of ∣ Φ 9 ∣ with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Fig. 3
Different graphical representations of Φ 9 with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Different graphical representations of ∣ Φ 25 ∣ with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Fig. 4
Different graphical representations of Φ 25 with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Different graphical representations of ∣ Φ 32 ∣ with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Fig. 5
Different graphical representations of Φ 32 with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Different graphical representations of ∣ Φ 33 ∣ with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .
Fig. 6
Different graphical representations of Φ 33 with n = 3 , m = 2 , ρ = 3 , Π = 1 , α = 0.5 , γ = 0.5 , Θ = 1 , β = 1 and s = 0.2 .

3

3 Bifurcations behavior and phase portraits

According to ideology of the planar dynamical systems, equilibrium point ( u q , z q ) is exclaimed as the saddle point if J < 0 , if J > 0 , then center and T 1 = 0 , a node if J > 0 and T 1 2 - 4 J > 0 while, when J = 0 then zero point and Poincaré index of ( u q , z q ) is zero. Where, J and T 1 exhibits the trace of coefficients matrix and also Jacobian matrix for any linearized system of (47). For the categorization of non-identical orbits in phase portraits of dynamical system (47), some symbols will be practised:

  • (1):

    Super non-linear periodic orbit is manifested by SNPO(e,s),

  • (2):

    Non-linear homoclinic orbit is manifested by NHO(e,s),

  • (3):

    Non-linear heteroclinic orbit is manifested by NHTO(e,s),

  • (4):

    Non-linear periodic orbit is manifested by NPO(e,s),

where ‘e’ exemplifies the equilibrium points and ‘s’ exhibits the seperatix layers enveloped by an orbit. Every one phase orbit is the closed non-self-intersecting curve on the phase plane. Phase portrait of dynamical systems is a certain class of such nested phase directions. Eq. (44) can be demonstrated as a system of non-linear dynamical equations:
(47)
du d ξ = z , dz d ξ = ( k 2 α + ω ) u K 2 ( γ + α ) + β u 3 K 2 ( γ + α ) .

The system (47) delineates the planar Hamiltonian kind. Hamiltonian functions prevailed by integrating (47):

(48)
H ( u , z ) = z 2 2 - ( k 2 α + ω ) u 2 2 K 2 ( γ + α ) - β u 4 4 K 2 ( γ + α ) = h .

From (48), validity can be attested as:

(49)
du d ξ = H z and dz d ξ = - H u .

As a system (47) is a planar Hamiltonian system and from (49), it can be concluded that system (47) is conservative and thus phase orbits expressed by the vector field of (47) will posses each traveling waves solution of Eq. (44) (for detail see Li et al., 2015 and references therein).

Level curves L h ( u , z ) in respect of energy level h can be defined in the following fashion: L h = { ( u , z ) R × R : H ( u , z ) = h } , where, H ( u , z ) is expounded in (48) and h is known as the energy level. In phase portraits against every energy level, h one can have an orbit. To investigate the relations between closed orbit and the energy level of the system (47), let us define:

(50)
E h ( u ) = h + ( k 2 α + ω ) u 2 2 K 2 ( γ + α ) + β u 4 4 K 2 ( γ + α ) .

From (48), one can easily find the following relation:

(51)
z = ± 2 h + ( k 2 α + ω ) u 2 K 2 ( γ + α ) + β u 4 2 K 2 ( γ + α ) , which means z 2 2 = E h ( u ) . Graphical illustration of (50) is given in Fig. 9(a-b). Here, all feasible phase trajectories for dynamical system (47) are presented and categorized.

System (47) has three equilibrium points: u 1 = ( 0 , 0 ) , u 2 = ( ( k 2 α + ω ) β , 0 ) , u 3 = ( - ( k 2 α + ω ) β , 0 ) .

The linearized system (47) reported as coefficient matrix at an equilibrium point ( u q , z q ) :

(52)
M = 0 1 ( k 2 α + ω ) K 2 ( γ + α ) 0 , while Jacobian for the system (47) is:
(53)
J = 0 1 ( k 2 α + ω ) K 2 ( γ + α ) + 3 β u 2 K 2 ( γ + α ) 0 .

It yields the following cases:

3.0.1

3.0.1 K 2 ( γ + α ) > 0 , ( - k 2 α - ω ) < 0 , C > 0 or K 2 ( γ + α ) < 0 , ( - k 2 α - ω ) > 0 , β < 0

There are three equilibrium points of system (47) u 1 , u 2 , and u 3 for this type. For this J ( u 1 ) < 0 , J ( u 2 ) > 0 , J ( u 3 ) > 0 while T 1 ( M ( u 2 ) ) = 0 and T 1 ( M ( u 3 ) ) = 0 . Above information helps to claim that u 1 is the saddle point and u 2 , u 3 are the center points (see Fig. 7(a)).

Phase portraits of nonlinear dynamical system.
Fig. 7
Phase portraits of nonlinear dynamical system.

For this case, the phase portraits of the nonlinear dynamical system (47) is presented in Fig. 7(a). This phase portrait encompasses a class of SNPO(3,1), where the family of SNPO(3,1) carries all included equilibrium points of the considered dynamical model. It also carries two families of NPO(1,0), which accommodates u 2 and u 3 . There is also a pair of NHO(1,0) at u 1 which carries u 2 and u 3 .

3.0.2

3.0.2 K 2 ( γ + α ) < 0 , ( - k 2 α - ω ) > 0 , β > 0

A single equilibrium point u 1 incorporated in the system (47), where J ( u 1 ) < 0 thus u 1 is a saddle point. (see Fig. 7(b)).

3.0.3

3.0.3 K 2 ( γ + α ) > 0 , ( - k 2 α - ω ) > 0 , C < 0 . or K 2 ( γ + α ) < 0 , ( - k 2 α - ω ) < 0 , β > 0

There are three equilibrium points of system (47) u 1 , u 2 , and u 3 for this type. For this J ( u 1 ) > 0 and T 1 ( M ( u 1 ) ) = 0 , so u 1 is the center point, while J ( u 2 ) > 0 and J ( u 3 ) > 0 with Poincarè index is zero thus u 2 , u 3 are cusp points, (see Fig. 8(a)). The phase portraits of the carried system of non-linear ODEs (47) is given in Fig. 8(a). This phase portrait carries a family of NPO(1,0) which envelops u 1 .

Phase portraits of nonlinear dynamical system.
Fig. 8
Phase portraits of nonlinear dynamical system.

3.0.4

3.0.4 K 2 ( γ + α ) > 0 , ( - k 2 α - ω ) > 0 , β > 0

Only single equilibrium point u 1 is assimilated in (47). For which J ( u 1 ) > 0 and T 1 ( M ( u 1 ) ) = 0 , thus u 1 is a center point. The phase portraits for this case is presented in Fig. 8(b), which shows that there is a family of NPO(1,0) which carries u 1 .

4

4 Conserved quantities

In this section, nontrivial conserved quantities (Nother et al. (1918) and Olver (1986)) are enumerated by using the method given by Anco and Bluman Anco and Bluman, 2002. They advocated a systematic approach to built non-trivial conservation laws.

Phase portraits of nonlinear dynamical system with different energy levels.
Fig. 9
Phase portraits of nonlinear dynamical system with different energy levels.

4.1

4.1 Multiplier approach

A system containing two partial differential equations of second order with three independent χ = ( t , x ) and two dependent variables Ψ = ( u , v ) and is denoted by

(54)
R 1 [ Ψ ] = F 1 χ , Ψ , Ψ χ , , Ψ χ χ , R 2 [ Ψ ] = F 2 χ , Ψ , Ψ χ , , Ψ χ χ , where Ψ χ and Ψ χ χ stand for the first and second order partial derivatives of the dependent variables in connection with independent variables, respectively. Let U = ( U 1 , U 2 ) represents the stream of arbitrary functions of independent variable χ , U χ and U χ χ etc. A set of multipliers (factors, characteristics) Λ = ( Λ 1 , Λ 2 ) yields a divergence expression for the system given in Eq. (54) if the identity
(55)
Λ 1 [ U ] R 1 [ Ψ ] + Λ 2 [ U ] R 2 [ Ψ ] = D χ C χ [ U ]
holds for arbitrary function U ( χ ) . In Eq. (55), T χ are called the conserved densities (fluxes) while D χ is the total derivative:
(56)
D χ = χ + Ψ χ Ψ χ + ·

If U = ( U 1 , U 2 ) is the solution of Eq. (54), form Eq. (55), one can derive the local conserved quantity by using the following equation

(57)
D χ T χ [ Ψ ] = 0 .

System (54) contains the set of multipliers for the conserved quantities if and only if following identity holds:

(58)
δ δ Ψ Λ 1 [ U ] R 1 [ Ψ ] + Λ 2 [ U ] R 2 [ Ψ ] = 0 , where δ δ U is said to be Euler operators and defined as:
(59)
δ δ Ψ = Ψ - D χ Ψ χ + ·

Eq. (58) leads to a set of over-diagnosed system of determining equations in expressions of multipliers Λ = ( Λ 1 , Λ 2 ) . The solution of the obtained determining equations with some computation further gives the conserved quantities. In this section, first-order nontrivial conserved quantities are computed by using the method given by Anco and Bluman Anco and Bluman, 2002. They advocated a systematic approach to built non-trivial conservation laws. In this method, multipliers Λ of precise order for undertook problem is required, additionally, which is used to get their analogous fluxes ϒ . Each class of multipliers and fluxes fabricates a local conservation law D ϒ = 0 holding for all solutions of considered differential equations.

4.2

4.2 Conserved quantities

In this section, conserved quantities of Eq. (1) are computed (Anco and Bluman, 2002; Cheviakov, 2007).

Eq. (44) with complex envelope:

(60)
Φ ( t , x ) = u ( t , x ) e iv ( t , x ) converts into a complex partial differential equation, after splitting into real and imaginary parts it yields:
(61)
β u 3 - uv t - α u ( v x ) 2 + ( α + γ ) u xx = 0 , u t + 2 α u x v x + α uv xx = 0 .

Substituting system (61) in Eq. (58) gives:

(62)
δ δ Ψ Λ 1 ( β u - uv t - α u ( v x ) 2 + ( α + γ ) u xx ) + Λ 2 ( u t + 2 α u x v x + α uv xx ) = 0 .

Equating the coefficients of derivatives of dependent variables concerning independent variables in Eq. (62), we get a linear homogeneous over-determined system of partial differential equations. After solving the obtained system of partial differential equations for Λ 1 = Λ 1 ( t , x , u , v , u t , v t , u x , v x ) and Λ 2 = Λ 2 ( t , x , u , v , u t , v t , u x , v x ) with the assistance of Maple, we will secure the following consequences:

(63)
Λ 1 = ( c 1 t + c 3 ) u x + c 2 u t , Λ 2 = c 1 v x t - c 1 x 2 α + c 2 v t + c 3 v x + c 4 u .

The next step is to find the fluxes by using the multiplier given in (63). For instance, the multipliers Λ 1 and Λ 2 for the constants c i give the following conservation laws:

  1. Λ 1 1 = tu x , Λ 2 1 = ( v x t - x 2 α ) u

    (64)
    T t 1 = u 2 4 α ( 2 α tv x - x ) , T x 1 = 1 4 β tu 4 - 2 tv t u 2 + 2 α tu 2 v x 2 + 2 α tu x 2 + 2 γ tu x 2 - 2 v x xu 2 .

  2. Λ 1 2 = u t , Λ 2 2 = v t u

    (65)
    T t 1 = 1 4 β u 4 - 2 α u 2 v x 2 - 2 u x 2 α - 2 u x 2 γ , T x 1 = α u 2 v t v x + α u t u x + γ u t u x .

  3. Λ 1 3 = u x , Λ 2 3 = v x u

    (66)
    T t 1 = u 2 v x 2 , T x 1 = 1 4 β u 4 - 2 v t u 2 + 2 α u 2 v x 2 + 2 u x 2 α + 2 u x 2 γ .

  4. Λ 1 4 = 0 , Λ 2 4 = u

    (67)
    T t 1 = u 2 2 , T x 1 = v x α u 2 .

4.3

4.3 Conclusion

To be brief, the new direct extended algebraic method (Rezazadeh, 2018) is applied to perceive the exact solutions of the resonant non-linear Schrödinger equation with Kerr law nonlinearity. The proffered technique gave a class of solutions that may be worthwhile for explaining certain physical phenomena accurately. Moreover, the physical composition of these solutions is described via their graphical presentation. Four portraits of a dynamical system (47) are obtained and the existence of the traveling wave solutions is discussed as well. It is observed that obtained solutions are new and not available in the literature. Further, all possible cases of the system parameters are considered by using the phase portraits, and the effect of different situations are shown in detail. Moreover, nontrivial, first order and new conserved quantities are given by using the multiplier approach.

In the future, the perturbed dynamical structure of the considered model can be investigated in light of chaos and sensitivity analysis.

5

5 Ethical standards

The authors state that this research complies with ethical standards. This research does not involve either human participants or animals.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. , , , , , , . Exact solutions, conservation laws, bifurcation of nonlinear and super nonlinear traveling waves for Sharma-Tasso-Olver equation. Nonlinear Dyn.. 2018;94:1791-1801.
    [Google Scholar]
  2. , , . Direct construction method for conservation laws of partial differential equations Part II: General treatment. Eur. J. Appl. math.. 2002;41:567-585.
    [Google Scholar]
  3. , , , , . Dark optical solitons and conservation laws to the resonance nonlinear Schrödingers equation with Kerr law nonlinearity. Optik. 2017;147:248-255.
    [Google Scholar]
  4. , , , . New Complex Hyperbolic Structures to the Lonngren-Wave Equation by Using Sine-Gordon Expansion Method. Appl. Math. Nonlinear Sci.. 2019;4(1):141-150.
    [Google Scholar]
  5. , . Über die Erhaltumgsatzeder Elektrodynamik. Math. Ann.. 1921;84:258-276.
    [Google Scholar]
  6. , , . Introduction to Non-Kerr Law Optical Solitons. Boca Raton, FL, USA: CRC Press; .
  7. Cattani, C., scalia, M., Laserra, E., Bochicchio, I., Nandi, K.K., 2013. Correct light deflection in Weyl conformal gravity. Phys. Rev. D 87 (4), 047503.
  8. , , , . Information and modeling in complexity. Math. Prob. Eng.. 2012;2012 868413
    [Google Scholar]
  9. , . GeM software package for computation of symmetries and conservation laws of differential equations. Comput. Phys. Commun.. 2007;176:48-61.
    [Google Scholar]
  10. Dubinov, A.E., Kolotkov, D. Yu, Sazonkin, M.A., 2012. Supernonlinear waves in plasma. Plasma Phys. Rept., 38 (10), 833–844.
  11. , , , . New Complex and Hyperbolic Forms for Ablowitz-Kaup-Newell-Segur Wave Equation with Fourth Order. Appl. Math. Nonlinear Sci.. 2019;4(1):105-112.
    [Google Scholar]
  12. , , , . Soliton solutions of the resonant nonlinear Schrödingers equation in optical fibers with time dependent coefficients by simplest equation approach. J. Mod. Opt.. 2013;60(19):1627-1636.
    [Google Scholar]
  13. , , , , . An efficient computational method for solving nonlinear stochastic It integral equations: application for stochastic problems in physics. J. Comput. Phys.. 2015;283:148-168.
    [Google Scholar]
  14. , . A new conservation theorem. J. Math. Anal. Appl.. 2007;333:311-328.
    [Google Scholar]
  15. , , . Noether-type symmetries and conservation laws via partial Lagrangian. Nonlinear Dyn.. 2006;45:367-383.
    [Google Scholar]
  16. , , , . On optimal system, exact solutions and conservation laws of the modified equal-width equation. Appl. Math. Nonlinear Sci.. 2018;3(2):409-418.
    [Google Scholar]
  17. , , . Travelling waves and conservation laws of a (2+1)-dimensional coupling system with Korteweg-de Vries equation. Appl. Math. Nonlinear Sci.. 2018;3(1):241-254.
    [Google Scholar]
  18. , , , . Effect of Hall current on the magneto hydrodynamic free convective flow between vertical walls with induced magnetic field. Eur. Phys. J. Plus. 2018;133:207.
    [Google Scholar]
  19. , , , . Influence of heat source/sink on MHD flow between vertical alternate conducting walls with Hall effect. Physica A. 2020;544 123562
    [Google Scholar]
  20. , , , , . Bifurcation analysis and solutions of a higher order nonlinear Schrödinger’s equation. Math. Prob. Eng.. 2015;2015 408586-408585
    [Google Scholar]
  21. , , , . Solutions and conservation laws of a generalized second extended (3+1)-dimensional Jimbo-Miwa equation. Appl. Math. Nonlinear Sci.. 2018;3(2):459-474.
    [Google Scholar]
  22. Nother, E., 1918. Invariante Variationsprobleme, Nacr. König. Gesell. Wissen., Göttingen, Math.-phys. Kl. Heft, 2, 235–257 (English translation in Transport Theory and Statistical Physics, 1(3), (1971) 186–207.
  23. , . Applications of Lie Groups to Differential Equations. New York: Springer-Verlag; .
  24. , . New soliton solutions of the complex Ginzburg-Landau equation with Kerr law nonlinearity. Optik. 2018;167:218-227.
    [Google Scholar]
  25. , , , . Wavelet analysis of the evolution of a solitary wave in a composite material. Int. Appl. Mech.. 2004;40:311-318.
    [Google Scholar]
  26. , , . Numerical simulation for fractional-order Bloch equation arising in nuclear magnetic resonance by using the Jacobi polynomials. Appl. Sci.. 2020;10(8):2850.
    [Google Scholar]
  27. Singh, H. Numerical simulation for fractional delay differential equations. Int. J. Dyn. Control doi: 10.1007/s40435-020-00671-6.
  28. , . Analysis for fractional dynamics of Ebola virus models. Chaos Solitons Fract.. 2020;138 109992
    [Google Scholar]
  29. , , , . Solving Non-Linear Fractional Variational Problems Using Jacobi Polynomials. Mathematics. 2019;7(3):224.
    [Google Scholar]
  30. , , , , . Legendre spectral method for the fractional Bratu problem. Math. Meth. Appl. Sci.. 2020;43(9):5941-5952.
    [Google Scholar]
  31. , , , . Magnetohydrodynamic flow in horizontal concentric cylinders. Inter. J. Indus. Math.. 2019;11:89-98.
    [Google Scholar]
  32. , , , . Bifurcation analysis and the travelling wave solutions of the Klein-Gordon-Zakharov equations. Pramana J. Phys.. 2013;80(1):41-59.
    [Google Scholar]
Show Sections