Home Insights from the nonlinear Schrödinger–Hirota equation with chromatic dispersion: Dynamics in fiber–optic communication
Article Open Access

Insights from the nonlinear Schrödinger–Hirota equation with chromatic dispersion: Dynamics in fiber–optic communication

  • Dean Chou , Hamood Ur Rehman , Muhammad Imran Asjad EMAIL logo and Ifrah Iqbal
Published/Copyright: April 29, 2025
Become an author with De Gruyter Brill

Abstract

Optical solitons have practical applications in communication systems as carriers of optical information. An advantage of an optical soliton is its ability to maintain its structure unchanged when interacting with other pulses. The purpose of this article is to strive for the optical soliton solutions of the nonlinear Schrödinger–Hirota equation incorporating chromatic dispersion, which is a governing model for understanding the dynamics of dispersive pulse propagation in optical fibers. Three distinct methodologies, namely, the 1 φ ( η ) , φ ( η ) φ ( η ) method, the new Kudryashov method, and the extended simple equation method, have been employed to derive solutions for this equation. Numerous optical soliton solutions including singular soliton, dark soliton, bright soliton, periodic singular soliton, dark–bright soliton, and dark-singular solution and singular periodic solutions are offered by this method. These techniques present an accurate and successful strategy for deriving exact solutions to various nonlinear partial differential equations. Graphical representations of some of the extracted solutions are shown to demonstrate the behavior of solitons. We also presented the eye diagram for the depiction of the signal to provide a clear visualization of the pulse shapes and their characteristics. The optical solitons produced in respect to this form have never been explored by the proposed techniques before, and the results have never been published. In this study, we perform a comprehensive perturbation analysis to explore the stability and dynamics of soliton solutions within the Schrödinger–Hirota equation.

1 Introduction

The exploration of optical soliton solutions to nonlinear partial differential equations (NPDEs) has emerged as a captivating and pivotal area of investigation within the realms of applied sciences and engineering [114]. Solitons, possessing numerous applications in science and engineering play a crucial role in various technological advancements. A paramount application of solitons lies in their ability to transmit digital information through optical fibers. The utilization of solitons in conveying information has become a cornerstone in the technological landscape [1522]. Additionally, solitons are subject to intensive examination in electromagnetics, serving as transverse electromagnetic waves passing between two strips of superconducting metal. In the field of nonlinear optics, optical solitons constitute a remarkable field of study including optical fibers, magneto-optics, optical couplers, crystals, meta-surfaces, and birefringent materials [2325]. The detailed study of optical solitons adds depth to our understanding of nonlinear optical phenomena and contributes in various optical technologies.

Recently, a significant focus in scientific research has been directed toward exploring mathematical models that elucidate the propagation of pulses in optical media. This heightened interest is primarily attributed to the rapid advancements in information transmission systems, especially over long distances. Among the widely employed nonlinear equations for characterizing pulse propagation in optical media is the nonlinear Schrödinger-Hirota equation (NLSHE) [2628]. The NLSE is an NPDE that has application in diverse physical scenarios incluiding fluid mechanics, hydrodynamic, and the depiction of surface gravity water wave evolution. Moreover, it serves as a fundamental model in various fields, including nonlinear optics, fluid dynamics, nuclear physics, mathematical finance, biochemistry, plasma physics, and superconductivity. Specifically, it describes phenomena such as solitary wave propagation in heat pulse propagation in solids, piezoelectric semiconductors in condensed matter physics, and more. Another variant that has garnered attention in several studies is the NLSHE. This equation accounts for third-order dispersion (TOD) and finds application in the analysis of pulse propagation processes within optical media [29].

Chromatic dispersion (CD) and nonlinear refractive index (NRI) represent essential components in the foundational principles governing soliton propagation within optical fibers [30]. CD arises due to the dependency of the phase and group velocity of light in a transparent medium on optical frequency. This dependence leads to variations in wavelength and propagation velocity [31]. While the impact of CD in fibers can be mitigated, it is an inherent phenomenon in the production of glass fibers. Managing CD becomes imperative for optimizing and minimizing the effects of nonlinear terms in wavelength-division multiplexing systems. Despite efforts to reset the effects of CD, achieving a perfect balance remains a significant challenge. CD introduces issues in fibers, manifesting as degraded signal quality, necessitating dispersion compensation. A highly effective approach to address CD concerns is the utilization of fiber gratings [32]. The NLSHE with CD is given as follows [33,34]:

(1) ι y t + a 1 y x x + b 1 y 2 y + ι ( c 1 y x x x + d 1 y 2 y x ) = 0 ,

where a 1 , b 1 , c 1 , and d 1 show the coefficients associated with group velocity dispersion (GVD), Kerr-law nonlinearity, TOD, and nonlinear dispersion terms, respectively.

The NLSHE with CD, as presented in our model, is applicable to dispersion-flattened fibers and waveguides exhibiting polarization-mode dispersion (PMD). The addition of Kerr nonlinearity, GVD, TOD, and nonlinear dispersion terms accurately investigates the behavior of optical pulses in these systems. Dispersion-flattened fibers are particularly relevant here due to the need for managing both second- and third-order dispersive effects, while the nonlinear dispersion term addresses the impact of PMD in birefringent fibers. This model can be extended to grating structures, and its primary applicability is in fiber optic communication systems where higher-order nonlinearities and dispersion significantly affects the soliton dynamics. The fibers considered in our model are step-index fibers. This choice is due to their distinct core and cladding refractive index profiles, which allow us to analyze the nonlinear dynamics and higher-order dispersive effects, such as GVD and TOD, more clearly. Step-index fibers, particularly in terms of dispersion-flattened designs, give a straightforward framework for studying soliton propagation, PMD, and other nonlinear effects, which are central to the model.

Numerous methodologies have been developed to extract soliton solutions from NPDEs, including, but not limited to, the following: Kudryashov’s method [35,36], extended hyperbolic function method [3741], Bernoulli sub-ODE method [42,43], Sardar sub-equation method [44,45], newly extended direct algebraic method [46,47], Jacobi’s elliptic function [48,49], mapping method [5053], and ϕ 6 -expansion method [54]. In this study, we have adopted three distinct approaches, namely, the 1 φ ( η ) , φ ( η ) φ ( η ) method [21], the new Kudryashov method [55], and the extended simple equation method [56], for the first time. Through these novel approaches, we have successfully extracted diverse solutions, including dark, bright, periodic singular, singular, combined dark–bright, and combined dark-singular optical solitons.

The remaining structure of this article is as follows: Section 2 gives mathematical analysis of the equation. Sections 3, 4, and 5 give the descriptions of applied methods with applications. The perturbation analysis for slowdown effect is given in Section 6. The discussion of this article is given in Section 7. Section 8 is of conclusion.

2 Mathematical analysis

Suppose

(2) y ( x , t ) = Y ( η ) e ι ξ , η = x v t , and ξ = α x + β t + ψ ,

where v , α , β , and ψ represent the velocity, frequency, wave number, and phase constant, respectively. By substituting (2) into (1), the real and imaginary parts are obtained as follows:

(3) ( α d 1 + b 1 ) Y 3 + ( α 3 c 1 a 1 α 2 β ) Y + ( 3 α c 1 + a 1 ) Y = 0 ,

(4) d 1 Y 2 Y + ( 3 α 2 c 1 2 a 1 α v ) Y + c 1 Y = 0 .

By integrating (4) and let integration of constant equals to zero

(5) d 1 Y 3 + ( 9 α 2 c 1 6 a 1 α 3 v ) Y + 3 c 1 Y = 0 .

Now, from (3) and (5), we obtained the following values of constants:

(6) d 1 = 3 b 1 c 1 a 1 , β = 8 α 3 c 1 2 + 8 a 1 α c 1 + 2 a 1 2 α + 3 α c 1 v + a 1 v c 1 .

Under these conditions, consider (3) or (5) as nonlinear ordinary differential equation (ODE).

3 Description of 1 φ ( η ) , φ ( η ) φ ( η ) method

Let the solution of (3) be

(7) Y ( η ) = w 0 + i = 1 m w i + v i φ ( η ) i φ ( η ) i ,

where w 0 , w i , and v i   ( i = 1 , 2 , , m ) are the constants. m can be obtained using the homogeneous balancing rule and φ ( η ) represents the following ODE:

(8) φ ( η ) 2 = φ ( η ) 2 ϱ ,

where

(9) φ ( η ) = c e η + ϱ 4 c e η .

Now, by substituting (7) along (8) into (3), the system of equations is attained, and by solving it, we obtain the values of constants.

3.1 Application of 1 φ ( η ) , φ ( η ) φ ( η ) method

Now, using the homogeneous balance rule, we obtain m = 1 .

(10) Y ( η ) = w 0 + w 1 + v 1 φ ( η ) φ ( η ) .

Now, by substituting Eq. (10) along Eq. (8) and Eq. (9) into (3) , the system of equations is attained, and by solving it, we obtain the following values of constants:

Set1:

w 0 = 0 , w 1 = v 1 ϱ ( 6 a 2 + 4 α 2 + 15 ) 6 a 2 4 α 2 + 3 , c 1 = 2 ( 2 a 2 β + 3 β ) α ( 2 α 2 + 3 ) , b 1 = 4 α 2 β 6 a 2 β 3 β 2 α 3 d 1 v 1 2 3 α d 1 v 1 2 ( 2 α 2 + 3 ) v 1 2 , a 1 = 2 β .

By substituting these values into Eq. (10), we have

y 1 ( x , t ) = v 1 4 c e η ϱ ( 6 a 2 + 4 α 2 + 15 ) 6 a 2 4 α 2 + 3 + c e η ϱ 4 c 2 e 2 η + ϱ e ι ξ .

By taking ϱ = ± 4 c 2 , we obtain

(11) y 1,1 ( x , t ) = v 1 ( 6 a 2 + 4 α 2 + 15 ) sech ( η ) 6 a 2 4 α 2 + 3 + tanh ( η ) e ι ξ .

(12) y 1,2 ( x , t ) = v 1 6 a 2 + 4 α 2 + 15 csch ( η ) 6 a 2 4 α 2 + 3 + coth ( η ) e ι ξ .

Set2:

v 1 = 0 , d 1 = a 1 b 1 3 ( a 1 a 2 + β ) , α = 3 ( a 1 a 2 + β ) a 1 , c 1 = a 1 3 2 3 3 ( a 1 a 2 + β ) .

By substituting these values into Eq. (10), we have

y 1 , * ( x , t ) = 4 c w 1 e η 4 c 2 e 2 η + ϱ + w 0 e ι ξ .

By taking ϱ = ± 4 c 2 , we obtain

(13) y 1,3 ( x , t ) = w 1 sech ( η ) 2 c + w 0 e ι ξ ,

(14) y 1,4 ( x , t ) = w 1 csch ( η ) 2 c + w 0 e ι ξ .

Set3:

w 0 = 0 , w 1 = 0 , c 1 = a 1 a 2 + 2 a 1 + β α ( α 2 + 6 ) , d 1 = 6 a 1 a 2 + 2 α 2 a 1 α 2 b 1 v 1 2 6 b 1 v 1 2 6 β ( α 3 + 6 α ) v 1 2 .

By substituting these values into Eq. (10), we have

y 1 , * * ( x , t ) = v 1 ( 4 d 2 e 2 η ϱ ) 4 c 2 e 2 η + ϱ e ι ξ .

By taking ϱ = ± 4 c 2 , we obtain

(15) y 1,5 ( x , t ) = ( v 1 tanh ( η ) ) e ι ξ ,

(16) y 1,6 ( x , t ) = ( v 1 coth ( η ) ) e ι ξ .

4 New Kudryashov’s method

Let the solution of (3) be

(17) Y ( η ) = i = 1 M h i χ ( η ) ,

where h i   ( i = 1 , 2 , , M ) are the constants. M can be obtained by using the homogeneous balancing rule, and χ ( η ) represents the following ODE:

(18) χ ( η ) 2 = δ 2 χ ( η ) 2 ( 1 λ χ ( η ) 2 ) ,

where

(19) χ ( η ) = 4 R 4 e δ η R 2 + e δ η λ .

Now, by substituting (17) along (18) into (3), the system of equations is attained, and by solving it, we obtain the values of constants.

4.1 Application of new Kudryashov’s method

As the balancing number M = 1 , so from (17), we have

(20) Y ( η ) = h 0 + h 1 χ ( η ) .

By substituting (20) along (18) into (3), the system of equations is attained, and by solving it, we obtain the values of constants:

h 0 = 0 , a 1 = β + α 3 ( c 1 ) + 3 α c 1 δ 2 a 2 δ 2 , b 1 = a 2 α h 1 2 d 1 + 6 a 2 α c 1 δ 2 λ 2 β δ 2 λ + α h 1 2 δ 2 d 1 2 α 3 c 1 δ 2 λ h 1 2 ( a 2 δ 2 ) .

By substituting these values of constants into (20), we have the following solutions:

(21) y 2 ( x , t ) = 4 h 1 R 4 R 2 e δ η + λ e δ ( η ) e ι ξ .

By taking λ = ± 4 R 2 , we have the following solutions:

(22) y 2,1 ( x , t ) = h 1 sech ( δ η ) 2 R e ι ξ ,

(23) y 2,2 ( x , t ) = h 1 csch ( δ η ) 2 R e ι ξ .

5 Extended simple equation method

Consider the solution of (3) is

(24) Y ( η ) = i = M M g i F ( η ) ,

where g i ( i = 1 , 2 , , M ) are the constants. M can be obtained by using the homogeneous balancing rule, and F ( η ) represents the following ODE:

(25) F ( η ) = p 3 F ( η ) 3 + p 2 F ( η ) 2 + p 1 F ( η ) + p 0 ,

where p i , i = 0 , 1, 2, 3 are the arbitrary constants and (25) have the following solutions:

Case 1 When p 0 = p 3 = 0 ,

F 1 ( η ) = p 1 exp ( η p 1 ) 1 p 2 exp ( η p 1 ) , p 1 > 0 , F 2 ( η ) = p 1 exp ( η p 1 ) 1 p 2 exp ( η p 1 ) , p 1 < 0 .

Case 2 When p 1 = p 3 = 0 ,

F 3 ( η ) = p 0 p 2 tan ( p 0 p 2 η ) p 2 , p 0 p 2 > 0 , F 4 ( η ) = p 0 p 2 tanh ( p 0 p 2 η ) p 2 , p 0 p 2 < 0 .

Case 3 When p 0 = p 2 = 0 ,

F 5 ( η ) = p 1 exp ( η p 1 ) 1 exp ( 2 η p 1 p 3 ) , p 1 > 0 , F 6 ( η ) = p 1 exp ( η p 1 ) exp ( 2 η p 1 p 3 ) , p 1 < 0 .

5.1 Application of the extended simple equation method

As homogeneous balancing rule gives M = 1 , so we acquired from (24):

(26) Y ( η ) = g 0 + g 1 F ( η ) + g 1 F ( η ) 1 .

By substituting (26) and (25) into (3), we obtain following sets of solutions:

Set 1 a 1 = 3 α c 1 , b 1 = α d 1 , β = α c 1 ( α 2 3 a 2 ) .

Using these values along solutions of (25), we obtained

Case 1 When p 0 = p 3 = 0 ,

y 3,1 ( η ) = g 1 ( e η p 1 p 2 ) p 1 + g 1 p 1 e η p 1 1 p 2 e η p 1 + g 0 e ι ξ , y 3,2 ( η ) = g 1 ( e η p 1 + p 2 ) p 1 + g 1 p 1 e η p 1 1 + p 2 e η p 1 + g 0 e ι ξ .

Case 2 When p 1 = p 3 = 0 ,

y 3,3 ( η ) = g 0 g 1 p 2 coth ( η p 0 p 2 ) + g 1 p 0 tanh ( η p 0 p 2 ) p 0 p 2 e ι ξ , y 3,4 ( η ) = g 0 + g 1 p 2 cot ( η p 0 p 2 ) + g 1 p 0 tan ( η p 0 p 2 ) p 0 p 2 e ι ξ .

Case 3 When p 0 = p 2 = 0 ,

y 3,5 ( η ) = g 1 1 e 2 η p 1 p 3 p 1 e η p 1 + g 1 p 1 e η p 1 1 e 2 η p 1 p 3 + g 0 e ι ξ , y 3,6 ( η ) = g 1 e 2 η p 1 p 3 p 1 ( e η p 1 ) + g 1 p 1 ( e η p 1 ) e 2 η p 1 p 3 + g 0 e ι ξ .

Set 2 g 1 = 0 , b 1 = α d 1 , a 1 = 3 α c 1 , β = α c 1 ( α 2 3 a 2 ) .

Using these values along solutions of (25), we obtained

Case 1 When p 0 = p 3 = 0 ,

y 3,7 ( η ) = g 1 e η p 1 g 1 p 2 + p 1 g 0 p 1 e ι ξ , y 3,8 ( η ) = g 1 e η p 1 + g 1 p 2 + p 1 g 0 p 1 e ι ξ .

Case 2 When p 1 = p 3 = 0 ,

y 3,9 ( η ) = g 0 + g 1 p 2 cot ( η p 0 p 2 ) p 0 p 2 e ι ξ , y 3,10 ( η ) = g 0 + g 1 p 2 coth ( η p 0 p 2 ) p 0 p 2 e ι ξ .

Case 3 When p 0 = p 2 = 0 ,

y 3,11 ( η ) = g 1 1 e 2 η p 1 p 3 p 1 e η p 1 + g 0 e ι ξ , y 3,12 ( η ) = g 1 e 2 η p 1 p 3 p 1 ( e η p 1 ) + g 0 e ι ξ .

Set 3 g 1 = 0 , b 1 = α d 1 , a 1 = 3 α c 1 , β = α c 1 ( α 2 3 a 2 ) .

Using these values along solutions of (25), we obtained

Case 1 When p 0 = p 3 = 0 ,

y 3,13 ( x , t ) = g 0 + g 1 p 1 exp ( η p 1 ) 1 p 2 exp ( η p 1 ) e ι ξ , y 3,14 ( x , t ) = g 0 + g 1 p 1 exp ( η p 1 ) 1 p 2 exp ( η p 1 ) e ι ξ .

Case 2 When p 1 = p 3 = 0 ,

y 3,15 ( x , t ) = g 0 + g 1 p 0 p 2 tan ( p 0 p 2 η ) p 2 e ι ξ , y 3,16 ( x , t ) = g 0 + g 1 p 0 p 2 tanh ( p 0 p 2 η ) p 2 e ι ξ .

Case 3 When p 0 = p 2 = 0 ,

y 3,17 ( x , t ) = g 0 + g 1 p 1 exp ( η p 1 ) 1 exp ( 2 η p 1 p 3 ) e ι ξ , y 3,18 ( x , t ) = g 0 + g 1 p 1 exp ( η p 1 ) exp ( 2 η p 1 p 3 ) e ι ξ .

6 Perturbative analysis

In this section, we extend our analysis of the soliton solution by incorporating the effects of TOD. The presence of TOD is known to induce both a slowdown in the soliton’s velocity and the emission of radiation in the form of small-amplitude dispersive waves. To obtain these physical effects, the soliton solution is perturbed by adding a small parameter, ε , which accounts for the deviations caused by TOD [57].

We begin with the unperturbed soliton solution, which is given by

(27) y 1,1 ( x , t ) = h 1 csch ( δ η ) 2 R e i ξ , y 2,1 ( x , t ) = h 1 sech ( δ η ) 2 R e i ξ ,

where η = x v t shows the soliton’s core moving at velocity v , h 1 presents the soliton amplitude, and R is a scaling parameter. The phase factor ξ = α x β t describes the wave number α and frequency β of the soliton. The sech ( δ η ) proposes the soliton’s localized shape, with δ controlling the soliton width.

To incorporate the effects of TOD, a small perturbation parameter ε is introduced to modify both the soliton’s velocity (slowdown effect) and induced radiation. The perturbed soliton solution is written as:

(28) y 1,1 ( x , t ) = h 1 csch ( δ η slowdown ) 2 R e i ξ slowdown + ε A e γ t e i ( κ x ω r t ) , y 2,1 ( x , t ) = h 1 sech ( δ η slowdown ) 2 R e i ξ slowdown + ε A e γ t e i ( κ x ω r t ) ,

where η slowdown = x ( v + ε t ) t introduces a time-dependent velocity correction, representing the soliton slowdown due to TOD. Similarly, the phase factor is modified as:

(29) ξ slowdown = α x ( β + ε t ) t .

The second term in (28), ε A e γ t e i ( κ x ω r t ) , represents the radiative component, where A is the amplitude of the emitted radiation, γ is a decay constant that controls how quickly the radiation diminishes over time, κ is the wavenumber of the radiative wave, and ω r is the frequency of the radiative wave. This term models the small-amplitude dispersive waves that are emitted by the soliton as it propagates and gradually loses energy due to the perturbation. The introduction of the perturbation term has two key effects: the soliton experiences a time-dependent reduction in its velocity due to TOD. This is formed by the perturbative shift in the soliton’s velocity and is given by ( v + ε t ) , where ε shows the magnitude of the slowdown. As time progresses, the soliton decelerates, as evidenced by the shift in the position of the soliton core. The perturbation also triggers the emission of radiation in the form of small dispersive waves. These waves, described by the oscillatory term ε A e γ t e i ( κ x ω r t ) , gradually radiate energy away from the soliton. The amplitude of these radiative waves decays over time due to the factor e γ t , and their wavenumber κ and frequency ω r characterize their oscillatory behavior. To visualize the effects of the perturbation, we numerically solved (28) for both the unperturbed and perturbed soliton solutions. The results are shown in Figure 8.

Figure 1 
               Plots of 
                     
                        
                        
                           
                              
                                 y
                              
                              
                                 1,1
                              
                           
                           
                              (
                              
                                 x
                                 ,
                                 t
                              
                              )
                           
                        
                        {y}_{\mathrm{1,1}}\left(x,t)
                     
                   by taking 
                     
                        
                        
                           α
                           =
                           1.5
                        
                        \alpha =1.5
                     
                  , 
                     
                        
                        
                           a
                           =
                           0.5
                           ,
                           β
                           =
                           1.3
                        
                        a=0.5,\beta =1.3
                     
                  , 
                     
                        
                        
                           
                              
                                 v
                              
                              
                                 1
                              
                           
                           =
                           0.5
                        
                        {v}_{1}=0.5
                     
                  , 
                     
                        
                        
                           v
                           =
                           1.5
                        
                        v=1.5
                     
                   and 
                     
                        
                        
                           ψ
                           =
                           0.5
                        
                        \psi =0.5
                     
                  .
Figure 1

Plots of y 1,1 ( x , t ) by taking α = 1.5 , a = 0.5 , β = 1.3 , v 1 = 0.5 , v = 1.5 and ψ = 0.5 .

In the study by Kaur et al. Figure 8, we observe that the perturbed soliton (red dashed line) shifts to the left compared to the original soliton (blue solid line), indicating the deceleration caused by the perturbation. Furthermore, small oscillations trailing the soliton illustrate the radiative waves that are emitted due to the perturbation.

7 Results and discussion

In the study by Kaur et al. [58], the Schrödinger–Hirota equation with variable coefficients and power-law nonlinearity was investigated, focusing on bright and dark optical solitons. Similarly, bright, dark, and singular soliton solutions of the perturbed Schrödinger–Hirota equation, incorporating spatio-temporal dispersion and Kerr nonlinearity, were developed in the study by Yildirim [59]. New soliton solutions of the time-fractional Schrödinger–Hirota equation were explored in the study by Zayed et al. [60], where various dispersive optical soliton solutions were derived. The dynamic behaviors and optical solitons within DWDM networks, modeled by the Schrödinger–Hirota equation, were examined in the study by Tang [61] using traveling wave reductions to a plane dynamical system. The existing literature, as referenced, establishes that while numerous aspects of the Schrödinger–Hirota equation have been explored, our work presents a novel contribution, further advancing the field by addressing gaps in previous studies and introducing new solutions within this framework.

This study delves into the exploration of soliton solutions for the NLSHE through the application of three distinct approaches. The obtained solutions encompass various types, providing a comprehensive understanding of the equation’s behavior. In this section, we present graphical representations of these solutions through 3D and 2D plots, showcasing absolute, real, and imaginary aspects. These visualizations serve as crucial tools for comprehending the intricate dynamics of waves and play a vital role in validating the accuracy of the obtained solutions. The eye diagram is presented for the signal, providing a clear visualization of the pulse shapes and their characteristics. The graphical representations delineate diverse wave patterns influenced by constants, with potential applications in marine engineering for the design of resilient structures and the control of wave propagation. The details of each figure are elucidated as follows:

Figure 1(a)–(c) presents the absolute value graph of y 1,1 ( x , t ) , revealing a distinctive combined dark–bright optical soliton. In Figure 1(d)–(f), the real value graphs of y 1,1 ( x , t ) illustrate the periodic patterns. Figure 1(g)–(i) portrays the imaginary value graphs of y 1,1 ( x , t ) , exhibiting periodic characteristics similar to the real value graphs but with differences in amplitude. These graphical representations correspond to the specified constants: α = 1.5 , a = 0.5 , β = 1.3 , v 1 = 0.5 , v = 1.5 and ψ = 0.5 . Figure 2(a)–(c) showcases the absolute value graph of y 1,4 ( x , t ) , unveiling a distinctive singular optical soliton.

Figure 2 
               Plots of 
                     
                        
                        
                           
                              
                                 y
                              
                              
                                 1,4
                              
                           
                           
                              (
                              
                                 x
                                 ,
                                 t
                              
                              )
                           
                        
                        {y}_{\mathrm{1,4}}\left(x,t)
                     
                   by taking 
                     
                        
                        
                           α
                           =
                           1.5
                        
                        \alpha =1.5
                     
                  , 
                     
                        
                        
                           a
                           =
                           0.5
                           ,
                           β
                           =
                           1.3
                        
                        a=0.5,\beta =1.3
                     
                  , 
                     
                        
                        
                           
                              
                                 w
                              
                              
                                 1
                              
                           
                           =
                           0.5
                        
                        {w}_{1}=0.5
                     
                  , 
                     
                        
                        
                           v
                           =
                           1.5
                        
                        v=1.5
                     
                  , 
                     
                        
                        
                           
                              
                                 w
                              
                              
                                 0
                              
                           
                           =
                           0.5
                        
                        {w}_{0}=0.5
                     
                  , 
                     
                        
                        
                           c
                           =
                           1
                        
                        c=1
                     
                  , and 
                     
                        
                        
                           ψ
                           =
                           0.5
                        
                        \psi =0.5
                     
                  .
Figure 2

Plots of y 1,4 ( x , t ) by taking α = 1.5 , a = 0.5 , β = 1.3 , w 1 = 0.5 , v = 1.5 , w 0 = 0.5 , c = 1 , and ψ = 0.5 .

Figure 3 
               Plots of 
                     
                        
                        
                           
                              
                                 y
                              
                              
                                 1,5
                              
                           
                           
                              (
                              
                                 x
                                 ,
                                 t
                              
                              )
                           
                        
                        {y}_{\mathrm{1,5}}\left(x,t)
                     
                   by taking 
                     
                        
                        
                           α
                           =
                           1.5
                        
                        \alpha =1.5
                     
                  , 
                     
                        
                        
                           a
                           =
                           0.5
                           ,
                           β
                           =
                           1.3
                        
                        a=0.5,\beta =1.3
                     
                  , 
                     
                        
                        
                           v
                           =
                           1.5
                        
                        v=1.5
                     
                  , 
                     
                        
                        
                           
                              
                                 v
                              
                              
                                 1
                              
                           
                           =
                           1.5
                        
                        {v}_{1}=1.5
                     
                  , and 
                     
                        
                        
                           ψ
                           =
                           0.5
                        
                        \psi =0.5
                     
                  .
Figure 3

Plots of y 1,5 ( x , t ) by taking α = 1.5 , a = 0.5 , β = 1.3 , v = 1.5 , v 1 = 1.5 , and ψ = 0.5 .

In Figure 2(d)–(f), the real value graphs of y 1,4 ( x , t ) illustrate the periodic-singular patterns. Figure 2(g)–(i) depict the imaginary value graphs of y 1,4 ( x , t ) , revealing periodic singular characteristics akin to the real value graphs but with variations in amplitude. These graphical representations correspond to the specified constants: α = 1.5 , a = 0.5 , β = 1.3 , w 1 = 0.5 , v = 1.5 , w 0 = 0.5 , c = 1 , and ψ = 0.5 .

Figure 3(a)–(c) showcases the absolute value graph of y 1,5 ( x , t ) , unveiling a distinctive dark optical soliton. In Figures 3(d)–(f), the real value graphs of y 1,5 ( x , t ) illustrate the periodic patterns. Figure 3(g)–(i) depicts the imaginary value graphs of y 1,5 ( x , t ) , revealing periodic characteristics same to the real value graphs but with variations in amplitude. These graphical representations correspond to the specified constants: α = 1.5 , a = 0.5 , β = 1.3 , v = 1.5 , v 1 = 1.5 , and ψ = 0.5 .

Figure 4 
               Plots of 
                     
                        
                        
                           
                              
                                 y
                              
                              
                                 2,1
                              
                           
                           
                              (
                              
                                 x
                                 ,
                                 t
                              
                              )
                           
                        
                        {y}_{\mathrm{2,1}}\left(x,t)
                     
                   by taking 
                     
                        
                        
                           α
                           =
                           1.5
                        
                        \alpha =1.5
                     
                  , 
                     
                        
                        
                           a
                           =
                           0.5
                           ,
                           β
                           =
                           1.3
                        
                        a=0.5,\beta =1.3
                     
                  , 
                     
                        
                        
                           v
                           =
                           1.5
                        
                        v=1.5
                     
                  , 
                     
                        
                        
                           R
                           =
                           1
                        
                        R=1
                     
                  , 
                     
                        
                        
                           
                              
                                 h
                              
                              
                                 1
                              
                           
                           =
                           1
                        
                        {h}_{1}=1
                     
                  , 
                     
                        
                        
                           δ
                           =
                           0.45
                        
                        \delta =0.45
                     
                  , and 
                     
                        
                        
                           ψ
                           =
                           0.5
                        
                        \psi =0.5
                     
                  .
Figure 4

Plots of y 2,1 ( x , t ) by taking α = 1.5 , a = 0.5 , β = 1.3 , v = 1.5 , R = 1 , h 1 = 1 , δ = 0.45 , and ψ = 0.5 .

Figure 4(a)–(c) showcases the absolute value graph of y 2,1 ( x , t ) , which represents the bright optical soliton. In Figure 4(d)–(f), the real value graphs of y 2,1 ( x , t ) illustrate the dark–bright optical soliton. Figure 4(g)–(i) depicts the imaginary value graphs of y 2,1 ( x , t ) , revealing dark–bright characteristics same to the real value graphs but with variations in amplitude. These graphical representations correspond to the specified constants: α = 1.5 , a = 0.5 , β = 1.3 , v = 1.5 , R = 1 , h 1 = 1 , δ = 0.45 , and ψ = 0.5 .

Figure 5 
               Plots of 
                     
                        
                        
                           
                              
                                 y
                              
                              
                                 3,3
                              
                           
                           
                              (
                              
                                 x
                                 ,
                                 t
                              
                              )
                           
                        
                        {y}_{\mathrm{3,3}}\left(x,t)
                     
                   by taking 
                     
                        
                        
                           α
                           =
                           1.5
                        
                        \alpha =1.5
                     
                  , 
                     
                        
                        
                           
                              
                                 g
                              
                              
                                 1
                              
                           
                           =
                           2
                        
                        {g}_{1}=2
                     
                  , 
                     
                        
                        
                           a
                           =
                           0.5
                           ,
                           β
                           =
                           1.3
                        
                        a=0.5,\beta =1.3
                     
                  , 
                     
                        
                        
                           v
                           =
                           1.5
                        
                        v=1.5
                     
                  , 
                     
                        
                        
                           
                              
                                 p
                              
                              
                                 0
                              
                           
                           =
                           0.1
                           ,
                           
                              
                                 p
                              
                              
                                 2
                              
                           
                           =
                           −
                           0.2
                        
                        {p}_{0}=0.1,{p}_{2}=-0.2
                     
                  , 
                     
                        
                        
                           
                              
                                 g
                              
                              
                                 −
                                 1
                              
                           
                           =
                           1
                           ,
                           
                              
                                 g
                              
                              
                                 0
                              
                           
                           =
                           0.5
                        
                        {g}_{-1}=1,{g}_{0}=0.5
                     
                  , and 
                     
                        
                        
                           ψ
                           =
                           0.5
                        
                        \psi =0.5
                     
                  .
Figure 5

Plots of y 3,3 ( x , t ) by taking α = 1.5 , g 1 = 2 , a = 0.5 , β = 1.3 , v = 1.5 , p 0 = 0.1 , p 2 = 0.2 , g 1 = 1 , g 0 = 0.5 , and ψ = 0.5 .

Figure 5(a)–(c) showcases the absolute value graph of y 3,3 ( x , t ) , which represents the combined dark–singular optical soliton. In Figure 5(d)–(f), the real value graphs of y 3,3 ( x , t ) illustrate the singular optical soliton. Figure 5(g)–(i) depicts the imaginary value graphs of y 3,3 ( x , t ) , revealing singular optical solitons similar to the real value graphs but with variations in amplitude. These graphical representations correspond to the specified constants: α = 1.5 , a = 0.5 , β = 1.3 , v = 1.5 , p 0 = 0.1 , p 2 = 0.2 , g 1 = 1 , g 0 = 0.5 , and ψ = 0.5 .

Figure 6 
               Plots of 
                     
                        
                        
                           
                              
                                 y
                              
                              
                                 3,5
                              
                           
                           
                              (
                              
                                 x
                                 ,
                                 t
                              
                              )
                           
                        
                        {y}_{\mathrm{3,5}}\left(x,t)
                     
                   by taking 
                     
                        
                        
                           α
                           =
                           0.1
                        
                        \alpha =0.1
                     
                  , 
                     
                        
                        
                           a
                           =
                           0.99
                        
                        a=0.99
                     
                  , 
                     
                        
                        
                           β
                           =
                           1
                           ,
                           
                              
                                 g
                              
                              
                                 −
                                 1
                              
                           
                           =
                           −
                           0.5
                           ,
                           
                              
                                 g
                              
                              
                                 1
                              
                           
                           =
                           0.5
                           ,
                           
                              
                                 p
                              
                              
                                 1
                              
                           
                           =
                           0.1
                           ,
                           
                              
                                 p
                              
                              
                                 3
                              
                           
                           =
                           −
                           0.2
                           ,
                           v
                           =
                           0.5
                        
                        \beta =1,{g}_{-1}=-0.5,{g}_{1}=0.5,{p}_{1}=0.1,{p}_{3}=-0.2,v=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 g
                              
                              
                                 0
                              
                           
                           =
                           1
                           .
                        
                        {g}_{0}=1.
                     
                  .
Figure 6

Plots of y 3,5 ( x , t ) by taking α = 0.1 , a = 0.99 , β = 1 , g 1 = 0.5 , g 1 = 0.5 , p 1 = 0.1 , p 3 = 0.2 , v = 0.5 , and g 0 = 1 . .

Figure 6(a)–(c) depicts the absolute value graph of y 3,5 ( x , t ) , which represents the singular optical soliton. In Figure 6(d)–(f), the real value graphs of y 3,5 ( x , t ) illustrate the singular optical soliton. Figures 6(g)–(i) depicts the imaginary value graphs of y 3,5 ( x , t ) , revealing singular optical solitons same to the real value graphs but with variations in amplitude under specified constants: α = 0.1 , a = 0.99 , β = 1 , g 1 = 0.5 , g 1 = 0.5 , p 1 = 0.1 , p 3 = 0.2 , v = 0.5 , and g 0 = 1 .

Figure 7 
               Plots of 
                     
                        
                        
                           
                              
                                 y
                              
                              
                                 3,15
                              
                           
                           
                              (
                              
                                 x
                                 ,
                                 t
                              
                              )
                           
                        
                        {y}_{\mathrm{3,15}}\left(x,t)
                     
                   by taking 
                     
                        
                        
                           α
                           =
                           0.1
                        
                        \alpha =0.1
                     
                  , 
                     
                        
                        
                           a
                           =
                           0.99
                        
                        a=0.99
                     
                  , 
                     
                        
                        
                           β
                           =
                           1
                           ,
                           
                              
                                 g
                              
                              
                                 1
                              
                           
                           =
                           −
                           0.5
                           ,
                           
                              
                                 p
                              
                              
                                 0
                              
                           
                           =
                           0.1
                           ,
                           
                              
                                 p
                              
                              
                                 2
                              
                           
                           =
                           0.2
                           ,
                           v
                           =
                           0.5
                        
                        \beta =1,{g}_{1}=-0.5,{p}_{0}=0.1,{p}_{2}=0.2,v=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 g
                              
                              
                                 0
                              
                           
                           =
                           1
                           .
                        
                        {g}_{0}=1.
Figure 7

Plots of y 3,15 ( x , t ) by taking α = 0.1 , a = 0.99 , β = 1 , g 1 = 0.5 , p 0 = 0.1 , p 2 = 0.2 , v = 0.5 , and g 0 = 1 .

Figure 7(a)–(c) illustrates the absolute value graph of y 3,15 ( x , t ) , which represents the periodic-singular optical soliton. In Figure 7(d)–(f), the real value graphs of y 3,15 ( x , t ) illustrate periodic singular optical soliton. Figure 7(g)–(i) depicts the imaginary value graphs of y 3,15 ( x , t ) , revealing periodic-singular optical solitons same to the real value graphs but with variations in amplitude and phase component under specified constants: α = 0.1 , a = 0.99 , β = 1 , g 1 = 0.5 , p 0 = 0.1 , p 2 = 0.2 , v = 0.5 , and g 0 = 1 .

Figure 8 
               Comparison of the original soliton (solid blue) and perturbed soliton (dashed red) showing both the slowdown effect and radiation emission. The perturbed soliton shows a backward shift due to the slowdown, and small oscillatory waves are emitted as radiation.
Figure 8

Comparison of the original soliton (solid blue) and perturbed soliton (dashed red) showing both the slowdown effect and radiation emission. The perturbed soliton shows a backward shift due to the slowdown, and small oscillatory waves are emitted as radiation.

The inclusion of TOD is crucial for accurately modeling soliton behavior in dispersive media. Without considering the effects of TOD, the analysis is incomplete, as it neglects the physically important phenomena of soliton slowdown and radiation. The perturbative approach used here captures these effects and aligns with the expected behavior of solitons in the presence of higher-order dispersion, as shown in Figure 8.

8 Conclusion

This study has delved into an exploration of the NLSHE, which incorporates CD and serves as a governing model for comprehending the dynamics of dispersive pulse propagation in optical fibers. Employing three distinct methodologies, namely, the 1 φ ( η ) , φ ( η ) φ ( η ) method, the new Kudryashov method, and the extended simple equation method, we have successfully derived solutions for this equation. Visualization of the results through 3D and 2D graphs, illustrating absolute, real, and imaginary values, provides a comprehensive understanding of soliton dynamics. These graphical representations underscore the existence of diverse soliton solutions, encompassing dark, bright, singular, periodic singular, dark–bright, and dark-singular solitons. The implications of these findings extend to the field of fiber–optic communication, where they significantly contribute to advancing our insights into the intricate dynamics of solitons. This research holds significant importance and may pave the way for further advancements and applications in the field of fiber–optic communication. By perturbing the soliton solution to include the effects of TOD, we have shown that solitons experience a gradual deceleration and emit small radiative waves. These effects are consistent with the underlying physics of dispersive systems, making the analysis more comprehensive and physically realistic.

Acknowledgments

The authors gratefully acknowledge the financial support provided by the National Science and Technology Council (NSTC), Taiwan, under Grant No. 112-2115-M-006-002.

  1. Funding information: This work was supported by the National Science and Technology Council (NSTC), Taiwan, under Grant No. 112-2115-M-006-002.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: All data generated or analyzed during this study are included in this published article.

References

[1] Asghari Y, Eslami M, Matinfar M, Rezazadeh H. Novel soliton solution of discrete nonlinear Schrödinger system in nonlinear optical fiber. Alexandr Eng J. 2024;90:7–16. 10.1016/j.aej.2024.01.024Search in Google Scholar

[2] Eslami M, Rezazadeh H. The first integral method for Wu-Zhang system with conformable time-fractional derivative. Calcolo. 2016;53:475–85. 10.1007/s10092-015-0158-8Search in Google Scholar

[3] Yang J. Computer data encryption system based on nonlinear partial differential equations. Mobile Inform Syst. 2022;2022:3395019. 10.1155/2022/3395019Search in Google Scholar

[4] Liang BL, Wang Y, Zhang SH, Guo QL, Wang SF, Fu GS, et al. Optical image processing by using a photorefractive spatial soliton waveguide. Phys Lett A. 2017;381(13):1207–12. 10.1016/j.physleta.2017.01.035Search in Google Scholar

[5] Borgese G, Vena S, Pantano P, Pace C, Bilotta E. Simulation, modeling, and analysis of soliton waves interaction and propagation in CNN transmission lines for innovative data communication and processing. Disc Dyn Nat Soc. 2015;2015:1207–12. 10.1155/2015/139238Search in Google Scholar

[6] Guo S, Fritsch AR, Greenberg C, Spielman IB, Zwolak JP. Machine-learning enhanced dark soliton detection in Bose-Einstein condensates. Machine Learn Sci Tech. 2021;2(3):035020. 10.1088/2632-2153/abed1eSearch in Google Scholar PubMed PubMed Central

[7] Al Qarni A, Bodaqah A, Mohammed A, Alshaery A, Bakodah H, Biswas A. Cubic-quartic optical solitons for Lakshmanan-Porsezian-Daniel equation by the improved adomian decomposition scheme. Ukr J Phys Opt. 2022;23(4):228–42. 10.3116/16091833/23/4/228/2022Search in Google Scholar

[8] AlQarni AA, Bodaqah AM, Mohammed ASHF, Alshaery AA, Bakodah HO, et al. Dark and singular cubic-quartic optical solitons with Lakshmanan-Porsezian-Daniel equation by the improved adomian decomposition scheme. Ukrain J Phys Opt. 2023;24(1):46–61. 10.3116/16091833/24/1/46/2023Search in Google Scholar

[9] Jawad AJM, Abu-AlShaeer MJ. Highly dispersive optical solitons with cubic law and cubic-quintic-septic law nonlinearities by two methods. Al-Rafidain J Eng Sci. 2023;1(1):1–8. 10.61268/sapgh524Search in Google Scholar

[10] Jawad AJM, Biswas A, Yildirim Y, Alshomrani AS. Dark-singular straddled optical solitons for the dispersive concatenation model with power-law of self-phase modulation by Tanh-Coth approach. Contemp Math. 2024;5(3):3198–214. 10.37256/cm.5320245013Search in Google Scholar

[11] Jawad AJM, Biswas A, Yildirim Y, Alshomrani AS. Optical solitons with differential group delay and inter-modal dispersion singlet. Contemp Math. 2024;5(1):1054–71. 10.37256/cm.5120244121Search in Google Scholar

[12] Jawad A, Biswas A. Solutions of resonant nonlinear Schrödingeras equation with exotic non-Kerr law nonlinearities. Al-Rafidain J Eng Sci. 2024:43–50. 10.61268/2bz73q95Search in Google Scholar

[13] Jawad AJM, Biswas A, Yiiildiiiriiim Y, Alshomrani AS. A fresh perspective on the concatenation model in optical fibers with kerr law of self-phase modulation. Engineering Science & Technology. 2024;5:195–208. 10.37256/est.5220244160Search in Google Scholar

[14] Jawad AJM, Yildirim Y, Biswas A, Ibraheem IK, Alshomrani AS. Highly dispersive optical solitons with differential group delay for Kerr law of self-phase modulation by Sardar sub-equation approach. Contemp Math. 2024:3839–57. 10.37256/cm.5320245185Search in Google Scholar

[15] Adem AR, Biswas A, Yildirim Y, Alshomrani AS. Revisitation of “implicit quiescent optical solitons with complex Ginzburg-Landau equation having nonlinear chromatic dispersion”: generalized temporal evolution. J Optics. 2024:1–10. 10.1007/s12596-024-01759-4Search in Google Scholar

[16] González-Gaxiola O, Biswas A, Belic MR. Optical soliton perturbation of Fokas-Lenells equation by the Laplace-Adomian decomposition algorithm. J Europ Opt Soc-Rapid Publ. 2019;15:1–9. 10.1186/s41476-019-0111-6Search in Google Scholar

[17] Zayed EM, Alurrfi KA, Alshbear RA. On application of the new mapping method to magneto-optic waveguides having Kudryashovas law of refractive index. Optik. 2023;287:171072. 10.1016/j.ijleo.2023.171072Search in Google Scholar

[18] González-Gaxiola O, Biswas A, Yildirim Y, Alshomrani AS. Bright optical solitons for the concatenation model with power-law nonlinearity: Laplace-Adomian decomposition. Contemporary Mathematics. 2023;4(4):1234–48. 10.37256/cm.4420233705Search in Google Scholar

[19] Topkara E, Milovic D, Sarma AK, Majid F, Biswas A. A study of optical solitons with Kerr and power law nonlinearities by He’s variational principle. J Europ Opt Soc-Rapid Publ. 2009;4:09050. 10.2971/jeos.2009.09050Search in Google Scholar

[20] Zayed EM, El-Horbaty M, Alngar ME, Shohib RM, Biswas A, Yiiildiiiriiim Y, et al. Dynamical system of optical soliton parameters by variational principle (super-Gaussian and super-sech pulses). J Europ Opt Soc-Rapid Publ. 2023;19(2):38. 10.1051/jeos/2023035Search in Google Scholar

[21] Arnous AH, Nofal TA, Biswas A, Yıldırım Y, Asiri A. Cubic-quartic optical solitons of the complex Ginzburg-Landau equation: a novel approach. Nonl Dyn. 2023;111(21):20201–16. 10.1007/s11071-023-08854-4Search in Google Scholar

[22] Jihad N, Abd Almuhsan M. Evaluation of impairment mitigations for optical fiber communications using dispersion compensation techniques. Rafidain J Eng Sci. 2023;1(1):81–92. 10.61268/0dat0751Search in Google Scholar

[23] Drazin PG, Johnson RS. Solitons: an introduction. vol. 2. Cambridge: Cambridge University Press; 1989. 10.1017/CBO9781139172059Search in Google Scholar

[24] Mirzazadeh M, Eslami M, Biswas A. Dispersive optical solitons by Kudryashovas method. Optik. 2014;125(23):6874–80. 10.1016/j.ijleo.2014.02.044Search in Google Scholar

[25] Nazarzadeh A, Eslami M, Mirzazadeh M. Exact solutions of some nonlinear partial differential equations using functional variable method. Pramana. 2013;81:225–36. 10.1007/s12043-013-0565-9Search in Google Scholar

[26] Ur Rehman H, Asjad Imran M, Bibi M, Riaz M, Akgül A. New soliton solutions of the 2D-chiral nonlinear Schrodinger equation using two integration schemes. Math Meth Appl Sci. 2021;44(7):5663–82. 10.1002/mma.7140Search in Google Scholar

[27] Sultan AM, Lu D, Arshad M, Rehman HU, Saleem MS. Soliton solutions of higher order dispersive cubic-quintic nonlinear Schrödinger equation and its applications. Chin J Phys. 2020;67:405–13. 10.1016/j.cjph.2019.10.003Search in Google Scholar

[28] Plastino AR, Tsallis C. Nonlinear Schroedinger equation in the presence of uniform acceleration. J Math Phys. 2013;54(4):041505. 10.1063/1.4798999Search in Google Scholar

[29] Ozdemir N, Secer A, Ozisik M, Bayram M. Perturbation of dispersive optical solitons with Schrödinger–Hirota equation with Kerr law and spatio-temporal dispersion. Optik. 2022;265:169545. 10.1016/j.ijleo.2022.169545Search in Google Scholar

[30] Xu T, Jacobsen G, Popov S, Li J, Vanin E, Wang K, et al. Chromatic dispersion compensation in coherent transmission system using digital filters. Optics express. 2010;18(15):16243–57. 10.1364/OE.18.016243Search in Google Scholar PubMed

[31] Udayakumar R, Khanaa V, Saravanan T. Chromatic dispersion compensation in optical fiber communication system and its simulation. Indian J Sci Tech. 2013;6(6):4762–6. 10.17485/ijst/2013/v6isp6.2Search in Google Scholar

[32] Dar AB, Jha RK. Chromatic dispersion compensation techniques and characterization of fiber Bragg grating for dispersion compensation. Opt Quant Electron. 2017;49:1–35. 10.1007/s11082-017-0944-4Search in Google Scholar

[33] Ekici M, Mirzazadeh M, Sonmezoglu A, Ullah MZ, Asma M, Zhou Q, et al. Dispersive optical solitons with Schrödinger–Hirota equation by extended trial equation method. Optik. 2017;136:451–61. 10.1016/j.ijleo.2017.02.042Search in Google Scholar

[34] Ozisik M, Onder I, Esen H, Cinar M, Ozdemir N, Secer A, et al. On the investigation of optical soliton solutions of cubic-quartic Fokas-Lenells and Schrödinger–Hirota equations. Optik. 2023;272:170389. 10.1016/j.ijleo.2022.170389Search in Google Scholar

[35] Rehman HU, Ullah N, Imran M. Highly dispersive optical solitons using Kudryashovas method. Optik. 2019;199:163349. 10.1016/j.ijleo.2019.163349Search in Google Scholar

[36] Ryabov PN, Sinelshchikov DI, Kochanov MB. Application of the Kudryashov method for finding exact solutions of the high order nonlinear evolution equations. Appl Math Comput. 2011;218(7):3965–72. 10.1016/j.amc.2011.09.027Search in Google Scholar

[37] Shi D, Rehman HU, Iqbal I, Vivas-Cortez M, Saleem MS, Zhang X. Analytical study of the dynamics in the double-chain model of DNA. Results Phys. 2023;52:106787.10.1016/j.rinp.2023.106787Search in Google Scholar

[38] Rehman HU, Awan AU, Tag-ElDin EM, Alhazmi SE, Yassen MF, Haider R. Extended hyperbolic function method for the (2+1)-dimensional nonlinear soliton equation. Results Phys. 2022;40:105802. 10.1016/j.rinp.2022.105802Search in Google Scholar

[39] Hong Z, Ji-Guang H, Wei-Tao W, Hong-Yong A. Applications of extended hyperbolic function method for quintic discrete nonlinear Schrödinger equation. Commun Theoret Phys. 2007;47(3):474. 10.1088/0253-6102/47/3/020Search in Google Scholar

[40] Ullah N, Asjad MI, Muhammad T, Akgül A. Analysis of power law non-linearity in solitonic solutions using extended hyperbolic function method. AIMS Math. 2022;7:18603–15.10.3934/math.20221023Search in Google Scholar

[41] Rehman HU, Imran MA, Ullah N, Akgül A. Exact solutions of (2+1)-dimensional Schrödingeras hyperbolic equation using different techniques. Numer Meth Partial Differ Equ. 2020;39:4575–94.10.1002/num.22644Search in Google Scholar

[42] Salam MA, Uddin MS, Dey P. Generalized Bernoulli Sub-ODE method and its applications. An Pure Appl Math. 2015;10(1):1–6. Search in Google Scholar

[43] Zheng B. A new Bernoulli sub-ODE method for constructing traveling wave solutions for two nonlinear equations with any order. UPB Scient Bullet Ser A. 2011;73(3):85–94. Search in Google Scholar

[44] Rehman HU, Iqbal I, Subhi Aiadi S, Mlaiki N, Saleem MS. Soliton solutions of Klein-Fock-Gordon equation using Sardar subequation method. Mathematics. 2022;10(18):3377. 10.3390/math10183377Search in Google Scholar

[45] Iqbal I, Rehman HU, Mirzazadeh M, Hashemi MS. Retrieval of optical solitons for nonlinear models with Kudryashov’s quintuple power law and dual-form nonlocal nonlinearity. Opt Quantum Electronics. 2023;55(7):588. 10.1007/s11082-023-04866-xSearch in Google Scholar

[46] Rehman HU, Ullah N, Imran M. Optical solitons of Biswas-Arshed equation in birefringent fibers using extended direct algebraic method. Optik. 2021;226:165378. 10.1016/j.ijleo.2020.165378Search in Google Scholar

[47] Kurt A, Tozar A, Tasbozan O. Applying the new extended direct algebraic method to solve the equation of obliquely interacting waves in shallow waters. J Ocean Univ China. 2020;19:772–80. 10.1007/s11802-020-4135-8Search in Google Scholar

[48] Liu S, Fu Z, Liu S, Zhao Q. Jacobi elliptic function expansion method and periodic wave solutions of nonlinear wave equations. Phys Lett A. 2001;289(1–2):69–74. 10.1016/S0375-9601(01)00580-1Search in Google Scholar

[49] Allahyani SA, Rehman HU, Awan AU, Tag-ElDin EM, Hassan MU. Diverse variety of exact solutions for nonlinear Gilson-Pickering equation. Symmetry. 2022;14(10):2151. 10.3390/sym14102151Search in Google Scholar

[50] Rehman HU, Saleem MS, Zubair M, Jafar S, Latif I. Optical solitons with Biswas-Arshed model using mapping method. Optik. 2019;194:163091. 10.1016/j.ijleo.2019.163091Search in Google Scholar

[51] Zeng X, Yong X. A new mapping method and its applications to nonlinear partial differential equations. Phys Lett A. 2008;372(44):6602–7. 10.1016/j.physleta.2008.09.025Search in Google Scholar

[52] Mohammed WW, Cesarano C. The soliton solutions for the (4+1)-dimensional stochastic Fokas equation. Math Meth Appl Sci. 2023;46(6):7589–97. 10.1002/mma.8986Search in Google Scholar

[53] Mohammed WW, Al-Askar FM, Cesarano C. The analytical solutions of the stochastic mKdV equation via the mapping method. Mathematics. 2022;10(22):4212. 10.3390/math10224212Search in Google Scholar

[54] Shahzad MU, Rehman HU, Awan AU, Zafar Z, Hassan AM, Iqbal I. Analysis of the exact solutions of nonlinear coupled Drinfeld-Sokolov-Wilson equation through phi6-model expansion method. Results Phys. 2023;52:106771. 10.1016/j.rinp.2023.106771Search in Google Scholar

[55] Rehman HU, Said GS, Amer A, Ashraf H, Tharwat M, Abdel-Aty M, et al. Unraveling the (4+1)-dimensional Davey-Stewartson-Kadomtsev-Petviashvili equation: Exploring soliton solutions via multiple techniques. Alexandr Eng J. 2024;90:17–23. 10.1016/j.aej.2024.01.058Search in Google Scholar

[56] Apeanti WO, Lu D, Zhang H, Yaro D, Akuamoah SW. Traveling wave solutions for complex nonlinear space-time fractional order (2+1)-dimensional Maccari dynamical system and Schrödinger equation with dual power law nonlinearity. SN Appl Sci. 2019;1:1–10. 10.1007/s42452-019-0537-xSearch in Google Scholar

[57] Biswas A, Porsezian K. Soliton perturbation theory for the modified nonlinear Schrödinger’s equation. Commun Nonl Sci Numer Simulat. 2007;12(6):886–903. 10.1016/j.cnsns.2005.11.006Search in Google Scholar

[58] Kaur L, Wazwaz AM. Bright-dark optical solitons for Schrödinger–Hirota equation with variable coefficients. Optik. 2019;179:479–84. 10.1016/j.ijleo.2018.09.035Search in Google Scholar

[59] Yildirim Y. Optical solitons to Schrödinger–Hirota equation in DWDM system with trial equation integration architecture. Optik. 2019;182:275–81. 10.1016/j.ijleo.2019.01.017Search in Google Scholar

[60] Zayed EM, Shohib RM, Biswas A, Ekici M, Alshomrani AS, Khan S, et al. Dispersive solitons in optical fibers and DWDM networks with Schrödinger–Hirota equation. Optik. 2019;199:163214. 10.1016/j.ijleo.2019.163214Search in Google Scholar

[61] Tang L. Bifurcations and dispersive optical solitons for the nonlinear Schrödinger–Hirota equation in DWDM networks. Optik. 2022;262:169276. 10.1016/j.ijleo.2022.169276Search in Google Scholar

Received: 2024-09-02
Revised: 2024-10-30
Accepted: 2024-11-05
Published Online: 2025-04-29

© 2025 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Research Articles
  2. Generalized (ψ,φ)-contraction to investigate Volterra integral inclusions and fractal fractional PDEs in super-metric space with numerical experiments
  3. Solitons in ultrasound imaging: Exploring applications and enhancements via the Westervelt equation
  4. Stochastic improved Simpson for solving nonlinear fractional-order systems using product integration rules
  5. Exploring dynamical features like bifurcation assessment, sensitivity visualization, and solitary wave solutions of the integrable Akbota equation
  6. Research on surface defect detection method and optimization of paper-plastic composite bag based on improved combined segmentation algorithm
  7. Impact the sulphur content in Iraqi crude oil on the mechanical properties and corrosion behaviour of carbon steel in various types of API 5L pipelines and ASTM 106 grade B
  8. Unravelling quiescent optical solitons: An exploration of the complex Ginzburg–Landau equation with nonlinear chromatic dispersion and self-phase modulation
  9. Perturbation-iteration approach for fractional-order logistic differential equations
  10. Variational formulations for the Euler and Navier–Stokes systems in fluid mechanics and related models
  11. Rotor response to unbalanced load and system performance considering variable bearing profile
  12. DeepFowl: Disease prediction from chicken excreta images using deep learning
  13. Channel flow of Ellis fluid due to cilia motion
  14. A case study of fractional-order varicella virus model to nonlinear dynamics strategy for control and prevalence
  15. Multi-point estimation weldment recognition and estimation of pose with data-driven robotics design
  16. Analysis of Hall current and nonuniform heating effects on magneto-convection between vertically aligned plates under the influence of electric and magnetic fields
  17. A comparative study on residual power series method and differential transform method through the time-fractional telegraph equation
  18. Insights from the nonlinear Schrödinger–Hirota equation with chromatic dispersion: Dynamics in fiber–optic communication
  19. Mathematical analysis of Jeffrey ferrofluid on stretching surface with the Darcy–Forchheimer model
  20. Exploring the interaction between lump, stripe and double-stripe, and periodic wave solutions of the Konopelchenko–Dubrovsky–Kaup–Kupershmidt system
  21. Computational investigation of tuberculosis and HIV/AIDS co-infection in fuzzy environment
  22. Signature verification by geometry and image processing
  23. Theoretical and numerical approach for quantifying sensitivity to system parameters of nonlinear systems
  24. Chaotic behaviors, stability, and solitary wave propagations of M-fractional LWE equation in magneto-electro-elastic circular rod
  25. Dynamic analysis and optimization of syphilis spread: Simulations, integrating treatment and public health interventions
  26. Visco-thermoelastic rectangular plate under uniform loading: A study of deflection
  27. Threshold dynamics and optimal control of an epidemiological smoking model
  28. Numerical computational model for an unsteady hybrid nanofluid flow in a porous medium past an MHD rotating sheet
  29. Regression prediction model of fabric brightness based on light and shadow reconstruction of layered images
  30. Dynamics and prevention of gemini virus infection in red chili crops studied with generalized fractional operator: Analysis and modeling
  31. Qualitative analysis on existence and stability of nonlinear fractional dynamic equations on time scales
  32. Fractional-order super-twisting sliding mode active disturbance rejection control for electro-hydraulic position servo systems
  33. Analytical exploration and parametric insights into optical solitons in magneto-optic waveguides: Advances in nonlinear dynamics for applied sciences
  34. Bifurcation dynamics and optical soliton structures in the nonlinear Schrödinger–Bopp–Podolsky system
  35. Review Article
  36. Haar wavelet collocation method for existence and numerical solutions of fourth-order integro-differential equations with bounded coefficients
  37. Special Issue: Nonlinear Analysis and Design of Communication Networks for IoT Applications - Part II
  38. Silicon-based all-optical wavelength converter for on-chip optical interconnection
  39. Research on a path-tracking control system of unmanned rollers based on an optimization algorithm and real-time feedback
  40. Analysis of the sports action recognition model based on the LSTM recurrent neural network
  41. Industrial robot trajectory error compensation based on enhanced transfer convolutional neural networks
  42. Research on IoT network performance prediction model of power grid warehouse based on nonlinear GA-BP neural network
  43. Interactive recommendation of social network communication between cities based on GNN and user preferences
  44. Application of improved P-BEM in time varying channel prediction in 5G high-speed mobile communication system
  45. Construction of a BIM smart building collaborative design model combining the Internet of Things
  46. Optimizing malicious website prediction: An advanced XGBoost-based machine learning model
  47. Economic operation analysis of the power grid combining communication network and distributed optimization algorithm
  48. Sports video temporal action detection technology based on an improved MSST algorithm
  49. Internet of things data security and privacy protection based on improved federated learning
  50. Enterprise power emission reduction technology based on the LSTM–SVM model
  51. Construction of multi-style face models based on artistic image generation algorithms
  52. Research and application of interactive digital twin monitoring system for photovoltaic power station based on global perception
  53. Special Issue: Decision and Control in Nonlinear Systems - Part II
  54. Animation video frame prediction based on ConvGRU fine-grained synthesis flow
  55. Application of GGNN inference propagation model for martial art intensity evaluation
  56. Benefit evaluation of building energy-saving renovation projects based on BWM weighting method
  57. Deep neural network application in real-time economic dispatch and frequency control of microgrids
  58. Real-time force/position control of soft growing robots: A data-driven model predictive approach
  59. Mechanical product design and manufacturing system based on CNN and server optimization algorithm
  60. Application of finite element analysis in the formal analysis of ancient architectural plaque section
  61. Research on territorial spatial planning based on data mining and geographic information visualization
  62. Fault diagnosis of agricultural sprinkler irrigation machinery equipment based on machine vision
  63. Closure technology of large span steel truss arch bridge with temporarily fixed edge supports
  64. Intelligent accounting question-answering robot based on a large language model and knowledge graph
  65. Analysis of manufacturing and retailer blockchain decision based on resource recyclability
  66. Flexible manufacturing workshop mechanical processing and product scheduling algorithm based on MES
  67. Exploration of indoor environment perception and design model based on virtual reality technology
  68. Tennis automatic ball-picking robot based on image object detection and positioning technology
  69. A new CNN deep learning model for computer-intelligent color matching
  70. Design of AR-based general computer technology experiment demonstration platform
  71. Indoor environment monitoring method based on the fusion of audio recognition and video patrol features
  72. Health condition prediction method of the computer numerical control machine tool parts by ensembling digital twins and improved LSTM networks
  73. Establishment of a green degree evaluation model for wall materials based on lifecycle
  74. Quantitative evaluation of college music teaching pronunciation based on nonlinear feature extraction
  75. Multi-index nonlinear robust virtual synchronous generator control method for microgrid inverters
  76. Manufacturing engineering production line scheduling management technology integrating availability constraints and heuristic rules
  77. Analysis of digital intelligent financial audit system based on improved BiLSTM neural network
  78. Attention community discovery model applied to complex network information analysis
  79. A neural collaborative filtering recommendation algorithm based on attention mechanism and contrastive learning
  80. Rehabilitation training method for motor dysfunction based on video stream matching
  81. Research on façade design for cold-region buildings based on artificial neural networks and parametric modeling techniques
  82. Intelligent implementation of muscle strain identification algorithm in Mi health exercise induced waist muscle strain
  83. Optimization design of urban rainwater and flood drainage system based on SWMM
  84. Improved GA for construction progress and cost management in construction projects
  85. Evaluation and prediction of SVM parameters in engineering cost based on random forest hybrid optimization
  86. Museum intelligent warning system based on wireless data module
  87. Optimization design and research of mechatronics based on torque motor control algorithm
  88. Special Issue: Nonlinear Engineering’s significance in Materials Science
  89. Experimental research on the degradation of chemical industrial wastewater by combined hydrodynamic cavitation based on nonlinear dynamic model
  90. Study on low-cycle fatigue life of nickel-based superalloy GH4586 at various temperatures
  91. Some results of solutions to neutral stochastic functional operator-differential equations
  92. Ultrasonic cavitation did not occur in high-pressure CO2 liquid
  93. Research on the performance of a novel type of cemented filler material for coal mine opening and filling
  94. Testing of recycled fine aggregate concrete’s mechanical properties using recycled fine aggregate concrete and research on technology for highway construction
  95. A modified fuzzy TOPSIS approach for the condition assessment of existing bridges
  96. Nonlinear structural and vibration analysis of straddle monorail pantograph under random excitations
  97. Achieving high efficiency and stability in blue OLEDs: Role of wide-gap hosts and emitter interactions
  98. Construction of teaching quality evaluation model of online dance teaching course based on improved PSO-BPNN
  99. Enhanced electrical conductivity and electromagnetic shielding properties of multi-component polymer/graphite nanocomposites prepared by solid-state shear milling
  100. Optimization of thermal characteristics of buried composite phase-change energy storage walls based on nonlinear engineering methods
  101. A higher-performance big data-based movie recommendation system
  102. Nonlinear impact of minimum wage on labor employment in China
  103. Nonlinear comprehensive evaluation method based on information entropy and discrimination optimization
  104. Application of numerical calculation methods in stability analysis of pile foundation under complex foundation conditions
  105. Research on the contribution of shale gas development and utilization in Sichuan Province to carbon peak based on the PSA process
  106. Characteristics of tight oil reservoirs and their impact on seepage flow from a nonlinear engineering perspective
  107. Nonlinear deformation decomposition and mode identification of plane structures via orthogonal theory
  108. Numerical simulation of damage mechanism in rock with cracks impacted by self-excited pulsed jet based on SPH-FEM coupling method: The perspective of nonlinear engineering and materials science
  109. Cross-scale modeling and collaborative optimization of ethanol-catalyzed coupling to produce C4 olefins: Nonlinear modeling and collaborative optimization strategies
  110. Unequal width T-node stress concentration factor analysis of stiffened rectangular steel pipe concrete
  111. Special Issue: Advances in Nonlinear Dynamics and Control
  112. Development of a cognitive blood glucose–insulin control strategy design for a nonlinear diabetic patient model
  113. Big data-based optimized model of building design in the context of rural revitalization
  114. Multi-UAV assisted air-to-ground data collection for ground sensors with unknown positions
  115. Design of urban and rural elderly care public areas integrating person-environment fit theory
  116. Application of lossless signal transmission technology in piano timbre recognition
  117. Application of improved GA in optimizing rural tourism routes
  118. Architectural animation generation system based on AL-GAN algorithm
  119. Advanced sentiment analysis in online shopping: Implementing LSTM models analyzing E-commerce user sentiments
  120. Intelligent recommendation algorithm for piano tracks based on the CNN model
  121. Visualization of large-scale user association feature data based on a nonlinear dimensionality reduction method
  122. Low-carbon economic optimization of microgrid clusters based on an energy interaction operation strategy
  123. Optimization effect of video data extraction and search based on Faster-RCNN hybrid model on intelligent information systems
  124. Construction of image segmentation system combining TC and swarm intelligence algorithm
  125. Particle swarm optimization and fuzzy C-means clustering algorithm for the adhesive layer defect detection
  126. Optimization of student learning status by instructional intervention decision-making techniques incorporating reinforcement learning
  127. Fuzzy model-based stabilization control and state estimation of nonlinear systems
  128. Optimization of distribution network scheduling based on BA and photovoltaic uncertainty
  129. Tai Chi movement segmentation and recognition on the grounds of multi-sensor data fusion and the DBSCAN algorithm
  130. Special Issue: Dynamic Engineering and Control Methods for the Nonlinear Systems - Part III
  131. Generalized numerical RKM method for solving sixth-order fractional partial differential equations
Downloaded on 14.11.2025 from https://www.degruyterbrill.com/document/doi/10.1515/nleng-2024-0058/html
Scroll to top button