Next Article in Journal
Shrinkage and Redundant Feature Elimination Network-Based Robust Image Zero-Watermarking
Previous Article in Journal
Optimization Design and Parameter Analysis of a Wheel with Array Magnets
Previous Article in Special Issue
Propagation and Parametric Amplification in Four-Wave Mixing Processes: Intramolecular Coupling and High-Order Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gap Solitons in Fiber Bragg Gratings Having Polynomial Law of Nonlinear Refractive Index and Cubic–Quartic Dispersive Reflectivity by Lie Symmetry

1
Department of Mathematics and Statistics, Central University of Punjab, Bathinda 151401, Punjab, India
2
Department of Mathematics and Physics, Grambling State University, Grambling, LA 71245, USA
3
Mathematical Modeling and Applied Computation (MMAC) Research Group, Department of Mathematics, King Abdulaziz University, Jeddah 21589, Saudi Arabia
4
Department of Applied Mathematics, National Research Nuclear University, 31 Kashirskoe Hwy, Moscow 115409, Russia
5
Department of Applied Sciences, Cross-Border Faculty of Humanities, Economics and Engineering, Dunarea de Jos University of Galati, 111 Domneasca Street, 800201 Galati, Romania
6
Department of Mathematics and Applied Mathematics, Sefako Makgatho Health Sciences University, Medunsa 0204, South Africa
7
Department of Computer Engineering, Biruni University, 34010 Istanbul, 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
9
Department of Computer Science and Information Technology, Faculty of Automation, Computers, Electrical Engineering and Electronics, Dunarea de Jos University of Galati, 47 Domneasca Street, 800008 Galati, Romania
10
Department of Mathematics and Statistics, Tshwane University of Technology, Pretoria 0008, South Africa
11
Department of Business Administration, Faculty of Economics and Business Administration, Dunarea de Jos University of Galati, 59-61 Nicolae Balcescu Street, 800001 Galati, Romania
12
Department of Industrial Engineering, University of Tabuk, Tabuk 45512, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(5), 963; https://doi.org/10.3390/sym15050963
Submission received: 7 February 2023 / Revised: 4 April 2023 / Accepted: 21 April 2023 / Published: 23 April 2023
(This article belongs to the Special Issue Nonlinear Optical Research)

Abstract

:
The current paper recovers cubic–quartic optical solitons in fiber Bragg gratings having polynomial law of nonlinear refractive index structures. Lie symmetry analysis is carried out, starting with the basic analysis. Then, it is followed through with improved Kudryashov and generalized Arnous schemes. The parameter constraints are also identified for the existence of such solitons. Numerical surface plots support the adopted applied analysis.

1. Introduction

The concept of cubic–quartic (CQ) solitons emerged about half a decade ago when it was realized that soliton sustainability becomes questionable with catastrophic consequences when chromatic dispersion (CD) tends to become depleted [1,2,3,4,5]. Therefore, third-order dispersion (3OD) and fourth-order dispersion (4OD) effects are included to replenish the low count of CD. Thus, CQ solitons were conceived [6,7,8,9,10]. Aside from this measure to counter low counts of CD, another approach is from an engineering marvel that introduces a grating structure in the internal walls of the fiber core, known as Bragg gratings after the engineer who first proposed this structure [11,12,13,14,15]. This grating structure also produces dispersive reflectivity, compensating for low dispersive effects. Thus, the two approaches would collectively ensure the safe transmission of solitons across intercontinental distances without any adverse effect.
The current paper models such a physical situation in fiber Bragg gratings with CQ dispersive effects. The source of self–phase modulation (SPM) comes from the polynomial law of nonlinear refractive index [3,4,11,14]. The model is therefore the coupled nonlinear Schrodinger’s equation (NLSE) that is addressed in the paper. A Lie symmetry analysis first reduces the corresponding partial differential equations (PDEs) into ordinary differential equations (ODEs) that are subsequently addressed using the improved Kudryashov scheme and the generalized Arnous method. This gave way to retrieving singular and bright optoelectronic wave fields with the model presented in the paper. The supporting surface plots for the bright solitons are presented. The details are addressed after a succinct picture of the model.

Governing Equations

Periodic structures called Bragg gratings can be created by etching into an optical fiber [16,17,18,19,20]. Bragg gratings have a unique property: they reflect light of a specific wavelength while allowing light of other wavelengths to pass through [21,22,23,24,25]. This is possible because the grating’s period and the fiber core’s refractive index determine the wavelength of the reflected light. To address the problem of limited bandwidth in fiber optic communication systems due to increasing data demands, Bragg gratings can serve as a solution by increasing the capacity of the optical fiber [26,27,28,29,30]. This capacity issue occurs when the available bandwidth of a system becomes depleted, and is often referred to as CD running low. By incorporating Bragg gratings, the system can be optimized to accommodate higher data rates while maintaining the integrity of the signal [31,32,33,34,35]. Wavelength division multiplexing (WDM) is a technique that can be used to increase the capacity of an optical fiber by transmitting multiple signals simultaneously over different wavelengths in the same fiber. To implement WDM, Bragg gratings can be used to create filters that are specific to certain wavelengths [36,37,38,39,40]. These filters can separate the different signals, allowing multiple channels to operate at different wavelengths within the same fiber. By using Bragg gratings in this way, the fiber’s capacity can be significantly increased, enabling it to accommodate a greater amount of data traffic [41,42,43,44,45,46]. Bragg gratings have another important application in the field of optical amplifiers, which are used to increase the strength of the signal in the fiber. By using Bragg gratings in conjunction with erbium-doped fiber amplifiers (EDFAs), it is possible to achieve wavelength-specific amplification. This feature enables multiple signals at different wavelengths to be amplified simultaneously, resulting in a significant increase in the overall capacity of systems. By implementing Bragg gratings in optical amplifiers, it is possible to achieve more efficient and effective signal amplification, making it possible to transmit data over longer distances with fewer errors. In the field of fiber optic communication, Bragg gratings have emerged as a versatile and powerful tool for increasing the capacity of the system. By enabling wavelength division multiplexing and providing wavelength-specific amplification in optical amplifiers, Bragg gratings offer a reliable and effective means of boosting the performance of fiber optic communication systems. With their ability to accommodate multiple signals at different wavelengths and amplify these signals with precision, Bragg gratings are an essential technology for meeting the growing demands of modern data transmission.
Bragg gratings have a multitude of applications in the context of chromatic dispersion (CD) running low in optical communication systems. Chromatic dispersion is a phenomenon that can cause signal distortion as different wavelengths of light travel at varying speeds through the fiber. Bragg gratings can be utilized to mitigate the effects of CD by creating wavelength-specific filters that compensate for the differences in propagation velocities. This approach can help to improve the quality and reliability of optical communication systems, particularly over long distances, by minimizing signal distortion and maintaining the integrity of the transmitted data. In addition to the applications mentioned earlier, Bragg gratings have proven to be a powerful tool for addressing CD-related issues in fiber optic communication. Furthermore, Bragg gratings have several other important applications and physical significance related to chromatic dispersion (CD) running low. These applications include the following: one key application of Bragg gratings in the context of chromatic dispersion (CD) running low is the compensation of CD-induced signal distortion. Bragg gratings can be designed to induce a precise phase shift in reflected light, which can offset the phase shift induced by CD. By doing so, the distortion of the optical signal can be reduced or eliminated, thereby improving the quality and reliability of data transmission. This capability is particularly important in high-capacity optical communication systems, where the effects of CD can become a significant limiting factor in the performance of the system. Bragg gratings are a versatile and essential tool in optical fiber communication engineering. They can be used to address a range of issues related to chromatic dispersion, including compensating for dispersion slope and mitigating the effects of CD-induced signal distortion. Additionally, Bragg gratings can provide sensing capabilities for monitoring a variety of environmental parameters, making them an important technology for a range of applications. With their ability to enhance the performance and reliability of fiber optic communication systems, Bragg gratings are a key component in the design and implementation of cutting-edge optical communication networks.
The mathematical model that governs CQ solitons in fiber Bragg gratings with the polynomial form of SPM [47] is structured as:
i q t + i a 1 r x x x + b 1 r x x x x + c 1 | q | 2 + d 1 | r | 2 q + ξ 1 | q | 4 + η 1 | q | 2 | r | 2 + ζ 1 | r | 4 q + l 1 | q | 6 + m 1 | q | 4 | r | 2 + n 1 | q | 2 | r | 4 + p 1 | r | 6 q + i α 1 q x + β 1 r = 0 ,
and
i r t + i a 2 q x x x + b 2 q x x x x + c 2 | r | 2 + d 2 | q | 2 r + ξ 2 | r | 4 + η 2 | q | 2 | r | 2 + ζ 2 | r | 4 r + l 2 | r | 6 + m 2 | r | 4 | q | 2 + n 2 | r | 2 | q | 4 + p 2 | q | 6 r + i α 2 r x + β 2 q = 0 .
For j = 1 , 2 , β j gives the detuning parameters, and α j yields the intermodal dispersion. d j comes from the cross–phase modulation (XPM) for cubic nonlinearity, while q ( x , t ) and r ( x , t ) are the wave profiles propagating forwards and backwards in sequence. Next, a j and b j indicate the coefficients of 3OD and 4OD reflectivity in sequence, while c j stems from the SPM. The XPM coefficients are given by m j , n j and p j , while the coefficients ξ j are for the SPM. Furthermore, from septic nonlinearity, the SPM coefficients are l j , while η j and ζ j arise from the XPM.

2. Lie Symmetry Analysis

The Norwegian mathematician Marius Sophus Lie (1842–1899) discovered the Lie group theory. According to Lie, the various approaches to solving differential equations precisely were simply various applications of a general approach to integration, known as the theory of transformation groups. Lie theory is the mathematical foundation for understanding the symmetry of differential equations. Lie’s introduction of the theory was so insightful that it allowed mathematicians to come up with many different methods and techniques for solving differential equations. Differential equations using symmetry analysis are becoming increasingly common in mathematics and the physical sciences, due to their effectiveness in solving complex problems. Tanwar et al. studied the propagation of some nonlinear dispersive water waves, such as waves in an infinitely narrow channel [48], waves with nonuniform velocities [49], waves in two horizontal directions in shallow water [50], the long wave equation with a large wavelength [51], and undulating bores in shallow water [52], by using Lie symmetry analysis. Some recent and very useful books devoted to Lie symmetry analysis are [53,54,55,56,57].
Here, Lie group analysis is used to derive the symmetry reduction of the system (1) and (2). Let us start by defining r ( x , t ) and q ( x , t ) as follows:
r ( x , t ) = u 2 ( x , t ) + ι v 2 ( x , t ) ,
and
q ( x , t ) = u 1 ( x , t ) + ι v 1 ( x , t ) ,
which divide the system (1) and (2) into the following systems
0 = v i t a i v j x x x + b i u j x x x x + c i u i 2 + v i 2 + d i u j 2 + v j 2 u i + ξ i u i 2 + v i 2 2 + η i u i 2 + v i 2 u j 2 + v j 2 + ζ i u j 2 + v j 2 2 u i + l i u i 2 + v i 2 3 + m i u i 2 + v i 2 2 u j 2 + v j 2 + n i u i 2 + v i 2 u j 2 + v j 2 2 + p i u j 2 + v j 2 3 u i α i v i x + β i u j ,
and
0 = u i t + a i u j x x x + b i v j x x x x + c i u i 2 + v i 2 + d i u j 2 + v j 2 v i + ξ i u i 2 + v i 2 2 + η i u i 2 + v i 2 u j 2 + v j 2 + ζ i u j 2 + v j 2 2 v i + l i u i 2 + v i 2 3 + m i u i 2 + v i 2 2 u j 2 + v j 2 + n i u i 2 + v i 2 u j 2 + v j 2 2 + p i u j 2 + v j 2 3 v i + α i u i x + β i v j ,
where j = 3 i and i = 1 , 2 . To derive symmetries for the system (5) and (6), the Lie group of point transformations is presented as below:
v 2 * = v 2 + ϵ ψ 2 ( x , t , u 1 , u 2 , v 1 , v 2 ) + O ( ϵ 2 ) , v 1 * = v 1 + ϵ ψ 1 ( x , t , u 1 , u 2 , v 1 , v 2 ) + O ( ϵ 2 ) , u 2 * = u 2 + ϵ ϕ 2 ( x , t , u 1 , u 2 , v 1 , v 2 ) + O ( ϵ 2 ) , u 1 * = u 1 + ϵ ϕ 1 ( x , t , u 1 , u 2 , v 1 , v 2 ) + O ( ϵ 2 ) , t * = t + ϵ τ ( x , t , u 1 , u 2 , v 1 , v 2 ) + O ( ϵ 2 ) , x * = x + ϵ ξ ( x , t , u 1 , u 2 , v 1 , v 2 ) + O ( ϵ 2 ) ,
where ϕ 2 , τ , ϕ 1 , ξ , ψ 2 and ψ 1 are infinitesimals that can be achieved by linearizing the invariance condition about ϵ = 0 . The vector field V connected to the Lie group of transformations (7) mentioned above is
V = ξ x + τ t + ϕ 1 u 1 + ϕ 2 u 2 + ψ 1 v 1 + ψ 2 v 2 .
The system of invariance conditions is generated by applying the transformations (7) to the system (5) and (6). The values of the appropriate extended infinitesimals should then be substituted, and the coefficient of the various powers of the differentials of u 1 , u 2 , v 1 , v 2 should be equalized to provide a system of determining equations. The following are the resultant infinitesimals that are acquired by resolving the system of determining equations:
ξ = C 2 , τ = C 1 , ϕ 1 = C 3 v 1 , ϕ 2 = C 3 v 2 , ψ 1 = C 3 u 1 , ψ 2 = C 3 u 2 ,
where C 1 , C 2 and C 3 are arbitrary constants. As a result, the infinitesimal generators derived from (9) are extracted as
V 1 = t , V 2 = x , V 3 = u 2 v 2 + u 1 v 1 v 2 u 2 v 1 u 1 .
The Lie algebra (10) comprises some additional Lie symmetry for the system (1) and (2). For example, if we consider system (1) and (2) with the non-zero b i , c i , d i , and all other parameters vanish, then the Lie algebra (10) admits the additional Lie symmetry
V 4 = 4 t t + x x 2 u 1 u 1 2 u 2 u 2 2 v 1 v 1 2 v 2 v 2 .
Consider the system (1) and (2) with the non-zero b i , ξ i , η i , ζ i ; all other parameters vanish. In this case, the system (1) and (2) admits the additional Lie symmetry as
V 4 = 4 t t + x x u 1 u 1 u 2 u 2 v 1 v 1 v 2 v 2 .
Consider the system (1) and (2) with the non-zero a i , c i , d i ; all other parameters vanish. In this case, the system (1) and (2) admits the additional Lie symmetry as
V 4 = 3 t t + x x 3 u 1 2 u 1 3 u 2 2 u 2 3 v 1 2 v 1 3 v 2 2 v 2 .
Consider the system (1) and (2) with the non-zero a i , ξ i , η i , ζ i ; all other parameters vanish. In this case, the system (1) and (2) admits the additional Lie symmetry as
V 4 = 3 t t + x x 3 u 1 4 u 1 3 u 2 4 u 2 3 v 1 4 v 1 3 v 2 4 v 2 .

Symmetry Reduction

In this subsection, we reduce the system (1) and (2) into coupled ordinary differential equations by considering the Lie algebra (10)–(14).
Case 1: For the vector field
V 3 + λ V 2 + μ V 1 = u 2 v 2 + u 1 v 1 v 2 u 2 v 1 u 1 + μ t + λ x ,
the corresponding similarity variables are recovered as
v 2 ( x , t ) = P 2 ( σ ) sin Q 2 ( σ ) , u 2 ( x , t ) = P 2 ( σ ) cos Q 2 ( σ ) , v 1 ( x , t ) = P 1 ( σ ) sin Q 1 ( σ ) , u 1 ( x , t ) = P 1 ( σ ) cos Q 1 ( σ ) , σ = μ x λ t ,
where λ and μ are non-zero constants, while P 1 and P 2 are new dependent variables. By utilizing (16), along with the relation
Q 1 ( σ ) = Q 2 ( σ ) = Q ( σ ) ,
the system (5) and (6) reduces to
0 = λ Q P i a i μ 3 Q 3 P j + 3 Q P j + 3 Q P j + Q P j + b i μ 4 Q 4 P j 6 Q 2 P j 12 Q Q P j 4 Q Q P j 3 Q 2 P j + P j + c i P i 2 + d i P j 2 P i + ξ i P i 4 + η i P i 2 P j 2 + ζ i P j 2 P i + l i P i 6 + m i P i 4 P j 2 + n i P i 2 P j 4 + p i P j 6 P i α i μ Q P i + β i P j ,
and
0 = λ P i + a i μ 3 3 Q 2 P j 3 Q Q P j + P j + b i μ 4 4 Q 3 P j 6 Q 2 Q P j + 4 Q P j + 6 Q P j + 4 Q P j + Q P j + α i μ P i ,
where j = 3 i and i = 1 , 2 . Now, from (19), comparing the coefficient of P j equal to zero, gives
Q ( σ ) = a i 4 μ b i σ + k 1 ,
where k 1 is the constant of integration. Next, let us consider
P 2 = ρ P 1 , λ = μ α 1 μ ρ a 1 3 8 b 1 2 = μ α 2 μ ρ a 2 3 8 b 2 2 , ( ρ 6 p 1 + ρ 4 n 1 + ρ 2 m 1 + l 1 ) = ρ ( ρ 6 l 2 + ρ 4 m 2 + ρ 2 n 2 + p 2 ) , ( ρ 4 ζ 1 + ρ 2 η 1 + ξ 1 ) = ρ ( ρ 4 ξ 2 + ρ 2 η 2 + ζ 2 ) , ( ρ 2 d 1 + c 1 ) = ( ρ 2 c 2 + d 2 ) , β 1 3 a 1 4 256 b 1 3 ρ a 1 4 μ b 1 ( λ μ α 1 ) = a 2 ρ 4 μ b 2 ( λ μ α 2 ) + β 2 3 a 2 4 256 b 2 3 , a 1 2 ρ b 1 = a 2 2 b 2 , b 1 ρ = b 2 .
Utilizing (20), along with the restrictions of (21), the system of Equations (18) and (19) leads to
S 1 P 1 + S 2 P 1 + S 3 P 1 7 + S 4 P 1 5 + S 5 P 1 3 + S 6 P 1 = 0 ,
where
S 1 = ρ b 1 μ 4 , S 2 = 3 ρ μ 2 a 1 2 8 b 1 , S 3 = ( ρ 6 p 1 + ρ 4 n 1 + ρ 2 m 1 + l 1 ) , S 4 = ( ρ 4 ζ 1 + ρ 2 η 1 + ξ 1 ) , S 5 = ( ρ 2 d 1 + c 1 ) , S 6 = β 1 3 a 1 4 256 b 1 3 ρ a 1 4 μ b 1 ( λ μ α 1 ) .
Now, by setting
P 1 ( σ ) = U ( σ ) 2 3 ,
Equation (22) changes to
81 S 3 U 8 + 81 S 4 U 20 3 + 81 S 5 U 16 3 + 81 S 6 U 4 + 54 S 1 U + S 2 U U 3 18 S 2 U 2 + 4 S 1 U U + 3 S 1 U 2 U 2 + 144 S 1 U U 2 U 56 S 1 U 4 = 0 .
For integrability aspects, one can set
S 4 = S 5 = 0 .
So, Equation (25) becomes
81 S 3 U 8 + 81 S 6 U 4 + 54 S 1 U + S 2 U U 3 18 S 2 U 2 + 4 S 1 U U + 3 S 1 U 2 U 2 + 144 S 1 U U 2 U 56 S 1 U 4 = 0 .
Case 2: For the vector field given by (11)
V 4 = 4 t t + x x 2 u 1 u 1 2 u 2 u 2 2 v 1 v 1 2 v 2 v 2 ,
the corresponding similarity variables are recovered as
v 2 ( x , t ) = Q 2 ( σ ) t , u 2 ( x , t ) = P 2 ( σ ) t , v 1 ( x , t ) = Q 1 ( σ ) t , u 1 ( x , t ) = P 1 ( σ ) t , σ = x 4 t ,
where P i and Q i are new dependent variables. By utilizing (29), the system (5) and (6) reduces to
0 = c i P i 3 + d i Q j 2 + c i Q i 2 + d i P j 2 P i + 1 2 Q i + 256 σ 3 b i P j + 1152 σ 2 b i P j + 816 σ b i P j + σ Q i + 24 b i P j ,
and
0 = c i P i 2 Q i 1 2 P i + c i Q i 3 + d i P j 2 + Q j 2 Q i + 256 σ 3 b i Q j + 1152 σ 2 b i Q j + 816 σ b i Q j σ P i + 24 b i Q j ,
where j = 3 i and i = 1 , 2 .
Case 3: For the vector field given by (12)
V 4 = 4 t t + x x u 1 u 1 u 2 u 2 v 1 v 1 v 2 v 2 ,
the corresponding similarity variables are recovered as
v 2 ( x , t ) = Q 2 ( σ ) x , u 2 ( x , t ) = P 2 ( σ ) x , v 1 ( x , t ) = Q 1 ( σ ) x , u 1 ( x , t ) = P 1 ( σ ) x , σ = x 4 t ,
where P i and Q i are new dependent variables. By utilizing (33), the system (5) and (6) reduces to
0 = ξ i P i 5 + 2 ξ i Q i 2 + η i Q j 2 + P j 2 P i 3 + ξ i Q i 4 + η i Q j 2 + P j 2 Q i + ζ i Q j 2 + P j 2 2 P i + 24 b i P j + 432 b i σ 2 P j + 896 b i σ 3 P j + 256 b i σ 4 P j + σ 2 Q i 24 b i σ P j ,
and
0 = η i Q j 2 + P j 2 Q i 3 + ξ i Q i 5 + 2 ξ i Q i 2 + η i Q j 2 + P j 2 Q i P i 2 + ζ i Q j 2 + P j 2 2 Q i + 24 b i Q j + ξ i P i 4 Q i + 256 b i σ 4 Q j 24 b i σ Q j σ 2 P i + 896 b i σ 3 Q j + 432 b i σ 2 Q j ,
where j = 3 i and i = 1 , 2 .
Case 4: For the vector field given by (13)
V 4 = 3 t t + x x 3 u 1 2 u 1 3 u 2 2 u 2 3 v 1 2 v 1 3 v 2 2 v 2 ,
the corresponding similarity variables are recovered as
v 2 ( x , t ) = Q 2 ( σ ) t , u 2 ( x , t ) = P 2 ( σ ) t , v 1 ( x , t ) = Q 1 ( σ ) t , u 1 ( x , t ) = P 1 ( σ ) t , σ = x 3 t ,
where P i and Q i are new dependent variables. By utilizing (37), the system (5) and (6) reduces to
0 = c i P i 3 + d i Q j 2 + c i Q i 2 + d i P j 2 P i + 1 2 Q i 54 σ a i Q j 27 σ 2 a i Q j + σ Q i 6 a i Q j ,
and
0 = P i 2 Q i c i P i 2 + Q i 3 c i + Q j 2 d i + P j 2 d i Q i + 54 σ a i P j + 27 σ 2 a i P j σ P i + 6 a i P j ,
where j = 3 i and i = 1 , 2 .
Case 5: For the vector field given by (14)
V 4 = 3 t t + x x 3 u 1 4 u 1 3 u 2 4 u 2 3 v 1 4 v 1 3 v 2 4 v 2 ,
the corresponding similarity variables are recovered as
v 2 ( x , t ) = Q 2 ( σ ) t 4 , u 2 ( x , t ) = P 2 ( σ ) t 4 , v 1 ( x , t ) = Q 1 ( σ ) t 4 , u 1 ( x , t ) = P 1 ( σ ) t 4 , σ = x 3 t ,
where P i and Q i are new dependent variables. By utilizing (41), the system (5) and (6) reduces to
0 = ξ i P i 5 + 2 ξ i Q i 2 + η i Q j 2 + P j 2 P i 3 + ξ i Q i 4 + η i Q j 2 + P j 2 Q i 2 + ζ i Q j 4 + 2 ζ i P j 2 Q j 2 + ζ i P j 4 P i + 1 4 Q i 54 σ a i Q j 27 σ 2 a i Q j + σ Q i 6 a i Q j ,
and
0 = ξ i Q i 5 1 4 P i + 2 ξ i Q i + η i Q j 2 + P j 2 Q i P i 2 + 2 ζ i P j 2 Q j 2 + ζ i P j 4 + ζ i Q j 4 Q i + ξ i P i 4 Q i + η i Q j 2 + P j 2 Q i 3 + 27 σ 2 a i P j + 54 σ a i P j + 6 a i P j σ P i ,
where j = 3 i and i = 1 , 2 .
It is difficult to express the exact solutions of coupled ordinary differential equations derived in the above cases (2)–(5), but one can express their solutions in the numerical forms.

3. Improved Kudryashov Method

Under the circumstances of the improved Kudryashov method [58], Equation (27) satisfies the explicit solution of the form
U ( σ ) = B 0 + B 1 R ( σ ) ,
where B 0 and B 1 are arbitrary constants, with B 1 0 . In (44), the function R ( σ ) fulfills the relation
R ( σ ) = δ ϱ R ( σ ) 3 R ( σ ) ,
and possesses the solution
R ( σ ) = 1 ϱ + k e 2 δ σ ,
where k is a real constant. By inserting (44) into (27) and then utilizing (45), we receive a polynomial of R ( σ ) . Then, by gathering the coefficients of different powers of derived polynomials equal to zero, one secures a system of equations that permits the following results:
B 0 = 0 , B 1 = ± 4 δ 3 22 ϱ 2 S 1 S 3 4 , S 2 = 68 δ 2 S 1 9 , S 6 = 256 δ 4 S 1 81 .
Therefore, the straddled soliton profiles are presented as below:
q ( x , t ) = ± 4 δ 3 22 ϱ 2 S 1 S 3 4 1 ϱ + k e 2 δ μ x μ α 1 μ ρ a 1 3 8 b 1 2 t 2 3 e i a 1 4 μ b 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t + k 1 ,
and
r ( x , t ) = ± ρ 4 δ 3 22 ϱ 2 S 1 S 3 4 1 ϱ + k e 2 δ μ x μ α 1 μ ρ a 1 3 8 b 1 2 t 2 3 e i a 1 4 μ b 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t + k 1 ,
with constraint S 1 S 3 < 0 .

4. Generalized Arnous Method

In this section, we extract the optical soliton solutions of the system (1) and (2) by using the generalized Arnous scheme [59] on Equation (27). According to this method, Equation (27) holds the solution
P ( σ ) = ϱ 0 + ϱ 1 + ν 1 Φ ( σ ) Φ ( σ ) ,
where ϱ 0 , ϱ 1 and ν 1 are arbitrary constants. Here Φ ( σ ) satisfies the following relation
Φ ( σ ) 2 = Φ ( σ ) 2 χ ln ( A ) 2 ,
with
Φ ( n ) ( σ ) = Φ ( σ ) ln ( A ) n 1 , n is odd , Φ ( σ ) ln ( A ) n , n is even , n 2
where A > 0 and A 1 . Equation (51) admits a solution as shown below:
Φ ( σ ) = χ 4 k ln ( A ) A σ + k ln ( A ) A σ ,
where χ and k are arbitrary parameters.
By inserting (50) along with (51) into Equation (22), we receive a polynomial of 1 Φ ( σ ) Φ ( σ ) Φ ( σ ) . Then, by gathering the coefficients of different powers of derived polynomial equal to zero, one obtains a system of equations that admits the following results:
ϱ 0 = 0 , ν 1 = 0 , ϱ 1 = ± 2 ln ( A ) 3 55 χ 2 S 1 S 3 4 , S 2 = 68 ln ( A ) 2 S 1 9 , S 6 = 256 ln ( A ) 4 S 1 81 .
Therefore, we reveal the straddled soliton profiles
q ( x , t ) = ± 2 ln ( A ) 3 55 χ 2 S 1 S 3 4 4 k ln ( A ) A μ x μ α 1 μ ρ a 1 3 8 b 1 2 t 4 k 2 ln ( A ) 2 A 2 μ x 2 μ α 1 μ ρ a 1 3 8 b 1 2 t + χ 2 3 e i a 1 4 μ b 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t + k 1 ,
and
r ( x , t ) = ± ρ 2 ln ( A ) 3 55 χ 2 S 1 S 3 4 4 k ln ( A ) A μ x μ α 1 μ ρ a 1 3 8 b 1 2 t 4 k 2 ln ( A ) 2 A 2 μ x 2 μ α 1 μ ρ a 1 3 8 b 1 2 t + χ 2 3 e i a 1 4 μ b 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t + k 1 ,
with constraint S 1 S 3 < 0 .
If we consider χ = 4 k 2 log ( A ) 2 in solutions (55) and (56), then we recover the bright optoelectronic wave fields
q ( x , t ) = ± 1 3 k 55 χ 2 S 1 S 3 4 sec h μ x μ α 1 μ ρ a 1 3 8 b 1 2 t ln A 2 3 e i a 1 4 μ b 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t + k 1 ,
and
r ( x , t ) = ± ρ 1 3 k 55 χ 2 S 1 S 3 4 sec h μ x μ α 1 μ ρ a 1 3 8 b 1 2 t ln A 2 3 e i a 1 4 μ b 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t + k 1 ,
while, if we consider χ = 4 k 2 log ( A ) 2 in solutions (55) and (56), then we recover the singular soliton profiles
q ( x , t ) = ± 1 3 k 55 χ 2 S 1 S 3 4 csc h μ x μ α 1 μ ρ a 1 3 8 b 1 2 t ln A 2 3 e i a 1 4 μ b 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t + k 1 ,
and
r ( x , t ) = ± ρ 1 3 k 55 χ 2 S 1 S 3 4 csc h μ x μ α 1 μ ρ a 1 3 8 b 1 2 t ln A 2 3 e i a 1 4 μ b 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t + k 1 .
The surface plots of solitons (57) and (58) are depicted in Figure 1. The parameter values chosen are: k = 1 , ρ = 2 , μ = 1 , p 1 = 1 , n 1 = 1 , m 1 = 1 , l 1 = 1 , α 1 = 1 , a 1 = 1 , b 1 = −1 and A = e . The surface plot of a bright soliton is a smooth, bell-shaped curve that remains relatively constant over time and space. The curve is localized in space and time, and its height corresponds to the amplitude of the soliton, while the x-coordinate represents the spatial dimension, and the t-coordinate represents time. Bright solitons are known for their ability to maintain their shape and amplitude over long distances, due to a balance between dispersion and nonlinearity. This stability means that the surface plot of a bright soliton shows only minor fluctuations or deformations, which reflects the regular and predictable behavior of the soliton. The symmetrical nature of the surface plot indicates the stability and regularity of the soliton. These properties make bright solitons ideal for use in applications such as fiber optic communications or Bose-Einstein condensates in ultracold atomic gases, where stability and control are crucial. The surface plot of a bright soliton is a visual representation of a stable and localized wave phenomenon that can propagate over long distances without losing its shape or amplitude.

5. Extended Tanh Technique

In this section, we derive optical soliton solutions of the system (1) and (2) by performing the extended tanh method [60] on Equation (27). The extended tanh scheme suggests the solution of Equation (27) as
U ( σ ) = A 0 + A 1 tan h ( m ξ ) + A 1 tan h 1 ( m ξ )
where A 0 , A 1 and A 1 are arbitrary constants. Plugging (61) into (27) yields a system of equations that satisfies the results:
A 0 = 0 , A 1 = A 1 , m = 3 S 2 272 S 1 , S 3 = ± 55 S 2 2 4624 S 1 A 1 4 , S 6 = 16 S 2 2 289 S 1 .
Hence, we achieve the singular optoelectronic wave fields
q ( x , t ) = ± A 1 4 1 / 3 csc h 2 m 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t 2 3 e i a 1 4 μ b 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t + k 1 ,
and
r ( x , t ) = ± ρ A 1 4 1 / 3 csc h 2 m 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t 2 3 e i a 1 4 μ b 1 μ x μ α 1 μ ρ a 1 3 8 b 1 2 t + k 1 ,
with the constraint S 1 S 2 < 0 , and (62). So, in this case, we have singular solitons.

6. Conclusions

The current paper is about retrieving bright and singular CQ solitons in fiber Bragg gratings with the polynomial form of SPM. These soliton solutions are of great importance in the study of fiber Bragg gratings. The main approach is Lie symmetry analysis, which reduced the PDEs to the corresponding ODEs that were addressed using a couple of newly proposed integration schemes, namely, the improved Kudryashov and generalized Arnous schemes. These new findings in the literature of solitons will lead to a number of new avenues to walk on. The conservation laws for the full-blown mathematical model are to be recovered, and this would naturally lead to quasistationary solitons and perturbation theory. In our analysis, it is not possible to come up with the conservation laws for Bragg gratings. This is due to the presence of dispersive reflectivity, as opposed to chromatic dispersion. Therefore, the multiplier approach fails to retrieve the conservation laws and fluxes. The failure attributed to the usage of the Lagrangian approach is similar.
The improved Kudryashov method and generalized Arnous method are both analytical techniques used to solve nonlinear partial differential equations (PDEs) that arise in various physical systems, including fiber Bragg gratings with polynomial law of nonlinear refractive index structures. The improved Kudryashov method is a powerful and efficient method for finding analytical solutions for nonlinear PDEs, particularly those with singularities. It has been used to find bright soliton solutions for various nonlinear systems, including fiber Bragg gratings. Bright solitons are characterized by a localized, finite amplitude that propagates without changing its shape. They are typically stable and can be used for applications such as signal processing and data transmission. On the other hand, the generalized Arnous method is an analytical method that is used to find singular soliton solutions for nonlinear PDEs. Singular solitons are characterized by a divergent amplitude at a particular point or points. They are less stable than bright solitons, but they can still be useful for certain applications in fields such as nonlinear optics and plasma physics. In our analysis, it seems that the improved Kudryashov method was able to find both bright and singular soliton solutions, while the generalized Arnous method was only able to find singular solitons. This suggests that the improved Kudryashov method may be a more powerful and versatile method for finding solutions to nonlinear PDEs in fiber Bragg gratings with polynomial law of nonlinear refractive index structures. However, the choice of method should depend on the specific requirements of the analysis and the properties of the solutions.

Author Contributions

Conceptualization, S.M. (Sandeep Malik) and S.K.; methodology, A.B.; software, Y.Y. and L.M.; writing—original draft preparation, S.M. (Simona Moldovanu); writing—review and editing, C.I. and S.P.M.; project administration, D.B. and 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. Bayram, M. Optical bullets with Biswas–Milovic equation having Kerr and parabolic laws of nonlinearity. Optik 2022, 270, 170046. [Google Scholar] [CrossRef]
  2. Belyaeva, T.; Serkin, V. Wave-particle duality of solitons and solitonic analog of the Ramsauer-Townsend effect. Eur. Phys. J. D 2012, 66, 1–9. [Google Scholar] [CrossRef]
  3. Biswas, A. Optical soliton cooling with polynomial law of nonlinear refractive index. J. Opt. 2020, 49, 580–583. [Google Scholar] [CrossRef]
  4. Gonzalez–Gaxiola, O.; Biswas, A.; Yildirim, Y.; Alshehri, H.M. Bright optical solitons with polynomial law of nonlinear refractive index by Adomian decomposition scheme. Proc. Est. Acad. Sci. 2022, 71, 213–220. [Google Scholar] [CrossRef]
  5. 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]
  6. Kudryashov, N.A. Rational solutions of equations associated with the second Painlevé equation. Regul. Chaotic Dyn. 2020, 25, 273–280. [Google Scholar] [CrossRef]
  7. Kudryashov, N.A.; Safonova, D.V. Painlevé analysis and traveling wave solutions of the sixth order differential equation with non-local nonlinearity. Optik 2021, 244, 167586. [Google Scholar] [CrossRef]
  8. Kudryashov, N.A.; Safonova, D.V.; Biswas, A. Painlevé analysis and a solution to the traveling wave reduction of the Radhakrishnan-Kundu-Lakshmanan equation. Regul. Chaotic Dyn. 2019, 24, 607–614. [Google Scholar] [CrossRef]
  9. Tang, L. Bifurcation analysis and multiple solitons in birefringent fibers with coupled Schrödinger-Hirota equation. Chaos Solitons Fractals 2022, 161, 112383. [Google Scholar] [CrossRef]
  10. 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]
  11. Triki, H.; Zhou, Q.; Biswas, A.; Liu, W.; Yıldırım, Y.; Alshehri, H.M.; Belic, M.R. Chirped optical solitons having polynomial law of nonlinear refractive index with self-steepening and nonlinear dispersion. Phys. Lett. A 2021, 417, 127698. [Google Scholar] [CrossRef]
  12. Wang, M.Y. Optical solitons with perturbed complex Ginzburg–Landau equation in Kerr and cubic–quintic–septic nonlinearity. Results Phys. 2022, 33, 105077. [Google Scholar] [CrossRef]
  13. Wazwaz, A.M.; Mehanna, M. Higher-order Sasa–Satsuma equation: Bright and dark optical solitons. Optik 2021, 243, 167421. [Google Scholar] [CrossRef]
  14. Zayed, E.M.; Alngar, M.E.; Biswas, A.; Ekici, M.; Guggilla, P.; Khan, S.; Alshehri, H.M.; Belic, M.R. Optical solitons in fiber Bragg gratings with polynomial law nonlinearity and cubic–quartic dispersive reflectivity. Opt. Spectrosc. 2022, 130, 28–34. [Google Scholar] [CrossRef]
  15. Zhou, Q. Influence of parameters of optical fibers on optical soliton interactions. Chin. Phys. Lett. 2022, 39, 010501. [Google Scholar] [CrossRef]
  16. Akter, A.; Islam, M.J.; Atai, J. Quiescent Gap Solitons in Coupled Nonuniform Bragg Gratings with Cubic-Quintic Nonlinearity. Appl. Sci. 2021, 11, 4833. [Google Scholar] [CrossRef]
  17. Islam, M.J.; Atai, J. Stability of moving Bragg solitons in a semilinear coupled system with cubic–quintic nonlinearity. J. Mod. Opt. 2021, 68, 365–373. [Google Scholar] [CrossRef]
  18. Islam, M.J.; Atai, J. Dynamics of colliding counterpropagating solitons in coupled Bragg gratings with cubic-quintic nonlinearity. J. Mod. Opt. 2019, 66, 1498–1505. [Google Scholar] [CrossRef]
  19. Ahmed, T.; Atai, J. Soliton-soliton dynamics in a dual-core system with separated nonlinearity and nonuniform Bragg grating. Nonlinear Dyn. 2019, 97, 1515–1523. [Google Scholar] [CrossRef]
  20. Islam, M.J.; Atai, J. Soliton–soliton interactions in a grating-assisted coupler with cubic-quintic nonlinearity. J. Mod. Opt. 2018, 65, 2153–2159. [Google Scholar] [CrossRef]
  21. Islam, M.J.; Atai, J. Stability of moving gap solitons in linearly coupled Bragg gratings with cubic–quintic nonlinearity. Nonlinear Dyn. 2018, 91, 2725–2733. [Google Scholar] [CrossRef]
  22. Ahmed, T.; Atai, J. Bragg solitons in systems with separated nonuniform Bragg grating and nonlinearity. Phys. Rev. E 2017, 96, 032222. [Google Scholar] [CrossRef] [PubMed]
  23. Sultana, J.; Islam, M.S.; Atai, J.; Islam, M.R.; Abbott, D. Near-zero dispersion flattened, low-loss porous-core waveguide design for terahertz signal transmission. Opt. Eng. 2017, 56, 076114. [Google Scholar] [CrossRef]
  24. Islam, M.S.; Sultana, J.; Atai, J.; Abbott, D.; Rana, S.; Islam, M.R. Ultra low-loss hybrid core porous fiber for broadband applications. Appl. Opt. 2017, 56, 1232–1237. [Google Scholar] [CrossRef] [PubMed]
  25. Jahirul Islam, M.; Atai, J. Stability of Bragg grating solitons in a semilinear dual-core system with cubic–quintic nonlinearity. Nonlinear Dyn. 2017, 87, 1693–1701. [Google Scholar] [CrossRef]
  26. Chowdhury, S.S.; Atai, J. Interaction dynamics of Bragg grating solitons in a semilinear dual-core system with dispersive reflectivity. J. Mod. Opt. 2016, 63, 2238–2245. [Google Scholar] [CrossRef]
  27. Cao, H.; Atai, J.; Zuo, J.; Yu, Y.; Gbadebo, A.; Xiong, B.; Hou, J.; Liang, P.; Gao, Y.; Shu, X. Simultaneous multichannel carrier-suppressed return-to-zero to non-return-to-zero format conversion using a fiber Bragg grating. Appl. Opt. 2015, 54, 6344–6350. [Google Scholar] [CrossRef]
  28. Cao, H.; Shu, X.; Atai, J.; Zuo, J.; Xiong, B.; Shen, F.; Cheng, J. Fiber Bragg grating based notch filter for bit-rate-transparent NRZ to PRZ format conversion with two-degree-of-freedom optimization. J. Opt. 2015, 17, 025702. [Google Scholar] [CrossRef]
  29. Cao, H.; Shu, X.; Atai, J.; Gbadebo, A.; Xiong, B.; Fan, T.; Tang, H.; Yang, W.; Yu, Y. Optimally-designed single fiber Bragg grating filter scheme for RZ-OOK/DPSK/DQPSK to NRZ-OOK/DPSK/DQPSK format conversion. Opt. Express 2014, 22, 30442–30460. [Google Scholar] [CrossRef] [PubMed]
  30. Islam, M.J.; Atai, J. Stability of gap solitons in dual-core Bragg gratings with cubic-quintic nonlinearity. Laser Phys. Lett. 2014, 12, 015401. [Google Scholar] [CrossRef]
  31. Cao, H.; Atai, J.; Yu, Y.; Xiong, B.; Zhou, Y.; Cai, J.; Shu, X. Carrier-suppressed return-to-zero to non-return-to-zero format conversion based on a single fiber Bragg grating with knife-shaped spectra. Appl. Opt. 2014, 53, 5649–5653. [Google Scholar] [CrossRef]
  32. Dasanayaka, S.; Atai, J. Moving Bragg grating solitons in a cubic–quintic nonlinear medium with dispersive reflectivity. Phys. Rev. E 2013, 88, 022921. [Google Scholar] [CrossRef]
  33. Cao, H.; Atai, J.; Shu, X.; Chen, G. Direct design of high channel-count fiber Bragg grating filters with low index modulation. Opt. Express 2012, 20, 12095–12110. [Google Scholar] [CrossRef]
  34. Baratali, B.; Atai, J. Gap solitons in dual-core Bragg gratings with dispersive reflectivity. J. Opt. 2012, 14, 065202. [Google Scholar] [CrossRef]
  35. Dasanayaka, S.; Atai, J. Interactions of solitons in Bragg gratings with dispersive reflectivity in a cubic-quintic medium. Phys. Rev. E 2011, 84, 026613. [Google Scholar] [CrossRef] [PubMed]
  36. Dasanayaka, S.; Atai, J. Stability of Bragg grating solitons in a cubic–quintic nonlinear medium with dispersive reflectivity. Phys. Lett. A 2010, 375, 225–229. [Google Scholar] [CrossRef]
  37. Carroll, S.S.; Atai, J. Collision dynamics of solitons in a multi-channel stabilized dispersion managed link. J. Opt. A: Pure Appl. Opt. 2009, 11, 085407. [Google Scholar] [CrossRef]
  38. Neill, D.R.; Atai, J. Gap solitons in a hollow optical fiber in the normal dispersion regime. Phys. Lett. A 2007, 367, 73–82. [Google Scholar] [CrossRef]
  39. Atai, J.; Malomed, B.A.; Merhasin, I.M. Stability and collisions of gap solitons in a model of a hollow optical fiber. Opt. Commun. 2006, 265, 342–348. [Google Scholar] [CrossRef]
  40. Neill, D.R.; Atai, J. Collision dynamics of gap solitons in Kerr media. Phys. Lett. A 2006, 353, 416–421. [Google Scholar] [CrossRef]
  41. Atai, J.; Malomed, B.A. Gap solitons in Bragg gratings with dispersive reflectivity. Phys. Lett. A 2005, 342, 404–412. [Google Scholar] [CrossRef]
  42. Atai, J.; Malomed, B.A. Stability and interactions of solitons in asymmetric dual-core optical waveguides. Opt. Commun. 2003, 221, 55–62. [Google Scholar] [CrossRef]
  43. Atai, J.; Malomed, B.A. Spatial solitons in a medium composed of self-focusing and self-defocusing layers. Phys. Lett. A 2002, 298, 140–148. [Google Scholar] [CrossRef]
  44. Nistazakis, H.; Frantzeskakis, D.; Atai, J.; Malomed, B.; Efremidis, N.; Hizanidis, K. Multichannel pulse dynamics in a stabilized Ginzburg-Landau system. Phys. Rev. E 2002, 65, 036605. [Google Scholar] [CrossRef] [PubMed]
  45. Atai, J.; Malomed, B.A. Solitary waves in systems with separated Bragg grating and nonlinearity. Phys. Rev. E 2001, 64, 066617. [Google Scholar] [CrossRef]
  46. Atai, J.; Malomed, B.A. Families of Bragg-grating solitons in a cubic–quintic medium. Phys. Lett. A 2001, 284, 247–252. [Google Scholar] [CrossRef]
  47. Arnous, A.H.; Zhou, Q.; Biswas, A.; Guggilla, P.; Khan, S.; Yıldırım, Y.; Alshomrani, A.S.; Alshehri, H.M. Optical solitons in fiber Bragg gratings with cubic–quartic dispersive reflectivity by enhanced Kudryashov’s approach. Phys. Lett. A 2022, 422, 127797. [Google Scholar] [CrossRef]
  48. Tanwar, D.V.; Kumar, M.; Tiwari, A.K. Lie symmetries, invariant solutions and phenomena dynamics of Boiti–Leon–Pempinelli system. Phys. Scr. 2022, 97, 075209. [Google Scholar] [CrossRef]
  49. Tanwar, D.V. Lie symmetry reductions and generalized exact solutions of Date–Jimbo–Kashiwara–Miwa equation. Chaos Solitons Fractals 2022, 162, 112414. [Google Scholar] [CrossRef]
  50. Tanwar, D.V.; Kumar, M. On Lie symmetries and invariant solutions of Broer–Kaup–Kupershmidt equation in shallow water of uniform depth. J. Ocean. Eng. Sci. 2022, in press. [CrossRef]
  51. Tanwar, D.V. Optimal system, symmetry reductions and group-invariant solutions of (2+ 1)-dimensional ZK-BBM equation. Phys. Scr. 2021, 96, 065215. [Google Scholar] [CrossRef]
  52. Tanwar, D.V.; Wazwaz, A.M. Lie symmetries and exact solutions of KdV–Burgers equation with dissipation in dusty plasma. Qual. Theory Dyn. Syst. 2022, 21, 164. [Google Scholar] [CrossRef]
  53. Bluman, G.W.; Cheviakov, A.F.; Anco, S.C. Applications of Symmetry Methods to Partial Differential Equations; Springer: Berlin/Heidelberg, Germany, 2010; Volume 168. [Google Scholar]
  54. Arrigo, D.J. Symmetry Analysis of Differential Equations: An Introduction; John Wiley & Sons: Hoboken, NJ, USA, 2015. [Google Scholar]
  55. Cherniha, R.; Serov, M.; Pliukhin, O. Nonlinear Reaction-Diffusion-Convection Equations: Lie and Conditional Symmetry, Exact Solutions and Their Applications; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  56. Hydon, P.E. Symmetry Methods for Differential Equations: A Beginner’s Guide; Number 22; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  57. Bluman, G.; Anco, S. Symmetry and Integration Methods for Differential Equations; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2008; Volume 154. [Google Scholar]
  58. Habiba, U.; Salam, M.A.; Hossain, M.B.; Datta, M. Solitary wave solutions of Chafee-Infante equation and (2 + 1)-dimensional breaking soliton equation by the improved Kudryashov method. Glob. J. Sci. Front. Res. 2019, 19, 1–9. [Google Scholar] [CrossRef]
  59. Malik, S.; Kumar, S. Pure-cubic optical soliton perturbation with full nonlinearity by a new generalized approach. Optik 2022, 258, 168865. [Google Scholar] [CrossRef]
  60. Wazwaz, A.M. The extended tanh method for new solitons solutions for many forms of the fifth-order KdV equations. Appl. Math. Comput. 2007, 184, 1002–1014. [Google Scholar] [CrossRef]
Figure 1. Profiles of bright solitons in fiber Bragg gratings.
Figure 1. Profiles of bright solitons in fiber Bragg gratings.
Symmetry 15 00963 g001
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

Malik, S.; Kumar, S.; Biswas, A.; Yıldırım, Y.; Moraru, L.; Moldovanu, S.; Iticescu, C.; Moshokoa, S.P.; Bibicu, D.; Alotaibi, A. Gap Solitons in Fiber Bragg Gratings Having Polynomial Law of Nonlinear Refractive Index and Cubic–Quartic Dispersive Reflectivity by Lie Symmetry. Symmetry 2023, 15, 963. https://doi.org/10.3390/sym15050963

AMA Style

Malik S, Kumar S, Biswas A, Yıldırım Y, Moraru L, Moldovanu S, Iticescu C, Moshokoa SP, Bibicu D, Alotaibi A. Gap Solitons in Fiber Bragg Gratings Having Polynomial Law of Nonlinear Refractive Index and Cubic–Quartic Dispersive Reflectivity by Lie Symmetry. Symmetry. 2023; 15(5):963. https://doi.org/10.3390/sym15050963

Chicago/Turabian Style

Malik, Sandeep, Sachin Kumar, Anjan Biswas, Yakup Yıldırım, Luminita Moraru, Simona Moldovanu, Catalina Iticescu, Seithuti P. Moshokoa, Dorin Bibicu, and Abdulaziz Alotaibi. 2023. "Gap Solitons in Fiber Bragg Gratings Having Polynomial Law of Nonlinear Refractive Index and Cubic–Quartic Dispersive Reflectivity by Lie Symmetry" Symmetry 15, no. 5: 963. https://doi.org/10.3390/sym15050963

APA Style

Malik, S., Kumar, S., Biswas, A., Yıldırım, Y., Moraru, L., Moldovanu, S., Iticescu, C., Moshokoa, S. P., Bibicu, D., & Alotaibi, A. (2023). Gap Solitons in Fiber Bragg Gratings Having Polynomial Law of Nonlinear Refractive Index and Cubic–Quartic Dispersive Reflectivity by Lie Symmetry. Symmetry, 15(5), 963. https://doi.org/10.3390/sym15050963

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