Abstract
This study investigates new optical solutions to a model with an integrable equation for wave packet envelopes. For this purpose, we have employed two reliable techniques involving the modified extended tanh function and the exponential rational function procedures. We have also given the 3D graphics of the obtained solutions.
1. Introduction
Numerous complex phenomena that are encountered in relativity, fluid dynamics, mathematical physics, economics, and optical fibers are modeled via nonlinear partial differential equations (NPDEs) [1]. Recently, new methods have been recommended for finding solutions to the NPDEs. Some notable studies have been discussed. Alquran has obtained periodic solutions and bell-shaped soliton solutions to the two-component evolutionary system of homogeneous KdV equations of the second order, the Gardner equation, Benjamin–Bona–Mahony equation, and the Cassama-Holm equation via the sine-cosine method [2]. Akbulut et al. have employed the exp -expansion procedure to the ill-posed Boussinesq equation and Zakharov–Kuznetsov–Benjamin–Bona–Mahony equation to obtain trigonometric, hyperbolic, and rational functions for these equations [3]. In 2015, Kaplan et al. applied the auxiliary equation procedure to the nonlinear Zoomeron equation, equal width wave equation, and coupled Higgs equation [4]. Ali et al. discussed cubic optical solitons in a polarization-preserving fiber described via nonlinear Schrödinger equation (NLSE) via the new extended direct algebraic method [5]. Yao et al. used the improved generalized Riccati equation mapping procedure to get traveling wave solutions for the higher-order Sasa-Satsuma equation [6]. Durur et al. applied the subequation procedure to find exact solutions for the KdV6 equation [7]. In 2021, Mirzazadeh et al. have employed the generalized Riccati equation mapping method [8]. Ma et al. have employed the first integral, the exp -expansion, and sine-Gordon expansion procedures to extract soliton solutions of the Hirota–Maccari system [9]. In 2016, Osman used the generalized unified method for the Zakharov–Kuznetsov equation [10]. Inc et al. utilized the functional variable procedure [11]. Leta et al. applied the bifurcation technique to the (2 + 1)-dimensional Bogoyavlenskii coupled system [12]. Rezazadeh et al. applied the extended rational sin-cos method to the chiral nonlinear Schrödinger equation [13].
In recent years, optical solitary waves, temporal and spatial solitons, have been discussed by many researchers theoretically and experimentally. For example, in 2018, Biswas et al. have founded optical solitons for the Lakshmanan–Porsezian–Daniel model via the modified simple equation procedure [14]. Manafian, in 2016, has founded optical soliton solutions for NLSE via the tan -expansion procedure [15]. Kudryashov, in 2020, has proposed a procedure to find highly-dispersive optical solitons of NPDEs [16].
This paper studiedwhere , , , , , and are parameters, and is a complex-valued function. It is known that this equation is the generalization of some famous equations. It is utilized to describe traveling wave solutions in optical fibers. This equation is reduced to NLSE when we choose , and it models the slowly varying envelope dynamics of a weakly nonlinear quasimonochromatic wave packet in dispersive media [17]. Also, by setting in equation (1), the cubic-quintic NLS (cqNLS) equation is obtained which is used in the nonlinear optics [18, 19]. cqNLS is used for the nuclear hydrodynamics with Skyrme forces [20], the optical pulse propagations in dielectric media of non-Kerr type [21], and description of the boson gas with two and three-body interactions [22]. For this reason, we can consider equation (1) as a generalization of the NLSE. It could be utilized for describing nonlinear processes in optics [23].
The paper consists of the following sections: we present the preliminary information in Section 2. Then, we give a description of utilized techniques, respectively, and The modified extended tanh-function procedure and the exponential rational function procedure in Section 2.2. The application of the given methods is given in Section 3. Finally, the conclusion of the whole research is given in Section 4.
2. Preliminary Information
In this section, we firstly introduce the preliminary information for reducing NPDE to an ordinary differential equation (ODE). For this purpose, we will consider the following NPDE:where inherits and its partial derivatives, and is a complex-valued function.
Wave transformation is given bywhere , and are the constants. If we apply transformation (3) to equation (2) and equate the real and imaginary parts to zero, we obtain two equation systems. Then, we solve these systems to find conditions for parameters and use the values of the parameters. Therefore, we find an ODE, which will be integrated with respect to possible times, as follows:where the partial derivative is given with respect to [24].
The modified extended tanh function and the exponential rational function procedures are applicable to this reduced ODE.
2.1. The Modified Extended Tanh-Function Procedure
We will assume the exact solution of equation (4) byto employ the modified extended tanh-function procedure. Here, are constants, and they will be calculated; and cannot be zero at the same time. In equation (5), the term is called the balancing term. This term can be found balancing the nonlinear terms and the highest order derivatives. is the solution of following auxiliary equation:where is given as a constant. The solutions of equation (6) are given in three cases as follows:(1)When , the hyperbolic solution is(2)When , the trigonometric solution is(3)When the rational solution is
Surrogating equation (5) into equation (4) along with equation (6), then accepting all the coefficients of is zero, we get an algebraic equation system for If we solve an obtained determining equation system with the Maple, we find a variety of exact solutions [25].
2.2. The Exponential Rational Function Procedure
We will assume the solution of equation (4) as follows [26, 27]:where are constants. Here, is the balancing number as the previous technique. We surrogate equation (10) into equation (4) and collect all terms with the same order of together, and we make into the left-hand side of equation (4) and another polynomial in . Then, we equate each coefficient of this polynomial to zero, and it gives a set of algebraic equations for unknown parameters. Finally, we solve the equation system to construct a variety of exact solutions for equation (4).
3. Application of the Given Methods
In the current section, we will see new traveling wave solutions of equation (1). For this reason, firstly, we give the mathematical analysis of the considered equation.
3.1. Mathematical Analysis
If we substitute equation (3) into equation (1) and then equate the real and imaginary parts to zero, we find
From the first equation in the above system, we obtain
Using (12) in the second equation in System (11), we get following reduced ODE:
We will separately apply the modified extended tanh-function procedure and the exponential rational function procedure to this reduced ODE to obtain different exact solutions.
3.2. Modified Extended Tanh-Function Procedure
Homogeneous balance is obtained as from equation (13); thus, the solution of equation (13) is given by
Surrogating equation (14) along with equation (6) into equation (13), then assuming the coefficients of as zero,is acquired. If we solve System (15), we obtain values of the constants.
3.2.1. Solution Family 1
When
When ,
When ,
3.2.2. Solution Family 2
When ,
When
When ,
3.2.3. Solution Family 3
When ,
When ,
When , we obtain the same solution as the third solution in solution family 2.
3.2.4. Solution Family 4
When ,
When
When , we obtain the same solution as the third solution in solution family 2.
3.2.5. Solution Family 5
When we obtainwhere .
When , we getwhere .
When , we findwhere .
3.3. The Exponential Rational Function Procedure
Since the balancing number is 1, according to the exponential rational function procedure, the solution is
If we surrogate equation (34) into equation (13) and collect all the terms with the same order of together, and we make into the left-hand side of equation (4) and another polynomial in . Then, we equate each coefficient to zero to get an algebraic equation set for unknown parameters as follows:
Thus, we obtain the following solutions: Solution 1: By substituting these values, we find Solution 2: By substituting these values, we find Solution 3: By substituting these values, we find
4. Conclusion
NPDEs have a broad field of applications in fundamental sciences such as computer science, biology, chemistry, mathematics, physics, and engineering sciences such as chemical engineering, electricity, mechanics, and materials. The exact solutions of NPDEs have attracted considerable observation from researchers. For this reason, we have studied modified extended tanh-function and the exponential rational function procedures to determine various types of exact solutions of an integrable equation for wave packet envelope. Also, we have obtained 3D graphics of the obtained solutions (Figures 1–8). The obtained exact solutions have been verified by substituting them back into the original equation using Maple 19. According to our point of view, our obtained results have the potential to be applicable in nonlinear dynamics and optics. The obtained solutions can be classified as periodic and soliton solutions. For future works, supernonlinear, quasiperiodic, chaotic wave, and solitonic solutions will be potentially obtained. Finally, despite using Maple 19 during the calculations, the authors have noted that wxMaxima can be chosen due to the fact that is free and open-source software.








Data Availability
The datasets used or analyzed during the current study are available from the corresponding author on reasonable request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.