Bifurcation, chaotic behavior, sensitivity analysis, and dynamical investigations of third-order Schrödinger equation using new auxiliary equation method

Abstract

This current study presents a precise analytical examination of the generalized third-order nonlinear Schrödinger equation through the application of the new auxiliary equation method. The approach provides several classes of exact solutions, such as V-shaped, dark soliton, periodic, kink, and anti-kink soliton solutions, which prove its effectiveness in solving higher-order nonlinear wave equations. The derived solutions are well depicted through 2D, contour, and 3D plots to show their spatial and temporal evolution features. A complete dynamical system analysis is carried out by Galilean transformation, showing the system behavior through accurate phase portraits and bifurcation diagrams. The analysis offers valuable information on stability of the solutions and transition processes amongst solution types. The system sensitivity analysis to parameters provides significant stability conditions for the solutions obtained. All the outcomes are derived by strict analytical means, and graphical plots are used to support the mathematical analysis. The study makes a contribution to the theoretical basis of nonlinear wave propagation and offers a sound framework for similar nonlinear evolution equations.

Introduction

The nonlinear Schrödinger equation (NLSE) is fundamental to the modeling of nonlinear wave phenomena in many different physical systems, especially in nonlinear optics, plasma physics, and fluid dynamics. Perhaps its most important application is in the description of the propagation of optical pulses in nonlinear media, for example, optical fibers, where it describes the interaction between dispersion and nonlinearity. In plasma physics¹, the NLSE controls the dynamics of wave packets in Langmuir and ion-acoustic waves. In hydrodynamics², it also models the evolution of deep-water wave trains. A distinctive feature of the NLSE is that it can support soliton solutions. Solitons are stable, localized wave structures that maintain their shape and velocity over time as a result of a fine balance between nonlinearity and dispersion. These solitons not only improved our theoretical knowledge of nonlinear systems but also resulted in practical applications, like distortion-free long-distance communication using fiber-optic cables.

Based on the fundamental applications of NLSE, it is extended by adding higher effects like third-order dispersion, self-steepening, and stimulated Raman scattering. The most significant advancement is the generalized third-order nonlinear Schrödinger equation (GTNSE)³. The higher effects are required in the transmission of ultrashort optical pulses in high-speed communications or in complex nonlinear materials with complicated dispersion. The GTNSE offers a more precise model for describing complex wave dynamics in nonlinear optics, plasma waves, and fluid dynamics where second-order models are inadequate. Similar to the NLSE, the GTNSE also possesses soliton solutions, which can take more complicated forms like dark, bright, kink-type, and rogue solitons based on how the higher-order terms interact. The soliton solutions of the GTNSE are crucial in understanding advanced phenomena like pulse compression, and stability of localized wave packets in nonlinear and dispersive materials. Due to its relevance in modeling real phenomena, a broad range of analytical and numerical methods has been developed to study its solutions.

Development of new analytical methods in recent years has significantly improved our capacity for finding the exact solutions of such nonlinear evolution equations. Some of them are modified , the first-integral method⁴, the -expansion method⁵, the improved -expansion method⁶, Sine–cosine method⁷ and the sub-ODE method⁸. In addition, the Hirota bilinear method^9,10, Modified kudryashov expansion method¹¹, rational sine and rational cosine methods^12,13, Kudryashov expansion¹⁴, extended tanh(coth)-expansion method¹⁵ and updated rational sine-cosine-function method have gained attention¹⁶. Overall, these methods highlights the importance of analytical method in soliton dynamics. Sardar sub equation method¹⁷ and new extended and modified rational expansion method¹⁸ have been used to successfully solve the GTNSE and derive various types of solutions ranging from bright, dark, and singular solitons to periodic and chaotic waveforms. Ahmad, et al.¹⁹ investigate the fractional GTNSE by using generalized riccati equation and mapping method. Badshah, et al.²⁰ explore the soliton soutions by using improved F-expansion method, extended hyperbolic function technique and Adomian decomposition method and applied modulation stability to verify the results.

The primary objective of this research is to obtain traveling wave solutions of the GTNSE by an effective and rigorous technique known as the new auxiliary equation method (NAEM). The technique enables us to simplify the original PDE to an easier ODE by appropriate transformations and the introduction of auxiliary functions, enabling us to obtain explicit analytical solutions. Using the technique, we obtain several exact solutions, including soliton, periodic, and rational-type waveforms, demonstrating the intricate math structure of the model behind. In addition to discovering exact solutions, we perform a careful analysis of the traveling wave system. Bifurcation analysis shows how system parameters affect the model’s performance and the shape of the solutions. Specifically, we investigate dynamics of the reduced ODE system, identify the transitions between the many dynamical states and determine equilibrium points, phase portrait and inspect conditions under which saddle-node, center, or spiral-type bifurcations occur, which improves our comprehension of the solutions’ structural characteristics. This work integrates the new auxiliary equation method and bifurcation analysis to further investigate the intricate behavior modelled by the GTNSE. Our study generates new exact solutions and offers informative analytical explanations for the broader area of nonlinear wave theory.

The paper is organized into sections. Governing equation is describe in Sect. 2. Section 3 covered algorithm of new auxiliary method in detail. In Sect. 4, method is implemented. Visual evaluation and discussion is explained in Sect. 5. Bifurcation analysis is done in Sect. 6. Chaotic behaviour and sensitivity analysis is examined in Sects. 7 and 8, respectively. Physical significance of study is describe in Sect. 9. Then conclusion is provided.

Governing equation

The GTNSE is a mathematical model, an extension of the original nonlinear Schrödinger equation, which is applied to model wave propagation in nonlinear media. The GTNSE²¹ is:

2.1

where the dispersive coefficient is indicated by parameters and , the coefficient of cubic nonlinearity is represented by parameter and is a complex function. Eq. (2.1) has been used to model ultrashort pulses in optical fibers. The above equation typically also contains the term of second-order derivative. But as soon as the derivative expression of third order is added, the second-order derivative term can be eliminated through a gauge transformation formation. If , then Eq. (2.1) is termed as modified Korteweg-de Vries or Hirota equation in complex form, which can be solved by the inverse scattering transform technique.

Algorithm of new auxiliary equation method

The procedures for using the auxiliary equation approach²² to find traveling wave solutions of nonlinear partial differential equations are described in this section.

Consider the general form of nonlinear partial differential equations (NLPDE), involving two independent variables x and t, expressed as follows:

3.1

where function is complex function. is the polynomial expression involving functions and their classical derivatives, where the maximum order of derivative and nonlinear terms are co-related. Detailed and concise flowchart of solutions procedure is given in Fig. 1.

• Step 1: Apply the wave transformation, given as follow:

  3.2

  Here represents soliton’s veocity, wave number and frequency, respectively. Phase constant is denoted by .

• Step 2: The transformation in step 1, converts NLPDE in Eq. (3.1) to nonlinear ordinary differential eqation :

  3.3

• Step 3: Now according to NAEM, the Eq. (3.3) has series solution of the form:

  3.4

  where are constants. satisfies the following equation:

  3.5

  where , and are real constant parameters. Depending on these constants, specific solutions of (3.5) are determined.

• Step 4 The positive integer N is obtained by balancing highest-order derivative with the highest-degree nonlinear term in Eq. (3.3).

• Step 5 Then, put Eq. (3.4) along with Eq. (3.5) in Eq. (3.3). An algebraic system is obtained by setting all like-terms of equal to zero. Solve the obtained system to find the values of unknown parameters. Then, put these parameters into Eq. (3.4) along with solutions of stated below : Following families are produced using general solution with :

  • Cluster 1: When , and , then

  • Cluster 2: When , and , then

  • Cluster 3: When and , then

  • Cluster 4: When and , then

  • Cluster 5: When and , then

  • Cluster 6: When and , then

  • Cluster 7: When , then

  • Cluster 8: When , and , then

  • Cluster 9: When and , then

  • Cluster 10: When , then

  • Cluster 11: When and , then

  • Cluster 12: If and , then

  • Cluster 13: If , then

  • Cluster 14: When , then

  • Cluster 15: When , then

  • Cluster 16: When , then

  • Cluster 17: When , then

  • Cluster 18: When , then

  • Cluster 19: When and , then

  • Cluster 20: When , then

  Fig. 1 [Images not available. See PDF.]

  Solution procedure for GTNSE via NAEM.

  Applications

  Travelling wave transformation for GTNSE

  Substitute transformation provided in Eqs. (3.2) to (2.1), then separate real part (Re) and imaginary part (Im) of resulting equation.

  4.1

  Now, by integrating imaginary part and setting constant of integration to zero, yields:

  4.2

  The Eq. (4.1) and Eq. (4.2) are equal under following constraints:

  Plug these constraints into Eq. (4.1), yields

  4.3

  Application of new auxiliary equation method

  Using balancing principle that equate the highest order derivative with the highest degree nonlinear term in Eq. (4.3), results in . Now, with help of Eq. (3.4), series solution of Eq. (4.3) becomes:

  4.4

  A nonlinear system in is obtained by substituting the above solution along with Eq. (3.5) into Eq. (4.3). Now, arrange all terms in power of . An algebraic system is obtained by setting each coefficients of generated polynomial equal to zero, that is Then, solve this system by using software Wolfram Mathematica, the following parameters are obtained:

  Using the families of solutions described in section 3, the analytical solutions of GTNSE are presented as follows:

  • Cluster 1: When , and , then

  • Cluster 2: When , and , then

    4.5

  • Cluster 3: When and , then

  • Cluster 4: When and , then

    4.6

    4.7

  • Cluster 5: When and , then

    4.8

  • Cluster 6: When and , then

    4.9

  • Cluster 7: When , then

  • Cluster 8: When , and , then

  • Cluster 9: When and , then

  • Cluster 10: When , then

  • Cluster 11: When along with , then

  • Cluster 12: When and , then

  • Cluster 13: When , then

  • Cluster 14: When , then

  • Cluster 15: When , then

  • Cluster 16: When , then

  • Cluster 17: When , then

  • Cluster 18: When , then

  • Cluster 19: When and , then

  • Cluster 20: When , then

  where and .

  Limiting cases and comparison

  In addition to broadening the scope of current analytical techniques, the analytical solutions of the GTNSE obtained using the NAEM also recover a number of classical solutions when the parameters are limited to certain values. By choosing properly the auxiliary parameters and in the equation,

  the general solution simplifies to classical waves forms and solitons described in literature. Specifically, when , and , the solution leads to or , corresponding to kink and dark solitons. These are same as standard solutions for GTNSE derived by using - method²³. Similarly, when , trigonometric solutions such as or arise, corresponding with periodic wave solutions obtained by the modified -expansion approach^24,25. In the case , rational solutions of the form are derived, which represent algebraic or rogue wave-type solutions^26,27. Moreover, the NAEM framework consolidates other current methodologies. As a result, the NAEM proves itself effective in generating a range of accurate solutions in addition to acting as a unified framework incorporating a wide range of standard approaches as definitive examples. Such consistency with standard results thus verifies the method’s potency and reliance in solving nonlinear evolution equations, thus raising its theoretical significance in the broader application area of soliton theory.

  Visual evaluations and discussion

  This section presents a detailed analysis of the traveling wave solutions derived by using the NAEM for GTNSE. This equation has extensive application in the modeling of nonlinear complex wave processes in nonlinear optics, plasma physics, hydrodynamics, and ultra-short pulse propagation. With the implementation of the novel method, a richness of exact wave solutions have been obtained successfully, such as dark solitons, V-type solitons, periodic waves, kink, and anti-kink wave profiles. This diversity of profiles indicates the success and effectiveness of the proposed method in generating physically meaningful solutions. In contrast to earlier published results, the findings here prove higher diversity and accuracy in the solutions achieved. The invariance and stability of such solutions under propagation attest to their usability within the modeling of actual nonlinear systems. To further clarify the dynamical behavior and space evolution of the solutions derived, graphical illustrations are presented in 2D, 3D, and contour plots. Contour plots specifically assist in the description of the wave intensity distribution and topological structure in various spatial regions.

  The results of applying NAEM is discussed below:

  • Dark soliton represented in Fig. 2, obtained from Eq. (4.5), create stable intensity dip with jumps in phase when nonlinearity compensates normal dispersion. These solutions have good propagation properties and can be used in optical switching applications and dark pulse generation.

  • Figure 3 illustrate anti-kink soliton solution provided in Eq. (4.6), with inverted phase profiles in contrast to kink solitons.

  • Figure 4 illustrate periodic solution, obtained from Eq. (4.7), demonstrate stable intensity modulation. Such solutions are useful for designing all-optical buffers and waveguide arrays.

  • V-shaped soliton as shown in Fig. 5, is found in Eq. (4.8) which are defined by their vertex-like intensity profiles with conjugate phase jumps. These solutions have uses in optical pulse compression and frequency comb generation.

  • Kink solitons in Fig. 6 arise in Eq. (4.9), with monotonic phase transitions between stable solutions. Such topological solitons are used to make optical logic gates.

  Parameters values for all solutions-type (dark, anti-kink, periodic, V-shaped and kink) is given in Table 1 for clarity. These visualizations and propagation stability provide an intuitive representation of the wave structures, phase transitions, and space distributions that are necessary for the analysis of the physical mechanism.

  Fig. 2 [Images not available. See PDF.]

  Dark soliton solution of absolute behaviour for Eq. (4.5) by employing the parametric values: (a) 3D plot illustrates the intensity dip’s progression along the x-t plane; (b) 2D profile demonstrate the dip’s propagation at (red line), (purple line) and (green line); (c) a contour plot shows the phase jump over the soliton core.

  Fig. 3 [Images not available. See PDF.]

  Absolute behaviour for Eq. (4.6) represent anti-kink soliton obtained by employing the parametric values: (a) 3D plot illustrates the transition profile along the x-t. plane; (b) 2D profiles demonstrate the phase shift at (red line), (purple line) and (green line); (c) Contour graph shows the kinetics of phase reversal throughout the propagation domain.

  Fig. 4 [Images not available. See PDF.]

  Periodic wave solution of real behaviour for Eq. (4.7) with the parametric values: (a) 3D plot illustrates the periodic wave along the x-t plane; (b) 2D profiles demonstrate wave at (red line), (purple line) and (green line); (c) Contour plot depicts phase matching within the periodic structure.

  Fig. 5 [Images not available. See PDF.]

  V-soliton solution of absolute behaviour for Eq. (4.8) with the parametric values: (a) 3D plot illustrates the propagation of V-shaped profile along the x-t plane; (b) 2D profiles at (red line), (purple line) and (green line); (c) Contour map depicts the combined effect of spectral broadening and phase shift dynamics.

  Fig. 6 [Images not available. See PDF.]

  Kink soliton solution of imaginary behaviour for Eq. (4.8) with the parametric values: (a) 3D plot illustrates the phase transition along the x-t plane; (b) 2D plot demonstrate phase shift at (red line), (purple line) and (green line); (c) Contour map illustrate the phase shift dynamics.

  Table 1

  Summary of parameters used for different solutions.

  Figures	Parameter values
  Figure 2
  Figure 3
  Figure 4
  Figure 5
  Figure 6

  Bifurcation analysis

  Bifurcation in dynamical systems occurs when small parameter changes result in large changes in the system’s behavior. This can involve the appearance of periodic patterns, chaotic dynamics or new stable phases. The bifurcation hypothesis describes these abrupt changes and even forecasts how systems evolve through various phases of their behavior. One well-documented real-world example is the way that traffic patterns on a road change. With more cars on the road, the system can transition suddenly from free-flowing traffic to dense stop-and-go waves. This bifurcation describes how a small parameter change (e.g., in car density) results in a large change in the system’s behavior^28,29.

  Let , then the planar dynamical system of equation Eq. (4.3) is:

  6.1

  With Hamiltonian function:

  ,

  where, h is Hamiltonian constant, and . It is important to note that the four parameter cases in the phase portraits are not arbitrarily selected. They arise naturally from the different sign combinations of the coefficients and in the reduced dynamical system Eq. (6.1). Each case determines whether the equilibrium points act as centers, saddles, or cuspidal points, thereby covering all qualitatively distinct dynamical behaviors of the system. Physically, these cases correspond to regimes governed by different balances between dispersion and nonlinearity. For instance, positive or negative values of and represent scenarios relevant to optical fibers, plasma waves, and fluid systems under different operating conditions. Hence, the four cases are both mathematically exhaustive and physically meaningful, providing a comprehensive view of the system’s dynamics.

  In hamiltonian, represent kinetic energy, and denotes potential energy. For equilibrium points, put , then Eq. (6.1) becomes:

  6.2

  Solving this system yields,

  The determinant of the Jacobian matrix of Eq. (6.1) is Consequently,

  • If , then represent saddle point.

  • If , then represent center point.

  • If , then represent cuspidal point.

  Changing the values of parameters results in four case:

  • Case 1: Figure 7a illustrates the case with the equilibrium points and (2, 0) under the parameters , where (0, 0) is saddle point and represent center points.

  • Case 2: Figure 7b illustrates the case with equilibrium points and (3, 0) under the parameters , where (0, 0) denotes center point and represent saddle points.

  • Case 3: Figure 7c illustrates the case with equilibrium points and under the parameters , where (0, 0) denotes saddle point and represent center-like points.

  • Case 4: Figure 7d illustrates case with equilibrium points and under the parameters , where (0, 0) is center point and represent saddle points.

  Fig. 7 [Images not available. See PDF.]

  Phase diagrams for planar system in Eq. (6.1): (a) illustrate the phase portrait for Case 1; (b) illustrate the phase portrait for Case 2; (c) illustrate the phase portrait for Case 3; (d) illustrate the phase portrait for Case 4.

  Chaotic behaviour

  Chaotic behavior is dynamics that are seen as random or unpredictable but are really governed by hidden irregular patterns. These systems are highly sensitive to initial conditions. Chaotic systems are, however, founded on deterministic rules that are transparent and comprehensible through certain geometries or structures. Chaos can be observed in natural and artificial systems such as optics, atmospheric flow, and fluid dynamics.

  Here, we study the chaotic behavior of the resulting dynamical system through the structures of instability. In this study, we have a careful comparison of the 2D and 3D phase portraits. To proceed with this study, we focus our attention to the second term included in Eq. (6.1).

  7.1

  where, and represents the frequency distribution and external disturbance, respectively.

  Now, consider the parameter values . We examine how the system is affected by frequency and perturbation . Figures 8, 9 and 10 show chaotic behaviour of the system (6.1) for different values of and . The detailed description of these figures are described in the captions. These findings verify the existence of chaotic behavior in the system with the provided parameters. Generally, these graphs demonstrate the system to be chaotic, with sensitive dependence on initial values and no steady or recurring action.

  Fig. 8 [Images not available. See PDF.]

  Chaotic nature of the system (6.1) for and .

  Fig. 9 [Images not available. See PDF.]

  Chaotic nature of the system (6.1) for and .

  Fig. 10 [Images not available. See PDF.]

  Chaotic behaviour of the system (6.1) for and .

  To provide a quantitative verification of chaos in addition to qualitative phase portraits and time series, we computed the largest Lyapunov exponent (LLE) of the reduced planar dynamical system using the standard Benettin–Wolf algorithm with RK4 integration. The parameters were selected according to the chaotic regimes presented in Section 7. Specifically, for and corresponding to and , the LLE values were found to be positive across the tested cases: for and ; for and ; and for and . The positivity of these exponents confirms the presence of sensitive dependence on initial conditions, thereby offering strong numerical evidence of chaotic dynamics in the GTNSE system.

  Sensitivity analysis

  The sensitivity analysis is carried out to examine the reaction of the dynamical system to small variations in its initial conditions. Precisely, when small changes in initial values bring about slight differences in the evolution of the system, it is classified as low sensitivity. However, when small initial perturbations cause substantial variations in system behavior, it indicates high sensitivity. In this analysis, various solution curves are seen for constant parameters: , to analyze how the stability changes of Eq. (6.1), when the initial conditions are slightly modified. In Fig. 11, two solution paths are plotted with red and with bue line. In Fig. 12, two orbits are shown again: in red and in green. In Fig. 13, two additional initial condition pairs are used: in green and in blue. In Fig. 14, three initial conditions are tested: in red and in green. Graphical inspection shows that even slight variations in result in visible alterations in solution behavior, demonstrating the system’s strong sensitivity to initial conditions. This is a significant observation in confirming the responsiveness and rich dynamics of the model, which can lead to unstable behaviors under a small change in inputs.

  Fig. 11 [Images not available. See PDF.]

  Sensitivity response by considering two conditions: = (0.7, 0.2) and = (0.7, 1.5), which are represented by lines in red and blue, respectively.

  Fig. 12 [Images not available. See PDF.]

  Sensitivity response by considering two conditions: = (0.7, 0.1) and = (0.9, 0.3), which are represented by lines in red and green, respectively.

  Fig. 13 [Images not available. See PDF.]

  Sensitivity response by considering two conditions: = (0.9, 0.9) and = (0.6, 0.3), which are represented by lines in green and blue, respectively.

  Fig. 14 [Images not available. See PDF.]

  Sensitivity profile under conditions: (0.5, 0.3), (0.3, 0.1) and (0.7, 0.8) represented by blue, green and red, respectively.

  In addition to graphical inspection, the sensitivity of the system can also be interpreted in a quantitative sense. Small perturbations in initial conditions lead to trajectories that diverge at an exponential rate, which is consistent with a positive largest Lyapunov exponent. This implies that the growth rate of separation between initially close trajectories serves as a measurable indicator of sensitivity. Hence, while our visual results illustrate the phenomenon effectively, they are also supported by the underlying dynamical system theory that confirms sensitivity through norms and growth rates.

  It is worth noting that the sensitivity outcomes are directly aligned with the stability conditions established in the bifurcation analysis. In particular, the bounded and regular trajectories observed under small perturbations correspond to the center-type equilibria, confirming local stability. Conversely, the unbounded divergence of nearby trajectories reflects the presence of saddle-type equilibria, which denote instability. This consistency between bifurcation-based classification and sensitivity responses provides a robust validation of the dynamical behavior of the system and emphasizes how small perturbations may trigger transitions predicted by the bifurcation framework.

  Physical interpretation

  The GTNSE plays a central role in the description of nonlinear wave propagation in a variety of physical systems. The mathematical formulation includes higher-order dispersion, self-phase modulation, self-steepening phenomena, and nonlocal nonlinear interactions. Its implications are present in fields including nonlinear optics, fluid dynamics, plasma physics, and Bose-Einstein condensates, where the interplay between dispersion and nonlinearity is the controlling influence on the subtle dynamics of wave propagation. The applications of terms involving in equation are describe below:

  • Temporal evolution term is a measure of how the wave envelope u(x, t) varies with time, which describes how the shape of the wave evolves with time. Application Seen in optical fiber communication systems, in plasma waves, and in the investigation of how water waves evolve with time

  • Third-order dispersion term refers to higher-order dispersion effects, which illustrate how various frequency components of the wave travel at various velocities. Extremely important for propagation of ultrashort pulses in optical fibers, creation of supercontinuum, and study of internal ocean waves where third-order dispersion is important.

  • Cubic nonlinearity (Kerr effect illustrates how self-phase modulation occurs due to the Kerr effect. The refractive index varies as intensity varies, causing phase shifts and contributing to the formation of solitons. Application finds in optical fiber communications (soliton pulses), Bose-Einstein condensates, water wave packets.

  • The term explains the self-steepening phenomenon, in which areas of high-intensity waves phase shift faster, leading to pulse distortion and steepening. The phenomenon is important for short, high-intensity pulses in nonlinear fiber optics and is significant for supercontinuum light sources.

  • defines excited Raman scattering or higher order nonlinear dispersion. Enables energy transfer between higher and lower frequencies (self-frequency shift), typically manifesting itself as the Raman effect or nonlinear spectral recoil. The phenomenon is important in fiber optics (soliton self-frequency shift, Raman scattering), acoustic phonon coupling, energy transfer in plasmas.

  The NAEM is an analytical tool for obtaining exact solutions to nonlinear evolution equations, which are essential in mathematics, physics, and engineering. The method is effective in producing many different exact solutions-kink, anti-kink, periodic, V-type, and dark soliton solutions-characterizing real-world wave phenomena. Kink and anti-kink soliton solutions are characterized by sudden jumps between two states. This renders them suitable for the description of magnetic domain walls, nerve impulses in biological systems, and phase boundaries in condensed matter physics. In nonlinear optics, for instance, a kink soliton can possess the sudden leading edge of an optical pulse when it makes a jump between two states of intensity. This is significant in signal processing and optical switching devices. Periodic solutions have recurring wave patterns. They are useful in naturally repeating systems, like the propagation of energy in crystal structures, sound waves in periodic media, and water wave dynamics. Periodic solutions are especially important in the case of repeating boundary optical fibers and in photonic crystals, where wave propagation is governed by periodicity in the refractive index. Temporal dark solitons have been experimentally demonstrated in optical fibers, where dark solitons propagate as ”holes” on a continuous-wave background. One experiment used 5.3ps dark solitons generated from 36ps wide pulses generated from an 850nm titanium-doped sapphire laser, and then propagating 1km of optical fiber. This is used in nonlinear guided wave optics and optical switching. Kink solitons have been utilized as useful applications in the form of polarization switches and optical logic gates in nonlinear optics. They are employedin modeling dislocation dispersion in crystals and wave propagation in plasmas. Their stability and topological nature have been investigated through experimental and theoretical studies, and in applications in fluid dynamics and quantum field theory. Anti-kink solitons, which are often studied alongside kink solitons in specific systems, especially Sine-Gordon-type systems, are topological solitons like particle-antiparticle wave objects in condensate matter and nonlinear wave systems. Experimental and theoretical work has studied their dynamic interactions (kink-antikink collisions). V-shaped dark solitons were investigated in relation to their propagation characteristics in inhomogeneous optical fibers, their behavior and shape being controllable to enable soliton control in optical communication systems. Periodic wave train or periodic soliton solutions have been generated in birefringent fibers for pulse shaping and nonlinear fiber optics. Experimental demonstration shows varied envelope soliton pulses periodic or multi-peak structures such as W-shaped solitons. It is important to note that the physical relevance of the obtained solutions is strongly dependent on parameter selection. Our bifurcation and sensitivity analyses demonstrate that while stable soliton, periodic, and kink-type structures exist in specific regimes, certain parameter combinations lead to instability or chaotic behavior. In these cases, the Hamiltonian structure of the unperturbed system may be destroyed, resulting in non-conservation of energy and loss of physical applicability. Thus, the derived analytical solutions remain valid in form, but their stability and physical significance are restricted to ranges where the balance between dispersion and nonlinearity is preserved.

  Conclusion

  This paper investigates the GTNSE comprehensively by a new powerful method known as the auxiliary equation method. The proposed method is found to be very efficient for analytically solving considered equation when compared with classical methods of solutions. In contrast to classical procedures that can have difficulty with the intricate higher-order terms and complicated nonlinear interactions involved in the considered equation, NAEM method systematically reduces the PDE into a more tractable ODE through wave transformations and simplifications. This makes possible the derivation of exact solutions, including V-shaped, kink, dark soliton, periodic, and anti-kink solitons, without losing the key nonlinear and dispersive features of the system. Each of them is obtained and explained carefully mathematically. These solutions significantly enrich the known type of solutions for this significant class of nonlinear wave equations, demonstrating how efficient the method is in handling complicated nonlinear evolution equations. The dynamical behaviors of the system were explored in depth with the help of an advanced analytical Galilean transformation. Hamiltonian formalism provided valuable information about the conserved quantities and stability properties of the system, and phase portrait construction revealed the rich dynamics at play, such as periodic, quasi-periodic, and chaotic phases. Bifurcation analysis charted the significant points of transition between solution types and provided valuable information on parameter dependence and solution stability. To gain further insight, we conducted a careful sensitivity analysis to see how the system responds to changes in parameters and initial conditions. This was found to reveal strong solution patterns and to indicate when chaotic behavior can arise. The mathematical formalism reached here presents new avenues of research, including its extension to higher-dimensional systems, investigation of coupled nonlinear systems, and investigation of other solution classes using higher-level tools. Physically, the findings provide insight into nonlinear wave propagation in optical fibers, plasma, and fluid systems where higher-order effects like dispersion, self-steepening, and Raman scattering are important. Not only does this work enlarge the solution space of proposed equation but also sets a strong platform for applications in photonics, high-speed communications, and other nonlinear physical systems. The results here not only provide a basis for improving our simple conceptual image of nonlinear wave equations but also give a solid theoretical foundation for the investigation of similar systems in mathematical physics. Future prospects include the extension of these methods to other physically interesting nonlinear evolution equations and investigation of potential applications in other fields of physics and engineering.

  Acknowledgements

  The authors would like to acknowledge the Deanship of Graduate Studies and Scientific Research, Taif University for funding this work.

  Author contributions

  Mehreen Fatima: methodology, investigation, writing–original draft. Muhammad Abbas: methodology, supervision, writing–original draft, writing–review & editing. Yagoub.A.S.Arko: visualization, formal analysis, investigation, writing–review & editing. Tahir Nazir: supervision, writing–original draft, writing–review & editing. Asnake Birhanu: formal analysis, investigation, writing–review & editing. Muhammad Zain Yousaf: formal analysis, investigation, writing–review & editing. All authors have read and agreed to the published version of the manuscript.

  Funding

  The author(s) received no external funding for this study.

  Data availability

  The data and materials used to support the findings of this study are included in this article.

  Declarations

  Competing interests

  The authors declare no competing interests.

  References

  1.HelalMASeadawyAR<article-title>Benjamin-Feir instability in nonlinear dispersive wavesComput. Math. Appl.2012641135573568299253410.1016/j.camwa.2012.09.0062.ZabuskyNJKruskalMD<article-title>Interaction of solitons in a collisionless plasma and the recurrence of initial statesPhys. Rev. Lett.19651562401965PhRvL..15..240Z10.1103/PhysRevLett.15.240
  3.LuDSeadawyARWangJArshadMFarooqU<article-title>Soliton solutions of the generalised third-order nonlinear Schrödinger equation by two mathematical methods and their stabilityPramana2019933442019Prama..93...44L10.1007/s12043-019-1804-54.BeheraS<article-title>Analysis of traveling wave solutions of two space-time nonlinear fractional differential equations by the first-integral methodMod. Phys. Lett. B2024380423502472024MPLB...3850247B468892510.1142/S0217984923502470
  5.BeheraSAljahdalyNH<article-title>Soliton solutions of nonlinear geophysical KdV equation via two analytical methodsInt. J. Theor. Phys.2024635107473991010.1007/s10773-024-05647-2
  6.BeheraS<article-title>Optical solitons for the Hirota–Ramani equation via improved -expansion methodMod. Phys. Lett. B2025390124504032025MPLB...3950403B483533110.1142/S02179849245040377.BeheraS<article-title>Multiple soliton solutions of some conformable fractional nonlinear models using Sine-Cosine methodOpt. Quantum Electron.2024567123510.1007/s11082-024-06403-w
  8.BeheraS<article-title>Computational and numerical analysis of the soliton solutions to the geophysical KdV equation using two robust analytical methodsPramana20259931082025Prama..99..108B10.1007/s12043-025-02971-y9.AlquranM<article-title>Investigating fluctuation varieties in the propagation of the perturbed KdV equation with time-dependent perturbation coefficientPartial Differ. Equ. Appl. Math.matics20251410120610.AlquranMAlhamiR<article-title>Analysis of lumps, single-stripe, breather-wave, and two-wave solutions to the generalized perturbed-KdV equation by means of Hirota’s bilinear methodNonlinear Dyn.202210931985199210.1007/s11071-022-07509-0
  11.AlquranM<article-title>Physical properties for bidirectional wave solutions to a generalized fifth-order equation with third-order time-dispersion termResults Phys.20212810457710.1016/j.rinp.2021.104577
12.AlquranM<article-title>Variation of the influence of Atangana-conformable time-derivative on various physical structures in the fractional KP-BBM modelInt. J. Theor. Phys.2024639225479459010.1007/s10773-024-05767-9

13.AlquranM<article-title>Necessary conditions for convex-periodic, elliptic-periodic, inclined-periodic, and rogue wave-solutions to exist for the multi-dispersions Schrödinger equationPhys. Scr.20249920252482024PhyS...99b5248A10.1088/1402-4896/ad1fba14.AlquranMAl-deiakehR<article-title>Lie-Backlund symmetry generators and a variety of novel periodic-soliton solutions to the complex-mode of modified Korteweg-de Vries equationQual. Theory Dyn. Syst.202423295470041710.1007/s12346-023-00953-215.JaradatIAlquranMAliMAl-deiakehR<article-title>Modeling synchronized propagation of two symmetric waves in a new two-mode extension of the -dimensional Chaffee-Infante modelInt. J. Theor. Phys.202564245486112410.1007/s10773-025-05914-w16.AlquranM<article-title>New interesting optical solutions to the quadratic–cubic Schrodinger equation by using the Kudryashov-expansion method and the updated rational sine–cosine functionsOpt. Quantum Electron.2022541066610.1007/s11082-022-04070-317.PanJRahmanMURafiullahU<article-title>Breather-like, singular, periodic, interaction of singular and periodic solitons, and a-periodic solitons of third-order nonlinear Schrödinger equation with an efficient algorithmEur. Phys. J. Plus2023138109122023EPJP..138..912P10.1140/epjp/s13360-023-04530-z18.NasreenNSeadawyARLuDArshadM<article-title>Optical fibers to model pulses of ultrashort via generalized third-order nonlinear Schrödinger equation by using extended and modified rational expansion methodJ. Nonlinear Opt. Phys. Mater.2024330423500582024JNOPM..3350058N10.1142/S021886352350058319.AhmedAIAlgolamMSCesaranoCRizkDGassemFMohammedWW<article-title>Dynamical behavior of the fractional generalized nonlinear Schrödinger equation of third-orderOpt. Quantum Electron.202456584310.1007/s11082-024-06626-x20.BadshahFTariqKUIncMKazmiSR<article-title>Solitons, stability analysis and modulation instability for the third order generalized nonlinear Schrödinger model in ultraspeed fibersOpt. Quantum Electron.20235512109410.1007/s11082-023-05411-621.FaridiWATipuGHRiazMBMostafaAMAlqahtaniSAMyrzakulovRUmurzakhovaZ<article-title>Analyzing optical soliton solutions in Kairat-X equation via new auxiliary equation methodOpt. Quantum Electron.2024568131710.1007/s11082-024-07197-722.KazmiSSJhangeerARazaNAlrebdiHIAbdel-AtyAHEleuchH<article-title>The analysis of bifurcation, quasi-periodic and solitons patterns to the new form of the generalized q-deformed Sinh-Gordon equationSymmetry202315713242023Symm...15.1324K10.3390/sym1507132423.BeheraSMohantySVirdiJPS<article-title>Analytical solutions and mathematical simulation of traveling wave solutions to fractional order nonlinear equationsPartial Differ. Equ. Appl. Math.2023810053510.1016/j.padiff.2023.10053524.BeheraS<article-title>Dynamical solutions and quadratic resonance of nonlinear perturbed Schrödinger equationFront. Appl. Math. Stat.2023810867662023maum.book.....B10.3389/fams.2022.108676625.BeheraSAljahdalyNH<article-title>Nonlinear evolution equations and their traveling wave solutions in fluid media by modified analytical methodPramana20239731302023Prama..97..130B10.1007/s12043-023-02602-426.SeadawyARArshadMLuD<article-title>The weakly nonlinear wave propagation of the generalized third-order nonlinear Schrödinger equation and its applicationsWaves Random Complex Media20223228198312022WRCM...32..819S439369510.1080/17455030.2020.180208527.MalikSKumarSNisarKSSaleelCA<article-title>Different analytical approaches for finding novel optical solitons with generalized third-order nonlinear Schrödinger equationResults Phys.20212910475510.1016/j.rinp.2021.10475528.LiZLyuJHussainE<article-title>Bifurcation, chaotic behaviors and solitary wave solutions for the fractional Twin-Core couplers with Kerr law non-linearitySci. Rep.2024141226162024NatSR..1422616L10.1038/s41598-024-74044-w393498161144305529.LiZZhaoS<article-title>Bifurcation, chaotic behavior and solitary wave solutions for the Akbota equationAIMS Math.2024982259022601477620610.3934/math.20241100

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Helal, MA; Seadawy, AR. Benjamin-Feir instability in nonlinear dispersive waves. Comput. Math. Appl.; 2012; 64, 11 pp. 3557-3568.2992534 [DOI: https://dx.doi.org/10.1016/j.camwa.2012.09.006]

2. Zabusky, NJ; Kruskal, MD. Interaction of solitons in a collisionless plasma and the recurrence of initial states. Phys. Rev. Lett.; 1965; 15, 6 240.1965PhRvL.15.240Z [DOI: https://dx.doi.org/10.1103/PhysRevLett.15.240]

3. Lu, D; Seadawy, AR; Wang, J; Arshad, M; Farooq, U. Soliton solutions of the generalised third-order nonlinear Schrödinger equation by two mathematical methods and their stability. Pramana; 2019; 93, 3 44.2019Prama.93..44L [DOI: https://dx.doi.org/10.1007/s12043-019-1804-5]

4. Behera, S. Analysis of traveling wave solutions of two space-time nonlinear fractional differential equations by the first-integral method. Mod. Phys. Lett. B; 2024; 38, 04 2350247.2024MPLB..3850247B4688925 [DOI: https://dx.doi.org/10.1142/S0217984923502470]

5. Behera, S; Aljahdaly, NH. Soliton solutions of nonlinear geophysical KdV equation via two analytical methods. Int. J. Theor. Phys.; 2024; 63, 5 107.4739910 [DOI: https://dx.doi.org/10.1007/s10773-024-05647-2]

6. Behera, S. <article-title>Optical solitons for the Hirota–Ramani equation via improved -expansion method

Mod. Phys. Lett. B; 2025; 39, 01 2450403.2025MPLB..3950403B4835331 [DOI: https://dx.doi.org/10.1142/S0217984924504037]

16. Alquran, M. New interesting optical solutions to the quadratic–cubic Schrodinger equation by using the Kudryashov-expansion method and the updated rational sine–cosine functions. Opt. Quantum Electron.; 2022; 54, 10 666. [DOI: https://dx.doi.org/10.1007/s11082-022-04070-3]

17. Pan, J; Rahman, MU; Rafiullah, U. Breather-like, singular, periodic, interaction of singular and periodic solitons, and a-periodic solitons of third-order nonlinear Schrödinger equation with an efficient algorithm. Eur. Phys. J. Plus; 2023; 138, 10 912.2023EPJP.138.912P [DOI: https://dx.doi.org/10.1140/epjp/s13360-023-04530-z]

18. Nasreen, N; Seadawy, AR; Lu, D; Arshad, M. Optical fibers to model pulses of ultrashort via generalized third-order nonlinear Schrödinger equation by using extended and modified rational expansion method. J. Nonlinear Opt. Phys. Mater.; 2024; 33, 04 2350058.2024JNOPM.3350058N [DOI: https://dx.doi.org/10.1142/S0218863523500583]

19. Ahmed, AI; Algolam, MS; Cesarano, C; Rizk, D; Gassem, F; Mohammed, WW. Dynamical behavior of the fractional generalized nonlinear Schrödinger equation of third-order. Opt. Quantum Electron.; 2024; 56, 5 843. [DOI: https://dx.doi.org/10.1007/s11082-024-06626-x]

20. Badshah, F; Tariq, KU; Inc, M; Kazmi, SR. Solitons, stability analysis and modulation instability for the third order generalized nonlinear Schrödinger model in ultraspeed fibers. Opt. Quantum Electron.; 2023; 55, 12 1094. [DOI: https://dx.doi.org/10.1007/s11082-023-05411-6]

21. Faridi, WA; Tipu, GH; Riaz, MB; Mostafa, AM; Alqahtani, SA; Myrzakulov, R; Umurzakhova, Z. Analyzing optical soliton solutions in Kairat-X equation via new auxiliary equation method. Opt. Quantum Electron.; 2024; 56, 8 1317. [DOI: https://dx.doi.org/10.1007/s11082-024-07197-7]

22. Kazmi, SS; Jhangeer, A; Raza, N; Alrebdi, HI; Abdel-Aty, AH; Eleuch, H. The analysis of bifurcation, quasi-periodic and solitons patterns to the new form of the generalized q-deformed Sinh-Gordon equation. Symmetry; 2023; 15, 7 1324.2023Symm..15.1324K [DOI: https://dx.doi.org/10.3390/sym15071324]

23. Behera, S; Mohanty, S; Virdi, JPS. Analytical solutions and mathematical simulation of traveling wave solutions to fractional order nonlinear equations. Partial Differ. Equ. Appl. Math.; 2023; 8, 100535. [DOI: https://dx.doi.org/10.1016/j.padiff.2023.100535]

24. Behera, S. Dynamical solutions and quadratic resonance of nonlinear perturbed Schrödinger equation. Front. Appl. Math. Stat.; 2023; 8, 1086766.2023maum.book...B [DOI: https://dx.doi.org/10.3389/fams.2022.1086766]

25. Behera, S; Aljahdaly, NH. Nonlinear evolution equations and their traveling wave solutions in fluid media by modified analytical method. Pramana; 2023; 97, 3 130.2023Prama.97.130B [DOI: https://dx.doi.org/10.1007/s12043-023-02602-4]

26. Seadawy, AR; Arshad, M; Lu, D. The weakly nonlinear wave propagation of the generalized third-order nonlinear Schrödinger equation and its applications. Waves Random Complex Media; 2022; 32, 2 pp. 819-831.2022WRCM..32.819S4393695 [DOI: https://dx.doi.org/10.1080/17455030.2020.1802085]

27. Malik, S; Kumar, S; Nisar, KS; Saleel, CA. Different analytical approaches for finding novel optical solitons with generalized third-order nonlinear Schrödinger equation. Results Phys.; 2021; 29, 104755. [DOI: https://dx.doi.org/10.1016/j.rinp.2021.104755]

28. Li, Z; Lyu, J; Hussain, E. Bifurcation, chaotic behaviors and solitary wave solutions for the fractional Twin-Core couplers with Kerr law non-linearity. Sci. Rep.; 2024; 14, 1 22616.2024NatSR.1422616L [DOI: https://dx.doi.org/10.1038/s41598-024-74044-w] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/39349816][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11443055]

29. Li, Z; Zhao, S. Bifurcation, chaotic behavior and solitary wave solutions for the Akbota equation. AIMS Math.; 2024; 9, 8 pp. 22590-22601.4776206 [DOI: https://dx.doi.org/10.3934/math.20241100]