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. Review Article
  31. Haar wavelet collocation method for existence and numerical solutions of fourth-order integro-differential equations with bounded coefficients
  32. Special Issue: Nonlinear Analysis and Design of Communication Networks for IoT Applications - Part II
  33. Silicon-based all-optical wavelength converter for on-chip optical interconnection
  34. Research on a path-tracking control system of unmanned rollers based on an optimization algorithm and real-time feedback
  35. Analysis of the sports action recognition model based on the LSTM recurrent neural network
  36. Industrial robot trajectory error compensation based on enhanced transfer convolutional neural networks
  37. Research on IoT network performance prediction model of power grid warehouse based on nonlinear GA-BP neural network
  38. Interactive recommendation of social network communication between cities based on GNN and user preferences
  39. Application of improved P-BEM in time varying channel prediction in 5G high-speed mobile communication system
  40. Construction of a BIM smart building collaborative design model combining the Internet of Things
  41. Optimizing malicious website prediction: An advanced XGBoost-based machine learning model
  42. Economic operation analysis of the power grid combining communication network and distributed optimization algorithm
  43. Sports video temporal action detection technology based on an improved MSST algorithm
  44. Internet of things data security and privacy protection based on improved federated learning
  45. Enterprise power emission reduction technology based on the LSTM–SVM model
  46. Construction of multi-style face models based on artistic image generation algorithms
  47. Special Issue: Decision and Control in Nonlinear Systems - Part II
  48. Animation video frame prediction based on ConvGRU fine-grained synthesis flow
  49. Application of GGNN inference propagation model for martial art intensity evaluation
  50. Benefit evaluation of building energy-saving renovation projects based on BWM weighting method
  51. Deep neural network application in real-time economic dispatch and frequency control of microgrids
  52. Real-time force/position control of soft growing robots: A data-driven model predictive approach
  53. Mechanical product design and manufacturing system based on CNN and server optimization algorithm
  54. Application of finite element analysis in the formal analysis of ancient architectural plaque section
  55. Research on territorial spatial planning based on data mining and geographic information visualization
  56. Fault diagnosis of agricultural sprinkler irrigation machinery equipment based on machine vision
  57. Closure technology of large span steel truss arch bridge with temporarily fixed edge supports
  58. Intelligent accounting question-answering robot based on a large language model and knowledge graph
  59. Analysis of manufacturing and retailer blockchain decision based on resource recyclability
  60. Flexible manufacturing workshop mechanical processing and product scheduling algorithm based on MES
  61. Exploration of indoor environment perception and design model based on virtual reality technology
  62. Tennis automatic ball-picking robot based on image object detection and positioning technology
  63. A new CNN deep learning model for computer-intelligent color matching
  64. Design of AR-based general computer technology experiment demonstration platform
  65. Indoor environment monitoring method based on the fusion of audio recognition and video patrol features
  66. Health condition prediction method of the computer numerical control machine tool parts by ensembling digital twins and improved LSTM networks
  67. Establishment of a green degree evaluation model for wall materials based on lifecycle
  68. Quantitative evaluation of college music teaching pronunciation based on nonlinear feature extraction
  69. Multi-index nonlinear robust virtual synchronous generator control method for microgrid inverters
  70. Manufacturing engineering production line scheduling management technology integrating availability constraints and heuristic rules
  71. Analysis of digital intelligent financial audit system based on improved BiLSTM neural network
  72. Attention community discovery model applied to complex network information analysis
  73. A neural collaborative filtering recommendation algorithm based on attention mechanism and contrastive learning
  74. Rehabilitation training method for motor dysfunction based on video stream matching
  75. Research on façade design for cold-region buildings based on artificial neural networks and parametric modeling techniques
  76. Intelligent implementation of muscle strain identification algorithm in Mi health exercise induced waist muscle strain
  77. Optimization design of urban rainwater and flood drainage system based on SWMM
  78. Improved GA for construction progress and cost management in construction projects
  79. Evaluation and prediction of SVM parameters in engineering cost based on random forest hybrid optimization
  80. Special Issue: Nonlinear Engineering’s significance in Materials Science
  81. Experimental research on the degradation of chemical industrial wastewater by combined hydrodynamic cavitation based on nonlinear dynamic model
  82. Study on low-cycle fatigue life of nickel-based superalloy GH4586 at various temperatures
  83. Some results of solutions to neutral stochastic functional operator-differential equations
  84. Ultrasonic cavitation did not occur in high-pressure CO2 liquid
  85. Research on the performance of a novel type of cemented filler material for coal mine opening and filling
  86. Testing of recycled fine aggregate concrete’s mechanical properties using recycled fine aggregate concrete and research on technology for highway construction
  87. A modified fuzzy TOPSIS approach for the condition assessment of existing bridges
  88. Nonlinear structural and vibration analysis of straddle monorail pantograph under random excitations
  89. Achieving high efficiency and stability in blue OLEDs: Role of wide-gap hosts and emitter interactions
  90. Construction of teaching quality evaluation model of online dance teaching course based on improved PSO-BPNN
  91. Enhanced electrical conductivity and electromagnetic shielding properties of multi-component polymer/graphite nanocomposites prepared by solid-state shear milling
  92. Optimization of thermal characteristics of buried composite phase-change energy storage walls based on nonlinear engineering methods
  93. A higher-performance big data-based movie recommendation system
  94. Nonlinear impact of minimum wage on labor employment in China
  95. Nonlinear comprehensive evaluation method based on information entropy and discrimination optimization
  96. Application of numerical calculation methods in stability analysis of pile foundation under complex foundation conditions
  97. Research on the contribution of shale gas development and utilization in Sichuan Province to carbon peak based on the PSA process
  98. Characteristics of tight oil reservoirs and their impact on seepage flow from a nonlinear engineering perspective
  99. Nonlinear deformation decomposition and mode identification of plane structures via orthogonal theory
  100. 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
  101. Cross-scale modeling and collaborative optimization of ethanol-catalyzed coupling to produce C4 olefins: Nonlinear modeling and collaborative optimization strategies
  102. Special Issue: Advances in Nonlinear Dynamics and Control
  103. Development of a cognitive blood glucose–insulin control strategy design for a nonlinear diabetic patient model
  104. Big data-based optimized model of building design in the context of rural revitalization
  105. Multi-UAV assisted air-to-ground data collection for ground sensors with unknown positions
  106. Design of urban and rural elderly care public areas integrating person-environment fit theory
  107. Application of lossless signal transmission technology in piano timbre recognition
  108. Application of improved GA in optimizing rural tourism routes
  109. Architectural animation generation system based on AL-GAN algorithm
  110. Advanced sentiment analysis in online shopping: Implementing LSTM models analyzing E-commerce user sentiments
  111. Intelligent recommendation algorithm for piano tracks based on the CNN model
  112. Visualization of large-scale user association feature data based on a nonlinear dimensionality reduction method
  113. Low-carbon economic optimization of microgrid clusters based on an energy interaction operation strategy
  114. Optimization effect of video data extraction and search based on Faster-RCNN hybrid model on intelligent information systems
  115. Construction of image segmentation system combining TC and swarm intelligence algorithm
  116. Particle swarm optimization and fuzzy C-means clustering algorithm for the adhesive layer defect detection
  117. Optimization of student learning status by instructional intervention decision-making techniques incorporating reinforcement learning
  118. Fuzzy model-based stabilization control and state estimation of nonlinear systems
  119. Optimization of distribution network scheduling based on BA and photovoltaic uncertainty
Downloaded on 23.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/nleng-2024-0058/html
Scroll to top button