Abundant New Optical Soliton Solutions to the Biswas–Milovic Equation with Sensitivity Analysis for Optimization

: This study extensively explores the Biswas–Milovic equation (BME) with Kerr and power law nonlinearity to extract the unique characteristics of optical soliton solutions. These optical soliton solutions have different applications in the field of precision in optical switching, applications in waveguide design, exploration of nonlinear optical effects, imaging precision, reduced intensity fluctuations, suitability for optical signal processing in optical physics, etc. Through the powerful ( G ′ / G , 1/ G )-expansion analytical method, a variety of soliton solutions are expressed in three distinct forms: trigonometric, hyperbolic, and rational expressions. Rigorous validation using Mathematica software ensures precision, while dynamic visual representations vividly portray various soliton patterns such as kink, anti-kink, singular soliton, hyperbolic, dark soliton, and periodic bright soliton solutions. Indeed, a sensitivity analysis was conducted to assess how changes in parameters affect the exact solutions, aiding in the understanding of system behavior and informing decision-making, especially in accurately designing or analyzing real-world optical phenomena. This investigation reveals the significant influence of parameters λ , τ , c , B , and K on the precise solutions in Kerr and power law nonlinearities within the BME. Notably, parameter λ exhibits consistently high sensitivity across all scenarios, while parameters τ and c demonstrate pronounced sensitivity in scenario III. The outcomes derived from this method are distinctive and carry significant implications for the dynamics of optical fibers and wave phenomena across various optical systems.


Introduction
Nonlinear equations (NLEs) are not just mathematical puzzles; instead, they are the hidden language of nature's most captivating phenomena.Unlike the straightforward lines of linear equations, NLEs capture the intricate dances of heat flow, light waves, and more, revealing the true complexity of our world.These equations unlock the secrets of solitons, solitary waves that travel through optical fibers, and hold the key to unraveling countless scientific mysteries.From the scorching heat of a furnace to the shimmering ripples of light pulses, NLEs provide a powerful lens to understand and harness these phenomena.
By deciphering their stories, mathematicians and engineers can build better fiber optic networks, design efficient heat transfer systems, and unlock the mysteries beyond the realm of light, shaping the future in remarkable ways.Furthermore, understanding and controlling soliton dynamics are crucial for various optical applications, encompassing optical communication, nonlinear optics, and ultrafast photonics.Optical solitons offer unique benefits, such as stable propagation over long distances and the preservation of pulse shape and integrity.As a result, accurate modeling of soliton dynamics is of paramount importance for advancing and optimizing optical devices and systems tailored for these specific applications [1][2][3][4].Within the realm of optical fiber communication, the significance of the BME is notable, as it stands as a valuable tool for modeling diverse soliton solutions.These solitons represent solitary waves capable of traversing long distances without altering their form, a characteristic crucial for transmitting information without distortion [5,6].Furthermore, their equations are integrable, indicating that they harbor an infinite number of conservation laws.This investigation delves into the interaction between parabolic and power law nonlinearities within the framework of the BME.Here is the form of the BME [5]: ( where M = M(x, y, t) correspond to the complex function, and In Equation (1), there are three components: the first signifies the overall evolution, the second term encapsulates the function C, representing the nonlinearity, and the third term encompasses the group velocity dispersion.Also, C is a real-valued algebraic function, and n serves as a constraint that extends from the nonlinear Schrödinger equation (NLSE) to the BME.The variables x, y, and t are independent spatial and temporal variables, while B, K, and τ are real quantities.
Generally, n is larger than or equal to 1, and once n equals 1, Equation (1) shortens to the (2 + 1)-dimensional NLSE structure of the BME.This investigation centers on Equation (1), which shows evident nonlinearity, with the hypothesis that n corresponds 1 is a basic NLSE form for optical systems, particularly optical fibers.

Among various methods, the G ′
G , 1 G -expansion method is a powerful analytical technique utilized in solving NLEs.It involves representing solutions as power series, where coefficients are figured out using two variables G ′ G and 1 G .These coefficients are obtained by substituting the series into the equation and matching coefficients of similar terms.This method is widely employed to reveal optical solutions effectively.Numerous researchers [42][43][44][45][46][47][48][49] have utilized this technique to reveal solutions for various NLEs.Throughout history, researchers have applied diverse mathematical models like the parabolic law, Kerr law, power law, and Kudryashov's quintuple power law to probe optical soliton solutions within the framework of the BME equation, employing various solving techniques.While earlier investigations primarily centered on revealing bright and dark soliton solutions, our study takes a different approach, aiming to investigate a wider array of soliton patterns.This includes kink, anti-kink, singular soliton, hyperbolic, dark soliton, and periodic bright soliton patterns.As of present, there has been no investigation of BME employing this technique with Kerr law and power law nonlinearity.This investigation aims to use the above-mentioned methodology to derive soliton solutions for the provided nonlinear equations.In addition to our primary investigation, we embarked on a comprehensive first-order sensitivity analysis, colloquially termed local sensitivity analysis.This analysis aims to elucidate the nuanced impacts of minor perturbations in each input parameter on the resultant outcomes within a specific operational context.This analytical endeavor involves the meticulous computation of sensitivities about the output concerning each input parameter, employing numerical methodologies such as finite differences.By scrutinizing these sensitivities, we tried to unveil the intricate interplay between individual input parameters and their consequential effects on the overarching solutions.Such insights are pivotal in augmenting our understanding of the complex dynamics of the BME's exact solutions.This paper is divided into several sections: (i) Methodology overview in Section 2. (ii) Discussion of parabolic and power law nonlinearity in Section 3. (iii) Application of the methodology to derive solutions from the equations in Section 4. (iv) Exploration of dynamic representations, showcasing soliton behaviors through 3D graphs, 2D graphs, and contours, along with sensitivity analysis of obtained solutions in Section 5. (v) Concluding remarks in Section 6. (vi) Finally, a list of relevant references.

A Short Overview of (G ′ /G, 1/G)-Expansion Technique
Here, we provide an in-depth explanation of the basic steps required to apply the (G ′ /G, 1/G) technique, which was first suggested by [50], to the analysis of NLEs.To aid in this analytical process, a careful formulation of an auxiliary linear ordinary differential equation of the following form is required: where ' ′ ' denotes first derivatives concerning φ, and the variables are: Equation (3) satisfies the subsequent correlations: The result of Equation ( 2) mentioned earlier fluctuates depending on λ, which falls into three separate scenarios [43,51]: Scenario I: For λ > 0. Equation (2) gives the following form of the general solution (GS): which generates: where is the arbitrary constant.Scenario II: For λ < 0 Equation (2) yields the succeeding form of the GS: subsequent in: where is the arbitrary constant.Scenario III: For λ = 0 Equation (2) results in the following form of the GS: which gives: where l 1 and l 2 are the arbitrary constant.Consider a general form of NLEs incorporating three independent variables (x, y, and t) as follows: where L denotes a polynomial function reliant on the variables encompassed in M and ∂x∂y and so on.To transform Equation (11), a new variable φ is introduced, governed by the following relation: M(x, y, t) = e iθ Z(φ), where θ = x + y + νt + σ 0 and φ = cx + dy + ϖt.Equation ( 12) portrays a wave characterized by a velocity of ϖ, phase described by θ, and a fixed starting point defined by the phase constant σ 0 .However, to completely grasp its behavior, three real values, c, d, and ν, remain unresolved and require determination.Now, Equation ( 11) is converted by using the transformation written in Equation ( 12) to a new form as follows: In this context, J is the new polynomial encapsulating Z along with its ordinary derivatives.
The GS of Equation ( 13) is provided by the following equation, derived through the previously discussed method as follows: Here, a 0 , a i , and b i (i = 1, 2, 3, . . . . . . . . ., N) are the arbitrary constant fulfilling the requirement a 2 N + b 2 N ̸ = 0, and the parameter N is defined as a positive balance number.Exact Solution Computation Workflow: First, applying the homogeneous balance principle to the transformed ODE provides the value of the balance number, then switching the value of N into Equation ( 14) and integrating the modified equation into Equation (13), with the utilization of Equations ( 6), (8) and (10), the left-hand side of Equation ( 13) transforms a polynomial containing the terms ℊ and F .In this polynomial, the degree accompanying F does not exceed 1, while the degree of ℊ ranges from 0 to any integer.Establishing a system of algebraic equations involves equating the coefficients of terms with corresponding powers within the polynomial to 0, incorporating parameters such as a i , b i , λ, η, and others.Employing Mathematica software 14 aids in solving these algebraic equations to determine the values of these parameters.
Subsequently, these values are substituted into Equation (13), enabling the extraction of soliton solutions represented by trigonometric functions (as in Equation ( 5)), hyperbolic functions (as in Equation ( 7)), and rational functions (as in Equation ( 9)).Through the utilization of Equation ( 12), the exact solution of Equation ( 1) is derived, encompassing trigonometric, hyperbolic, and rational functions.

Implementation of Nonlinearities Laws
In this section, we analyzed Equation (1) using both the Ker and power law nonlinearity.The implementation of these laws on the BME is detailed in Sections 3.1 and 3.2 below.

Ker Law
For an exploration of Equation ( 1) under Kerr law nonlinearity, specifically when C(Θ) = Θ [52,53], the resulting NLSE is formulated as follows: By inserting Equation ( 12) into Equation ( 15), we obtain the following equations from the imaginary and real components: and Given that Z(φ) is non-zero and shows second-order derivatives, Equation ( 16) imposes the following restriction on the velocity: Hence, Equation ( 17) signifies the NODE (nonlinear ordinary differential equation) form of Equation (15) with Kerr law nonlinearity, subject to the constraint outlined in Equation (18).

Power Law
In examining Equation (1) under power law nonlinearity, specifically when C(Θ) = Θ n [52,53], the resulting NLSE is formulated as follows: In this scenario, for stability, it is imperative that 0 < n < 2. By inserting Equation (12) into Equation ( 19) and organizing the resulting expression, the following equations are derived: and Likewise, Equation ( 21) offers the following restriction on the velocity: Therefore, Equation ( 21) signifies the NODE form of Equation ( 19) with power law nonlinearity, subject to the constraint outlined in Equation (22).

Method's Application with Nonlinearities Laws
In this segment, we employ the methodologies elucidated in Section 2 to derive the optical soliton solution for Equation (1).

Application for Kerr Law
Through the application of the balancing principle in Equation ( 17) provides N = 1.Therefore, Equation ( 14) transforms into the following form: In this equation, the coefficients a 0 , a 1 and b 1 are constants that need to be decided.Scenario I. λ > 0 (For trigonometric) To obtain the solutions for Equation ( 17), we use Equations ( 4) and ( 6) and replace Equation (17) with Equation (23).With this technique, a polynomial with ℊ and F is created from the left side of Equation ( 17).This polynomial's coefficients, a 0 , a 1 and b 1 , are involved in a system of algebraic equations when each coefficient is set to zero.Solving this system of algebraic equations yields the follows arbitrary constants: , where τ = 1 2 2c − 4K + c 2 Kλ + d 2 Kλ and K < 0. Now, using these calculated values in Equation ( 14), we obtain the following solution of Equation ( 17): where φ = cx + dy + 2(c + d)Kt.
Reverting to its first form with the aid of Equation ( 12), we obtain the following form: where To derive the hyperbolic solution, apply the same procedure, resulting in: Now, using these calculated values in Equation ( 14), the following solution of Equation (17): where With the help of Equation ( 12), obtain solution of main PDEs: where When l 1 is confirmed to be non-zero and both η and l 2 are set to zero, Equation ( 27) takes on the form of the following soliton solution: where φ = cx + dy + 2(c + d)Kt and θ = x + y + νt + σ 0 .Scenario III.For λ = 0 (For rational) To achieve the rational solution, follow the same procedure outlined in scenario I, resulting in: , where K = c−τ 2 .Now, incorporating these computed values into Equation ( 14), we derive the following solution for Equation (17): where Inserting the transformation variables defined in Equation ( 12) into the above equation results in the following generalized solution for Equation (1): where When l 1 is confirmed to be non-zero and both η and l 2 are set to zero, Equation (30) takes on the form of the following soliton solution: where φ = cx + dy + 2(c + d)Kt and θ = x + y + νt + σ 0 .

Application for Power Law
Operating the balancing principle in Equation ( 21), we deduce that N = 1/n.Consequently, the following relation can be defined: Thus, employing the transformation, Equation (38) undergoes the following form: Applying the balancing principle once more and figuring out N = 1, the solution adopts the following structure: Now, employing a similar procedure as described in Section 4.1 for all scenarios, we derive the following solutions for Equation (33), then Equation ( 21), and finally Equation (1): Scenario I. λ > 0 (For trigonometric) The coefficients obtained from this case are as follows: Now, utilizing these computed values in Equation ( 34), we derive the following solution for Equation (33): where φ = cx + dy + 2(c + d)Kt.
Returning to its original form with the aid of Equation ( 12), we obtain the following equation: where When l 2 is confirmed to be non-zero and both η and l 1 are set to zero, Equation (36) takes on the form of the following soliton solution: where φ = cx + dy + 2(c + d)Kt and θ = x + y + νt + σ 0 .When l 1 is confirmed to be non-zero and both η and l 2 are set to zero, Equation (36) takes on the form of the following soliton solution: Plugging in the transformation variables outlined in Equation ( 12) into the given equation produces the following comprehensive solution for Equation (1): where φ = cx + dy + 2(c + d)Kt and θ = x + y + νt + σ 0 .

Graphs and Meanings
Harnessing the power of Mathematica, an advanced mathematical computing tool, we explored the intricate graphical patterns shown by the BME with Kerr law nonlinearity and BME with Power law nonlinearity.Our presentation encompasses a variety of visual representations, including 3D renderings, 2D graphical displays, and contour plots.These visualizations cover a broad range of parameter values for each relevant variable, allowing us to gain a thorough insight into the graphical behavior of these nonlinear equations across diverse parameter space.
To ensure clarity and conciseness, we opted for three representative solution sets for each case (BME with Kerr law nonlinearity and BME with power law nonlinearity) from our comprehensive results for visual representation, standardizing only three graphs, including 3D renderings, 2D graphical displays, and contour plots, for each nonlinearity.The specific constants linked to each graph are detailed in the corresponding figure captions.In 2D graphs, we consolidated multiple solutions within a single figure by varying the parameter t.

Visualization of the BME with Kerr Law Nonlinearity
The depiction in Figure 1, originating from Equation (25), vividly captures the essence of periodic bright solitons.These solitons, known for their periodic behavior and stable propagation, hold immense utility across various domains such as optical communications, nonlinear optics, and fiber optics.Their unique ability to convey information over extended distances without distortion is highly advantageous.Moreover, they serve as invaluable tools for studying wave dynamics and nonlinear phenomena, shedding light on the intricate behavior of light pulses within optical fibers and other waveguides.
The image portrayed in Figure 2, stemming from Equation (28), vividly illustrates the essence of kink-type solitons with singular characteristics.These solitons are recognized by their localized wave packets with a kink-shaped profile, showcasing distinctive behavior during optical fiber propagation.Renowned for their singular attributes, they hold significance across diverse optical fiber applications, spanning signal processing, nonlinear optics, and optical communications.Their exceptional features make them invaluable for enhancing the stability and efficiency of optical signal manipulation and transmission within fiber optic networks.
The depiction in Figure 1, originating from Equation (25), vividly captures the essence of periodic bright solitons.These solitons, known for their periodic behavior and stable propagation, hold immense utility across various domains such as optical communications, nonlinear optics, and fiber optics.Their unique ability to convey information over extended distances without distortion is highly advantageous.Moreover, they serve as invaluable tools for studying wave dynamics and nonlinear phenomena, shedding light on the intricate behavior of light pulses within optical fibers and other waveguides.The image portrayed in Figure 2, stemming from Equation ( 28), vividly illustrates the essence of kink-type solitons with singular characteristics.These solitons are recognized by their localized wave packets with a kink-shaped profile, showcasing distinctive behavior during optical fiber propagation.Renowned for their singular attributes, they hold significance across diverse optical fiber applications, spanning signal processing, nonlinear optics, and optical communications.Their exceptional features make them invaluable for enhancing the stability and efficiency of optical signal manipulation and transmission within fiber optic networks.Figure 3, originating from Equation ( 30), showcases dark-type solitons.These solitons, characterized by localized waveforms with reduced intensities, are essential for various optical communication tasks such as signal amplification, pulse compression, and signal regeneration within fiber optic networks.Their unique ability to keep waveform integrity during propagation through optical fibers ensures reliable and efficient data transmission over long distances in optical communication systems.Geometrically, a dark soliton is a depressed section within a wave, where the amplitude falls below the neighboring crests, creating a noticeable trough.This imagery mirrors a dim, less intense area   30), showcases dark-type solitons.These solitons, characterized by localized waveforms with reduced intensities, are essential for various optical communication tasks such as signal amplification, pulse compression, and signal regeneration within fiber optic networks.Their unique ability to keep waveform integrity during propagation through optical fibers ensures reliable and efficient data transmission over long distances in optical communication systems.Geometrically, a dark soliton is a depressed section within a wave, where the amplitude falls below the neighboring crests, creating a noticeable trough.This imagery mirrors a dim, less intense area amid the brighter wave peaks, showcasing the soliton's unique trait of maintaining its shape and conserving energy while moving through the medium.

Sensitivity Analysis:
Sensitivity analysis involves scrutinizing how variations in input parameters impact the output of a model or system within a localized region of the parameter space.In recent studies, researchers [54,55] have utilized this analysis to explore how adjustments in input parameters influence the output of NLEs.In this section, we conduct a local sensitivity analysis of our derived solution for each scenario (for  = 0,  < 0, and  > 0). Figure 4a illustrates the sensitivity of each parameter in Equation (25) for Scenario I. We note that parameters (, , ,  1 , Κ) show negative sensitivity, whereas  demonstrates notably high negative sensitivity.On the other hand, the remaining parameters (ℬ, ,  2 ) exhibit

Sensitivity Analysis
Sensitivity analysis involves scrutinizing how variations in input parameters impact the output of a model or system within a localized region of the parameter space.In recent studies, researchers [54,55] have utilized this analysis to explore how adjustments in input parameters influence the output of NLEs.In this section, we conduct a local sensitivity analysis of our derived solution for each scenario (for λ = 0, λ < 0, and λ > 0). Figure 4a illustrates the sensitivity of each parameter in Equation (25) for Scenario I. We note that parameters (λ, c, d, l 1 , K) show negative sensitivity, whereas λ demonstrates notably high negative sensitivity.On the other hand, the remaining parameters (B, η, l 2 ) exhibit positive sensitivity, although not as pronounced as λ.Additionally, it is apparent that l 2 exhibits almost negligible sensitivity compared to others.For Scenario II (λ < 0), the results remain largely consistent, with one notable difference: the parameter K shows slightly higher sensitivity compared to λ, as depicted in Figure 4b.In the case of a rational solution (λ = 0), the most sensitive parameters are c and τ, but they exhibit opposite directions of sensitivity, whereas other parameters are nearly negligible, as illustrated in Figure 4c.

The Graphical Representation of the BME with Power Law Nonlinearity
Figure 5, originating from Equation (38), depicts a singular soliton type.This variant of solitons within optical fibers enables the transmission of information over extensive distances with minimal distortion, representing notable advancements in telecommunications and optical signal processing.Renowned for their ability to preserve shape and velocity during propagation, these solitary waves play a critical role in elevating data transmission rates and mitigating signal degradation in fiber optic communication systems.

The Graphical Representation of the BME with Power Law Nonlinearity
Figure 5, originating from Equation (38), depicts a singular soliton type.This variant of solitons within optical fibers enables the transmission of information over extensive distances with minimal distortion, representing notable advancements in telecommunications and optical signal processing.Renowned for their ability to preserve shape and velocity during propagation, these solitary waves play a critical role in elevating data transmission rates and mitigating signal degradation in fiber optic communication systems.
The anti-kink soliton solution, derived from Equation (40) and depicted in Figure 6, is vital in optical fibers.It effectively counters kinks or distortions, preserving signal integrity and ensuring reliable long-distance communication.This is crucial in high-speed networks and fiber optic systems where signal degradation is detrimental.Additionally, anti-kink solitons enhance overall efficiency and stability, supporting the advancement of optical fiber-based technologies and high-capacity data transmission networks.The anti-kink soliton solution, derived from Equation (40) and depicted in Figure 6, is vital in optical fibers.It effectively counters kinks or distortions, preserving signal integrity and ensuring reliable long-distance communication.This is crucial in high-speed networks and fiber optic systems where signal degradation is detrimental.Additionally, anti-kink solitons enhance overall efficiency and stability, supporting the advancement of optical fiber-based technologies and high-capacity data transmission networks.The solution derived from Equation (43) and illustrated in Figure 7 is of the hyperbolic type.Within optical fibers, this type of solution ensures stable signal transmission by preserving shape and velocity, minimizing distortion across long distances.Offering high energy and bandwidth, hyperbolic solutions facilitate efficient signal amplification and support high-capacity data transmission.Their robust characteristics enable extended transmission distances and improved signal quality.Moreover, hyperbolic solitons enable nonlinear effects such as soliton self-frequency shifting and supercontinuum generation, broadening their applications in nonlinear optics and photonics.Overall, they significantly elevate the performance and capabilities of optical fiber communication networks.The solution derived from Equation ( 43) and illustrated in Figure 7 is of the hyperbolic type.Within optical fibers, this type of solution ensures stable signal transmission by preserving shape and velocity, minimizing distortion across long distances.Offering high energy and bandwidth, hyperbolic solutions facilitate efficient signal amplification and support high-capacity data transmission.Their robust characteristics enable extended transmission distances and improved signal quality.Moreover, hyperbolic solitons enable nonlinear effects such as soliton self-frequency shifting and supercontinuum generation, broadening their applications in nonlinear optics and photonics.Overall, they significantly elevate the performance and capabilities of optical fiber communication networks.

Sensitivity Analysis:
Similarly, we performed a sensitivity analysis on our derived solution for each scenario of  = 0,  < 0, and  > 0 for the power law nonlinearity of BME.In Figure 8a, we depict the sensitivity of each parameter in Equation ( 36) for Scenario I. Notably, parameters (, , , , , Κ) display negative sensitivity, with λ showing particularly high negative sensitivity.Conversely, the remaining parameters (ℬ,  1 ,  2 ) exhibit positive sensitivity, albeit not as pronounced as λ.Moreover,  1 demonstrates almost negligible sensitivity compared to others.For Scenario II ( < 0), the results remain largely consistent, with a notable difference: the parameter ,  shows positive sensitivity, as depicted in Figure 8b.In the case of a rational solution ( = 0), the most sensitive parameters are c and τ, but they exhibit opposite directions of sensitivity, whereas other parameters are nearly negligible, as illustrated in Figure 8c, which is similar to the Kerr law nonlinearity of BME sensitivity.

Sensitivity Analysis
Similarly, we performed a sensitivity analysis on our derived solution for each scenario of λ = 0, λ < 0, and λ > 0 for the power law nonlinearity of BME.In Figure 8a, we depict the sensitivity of each parameter in Equation ( 36) for Scenario I. Notably, parameters (λ, c, d, η, τ, K) display negative sensitivity, with λ showing particularly high negative sensitivity.Conversely, the remaining parameters (B, l 1 , l 2 ) exhibit positive sensitivity, albeit not as pronounced as λ.Moreover, l 1 demonstrates almost negligible sensitivity compared to others.For Scenario II (λ < 0), the results remain largely consistent, with a notable difference: the parameter c, d shows positive sensitivity, as depicted in Figure 8b.In the case of a rational solution (λ = 0), the most sensitive parameters are c and τ, but they exhibit opposite directions of sensitivity, whereas other parameters are nearly negligible, as illustrated in Figure 8c, which is similar to the Kerr law nonlinearity of BME sensitivity.

Conclusions
The ( / ,1/ ) method was applied, and several precise optical soliton solutions were obtained for the Kerr and power law nonlinearities of BME.We used an analytical approach to generate several types of nonlinear wave patterns inside these nonlinearities, and different exact solutions in rational, hyperbolic, and trigonometric forms were revealed.Additionally, a sensitivity analysis was conducted to identify parameters that were more sensitive in both the negative and positive senses.The analysis underscores the profound influence of parameters , , , ℬ, and Κ on the precise solutions derived across both Kerr law and power law nonlinearities within the context of BME.Notably, parameter  emerges as particularly critical across all scenarios, exhibiting heightened sensitivity.Moreover, parameters  and  demonstrate heightened sensitivity specifically within scenario III.These findings underscore the pivotal role of these parameters in shaping the outcomes of the studied phenomena, offering valuable insights into the nuanced interplay of resultant solutions.Validation through Mathematica software 14 ensured precision, while dynamic visual representations illustrated soliton solutions

Conclusions
The (G ′ /G, 1/G) method was applied, and several precise optical soliton solutions were obtained for the Kerr and power law nonlinearities of BME.We used an analytical approach to generate several types of nonlinear wave patterns inside these nonlinearities, and different exact solutions in rational, hyperbolic, and trigonometric forms were revealed.Additionally, a sensitivity analysis was conducted to identify parameters that were more sensitive in both the negative and positive senses.The analysis underscores the profound influence of parameters λ, τ, c, B, and K on the precise solutions derived across both Kerr law and power law nonlinearities within the context of BME.Notably, parameter λ emerges as particularly critical across all scenarios, exhibiting heightened sensitivity.Moreover, parameters τ and c demonstrate heightened sensitivity specifically within scenario III.These findings underscore the pivotal role of these parameters in shaping the outcomes of the studied phenomena, offering valuable insights into the nuanced interplay of resultant solutions.Validation through Mathematica software 14 ensured precision, while dynamic visual representations illustrated soliton solutions exhibiting different patterns, including bright solitons, dark solitons, singular solitons, king solitons, anti-kink, and parabolic types.These soliton solutions offer advantages such as energy localization, propagation stability, predictable interactions, precision in optical switching, applications in waveguide design, exploration of nonlinear optical effects, imaging precision, reduced intensity fluctuations, and suitability for optical signal processing in optical physics.The outcomes of this study are novel and hold promise for various forms of BME and related research in optical physics.Moreover, the demonstrated effectiveness and robustness of the applied method in obtaining diverse solitons within different forms of BME further underscore its significance.

Figure 3 ,
Figure3, originating from Equation(30), showcases dark-type solitons.These solitons, characterized by localized waveforms with reduced intensities, are essential for various optical communication tasks such as signal amplification, pulse compression, and signal regeneration within fiber optic networks.Their unique ability to keep waveform integrity during propagation through optical fibers ensures reliable and efficient data transmission over long distances in optical communication systems.Geometrically, a dark soliton is a depressed section within a wave, where the amplitude falls below the neighboring crests, creating a noticeable trough.This imagery mirrors a dim, less intense area amid the brighter wave peaks, showcasing the soliton's unique trait of maintaining its shape and conserving energy while moving through the medium.

Mathematics 2024 ,
12,  x FOR PEER REVIEW 14 of 21 directions of sensitivity, whereas other parameters are nearly negligible, as illustrated in Figure4c.