Novel dynamics of the Fokas-Lenells model in Birefringent fibers applying different integration algorithms

Optical soliton disturbance is one of the most discernible bounds of exploration in the telecommunication industry [1,2,3,4,5,6]. The realm of studying soliton transmission in optical fibers has flourished in the last few decades due to its rich prominence of occurrences explaining these events [7,8,9,10,11]. A soliton refers to a self-localized out-come of a non-linear evolutionary equation narrating the progression of a non-linear dynamical process with an unbounded number of degrees of freedom [12,13,14,15,16,17]. Solitary wave and soliton solutions can be distinguished by their unique interaction dynamics. Solitons possess the remarkable attribute of retaining their shape unchanged during interactions, leading to a mere phase shift as the outcome. On the other hand, the presence of solitary wave results indicates a precise equilibrium between dispersion and nonlinearity, often necessitating specific conditions that cannot be universally established [18,19,20]. If optical losses are not present, the fundamental equation governing the transmission of optical pulses is in mono-mode fiber. It is described by the renowned non-linear Schrödinger’s (NLS) equation [21,22,23]. Apart from its significance in non-linear optical phenomena, the non-linear Schrödinger’s equation also emerges in diverse domains such as Bose-Einstein condensates, non-linear quantum field theory, and water waves biomolecular kinetics. For the NLS equation, the inquisition of exact soliton solution and its augmentation is one of the most emergent tasks in non-linear science. These applications offer a wide range of uses in the development of fiber optic amplifiers, communication links based on solitary waves, and more. In recent years, a multitude of adaptations to the NLS equation have emerged, catering to distinct physical scenarios and specific requirements. Examples comprise the cubic-quintic NLS equation, the Fokas-Lenells (FL) system, the cubic-quintic-septic NLS equation, and so on. By employing the bi-Hamiltonian method, the FL model, which represents an integrable extension of the NLS equation, was first obtained [24].

The FL equation holds significant importance in explaining optical rogue waves and the behavior of pulse transmission in optical fibers [25]. Moreover, it is of great significance in analyzing the behavior of electromagnetic vibrations in birefringent media, especially when subjected to differential group delay caused by manufacturing flaws and other imperfections in fibers.

The FL equation has become widely recognized since its introduction. It has been investigated by many academics. The soliton outcomes for the FL system have been acquired via the generalized exponential function technique [26], the unified method [27], the mapping scheme [28], the extended trial function procedure [29], the trial equation strategy [30], and the complex envelope function hypothesis [31]. The bright and dark soliton solutions of the FL structure are also obtained utilizing a direct scheme by Matsuno [32, 33]. Optical soliton outcomes to the FL equation have recently been achieved using the φ⁶-model expansion approach [34]. A specific analytical depiction of the rogue fluctuations of the FL equation has been presented recently [35]. Optical soliton solutions of the mentioned model is obtained by the modified Kudryashov’s, the exp(ϕ (ξ)) function, and the Sine-Gordon expansion techniques [36].

Many researchers are dedicated to advancing the field of nonlinear models, exploring diverse perspectives in differential forms. Ma and Lee applied the transformed rational function procedure technique to extract exact solutions to the (3+1) dimensional Jimbo-Miwa model [37]. Ma and Fuchssteiner have used the Cole-Hopf transformation to observe the solutions of Kolmogorov-Petrovskii-Piskunov equation [38]. The (1/G′)-expansion technique was first presented by Yokus [39]. In this method, he used the second-order Riccati equation G″ + λ G′ + γ = 0, with arbitrary constants λ and γ. Then, many researchers utilized this technique to explore different nonlinear evolution equations (NLEEs) in order to derive solutions representing traveling waves. Relevant references for these studies can be found in [40]. The direct algebraic technique was initially introduced by Zhang in 2009 [41]. In this approach, he employed the first-order Riccati equation S′(ξ) = S²(ξ) + p, where p is an arbitrary constant. Subsequently, many researchers utilized this method to investigate various NLEEs and obtain solutions depicting traveling waves. The pertinent references for these studies are cited in [42, 43].

The construction of this work is outlined as follows: In Section 2, we present the model studied. In Section 3, we provide the ODE conversion of the FL model. Moving to Section 4, we explain the application of the direct algebraic method [41,42,43] and its solution procedure. Section 5 delves into the modified rational sine-cosine technique [44] and its practical implementation. The utilization of the (1/G′) expansion scheme [45] is elaborated in Section 6, accompanied by its application. Section 7 is dedicated to offering a clear understanding of the graphs through physical explanations. The comparison of outcomes and the novelty they bring forth are discussed in Section 8. Finally, Section 9 encompasses our concluding remarks on the study.

The FL model in its dimensionless form is introduced by [25,26,27,28,29,30,31,32,33,34,35,36]: (1) iψt+P1ψxx+P2ψxt+|ψ|2(qψ+iσψx)=i[kψx+b(|ψ|2nψ)x+c(|ψ|2n)xψ]. i{\psi _t} + {P_1}{\psi _{xx}} + {P_2}{\psi _{xt}} + |\psi {|^2}(q\psi + i\sigma {\psi _x}) = i[k{\psi _x} + b{(|\psi {|^{2n}}\psi )_x} + c{(|\psi {|^{2n}})_x}\psi ]. In the above equation ψ(x, t) corresponds to the complex wave amplitude. The value of P₁ is group velocity dispersion, the factor of P₂ denotes Spatiotemporal dispersion. The coefficient of q depicts the nonlinear component. Here, k corresponds to the inter-model dispersion when c is linear dispersion, and b comes from the self-steepening effect. Finally, n is form the full non-linearity exponent.

To solve the equation (1), the solution structure can be considered as [46] (2) ψ(ξ)=Q(ξ)eiθ, \psi (\xi ) = Q(\xi ){e^{i\theta }}, where ξ = x − mt for soliton velocity m with the soliton magnitude element Q and the phase component θ can be written as θ = −rx + lt + ɛ in which wave number l, frequency r, and phase value ɛ.

Substituting the value of equation (2) into equation (1) and parting the imaginary and real components.

The real part is (3) (P1−P2m)Q″+(P2lr−P1r2−kr−l)Q+(q+rσ)Q3−rQ[(2n+1)b+2nc]Q2n=0, ({P_1} - {P_2}m)Q'' + ({P_2}lr - {P_1}{r^2} - kr - l)Q + (q + r\sigma ){Q^3} - rQ[(2n + 1)b + 2nc]{Q^{2n}} = 0, and the imaginary part is (4) m+2rP1−P2(mr+l)−σQ2+k+[(2n+1)b+2cn]Q2n=0, m + 2r{P_1} - {P_2}(mr + l) - \sigma {Q^2} + k + [(2n + 1)b + 2cn]{Q^{2n}} = 0, (5) or,[(2n+1)b+2cn]Q2n=−m−2rP1+P2(mr+l)+σQ2−k. or,[(2n + 1)b + 2cn]{Q^{2n}} = - m - 2r{P_1} + {P_2}(mr + l) + \sigma {Q^2} - k. Putting [(2n + 1)b + 2cn]Q²ⁿ = 0, σ = 0, then we have follows b=−2cn2n+1, m=2rP1−P2l+kP2r−1. b = \frac{{ - 2cn}}{{2n + 1}},\,m = \frac{{2r{P_1} - {P_2}l + k}}{{{P_2}r - 1}}. Now from equation (3), we have, (6) (P1−P2m)Q″+(P1r2−l+rm(1−P2r))Q+qQ3=0. ({P_1} - {P_2}m)Q'' + ({P_1}{r^2} - l + rm(1 - {P_2}r))Q + q{Q^3} = 0.

Let the trial solution of the FL model be, (7) Q(ξ)=∑k=0NλkS(ξ), Q(\xi ) = \sum\limits_{k = 0}^N {\lambda _k}S(\xi ), where λ_k (k = 0, 1, 2, ⋯ , N) are real quantities and S(ξ) gratifies the next Riccati equation as, (8) S'(ξ)=S2(ξ)+p. {S^\prime}(\xi ) = {S^2}(\xi ) + p. Solutions of equation (8) are introduced in [41,42,43] and can be written as

for p < 0, (9) S(ξ)=−−ptanh(−pξ),−−pcoth(−pξ), S(\xi ) = \left\{ {\begin{array}{*{20}{c}}{ - \sqrt { - p} \tanh (\sqrt { - p} \xi ),}\\{ - \sqrt { - p} \coth (\sqrt { - p} \xi ),}\end{array}} \right.

for p > 0, (10) S(ξ)=ptan(p ξ),−pcot(p ξ), S(\xi ) = \left\{ {\begin{array}{*{20}{c}}{\sqrt p \tan (\sqrt p \,\xi ),}\\{ - \sqrt p \cot (\sqrt p \,\xi ),}\end{array}} \right.

for p = 0, (11) S(ξ)=−1ξ. S(\xi ) = - \frac{1}{\xi }.

By balancing between Q³ and Q″ in equation (6) one can easily obtained N = 1. For N = 1 equation equation (7) becomes (12) Q(ξ)=λ0+λ1S(ξ). Q(\xi ) = {\lambda _0} + {\lambda _1}S(\xi ). Now combining equation (12) and equation (8) and putting into equation (6), we get the solution set as, (13) λ0=0,k=−qr2P22λ12+2qpP22λ12−2qrP2λ12+qλ12+2P12P2,l=−qr2P2λ12+2qpP2λ12−qrλ12−2rP12P2,m=rP1+12qr2P2λ12+qpP2λ12−12qrλ12+krp2−1, \left\{ {\begin{array}{*{20}{l}}{{\lambda _0} = 0,k = - \frac{{q{r^2}P_2^2\lambda _1^2 + 2qpP_2^2\lambda _1^2 - 2qr{P_2}\lambda _1^2 + q\lambda _1^2 + 2{P_1}}}{{2{P_2}}},l = - \frac{{q{r^2}{P_2}\lambda _1^2 + 2qp{P_2}\lambda _1^2 - qr\lambda _1^2 - 2r{P_1}}}{{2{P_2}}},}\\{m = \frac{{r{P_1} + \frac{1}{2}q{r^2}{P_2}\lambda _1^2 + qp{P_2}\lambda _1^2 - \frac{1}{2}qr\lambda _1^2 + k}}{{r{p_2} - 1}},}\end{array}} \right. where q, r, p, P₁, P₂ and λ₁ are arbitrary constants.

Applying equations (8),(9),(10),(11),(12),(13) into equation (2), we get five types of exact solutions as below, ψ1(x,t)=(−−ptanh−p(−(rP1+12qr2P2λ12+qpP2λ12−12qrλ12+k)trp2−1+x))λ1eiθ,ψ2(x,t)=(−−pcoth−p(−(rP1+12qr2P2λ12+qpP2λ12−12qrλ12+k)trp2−1+x))λ1eiθ,ψ3(x,t)=(ptanp(−(rP1+12qr2P2λ12+qpP2λ12−12qrλ12+k)trp2−1+x))λ1eiθ,ψ4(x,t)=(−pcotp(−(rP1+12qr2P2λ12+qpP2λ12−12qrλ12+k)trp2−1+x))λ1eiθ,ψ5(x,t)=−λ1eiθ−rP1+12qr2P2λ12+qpP2λ12−12qrλ12+k)trp2−1+x, \begin{array}{*{20}{c}}{{\psi _1}(x,t) = ( - \sqrt { - p} \tanh \sqrt { - p} ( - \frac{{(r{P_1} + \frac{1}{2}q{r^2}{P_2}\lambda _1^2 + qp{P_2}\lambda _1^2 - \frac{1}{2}qr\lambda _1^2 + k)t}}{{r{p_2} - 1}} + x)){\lambda _1}{e^{i\theta }},}\\{{\psi _2}(x,t) = ( - \sqrt { - p} \coth \sqrt { - p} ( - \frac{{(r{P_1} + \frac{1}{2}q{r^2}{P_2}\lambda _1^2 + qp{P_2}\lambda _1^2 - \frac{1}{2}qr\lambda _1^2 + k)t}}{{r{p_2} - 1}} + x)){\lambda _1}{e^{i\theta }},}\\{{\psi _3}(x,t) = (\sqrt p \tan \sqrt p ( - \frac{{(r{P_1} + \frac{1}{2}q{r^2}{P_2}\lambda _1^2 + qp{P_2}\lambda _1^2 - \frac{1}{2}qr\lambda _1^2 + k)t}}{{r{p_2} - 1}} + x)){\lambda _1}{e^{i\theta }},}\\{{\psi _4}(x,t) = ( - \sqrt p \cot \sqrt p ( - \frac{{(r{P_1} + \frac{1}{2}q{r^2}{P_2}\lambda _1^2 + qp{P_2}\lambda _1^2 - \frac{1}{2}qr\lambda _1^2 + k)t}}{{r{p_2} - 1}} + x)){\lambda _1}{e^{i\theta }},}\\{{\psi _5}(x,t) = \frac{{ - {\lambda _1}{e^{i\theta }}}}{{ - \frac{{r{P_1} + \frac{1}{2}q{r^2}{P_2}\lambda _1^2 + qp{P_2}\lambda _1^2 - \frac{1}{2}qr\lambda _1^2 + k)t}}{{r{p_2} - 1}} + x}},}\end{array} where θ=−rx−12t(qr2P2λ12+2qpP2λ12−qrλ12−2rP1)P2+ɛ. \theta = - rx - \frac{1}{2}\frac{{t(q{r^2}{P_2}\lambda _1^2 + 2qp{P_2}\lambda _1^2 - qr\lambda _1^2 - 2r{P_1})}}{{{P_2}}} + \varepsilon . .

It is well known that the mentioned technique is a particular case of the transformed rational function procedure [37], and the adopted auxiliary equation (8) is a Riccati equation. The general solution for this equation was introduced in Reference [38].

The modified rational sine-cosine scheme presupposes that the resolution of equation (6) adopts the format [44], (14) Q(ξ)=Lsin(ξ)+1M+Ncos(ξ). Q(\xi ) = \frac{{L\sin (\xi ) + 1}}{{M + N\cos (\xi )}}. To find out the parameter’s value L, M and N, we insert equation (14) into equation (6) and gather the co-efficient of the fundamental functions, sin(ξ), cos(ξ), sin(ξ) cos(ξ), sin² (ξ), sin³ (ξ), and sin² (ξ) cos(ξ).

Next, we equate each obtained co-efficient to zero as the original equation is homogeneous. The process of solving these equations results in a set of solutions that form a family as (15) L=1,M=0,l=2r3P1P2+2r2kP2−2r2P1+rP1P2−2rk+kP2+P12r2P22−4rP2+P22+2,q=−N2(kP2+P1)2r2P22−4rP2+P22+2,m=2rP1−P2(2r3P1P2+2r2kP2−2r2P1+rP1P2−2rk+kP2+P1)2r2P22−4rP2+P22+2+kP2r−1, \left\{ {\begin{array}{*{20}{l}}{L = 1,M = 0,l = \frac{{2{r^3}{P_1}{P_2} + 2{r^2}k{P_2} - 2{r^2}{P_1} + r{P_1}{P_2} - 2rk + k{P_2} + {P_1}}}{{2{r^2}P_2^2 - 4r{P_2} + P_2^2 + 2}},q = - \frac{{{N^2}(k{P_2} + {P_1})}}{{2{r^2}P_2^2 - 4r{P_2} + P_2^2 + 2}},}\\{m = \frac{{2r{P_1} - {P_2}\frac{{(2{r^3}{P_1}{P_2} + 2{r^2}k{P_2} - 2{r^2}{P_1} + r{P_1}{P_2} - 2rk + k{P_2} + {P_1})}}{{2{r^2}P_2^2 - 4r{P_2} + P_2^2 + 2}} + k}}{{{P_2}r - 1}},}\end{array}} \right. where k, r, P₁, P₂ and N are arbitrary constants.

The exact solution gives the result, ψ6(x,t)=1+sin(−(2rp1−(2r3P1P2+2r2kP2−2r2P1+rP1P2−2rk+kP2+P1)P22r2P22−4rP2+P22+2+k)trP2−1+x)Ncos(−(2rp1−(2r3P1P2+2r2kP2−2r2P1+rP1P2−2rk+kP2+P1)P22r2P22−4rP2+P22+2+k)trP2−1+x)eiθ, {\psi _6}(x,t) = \frac{{1 + \sin ( - \frac{{(2r{p_1} - \frac{{(2{r^3}{P_1}{P_2} + 2{r^2}k{P_2} - 2{r^2}{P_1} + r{P_1}{P_2} - 2rk + k{P_2} + {P_1}){P_2}}}{{2{r^2}P_2^2 - 4r{P_2} + P_2^2 + 2}} + k)t}}{{r{P_2} - 1}} + x)}}{{N\cos ( - \frac{{(2r{p_1} - \frac{{(2{r^3}{P_1}{P_2} + 2{r^2}k{P_2} - 2{r^2}{P_1} + r{P_1}{P_2} - 2rk + k{P_2} + {P_1}){P_2}}}{{2{r^2}P_2^2 - 4r{P_2} + P_2^2 + 2}} + k)t}}{{r{P_2} - 1}} + x)}}{e^{i\theta }}, where θ=−rx+2r3P1P2+2r2kP2−2r2P1+rP1P2−2rk+kP2+P1)t2r2P22−4rP2+P22+2+ɛ \theta = - rx + \frac{{2{r^3}{P_1}{P_2} + 2{r^2}k{P_2} - 2{r^2}{P_1} + r{P_1}{P_2} - 2rk + k{P_2} + {P_1})t}}{{2{r^2}P_2^2 - 4r{P_2} + P_2^2 + 2}} + \varepsilon .

Let’s apply the fundamental assumption of equation (1) with the auxiliary form, (16) Q(ξ)=∑i=0Kβi1G′(ξ)i, Q(\xi ) = \sum\limits_{i = 0}^K {\beta _i}{\left( {\frac{1}{{G'(\xi )}}} \right)^i}, where β_k (k = 0, 1, 2, ⋯ , N) are real quantities and G(ξ) gratifies the next Riccati equation as, (17) G″(ξ)+λG′(ξ)+γ=0, G''(\xi ) + \lambda G'(\xi ) + \gamma = 0, whilst γ and λ are scalars. Solutions of equation (17) are introduced in [45] and can be written as (18) 1G′(ξ)=1−γλ+Bcosh(ξλ)−Bsinh(ξλ), \frac{1}{{G'(\xi )}} = \frac{1}{{ - \frac{\gamma }{\lambda } + B\cosh (\xi \lambda ) - B\sinh (\xi \lambda )}}, where B is a real constant. We can notice that for integer K = 1, there is a balance between Q³ and Q″ in equation (16). Then, the initial hypothesis in equation (16) is written as: (19) Q(ξ)=β0+β11G′(ξ). Q(\xi ) = {\beta _0} + {\beta _1}\left( {\frac{1}{{G'(\xi )}}} \right). Putting equation (16) and equation (17) into equation (6) and assigning each coefficient of 1, 1G′(ξ) \frac{1}{{G'(\xi )}} , 1(G′(ξ))2 \frac{1}{{{{(G'(\xi ))}^2}}} , and 1(G′(ξ))3 \frac{1}{{{{(G'(\xi ))}^3}}} to 0 gives the results, (20) β0=λβ12γ,k=−2qr2p22β12−qλ2p22β12−4qrP2β12+2qβ12+4γ2P14γ2P2,m=2rP1+2qr2P2β12−qλ2P2β12−2qrβ12−4rγ2P14γ2+kP2r−1,l=−2qr2P2β12−qλ2P2β12−2qrβ12−4rγ2P14γ2P2, \left\{ {\begin{array}{*{20}{l}}{{\beta _0} = \frac{{\lambda {\beta _1}}}{{2\gamma }},k = - \frac{{2q{r^2}p_2^2\beta _1^2 - q{\lambda ^2}p_2^2\beta _1^2 - 4qr{P_2}\beta _1^2 + 2q\beta _1^2 + 4{\gamma ^2}{P_1}}}{{4{\gamma ^2}{P_2}}},}\\{m = \frac{{2r{P_1} + \frac{{2q{r^2}{P_2}\beta _1^2 - q{\lambda ^2}{P_2}\beta _1^2 - 2qr\beta _1^2 - 4r{\gamma ^2}{P_1}}}{{4{\gamma ^2}}} + k}}{{{P_2}r - 1}},l = - \frac{{2q{r^2}{P_2}\beta _1^2 - q{\lambda ^2}{P_2}\beta _1^2 - 2qr\beta _1^2 - 4r{\gamma ^2}{P_1}}}{{4{\gamma ^2}{P_2}}},}\end{array}} \right. where q, r, P₁, P₂, λ, γ, and β₁ are arbitrary constants.

Applying equation (18), (19), (20) into equation (2), we get the exact result of the suggested model, Ψ7(x,t)=(λβ12γ+β1−γλ+Bcosh(ξλ)−Bsinh(ξλ))eiθ, {\Psi _7}(x,t) = (\frac{{\lambda {\beta _1}}}{{2\gamma }} + \frac{{{\beta _1}}}{{ - \frac{\gamma }{\lambda } + B\cosh (\xi \lambda ) - B\sinh (\xi \lambda )}}){e^{i\theta }}, where ξ=x−(2rP1+2qr2P2β12−qλ2P2β12−2qrβ12−4rγ2P14γ2+kP2r−1)t \xi = x - (\frac{{2r{P_1} + \frac{{2q{r^2}{P_2}\beta _1^2 - q{\lambda ^2}{P_2}\beta _1^2 - 2qr\beta _1^2 - 4r{\gamma ^2}{P_1}}}{{4{\gamma ^2}}} + k}}{{{P_2}r - 1}})t and θ=−rx−142qr2P2β12−qλ2P2β12−2qrβ12−4rγ2P1tγ2P2+ɛ \theta = - rx - \frac{1}{4}\left( {2q{r^2}{P_2}\beta _1^2 - q{\lambda ^2}{P_2}\beta _1^2 - 2qr\beta _1^2 - 4r{\gamma ^2}{P_1}} \right)\frac{t}{{{\gamma ^2}{P_2}}} + \varepsilon .

It is well known that the mentioned technique is a particular case of the transformed rational function procedure [37], and the adopted auxiliary equation (17) implies that 1G′ \frac{1}{{G'}} satisfies another Riccati equation. Therefore, all presented solutions here can be obtained from a general class of solutions to the general Riccati equation introduced in [38].

In the context of our analysis, we have explored the solution of ψ₁. The specific case of ψ₁ depicted in Fig.1(a–f). From the demonstration, we have noted distinctive characteristics associated with different components of ψ₁. Re(ψ₁) and Im(Ψ₁) reveal double periodic waves, while |ψ₁| gives a dark soliton solution. Solutions ψ₂ and ψ₅ display an identical wave pattern, which is depicted in Fig.2(a–f). ψ₅ is demonstrated in this instance. From the demonstration, distinctive characteristics associated with different components of ψ₅ have been noted. Re(ψ₅) and Im(ψ₅) reveal a singular bright-dark breather waves, whereas |ψ₅| gives a bright breather wave solution. The identical pattern is observed for solutions ψ₃, ψ₄, and ψ₆, which is illustrated in Fig.3(a–f) specifically for solution ψ₆. From the demonstration, distinctive characteristics associated with different components of ψ₆ have been noted. Re(ψ₆) and Im(ψ₆) reveal a singular multiple bright-dark breather waves, whereas |ψ₆| represents multiple bright breather waves solution. The specific case of ψ₇ depicted in Fig.4(a–f) is being discussed. From the demonstration, distinctive characteristics associated with different components of ψ₇ have been noted. Re(ψ₇) and Im(ψ₇) reveal a periodic breather wave, whereas |ψ₇| gives a bright breather waves solution.

Fig. 1

Fig. 2

Fig. 3

Fig. 4

This section outlines the uniqueness of the findings and the paper’s significant contribution to the field. Our research direction was established through a comparative assessment of our outcomes with those of recently published investigations [29,30,31,32,33,34,35,36]. A diverse array of analytical methodologies has been employed in the pursuit of soliton solutions to FL equation. The extended trial function technique has been applied to discover soliton solutions characterized as bright, dark, and singular for the FL equation [29]. A group of chirped solutions including bright, dark, and kink solutions is also obtained [30]. By using the complex envelop function ansatz, a class of exact solitary wave solution is found [31]. Matsuno employed a direct approach to successfully derive both bright and dark soliton solutions for the FL equation [32,33]. Optical soliton outcomes in the context of the FL model have been obtained by implementing the φ⁶-model expansion approach [34]. Recent developments include the presentation of a specific analytical depiction of rogue oscillations within the FL equation [35]. The methodologies involving the Modified Kudryashov’s, the exp(ϕ (ξ)) function, and the sine-Gordon expansion have been harnessed to procure new soliton solutions for the FL equation [36]. We observe that they found bright, dark soliton solutions, kink, and breather waves. However, our study expands the scope, revealing additional phenomena such as dark and bright periodic waves, dark soliton solutions, single and multiple bright breather waves, and single and multiple bright-dark breather waves with singular solutions. Importantly, several of our findings present distinct variations from those previously documented [29,30,31,32,33,34,35,36].

This study has presented the direct algebraic, the modified rational sine-cosine, and the (1/G′)-expansion schemes to integrate the FL model to retrieve optical solitons that can be used in birefringent fibers. These procedures allow us to simplify the proposed framework into an algebraic methodology that can be effectively solved using various computational packages. Through these techniques, we have been able to derive double periodic waves, bright soliton, dark soliton, single and multiple bright breather waves, and periodic breather waves of this model through symbolic computation. This underscores the reliability and utility of our suggested approaches, which can be readily applied to various complicated systems that arise in nonlinear science. These solutions are set to have a crucial impact on the progression of subsequent studies involving the model.