Next Article in Journal
Generic Model of Max Heteroassociative Memory Robust to Acquisition Noise
Next Article in Special Issue
Hamiltonians of the Generalized Nonlinear Schrödinger Equations
Previous Article in Journal
ConvGRU-RMWP: A Regional Multi-Step Model for Wave Height Prediction
Previous Article in Special Issue
On Construction of Partially Dimension-Reduced Approximations for Nonstationary Nonlocal Problems of a Parabolic Type
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optical Solitons for the Concatenation Model with Differential Group Delay: Undetermined Coefficients

by
Anjan Biswas
1,2,3,4,5,
Jose Vega-Guzman
6,
Yakup Yıldırım
7,
Luminita Moraru
8,*,
Catalina Iticescu
8 and
Abdulah A. Alghamdi
2
1
Department of Mathematics and Physics, Grambling State University, Grambling, LA 71245, USA
2
Mathematical Modeling and Applied Computation (MMAC) Research Group, Department of Mathematics, King Abdulaziz University, Jeddah 21589, Saudi Arabia
3
Department of Applied Mathematics, National Research Nuclear University, 31 Kashirskoe Hwy, Moscow 115409, Russia
4
Department of Applied Sciences, Cross-Border Faculty of Humanities, Economics and Engineering, Dunarea de Jos University of Galati, 111 Domneasca Street, 800201 Galati, Romania
5
Department of Mathematics and Applied Mathematics, Sefako Makgatho Health Sciences University, Medunsa 0204, South Africa
6
Department of Mathematics, Lamar University, Beaumont, TX 77710, USA
7
Department of Computer Engineering, Biruni University, Istanbul 34010, Turkey
8
Department of Chemistry, Physics and Environment, Faculty of Sciences and Environment, Dunarea de Jos University of Galati, 47 Domneasca Street, 800008 Galati, Romania
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(9), 2012; https://doi.org/10.3390/math11092012
Submission received: 30 March 2023 / Revised: 23 April 2023 / Accepted: 23 April 2023 / Published: 24 April 2023

Abstract

:
In the current study, the concatenation model of birefringent fibers is explored for the first time, and we present optical soliton solutions to the model. The integration algorithm used to achieve this retrieval is the method of undetermined coefficients, which yields a wide range of soliton solutions. The parameter constraints arise naturally during the derivation of the soliton solutions, which are essential for such solitons to exist.

1. Introduction

While there exists an abundance of models that govern the propagation of solitons across transcontinental and transoceanic distances through an optical fiber, an interesting equation that was proposed about a decade ago in 2014 by Ankiewicz et al. is the concatenation model [1,2]. This model is a combination of three well-known equations that describe the dynamics of soliton transmission across intercontinental distances. They are the nonlinear Schrödinger Equation (NLSE), the Lakshmanan–Porsezian–Daniel (LPD) model, and the Sasa–Satsuma Equation (SSE). These three models are concatenated in sequence, hence the name.
Recently, the concatenation model has gained further importance and has attracted considerable attention from a wide range of perspectives. The soliton solutions were derived through the use of undetermined coefficients. The conservation laws were also enumerated after they were identified by the multiplier approach. The model was also addressed numerically with the Laplace–Adomian decomposition scheme. The trial equation approach was also implemented, and Painleve analysis was carried out [3,4,5,6,7]. These studies were all conducted for the scalar version of the model. It is now time to turn the page and move on.
The current paper addresses the concatenation model in birefringent fibers, again with the help of undetermined coefficients. A full spectrum of solitons is thus made available, along with the solvability conditions, which are listed as parameter constraints. After providing a brief introduction to the model in birefringent fibers, in the rest of the paper, we explain the derivation of these solitons in detail using the integration algorithm. The results are then presented with full transparency.

Governing Model

The concatenation model for polarization-preserving fibers can be stated as follows:
i u t + a u x x + b u 2 u + c 1 σ 1 u x x x x + σ 2 u x 2 u * + σ 3 u x 2 u + σ 4 u 2 u x x + σ 5 u 2 u x x * + σ 6 u 4 u + i c 2 σ 7 u x x x + σ 8 u 2 u x + σ 9 u 2 u x * = 0 .
Equation (1) is the concatenation model that is formed by joining three well-studied models in fiber optics, namely NLSE, the LPD equation, and the SSE. It must be noted that when c 1 = c 2 = 0 , we have the NLSE; while if c 1 = 0 , we recover the SSE; and when c 2 = 0 , we obtain the LPD model. Thus, Equation (1), as it stands, is the model that is generated by concatenating the three globally familiar models from nonlinear fiber optics.
For birefringent fibers, Equation (1) is split into two components as:
i q t + a 1 q x x + b 1 q 2 + c 1 r 2 q + c 11 σ 11 q x x x x + α 1 q x 2 + β 1 r x 2 q * + γ 1 q x 2 + λ 1 r x 2 q + δ 1 q 2 + ζ 1 r 2 q x x + μ 1 q 2 + ρ 1 r 2 q x x * + f 1 q 4 + g 1 q 2 r 2 + h 1 r 4 q + i c 21 σ 71 q x x x + η 1 q 2 + θ 1 r 2 q x + ε 1 q 2 + τ 1 r 2 q x * = 0 ,
and
i r t + a 2 r x x + b 2 r 2 + c 2 q 2 r + c 12 σ 12 r x x x x + α 2 r x 2 + β 2 q x 2 r * + γ 2 r x 2 + λ 2 q x 2 r + δ 2 r 2 + ζ 2 q 2 r x x + μ 2 r 2 + ρ 2 q 2 r x x * + f 2 r 4 + g 2 r 2 q 2 + h 2 q 4 r + i c 22 σ 72 r x x x + η 2 r 2 + θ 2 q 2 r x + ε 2 r 2 + τ 2 q 2 r x * = 0 .
Equations (2) and (3) represent the concatenation model split into two components for a birefringent fiber. For j = 1 , 2 , a j represents the chromatic dispersion along the two components, while b j accounts for self-phase modulation, and c j represents the cross-phase modulation. Then, σ 1 j represents the fourth-order dispersions along the two components. Next up, α j , β j , γ j , λ j , δ j , ζ j , μ j , ρ j , f j , g j , and h j are the respective split-ups of the coefficients σ 2 to σ 6 from the LPD model, along the two components for a birefringent fiber. σ 7 j represents the coefficients of the fourth-order dispersion along the components, while in Equation (1), this effect is represented by the coefficient of σ 7 . Finally, η j , θ j , ϵ j , and τ j are the components of soliton self-frequency shift along the two components of a birefringent fiber, which are designated by σ 8 and σ 9 in the SSE part.
The concatenation model of (2) and (3) describes the propagation of solitons in birefringent optical fibers. It is composed of three main equations, namely the NLSE, the LPD model, and the SSE. These three equations are concatenated or combined to form a more comprehensive model that can describe the behavior of solitons in a wide range of situations. The NLSE is a fundamental equation in nonlinear optics and describes the propagation of light in a medium with a nonlinear refractive index. It is used to model the propagation of solitons in fiber optic communication systems. The LPD model is a modified version of the NLSE that includes additional terms to account for the effects of birefringence, which is the property of a medium that causes light to split into two polarized components. This model is used to study the behavior of solitons in birefringent fibers. The SSE is another important equation in soliton theory that describes the propagation of solitons in dispersive media. It is a generalization of the NLSE and is used to study the behavior of solitons in media with higher-order dispersion.
By combining these three equations into a single model, the concatenation model provides a more complete and accurate description of the behavior of solitons in birefringent optical fibers. It has many applications in fiber optic communication systems, where solitons are used to transmit information over long distances with minimal distortion.
It must be noted that the preliminary model was introduced about a decade ago [1,2]. It is not yet known what the exact form of an optical fiber would be or where such a form would be applicable. With a conjecture that it would make sense for erbium-doped fiber, in the present study, the governing model (1) is split and considered with differential group delay for the first time. This is proposed as a purely analytical model before any laboratory testing.
The integration of the coupled system (2) and (3) requires the following hypothesis to be applied:
q ( x , t ) = P 1 ( x , t ) e i κ x + ω t + Θ ,
and
r ( x , t ) = P 2 ( x , t ) e i κ x + ω t + Θ .
On the basis of this hypothesis, the waveform is described by the function P ( x , t ) , which is unique for each type of soliton; Θ is a phase constant; the wave number is represented by ω ; and the soliton frequency is given by κ . Inserting these hypotheses into (2) and (3) leads to the real part identity
[ ω + κ 2 a l + κ 3 ( c 2 l σ 7 l κ c 1 l σ 1 l ) ] P l + [ b l κ c 2 l ( ε l η l ) κ 2 c 1 l ( α l γ l + δ l + μ l ) ] P l 3 + c 1 l f l P l 5 + c 1 l g l P l 3 P 2 2 + c 1 l h l P l P l ¯ 4 + [ c l κ ( κ c 1 l ( β l + ξ l + ρ l λ l ) + c 2 l ( τ l θ l ) ) ] P l P l ¯ 2 + c 1 l ( α l + γ l ) P l P l x 2 + c 1 l ( β l + λ l ) P l P l ¯ x 2 + c 1 l [ ( δ l + μ l ) P l 2 + ( ξ l + ρ l ) P l ¯ 2 ] 2 P l x 2 + [ a l + 3 κ ( c 2 l σ 7 l 2 κ c 1 l σ 1 l ) ] 2 P l x 2 + c 1 l σ 1 l 4 P l x 4 = 0 ,
while the imaginary counterpart reduces to
P l t [ 2 κ a l ( 4 k c 1 l σ 1 l 3 c 2 l σ 7 l ) κ 2 ] P l x + [ c 2 l ( ε l + η l ) 2 κ c 1 l ( α l + δ l μ l ) ] P l 2 P l x + [ c 2 l ( θ l + τ l ) 2 κ c 1 l ( ξ l ρ l ) ] P l ¯ 2 P l x 2 κ β l c 1 l P l P l ¯ P l ¯ x + ( c 2 l σ 7 l 4 κ c 1 l σ 1 l ) P l x 3 = 0 ,
where l ¯ = 3 l and l = 1 , 2 . Equation (7) shows that the speed can be retrieved as
v = 2 κ ( a l + c 2 l σ 7 l κ ) ,
after considering
c 2 l ( ε l + η l ) = 2 κ c 1 l ( α l + δ l μ l ) ,
c 2 l ( θ l + τ l ) = 2 κ c 1 l ( ξ ρ l ) ,
β l = 0 ,
c 2 l σ 7 l = 4 κ c 1 l σ 1 l .
Additionally, comparing the two values of the speed yields
a 1 = a 2 = a , c 21 = c 22 = c 7 , σ 71 = σ 72 = σ 7 .
Thus, as a direct consequence, the speed is rewritten as
v = 2 κ ( a + c 7 σ 7 κ ) ,
while in view of constraints (9)–(12), the real part of Equation (6) is recast as
( ω + a κ 2 + 3 c 1 l σ 1 l κ 4 ) P l + [ b l κ c 7 ( ε l η l ) c 1 l κ 2 ( α l γ l + δ l + μ l ) ] P l 3 + c 1 l f l P l 5 + c 1 l g l P l 3 P l ¯ 2 + c 1 l h l P l P l ¯ 4 + [ c l κ ( κ c 1 l ( ξ l + ρ l λ l ) + c 7 ( τ l θ l ) ) ] P l P l ¯ 2 + c 1 l ( α l + γ l ) P l P l x 2 + c 1 l λ l P l P l ¯ x 2 + c 1 l [ ( δ l + μ l ) P l 2 + ( ξ l + ρ l ) P l ¯ 2 ] 2 P l x 2 + ( a + 6 c 1 l σ 1 l κ 2 ) 2 P l x 2 + c 1 l σ 1 l 4 P l x 4 = 0 .
It must be noted from (8) that the speeds of the solitons along the two components of a birefringent fiber appear to be different. This is because the parameters and coefficients are chosen generally. Then, equating the two components of the soliton velocity yields a number of parameter constraints and relations given by (13), leading to the recovery of the single speed at which the soliton travels along the two components, as expressed by (14).
Throughout the next section, we will explore the integrability of the real part Equation (15) from the point of view of four different types of soliton solutions.

2. Soliton Solutions

Soliton solutions are a type of mathematical solution that describes a self-reinforcing solitary wave that maintains its shape and velocity over a significant distance. These solutions arise in various fields, including physics, engineering, and applied mathematics, where they have important applications in areas such as optical fiber communications, plasma physics, and water waves. The study of solitons began in the 19th century when Scottish engineer John Scott Russell observed a solitary wave on a canal in Scotland. It was not until the 1960s, however, that the theory of solitons began to take shape with the development of the NLSE, which is a mathematical model that describes the behavior of solitons. Today, researchers continue to study soliton solutions in various settings, developing new models and methods for their analysis and application. The study of solitons remains an active area of research, with many new discoveries and applications expected in the years to come.
We now proceed to secure optical solitons with the concatenation model with differential group delay ((2) and (3)). Four different solitons are explored in the next four subsections by studying the integrability of the real part Equation (15) according to the corresponding waveform ( P l ( x , t ) ).
The principle of undetermined coefficients is fairly simple. With an appropriate hypothesis of the soliton solutions that can be deciphered from the scalar version of the model, namely (1), the governing model is simplified after inserting this hypothesis into the two components. The real and imaginary components are separated. The imaginary component gives way to the velocity of the soliton, along with a few parameter constraints. The real part gives a relation from which the coefficients of the linearly independent functions are set to zero. This yields the relation between the soliton amplitude and width, along with a few more parameter constraints. Thus, a complete picture of the soliton, along with its essential parameter relations, is derived. It is applicable to bright, dark, and singular solitons. A transparent shortcoming of this approach is that it fails to recover soliton radiation and cannot be applied to derive multiple soliton solutions.

2.1. Bright Solitons

Bright solitons are a type of soliton solution that describes a self-reinforcing wave that has a maximum amplitude at its center. They are called “bright” solitons because they can be easily observed in optical systems, where they appear as bright pulses of light. Bright soliton solutions are important in many areas of physics, including fiber optics, plasma physics, and Bose–Einstein condensates. Bright soliton solutions are typically described by the NLSE, which is a mathematical model that takes into account the nonlinearity of the wave equation. These solutions are characterized by their ability to maintain their shape and velocity over long distances, and they are highly stable, making them ideal for use in telecommunications and other applications. The study of bright soliton solutions is an active area of research, with many new models and methods for their analysis and application. Researchers continue to explore the properties of these solutions in various settings, and new discoveries are expected to emerge in the years to come. Bright soliton solutions are a fascinating and important area of study in modern physics and applied mathematics.
The first case to be addressed is the bright soliton, this waveform of which typically has a structure of the form,
P l ( x , t ) = A l sech p l φ , φ = B ( x v t ) ,
where the soliton amplitudes A l ( l = 1 , 2 ) have a width of B, which implies that there are two soliton amplitudes with a specified width. The value of the unknown exponent ( p l ) can be determined using the principle of balance. By plugging the expression (16) into Equation (15), we obtain a simplified form of the equation,
[ p l 2 ( a l + 6 c 1 l σ 1 l κ 2 ) B 2 + p l 4 c 1 l σ 1 l B 2 ( ω + a κ 2 + 3 c 1 l σ 1 l κ 4 ) ] A l sech p l φ + [ ( b l κ c 7 ( ε l η l ) c 1 l κ 2 ( α l γ l + δ l + μ l ) ) + p l 2 c 1 l ( α l + γ l + δ l + μ l ) B 2 ] A l 3 sech 3 p l φ + f l c 1 l A l 5 sech 5 p l φ p p ( 1 + p l ) [ ( a + 6 c 1 l σ 1 l κ 2 ) + 2 c 1 l σ 1 l ( 2 + p l ( 2 + p l ) ) B 2 ] A l B 2 sech 2 + p l φ + p l ( 1 + p l ) ( 2 + p l ) ( 3 + p l ) c 1 l σ 1 l A l B 4 sech 4 + p l φ p l [ p l c 1 l ( α l + γ l ) + ( 1 + p l ) c 1 l ( δ l + μ l ) ] A l 3 B 2 sech 2 + 3 p l φ + { [ c l κ ( κ c 1 l ( ξ l + ρ l λ l ) + c 7 ( τ l θ l ) ) ] + p l 2 c 1 l ( ξ l + ρ l ) B 2 + p l ¯ 2 λ l c 1 l B 2 } A l A l ¯ 2 sech p l + 2 p l ¯ φ + g l c 1 l A l 3 A l ¯ 2 sech 3 p l + 2 p l ¯ φ + h l c 1 l A l A l ¯ 4 sech p l + 4 p l ¯ φ [ p l ( 1 + p l ) c 1 l ( ξ l + ρ l ) + p l ¯ 2 λ l c 1 l ] A l A l ¯ 2 B 2 sech 2 + p l + 2 p l ¯ φ = 0 ,
where, as previously defined, l ¯ = 3 l and l = 1 , 2 . On the last identity (17), balancing nonlinearity and second-order dispersion leads to the comparison of the exponents of sech 3 p l φ and sech p l + 2 φ , from which,
p l = 1 ,
for l = 1 , 2 . The substitution of (18) into (17) allows us to set the coefficients of the functions sech j φ to zero for j = 1 , 3 , 5 . This provides us a wave number that can be expressed as follows:
ω = a B 2 κ 2 + c 1 l σ 1 l B 4 + 6 κ 2 B 2 3 κ 4 ,
with inverse width described by
B = ( Z 4 + Z 5 ) A l 2 2 Z 7 + ( c 1 l λ l + Z 6 ) A l ¯ 2 40 c 1 l σ 1 l ± [ ( Z 4 + Z 5 ) A l 2 2 Z 7 + ( c 1 l λ l + Z 6 ) A l ¯ 2 ] 2 + 80 c 1 l σ 1 l ( Z 2 A l 2 + Z 3 A l ¯ 2 ) 40 c 1 l σ 1 l 1 2 ,
provided the radicand is non-negative and the identity
h l c 1 l A l ¯ 4 + f l c 1 l A l 4 + 24 c 1 l σ 1 l B 4 = Z 4 + 2 Z 5 B 2 g l c 1 l A l ¯ 2 A l 2 + 2 Z 6 + c 1 l λ l A l ¯ 2 B 2 .
For simplicity, we adopted the following notations:
Z 1 = ( ω + a κ 2 + 3 c 1 l σ 1 l κ 4 ) ,
Z 2 = [ b l κ c 7 ( ε l η l ) c 1 l κ 2 ( α l γ l + δ l + μ l ) ] ,
Z 3 = [ c l κ ( κ c 1 l ( ξ l + ρ l λ l ) + c 7 ( τ l θ l ) ) ] ,
Z 4 = c 1 l ( α l + γ l ) ,
Z 5 = c 1 l ( δ l + μ l ) ,
Z 6 = c 1 l ( ξ l + ρ l ) ,
Z 7 = ( a + 6 c 1 l σ 1 l κ 2 ) .
Hence, we can derive the optical soliton solution for the concatenated model ((2) and (3)) as follows:
q ( x , t ) = A 1 sech [ B ( x v t ) ] e i κ x + ω t + Θ ,
and
r ( x , t ) = A 2 sech [ B ( x v t ) ] e i κ x + ω t + Θ .
Above, we discussed the associated constraints for the wave numbers, speed, inverse width, and amplitudes. Plots of the bright solitons described by Equations (29) and (30) with A 1 = 1 , A 2 = 1 , κ = 1 , a = 1 , c 7 = 1 , σ 7 = 1 , c 1 l = 1 , σ 1 l = 1 , λ l = 1 , α l = 1 , γ l = 1 , δ l = 1 , μ l = 1 , ξ l = 1 , ρ l = 1 , b l = 1 , ε l = 1 , η l = 1 , c l = 1 , τ l = 1 , and θ l = 1 are illustrated in Figure 1. Bright solitons are a type of soliton solution that describes a self-reinforcing wave that has a maximum amplitude at its center.

2.2. Dark Solitons

Dark solitons are another type of soliton solution that describes a self-reinforcing wave with a minimum amplitude at its center. They are called “dark” solitons because they appear as a localized dip or dark region in the waveform. Dark soliton solutions arise in many physical systems, including nonlinear optics, Bose–Einstein condensates, and ocean waves. Unlike bright solitons, dark solitons are typically unstable, meaning they eventually break up or decay over long distances. However, they are still of great interest to researchers because of their unique properties and potential applications. For example, they can be used to manipulate and control light in optical systems, and they can provide insight into the behavior of complex systems. The study of dark soliton solutions is an active area of research, with many new models and methods for their analysis and application. Researchers continue to explore the properties of these solutions in various settings, and new discoveries are expected to emerge in the years to come. Dark soliton solutions are an important area of study in modern physics and applied mathematics, with potential applications in a wide range of fields.
The second case to be considered is dark solitons, the wave form assumption ( P l ( x , t ) ) of which is
P l ( x , t ) = A l tanh p l φ , φ = B ( x v t ) .
Here, the soliton speed is represented by v, and A l and B are considered free parameters. The parameter p l is determined using the principle of balancing. By plugging (31) into (15), we arrive at the simplified expression:
{ ( ω + a κ 2 + 3 c 1 l σ 1 l κ 4 ) + 2 p l 2 [ ( a + 6 c 1 l σ 1 l κ 2 ) c 1 l σ 1 l ( 5 + 3 p l 2 ) B 2 ] B 2 } A l tanh p l φ + { [ b l κ c 7 ( ε l η l ) c 1 l κ 2 ( α l γ l + δ l + μ l ) ] 2 p l 2 c 1 l ( α l + γ l + δ l + μ l ) B 2 } A l 3 tanh 3 p l φ + f l c 1 l A l 5 tanh 5 p l φ + p l ( 1 + p l ) [ ( a + 6 c 1 l σ 1 l κ 2 ) 4 c 1 l σ 1 l ( 2 + p l ( 2 + p l ) ) B 2 ] A l B 2 tanh 2 + p l φ + c 1 l σ 1 l p l ( 1 + p l ) ( 2 + p l ) ( 3 + p l ) A l B 4 tanh 4 + p l φ + p l [ p l c 1 l ( α l + γ l ) + ( 1 + p l ) c 1 l ( δ l + μ l ) ] A l 3 B 2 tanh 2 + 3 p l φ p p p p + { [ c l κ ( κ c 1 l ( ξ l + ρ l λ l ) + c 7 ( τ l θ l ) ) ] ( p l 2 c 1 l ( ξ l + ρ l ) + λ l p l ¯ 2 c 1 l ) B 2 } A l A l ¯ 2 tanh p l + 2 p l ¯ φ + h l c 1 l A l A l ¯ 4 tanh p l + 4 p l ¯ φ + g l c 1 l A l 3 A l ¯ 2 tanh 3 p l + 2 p l ¯ φ p l ( 1 p l ) [ ( a + 6 c 1 l σ 1 l κ 2 ) 4 c 1 l σ 1 l ( 2 p l ( 2 p l ) ) B 2 ] A l B 2 tanh p l 2 φ p l ( 1 p l ) ( 2 p l ) ( 3 p l ) c 1 l σ 1 l A l B 4 tanh p l 4 φ + p l [ p l c 1 l ( α l + γ l ) ( 1 p l ) c 1 l ( δ l + μ l ) ] A l 3 B 2 tanh 3 p l 2 φ + [ p l ( p l 1 ) c 1 l ( ξ l + ρ l ) + p l ¯ 2 λ l c 1 l ] A l A l ¯ 2 B 2 tanh p l + 2 p l ¯ 2 φ + [ p l ( 1 + p l ) c 1 l ( ξ l + ρ l ) + p l ¯ 2 λ l c 1 l ] A l A l ¯ 2 B 2 tanh p l + 2 p l ¯ + 2 φ = 0 ,
where, as in the previous section, l ¯ = 3 l and l = 1 , 2 . By inspection, the parameter p l in this case can be verified to result in (18) for both l = 1 and l = 2 . Substituting the obtained value of p l into Equation (32) and setting the coefficients of the functions tanh j φ with j = 1 , 3 , 5 to zero results in the following expression for the wave number:
ω = a κ 2 + 2 B 2 + c 1 l α l + γ l A l 2 + λ l A l ¯ 2 B 2 + c 1 l σ 1 l 16 B 4 12 κ 2 B 2 3 κ 4 ,
and the free parameter,
B = Z 7 ( Z 4 + Z 5 ) A l 2 ( Z 6 + c 1 l λ l ) A l ¯ 2 40 c 1 l σ 1 l ± [ ( Z 4 + Z 5 ) A l 2 Z 7 + ( c 1 l λ l + Z 6 ) A l ¯ 2 ] 2 + 40 c 1 l σ 1 l ( Z 2 A l 2 + Z 3 A l ¯ 2 ) 40 c 1 l σ 1 l 1 2 ,
along with the solvability condition,
c 1 l f l A l 4 + c 1 l h l A l ¯ 4 + Z 4 + 2 Z 5 B 2 + c 1 l g l A l ¯ 2 A l 2 + 2 Z 6 + c 1 l λ l A l ¯ 2 B 2 + 24 c 1 l σ 1 l B 4 = 0 .
For convenience, we adopted the notations previously given on (22)–(28). It is clear from (34) that
c 1 l σ 1 l 0 ,
and that
Z 4 + Z 5 A l 2 Z 7 + c 1 l λ l + Z 6 A l ¯ 2 2 + 40 c 1 l σ 1 l Z 2 A l 2 + Z 3 A l ¯ 2 0 ,
along with the fact that the entire radicand on (34) must be positive. Consequently, the solution for the dark soliton in the considered system can be expressed as follows:
q ( x , t ) = A 1 tanh [ B ( x v t ) ] e i κ x + ω t + Θ ,
and
r ( x , t ) = A 2 tanh [ B ( x v t ) ] e i κ x + ω t + Θ .
Above, described the relationship between the respective free parameters, along with the speed (14) and wave number (33), subject to the constraints discussed above. Plots of the dark solitons described by Equations (38) and (39) with A 1 = 1 , A 2 = 1 , κ = 1 , a = 1 , c 7 = 1 , σ 7 = 1 , c 1 l = 1 , σ 1 l = 1 , λ l = 1 , α l = 1 , γ l = 1 , δ l = 1 , μ l = 1 , ξ l = 1 , ρ l = 1 , b l = 1 , ε l = 1 , η l = 1 , c l = 1 , τ l = 1 , and θ l = 1 are demonstrated in Figure 2. Dark solitons are a type of soliton solution that describes a self-reinforcing wave with a minimum amplitude at its center.

2.3. Singular Solitons

Singular solitons are a type of soliton solution that exhibits a singularity at the center of the wave profile. They are also referred to as peakons because they resemble a peaked function rather than a smooth wave. Singular solitons arise in various fields, including fluid dynamics, nonlinear optics, and plasma physics. Unlike bright and dark solitons, which have a finite energy and amplitude, singular solitons have infinite amplitude at the center of the wave profile. As a result, these solutions are typically highly unstable and cannot propagate over long distances. However, they are still of great interest to researchers because of their unique properties and potential applications. The study of singular soliton solutions is an active area of research, with many new models and methods for their analysis and application. Researchers continue to explore the properties of these solutions in various settings, and new discoveries are expected to emerge in the years to come. Singular soliton solutions are a fascinating and important area of study in modern physics and applied mathematics, with potential applications in fields such as shock wave theory and signal processing.

2.3.1. Singular Solitons (Type I)

In this case, we need to make a certain assumption regarding the waveform, which is:
P l ( x , t ) = A l csch p l φ , φ = B ( x v t ) .
Let us assume that A l and B are independent variables that are considered as free parameters. By substituting this assumption into (15), we obtain the following result:
{ p l 2 ( a + 6 c 1 l σ 1 l κ 2 ) B 2 + p l 4 c 1 l σ 1 l B 4 ( ω + a κ 2 + 3 c 1 l σ 1 l κ 4 ) } A l csch p l φ + { [ b l κ c 7 ( ε l η l ) c 1 l κ 2 ( α l γ l + δ l + μ l ) ] + p l 2 c 1 l ( α l + γ l + δ l + μ l ) B 2 } A l 3 csch 3 p l φ + f l c 1 l A l 5 csch 5 p l φ + p l ( 1 + p l ) [ ( a + 6 c 1 l σ 1 l κ 2 ) + 2 ( 2 + p l ( 2 + p l ) ) c 1 l σ 1 l B 2 ] A l B 2 csch 2 + p l φ + p l ( 1 + p l ) ( 2 + p l ) ( 3 + p l ) c 1 l σ 1 l A l B 4 csch 4 + p l φ + p l [ p l c 1 l ( α l + γ l ) + ( 1 + p l ) c 1 l ( δ l + μ l ) ] A l 3 B 2 csch 2 + 3 p l φ + { [ c l κ ( κ c 1 l ( ξ l + ρ l λ l ) + c 7 ( τ l θ l ) ) ] + ( p l ¯ 2 λ l c 1 l + p l 2 c 1 l ( ξ l + ρ l ) ) B 2 } A l A l ¯ 2 csch p l + 2 p l ¯ φ + g l c 1 l A l 3 A l ¯ 2 csch 3 p l + 2 p l ¯ φ + h l c 1 l A l A l ¯ 4 csch p l + 4 p l ¯ φ + [ p l ¯ 2 λ l c 1 l + p l ( 1 + p l ) c 1 l ( ξ l + ρ l ) ] A l A l ¯ 2 B 2 csch 2 + p l + 2 p l ¯ φ = 0 ,
after simplification, where l ¯ = 3 l and l = 1 , 2 . Balancing the cubic nonlinearity and second-order dispersion in the final equation results in the determination of p l as presented in Equation (18). We then substitute the obtained value of p l and, following the same procedure as for dark and bright solitons, gather the coefficients of csch j φ for j = 1 , 3 , 5 , which results in the same wave number as that of the dark soliton solution given by Equation (19). The free parameter (B) is also determined by
B = ( Z 4 + Z 5 ) A l 2 Z 7 ( c 1 l λ l + Z 6 ) A l ¯ 2 40 c 1 l σ 1 l ± ± [ ( Z 4 + Z 5 ) A l 2 + 2 Z 7 + ( c 1 l λ l + Z 6 ) A l ¯ 2 ] 2 80 c 1 l σ 1 l ( Z 2 A l 2 + Z 3 A l ¯ 2 ) 40 c 1 l σ 1 l 1 2 ,
while the identity equation is the same as in the dark soliton case (35). Here, we again adopted the expressions given in (22)–(28). Clearly, the constraint (36) also applies here, and (37) slightly changes to
Z 4 + Z 5 A l 2 + 2 Z 7 + c 1 l λ l + Z 6 A l ¯ 2 2 80 c 1 l σ 1 l Z 2 A l 2 + Z 3 A l ¯ 2 0 .
As a result, the singular soliton solutions for the system described by Equations (2) and (3) are introduced as below
q ( x , t ) = A 1 csch [ B ( x v t ) ] e i κ x + ω t + Θ ,
and
r ( x , t ) = A 2 csch [ B ( x v t ) ] e i κ x + ω t + Θ .
As previously discussed, the relevant parameters and their corresponding constraints have been outlined above. Plots of the singular solitons described by Equations (44) and (45) with A 1 = 1 , A 2 = 1 , κ = 1 , a = 1 , c 7 = 1 , σ 7 = 1 , c 1 l = 1 , σ 1 l = 1 , λ l = 1 , α l = 1 , γ l = 1 , δ l = 1 , μ l = 1 , ξ l = 1 , ρ l = 1 , b l = 1 , ε l = 1 , η l = 1 , c l = 1 , τ l = 1 , and θ l = 1 are depicted in Figure 3. Singular solitons are a type of soliton solution that exhibits a singularity at the center of the wave profile.

2.3.2. Singular Solitons (Type II)

In this case, the assumption for the waveform portion ( P l ( x , t ) ) is
P l ( x , t ) = A l coth p l φ , φ = B ( x v t ) .
Similarly, the parameters A l and B are considered to be free parameters. Substituting them into Equation (15) yields the following relation
{ ( ω + a κ 2 + 3 c 1 l σ 1 l κ 4 ) + 2 p l 2 [ ( a + 6 c 1 l σ 1 l κ 2 ) ( 5 + 3 p l 2 ) B 2 ] B 2 } A l coth p l φ + { [ b l κ c 7 ( ε l η l ) c 1 l κ 2 ( α l γ l + δ l + μ l ) ] 2 p l 2 c 1 l ( α l + γ l + δ l + μ l ) B 2 } A l 3 coth 3 p l φ + p l ( 1 + p l ) [ ( a + 6 c 1 l σ 1 l κ 2 ) 4 ( 2 + p l ( 2 + p l ) ) c 1 l σ 1 l B 2 ] A l B 2 coth 2 + p l φ + f l c 1 l A l 5 coth 5 p l φ + p l ( 1 + p l ) ( 2 + p l ) ( 3 + p l ) c 1 l σ 1 l A l B 4 coth 4 + p l φ + p l [ p l c 1 l ( α l + γ l ) + ( 1 + p l ) c 1 l ( δ l + μ l ) ] A l 3 B 2 coth 2 + 3 p l φ p l ( 1 p l ) [ ( a + 6 c 1 l σ 1 l κ 2 ) 4 ( 2 p l ( 2 p l ) ) c 1 l σ 1 l B 2 ] A l B 2 coth 2 + p l φ p l ( 1 p l ) ( 2 p l ) ( 3 p l ) c 1 l σ 1 l A l B 4 coth 4 + p l φ + p l [ p l c 1 l ( α l + γ l ) ( 1 p l ) c 1 l ( δ l + μ l ) ] A l 3 B 2 coth 2 + 3 p l φ + h l c 1 l A l A l ¯ 4 coth p l + 4 p l ¯ φ + g l c 1 l A l 3 A l ¯ 2 coth 3 p l + 2 p l ¯ φ + [ p l ¯ 2 λ l c 1 l + p l ( 1 + p l ) c 1 l ( ξ l + ρ l ) ] A l A l ¯ 2 B 2 coth 2 + p l + 2 p l ¯ φ + [ p l ¯ 2 λ l c 1 l p l ( 1 p l ) c 1 l ( ξ l + ρ l ) ] A l A l ¯ 2 B 2 coth 2 + p l + 2 p l ¯ φ + { [ c l κ ( κ c 1 l ( ξ l + ρ l λ l ) + c 7 ( τ l θ l ) ) ] 2 [ p l ¯ 2 λ l c 1 l + p l 2 c 1 l ( ξ l + ρ l ) ] B 2 } A l A l ¯ 2 coth p l + 2 p l ¯ φ = 0 .
In this case, the balancing algorithm and coefficients of the standalone elements, specifically coth 2 + p l φ and coth 4 + p l φ , provide us the value of p l as presented in Equation (18). Upon substituting the obtained value of p l from Equation (18) into Equation (47), the resulting solutions and corresponding solvability conditions are identical to those obtained for the case of a dark soliton, as shown in Equations (33)–(37). Consequently, the concatenation coupled model given by Equations (2) and (3) possesses type-II singular soliton solutions, which are introduced below
q ( x , t ) = A 1 coth [ B ( x v t ) ] e i κ x + ω t + Θ ,
and
r ( x , t ) = A 2 coth [ B ( x v t ) ] e i κ x + ω t + Θ ,
where all the relevant parameters and their corresponding constraints for these solitons are identical to those for dark solitons. Plots of the singular solitons described by Equations (48) and (49) with A 1 = 1 , A 2 = 1 , κ = 1 , a = 1 , c 7 = 1 , σ 7 = 1 , c 1 l = 1 , σ 1 l = 1 , λ l = 1 , α l = 1 , γ l = 1 , δ l = 1 , μ l = 1 , ξ l = 1 , ρ l = 1 , b l = 1 , ε l = 1 , η l = 1 , c l = 1 , τ l = 1 , and θ l = 1 are presented in Figure 4. Soliton solutions that display a singularity at the center of the wave profile are referred to as singular solitons.

3. Conclusions

In the current, we paper addressed optical solitons in a birefringent fiber modeled by the concatenation-type governing equation. Three standard equations are concatenated to formulate the model. They are the most familiar NLSE, the LPD model, and the SSE. The model of study was integrated using the principle of undetermined coefficients. This resulted in the retrieval of bright, dark, and singular (both kinds) one-soliton solutions to the model, along with the applicable constraint conditions that naturally emerged during the course of the derivation of the soliton solution.
While this is a simple approach, it has its merits and shortcomings. It is simple in the sense that after picking an appropriate hypothesis and inserting it into the equation, the parameter relations and constraints naturally emerge. The imaginary part relation gives the soliton velocity and additional constraint relations that must remain valid for the solitons to exist. One of the shortcomings is that it fails to recover multisoliton solutions, unlike inverse scattering transform (IST) or the Hirota’s bilinear approach. Another disadvantage is it is unable to go past the discrete regime and provide an answer for the continuous regime. In other words, the algorithm does not yield soliton radiation. In order to recover soliton radiation, one must implement additional approaches such as IST or, beyond all order asymptotics, even applying the theory of unfoldings. Another disadvantage is that the method fails to obtain additional soliton solutions such as the straddled solitons. However, these can be recovered with the aid of Kudryashov’s approach.
With the fundamental results in place, various avenues are opened, presenting opportunities for future studies. First, an immediate burning question would be, “What are the conservation laws?”. Such laws will be investigated in a separate paper, with conserved densities derived using the multipliers approach and listed after computing the conserved quantities from the bright one-soliton solution derived in this paper.
The model proposed in this paper involves differential group delay, while in a previous report, the same model was investigated but for the scalar case [3]. Dispersion-flattened fibers can be explored in future work. This would provide a newer perspective to extract the soliton solutions and provide a new perspective of possible conservation laws.
The inclusion of perturbation terms is our next target. Thus, the perturbed concatenation model will be analyzed for its integrability. The availability of conservation laws would mean that one can compute the adiabatic parameter dynamics of the soliton parameters via the soliton perturbation theory. One can also implement the moment method or the collective variables approach to recover this kind of dynamical system.
The integrability of the perturbed model also needs to be studied to recover the soliton solution when the concatenation model includes perturbation terms. This depends on two situations. If the perturbation terms are weak, with quasimonochromaticity, one must integrate the perturbed concatenation model approximately with the use of multiscale perturbation analysis, irrespective of the type of perturbation term, be it Hamiltonian or non-Hamiltonian. Thus, a quasistationary soliton solution would be revealed, along with the resonance conditions. However, if the perturbation terms are strong but of Hamiltonian type, several integration schemes are available to integrate the perturbed concatenation model. A few such schemes are Kudryashov’s approach, Jacobi’s elliptic function scheme, and the sin-Gordon equation scheme, just to mention a few. In the final situation, when the perturbation terms are strong but non-Hamiltonian type, a slick scheme needs to be implemented that considers the model with time-dependent coefficients and applies an integration scheme possibly including the method of undetermined coefficients. This would naturally apply constraints on the parameter and the time-dependent coefficients. Some of the non-Hamiltonian perturbations to which this scheme can be successfully applied are the linear attenuation effect and Raman scattering, among others. All of these possibilities are subject to one exception. If the perturbation term(s) contain(s) maximum intensity, all of the above-mentioned integration approaches fail miserably. In that case, one must resort to approximate integration with the aid of the semi-inverse variational principle. In this case, it is only the bright one-soliton solution that can be recovered, and the retrieval of dark or singular solitons is not possible because the stationary integral is rendered divergent.
With adiabatic parameter dynamics in place, which can be recovered with the use of soliton perturbation theory, one can address the suppression of intrachannel collision of solitons with the aid of quasiparticle theory, which is yet to be developed for the concatenation model.
Then, the application of the model to additional devices apart from optical fibers can be considered. This would involve the implementation of the model for magneto-optic wave guides, Bragg grating fibers, and optical couplers, as well as for the study the concatenation model in optical metamaterials. Thus, gap solitons and magneto-optic solitons would be subsequently recovered and reported. Moreover, the governing model should addressed with fractional temporal evolution rather than linear temporal evolution as in the current paper and most previous papers. The consideration of fractional temporal evolution would lead to the slow evolution of solitons, meaning that one can successfully control the Internet bottleneck effect, which is a growing problem in the telecommunication industry. Another viable approach to this problem is to introduce the spatiotemporal dispersion (STD), in addition to pre-existing chromatic dispersion (CD). This combination if CD and STD can also be used as an engineering marvel to mitigate the Internet bottleneck effect.
One of the viable aspects to further consider with respect to this concatenation model is the nonlinearity of the CD. This can occur due to rough handling of fibers, as well as the random injection of pulses at the initial end of the fiber. When the CD becomes nonlinear, soliton transmission across intercontinental distance is stalled, and this leads to the formation of quiescent solitons—a very unwanted feature. While this aspect has already been studied for the scalar case, it must be addressed for fibers with differential group delay and dispersion-flattened fibers.
Apart from the various considerations for the analytical approaches listed thus far, one must additionally consider a model with the use of numerical algorithms. Some of the most applicable numerical schemes are the Adomian decomposition method, the improved Adomian decomposition approach, the Laplace–Adomian decomposition scheme, and the variational iteration method. Several other numerical algorithms that can be implemented to address the problem include the finite difference method, the finite element method, the finite volume scheme, and spectral approaches [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35]. Thus, much work remains to be done in this field.

Author Contributions

Conceptualization, A.B.; methodology, J.V.-G.; software, Y.Y.; writing—original draft preparation, L.M.; writing—review and editing, C.I.; project administration, A.A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the anonymous referees whose comments helped to improve this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ankiewicz, A.; Akhmediev, N. Higher-order integrable evolution equation and its soliton solutions. Phys. Lett. A 2014, 378, 358–361. [Google Scholar] [CrossRef]
  2. Ankiewicz, A.; Wang, Y.; Wabnitz, S.; Akhmediev, N. Extended nonlinear Schrödinger equation with higher-order odd and even terms and its rogue wave solutions. Phys. Rev. E 2014, 89, 012907. [Google Scholar] [CrossRef] [PubMed]
  3. Biswas, A.; Vega-Guzman, J.; Kara, A.H.; Khan, S.; Triki, H.; González-Gaxiola, O.; Moraru, L.; Georgescu, P.L. Optical solitons and conservation laws for the concatenation model: Undetermined coefficients and multipliers approach. Universe 2023, 9, 15. [Google Scholar] [CrossRef]
  4. Triki, H.; Sun, Y.; Zhou, Q.; Biswas, A.; Yıldırım, Y.; Alshehri, H.M. Dark solitary pulses and moving fronts in an optical medium with the higher-order dispersive and nonlinear effects. Chaos Solitons Fractals 2022, 164, 112622. [Google Scholar] [CrossRef]
  5. Wang, M.Y.; Biswas, A.; Yıldırım, Y.; Moraru, L.; Moldovanu, S.; Alshehri, H.M. Optical Solitons for a Concatenation Model by Trial Equation Approach. Electronics 2023, 12, 19. [Google Scholar] [CrossRef]
  6. Yıldırım, Y.; Biswas, A.; Moraru, L.; Alghamdi, A.A. Quiescent Optical Solitons for the Concatenation Model with Nonlinear Chromatic Dispersion. Mathematics 2023, 11, 1709. [Google Scholar] [CrossRef]
  7. Kudryashov, N.A.; Biswas, A.; Borodina, A.G.; Yıldırım, Y.; Alshehri, H.M. Painlevé analysis and optical solitons for a concatenated model. Optik 2023, 272, 170255. [Google Scholar] [CrossRef]
  8. Kudryashov, N.A. First integrals and general solution of the traveling wave reduction for Schrödinger equation with anti-cubic nonlinearity. Optik 2019, 185, 665–671. [Google Scholar] [CrossRef]
  9. Kudryashov, N.A. Mathematical model of propagation pulse in optical fiber with power nonlinearities. Optik 2020, 212, 164750. [Google Scholar] [CrossRef]
  10. Kudryashov, N.A. A generalized model for description of propagation pulses in optical fiber. Optik 2019, 189, 42–52. [Google Scholar] [CrossRef]
  11. Kudryashov, N.A.; Antonova, E.V. Solitary waves of equation for propagation pulse with power nonlinearities. Optik 2020, 217, 164881. [Google Scholar] [CrossRef]
  12. Kudryashov, N.A. Optical solitons of the model with arbitrary refractive index. Optik 2020, 224, 165767. [Google Scholar] [CrossRef]
  13. Kudryashov, N.A. Highly dispersive optical solitons of equation with various polynomial nonlinearity law. Chaos Solitons Fractals 2020, 140, 110202. [Google Scholar] [CrossRef]
  14. Kudryashov, N.A. Optical solitons of model with integrable equation for wave packet envelope. Chaos Solitons Fractals 2020, 141, 110325. [Google Scholar] [CrossRef]
  15. Kudryashov, N.A. Periodic and solitary waves in optical fiber Bragg gratings with dispersive reflectivity. Chin. J. Phys. 2020, 66, 401–405. [Google Scholar] [CrossRef]
  16. Kudryashov, N.A. Solitary waves of the generalized Sasa–Satsuma equation with arbitrary refractive index. Optik 2021, 232, 166540. [Google Scholar] [CrossRef]
  17. Kudryashov, N.A. Revised results of Khalida Bibi on the Radhakrishnan–Kundu–Lakshmanan equation. Optik 2021, 240, 166898. [Google Scholar] [CrossRef]
  18. Kudryashov, N.A. The Lakshmanan–Porsezian–Daniel model with arbitrary refractive index and its solution. Optik 2021, 241, 167043. [Google Scholar] [CrossRef]
  19. Vega-Guzman, J.; Biswas, A.; Mahmood, M.F.; Zhou, Q.; Moshokoa, S.P.; Belic, M. Optical solitons with polarization mode dispersion for Lakshmanan–Porsezian–Daniel model by the method of undetermined coefficients. Optik 2018, 171, 114–119. [Google Scholar] [CrossRef]
  20. Yıldırım, Y.; Topkara, E.; Biswas, A.; Triki, H.; Ekici, M.; Guggilla, P.; Khan, S.; Belic, M.R. Cubic–quartic optical soliton perturbation with Lakshmanan–Porsezian–Daniel model by sine-Gordon equation approach. J. Opt. 2021, 50, 322–329. [Google Scholar] [CrossRef]
  21. Columbu, A.; Frassu, S.; Viglialoro, G. Refined criteria toward boundedness in an attraction–repulsion chemotaxis system with nonlinear productions. Appl. Anal. 2023. [Google Scholar] [CrossRef]
  22. Frassu, S.; Li, T.; Viglialoro, G. Improvements and generalizations of results concerning attraction-repulsion chemotaxis models. Math. Methods Appl. Sci. 2022, 45, 11067–11078. [Google Scholar] [CrossRef]
  23. Triki, H.; Pan, A.; Zhou, Q. Pure-quartic solitons in presence of weak nonlocality. Phys. Lett. A 2023, 459, 128608. [Google Scholar] [CrossRef]
  24. Wang, H.; Li, X.; Zhou, Q.; Liu, W. Dynamics and spectral analysis of optical rogue waves for a coupled nonlinear Schrödinger equation applicable to pulse propagation in isotropic media. Chaos Solitons Fractals 2023, 166, 112924. [Google Scholar] [CrossRef]
  25. Zhou, Q.; Sun, Y.; Triki, H.; Zhong, Y.; Zeng, Z.; Mirzazadeh, M. Study on propagation properties of one-soliton in a multimode fiber with higher-order effects. Results Phys. 2022, 41, 105898. [Google Scholar] [CrossRef]
  26. Wang, T.Y.; Zhou, Q.; Liu, W.J. Soliton fusion and fission for the high-order coupled nonlinear Schrödinger system in fiber lasers. Chin. Phys. B 2022, 31, 020501. [Google Scholar] [CrossRef]
  27. Kudryashov, N.A. Mathematical model with unrestricted dispersion and polynomial nonlinearity. Appl. Math. Lett. 2023, 138, 108519. [Google Scholar] [CrossRef]
  28. Kudryashov, N.A. Embedded Solitons of the Generalized Nonlinear Schrödinger Equation with High Dispersion. Regul. Chaotic Dyn. 2022, 27, 680–696. [Google Scholar] [CrossRef]
  29. Kudryashov, N.A. Optical solitons of the generalized nonlinear Schrödinger equation with Kerr nonlinearity and dispersion of unrestricted order. Mathematics 2022, 10, 3409. [Google Scholar] [CrossRef]
  30. Kudryashov, N.A. First integrals and general solution of the complex Ginzburg-Landau equation. Appl. Math. Comput. 2020, 386, 125407. [Google Scholar] [CrossRef]
  31. Kudryashov, N.A. On Integrability of the FitzHugh Rinzel Model. Russ. J. Nonlinear Dyn. 2019, 15, 13–19. [Google Scholar] [CrossRef]
  32. Ekici, M. Stationary optical solitons with Kudryashov’s quintuple power law nonlinearity by extended Jacobi’s elliptic function expansion. J. Nonlinear Opt. Phys. Mater. 2023, 32, 2350008. [Google Scholar] [CrossRef]
  33. Ekici, M. Stationary optical solitons with complex Ginzburg–Landau equation having nonlinear chromatic dispersion and Kudryashov’s refractive index structures. Phys. Lett. A 2022, 440, 128146. [Google Scholar] [CrossRef]
  34. Ekici, M. Kinky breathers, W-shaped and multi-peak soliton interactions for Kudryashov’s quintuple power-law coupled with dual form of non-local refractive index structure. Chaos Solitons Fractals 2022, 159, 112172. [Google Scholar] [CrossRef]
  35. Ekici, M. Optical solitons with Kudryashov’s quintuple power–law coupled with dual form of non–local law of refractive index with extended Jacobi’s elliptic function. Opt. Quantum Electron. 2022, 54, 279. [Google Scholar] [CrossRef]
Figure 1. Profiles of a bright 1-soliton along the two components of a birefringent fiber after neglecting the radiative effects.
Figure 1. Profiles of a bright 1-soliton along the two components of a birefringent fiber after neglecting the radiative effects.
Mathematics 11 02012 g001
Figure 2. Profiles of a dark 1-soliton along the two components of a birefringent fiber after neglecting the radiative effects.
Figure 2. Profiles of a dark 1-soliton along the two components of a birefringent fiber after neglecting the radiative effects.
Mathematics 11 02012 g002
Figure 3. Profiles of a singular 1-soliton along the two components of a birefringent fiber with a clear visibility of the blow-up effects.
Figure 3. Profiles of a singular 1-soliton along the two components of a birefringent fiber with a clear visibility of the blow-up effects.
Mathematics 11 02012 g003
Figure 4. Characteristics of singular solitons.
Figure 4. Characteristics of singular solitons.
Mathematics 11 02012 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Biswas, A.; Vega-Guzman, J.; Yıldırım, Y.; Moraru, L.; Iticescu, C.; Alghamdi, A.A. Optical Solitons for the Concatenation Model with Differential Group Delay: Undetermined Coefficients. Mathematics 2023, 11, 2012. https://doi.org/10.3390/math11092012

AMA Style

Biswas A, Vega-Guzman J, Yıldırım Y, Moraru L, Iticescu C, Alghamdi AA. Optical Solitons for the Concatenation Model with Differential Group Delay: Undetermined Coefficients. Mathematics. 2023; 11(9):2012. https://doi.org/10.3390/math11092012

Chicago/Turabian Style

Biswas, Anjan, Jose Vega-Guzman, Yakup Yıldırım, Luminita Moraru, Catalina Iticescu, and Abdulah A. Alghamdi. 2023. "Optical Solitons for the Concatenation Model with Differential Group Delay: Undetermined Coefficients" Mathematics 11, no. 9: 2012. https://doi.org/10.3390/math11092012

APA Style

Biswas, A., Vega-Guzman, J., Yıldırım, Y., Moraru, L., Iticescu, C., & Alghamdi, A. A. (2023). Optical Solitons for the Concatenation Model with Differential Group Delay: Undetermined Coefficients. Mathematics, 11(9), 2012. https://doi.org/10.3390/math11092012

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop