A study of the wave dynamics of the space–time fractional nonlinear evolution equations of beta derivative using the improved Bernoulli sub-equation function approach

The space–time fractional nonlinear Klein-Gordon and modified regularized long-wave equations explain the dynamics of spinless ions and relativistic electrons in atom theory, long-wave dynamics in the ocean, like tsunamis and tidal waves, shallow water waves in coastal sea areas, and also modeling several nonlinear optical phenomena. In this study, the improved Bernoulli sub-equation function method has been used to generate some new and more universal closed-form traveling wave solutions of those equations in the sense of beta-derivative. Using the fractional complex wave transformation, the equations are converted into nonlinear differential equations. The achieved outcomes are further inclusive of successfully dealing with the aforementioned models. Some projecting solitons waveforms, including, kink, singular soliton, bell shape, anti-bell shape, and other types of solutions are displayed through a three-dimensional plotline, a plot of contour, and a 2D plot for definite parametric values. It is significant to note that all obtained solutions are verified as accurate by substituting the original equation in each case using the computational software, Maple. Additionally, the results have been compared with other existing results in the literature to show their uniqueness. The proposed technique is effective, computationally attractive, and trustworthy to establish more generalized wave solutions.


Atangana beta-derivative
For modeling intricate systems and processes, fractional derivatives offer a more flexible and precise tool.Classical integer-order derivatives are unable to effectively capture the non-integer-order behavior of many realworld phenomena.Such systems can be more effectively represented using fractional calculus.The definition of a fractional derivative, which has been studied on by several academics, includes the Jumarie-modified Riemann-Liouville, conformable, and Caputo fractional derivative [36][37][38] .The beta-derivative is a newly proposed fractional derivative introduced by Atangana et al. 39 as follows: The abovementioned definition (1) satisfies the following characteristics: Based on previously mentioned features, the beta-derivative can easily convert the NLFPDEs into an ordinary differential equation that is non-linear with an integer order.Many researchers use the beta-derivative in various (1) (2) T σ a 1 f (t ) + a 2 g(t ) = a 1 T σ f (t ) + a 2 T σ g(t ) , for alla 1 , a 2 ǫR, T σ ( fog)(t www.nature.com/scientificreports/physical applications due to its flexibility, well-posedness of mathematical properties, and ability to generalize classical derivative to fractional order making it an effective tool in a wide range of scientific and engineering problems.Among all the definitions of fractional derivative that have been introduced so far, the Atangana beta derivative is more reliable.

Outline of the IBSEF method
The BSEFM has been extended to derive the IBSEF method 21 , which will be detailed in this section.
Step I. Let's consider the following equation for the fractional differential equation: where u = u(x, t) is a function which is not defined, and γ ∈ (0, 1] is the order of beta-derivative. Using the following wave transformation where m and k are the wave number and velocity respectively and 0 < γ ≤ 1 indicates fractional-order.It can be converted Eq. ( 8) into nonlinear ordinary differential equation (NLODE) based on the above transformation (9) as: Step II.According to the IBSEF method, the trial solution in Eq. ( 10) can be expressed as In accordance with the Bernoulli theory, we can assume the general form of Bernoulli differential equation for H ′ as underneath: where H = H(η) is the solution of Eq. (12)  Substituting the solution (11) into the Eq. ( 10) and making use of the Eq. ( 12), it gives the following polynomial �(H) of H: To determine the values of l, m and N , where l and m both are unknowns, we use homogeneous balancing technique of the largest order linear term with the highest order nonlinear term.
Step III.By setting the coefficients of �(H) to zero, we can obtain a system of equations: We must solve the system in order to find the values of p 0 , p 1 , p 2 , . . ., p d and q 0 , q 1 , q 2 , . . ., q c .
Step IV.The following couple of outcomes that we obtain from solving nonlinear BDE reliant on the values of the parameters α and β:

Advantages:
The approach has advantages over other methods, such as the exp-function method, tanhfunction method, basic G ′ /G -expansion method, etc.The suggested approach offers further exact traveling wave solutions with extra free parameters.These precise solutions are essential for revealing the underlying principles of physical phenomena.When applicable, the IBSEF method can yield closed-form solutions to space-time fractional equations, which provide a clear mathematical expression for the behavior of the system.On the other hand, when dealing with reduced differential equations of third order or below, using symbolic computation software like Maple substantially enhances the possibility of obtaining useful solutions to the related algebraic equations.However, it gets more difficult to guarantee the existence of solutions for the resultant algebraic equations in general as the order of the equations rises.The capacity of the recommended approach to accommodate a greater number of arbitrary constants than other existing methods makes them particularly useful in such circumstances.This strategy can solve problems where other approaches could fail thanks to this quality.

Limitations:
The method has a limitation that should be taken into account despite its many benefits.It occasionally creates solutions that are disguised versions of well-known solutions that may be found using other

Analysis of closed form solutions
By using the definition of beta derivative, we determined some advanced, broad-ranging and further inclusive closed form stable soliton solutions to the space-time fractional NLKG equation and mRLW equation with the help of IBSEF method.

The space-time fractional nonlinear Klein-Gordon equation
In 1926, the famous physicists Klein 40 and Gordon 41 first introduced the space-time fractional NLKG equation.Now consider this equation as follows 42 : where the constants a , ν and b are non-zero coefficients depending upon physical conditions of the system.These constants will be evaluated afterwards through the computational package Maple which will represent numerous phenomena of the solutions.Also, γ denotes a derivative whose order is fractional.
Applying the wave transformation where m be the number of waves, k be the wave speed, and η represents wave transformation.Substituting Eq. ( 18) into Eq.( 17), then Eq. ( 17) can be transformed to an integer order NLODE as follows: For the maximum order linear term V (ω) , the exponent of V is ωN + d − c − 1 and for the maximum degree nonlinear term V n whose exponent is n(d − c) + 1 .Since, soliton solutions are being sought, and thus the exponent of the linear term is identical to that of the nonlinear term.After balancing 43,44 these two terms we obtain ωN 19), we observe that V 3 (n = 3) is the nonlinear term of highest degree as well as the linear term of maximum order is V ′′ (ω = 2) .Now, using these values of ω and n we get the relationship shown below: Here N and c are random parameters, selecting c = 1 and N = 3 , which leads to d = 3 .Therefore, the general solution is provided by the trial solution to Eq. (19) in the following way: where H will be calculated from the IBSEF method.Now, to determine V ′ (η) and V ′′ (η) , we can write as follows: and 23) are used in Eq. ( 19) to produce a system of algebraic equations based on the coefficients of H polynomial.By solving the system of algebraic equations with the use of Maple 13, the constant values may be determined.Set 1.
Case 1a: When α = β By using the parameter values accumulated in (24) alongside (15) into the solution (20), we achieve the exponential function solution (EFS) to the space-time fractional NLKG equation written as: where the constants a , b , m , k are non-zero.

Case 1b: When α = β
Using the values of the coefficients accumulated in ( 24) alongside ( 16) into the solution (20), we achieve the hyperbolic trigonometric function solution (HTFS) to the space-time fractional NLKG equation given below: where the constants a , b , m , k are non-zero.
If we take the Eq. ( 27) coefficients along with Eq. ( 15) in Eq. ( 20), we find the EFS to the NLKG equation given as following: Case 2b: When α = β If we take the Eq. ( 27) coefficients along with Eq. ( 16) in Eq. ( 20), we find the HTFS shown as: Case 3a: When α = β.Putting Eq. ( 30) along with Eq. ( 15) in Eq. ( 20), we acquire the EFS to the NLKG equation written as: Vol:.( 1234567890)  30) along with Eq. ( 16) in Eq. ( 20), we acquire the HTFS to the NLKG equation given below: All the general solutions to the space-time fractional NLKG equation presented here are valid and effective, and interestingly, all of these solutions have not been previously investigated.

The space-time fractional modified regularized long-wave equation
The space-time fractional mRLW equation was initially deduced by Benjamin 45 in 1972 that explains generally the one-way propagation of long waves in certain nonlinear dispersive systems.This equation is regarded as a substitute for the Korteweg-de Vries equation, which is designed to describe numerous physical phenomena like plasma waves and shallow waters in coastal oceans.The space-time fractional mRLW equation is taken into consideration as 32 : where γ is a fractional derivative, and ζ , ε, and µ are arbitrary constants depend on the system's physical proper- ties.Following that, these constants will be assessed using the computational tool Maple, which will represent a variety of phenomena of the solutions.Using the complex fractional wave variable where m be the number of waves, k be the wave speed, and η represents wave transformation.Then the space-time fractional mRLW Eq. ( 24) is transformed into an integer order NLODE as follows: After integrating Eq. ( 25) with respect to η once then putting integrating constant = 0 , we get the equation of following form: Rewrite Eq. ( 26) as As the relationship among d , c , and N is covered in detail in section "The space-time fractional nonlinear Klein-Gordon equation", so it is not extensively outlined here.From Eq. ( 27), we observe that V ′′ is the largest order derivative term and V 3 is the nonlinear term of height order.Consequently, a relationship for d , c , and N is obtained by balancing the above-mentioned two maximum order terms in Eq. ( 27), which can be expressed in this way: and putting these values into Eq.( 11) which provides the trial solution of Eq. ( 27) in the ensuing structure: where H is the root of Eq. (12)  Setting the solution (38), alongside with ( 12) into (37), yields a polynomial in H .A set of equations for the parameters can be obtained, and by solving the system of algebraic equations using Maple 13, the constant values can be determined.
Category 1a: First, we assume α = β , because the solutions to the IBSEF approach improved Bernoulli equa- tion depend on these two variables.

Results discussion and physical explanation
The attained travelling wave solutions are represented in three types of diagrams namely a three-dimensional plotline, a plot of contour, and a 2D plot for diverse values of the arbitrary constants with the help of Mathematica.The importance of graphical representations is that they can aid in improving our comprehension of solutions and analysis of complex data 46 .The diagrams depict a range of solutions namely kink wave solution, singular soliton solution, periodic soliton wave solution, bell shape, anti-bell shape, and single kink type.Now, we have shown a graphical illustration of the proposed equations.The acquired results of the both selected equations for various parameter values are represented graphically in this section.For the values α = 3, β = 2, a = 1, b = 1 , solution V 1 1 (x, t) represents singular-kink shape wave featuring infinite tails in Fig. 1 within the duration 0 < x < 50 and 0 < t < 75 .The achieved solutions which represents in V 1 2 (x, t) and V 1 4 (x, t) indicates soliton shape wave within the duration 0 < x, t < 2 , 0 < x < 10 and 0 < t < 15 with the same parametric values are displayed in Figs. 2 and 3.The other solution V 1 3 (x, t) embodies in Fig. 4 for the values of α = 3, β = 2, a = 1, b = 1 within the interval 0 < x < 0.005 and 0 < t < 1.5 signifies (48) one-sided-kink wave.A kink wave is a type of wave that either travels from one asymptotic position to another or rises and remains constant at infinity.It should be noted that the solution represents in V 1 5 (x, t) and V 1 6 (x, t) deliver the same types of soliton solutions but for state forwardness we discard from here.
In the space-time fractional mRLW equation the attained solutions V 2 2 (x, t) and V 2 4 (x, t) are illustrated by the form of periodic soliton wave solutions designed for the ideals α   the interval 0 < x < 50 and 0 < t < 20 , 0 < x < 90 and 0 < t < 70 respectively are symbolized by Fig. 5  and 6.The profile of Fig. 7 in solution V 2 5 (x, t) shows the wave shape namely compacton for the values α = 3, β = 2, ζ = 1, ε = 1, µ = 1 in the interval 0 < x < 5.5 and 0 < t < 12.5 .A compacton is a special type of soliton that has a compact spatial profile, meaning it is limited to a finite region.Remarkably, compacton maintain their coherent shape after collisions, due to their soliton-like properties.The solution V    cropped here in Fig. 8 within the range 0 < x < 0.1 , 0 < t < 2 .In addition, the solution V 2 7 (x, t) illustrates in Fig. 9 depicts the compacton wave solution for α = 3, β = 2, ζ = 1, ε = 1, µ = 1 in the duration 0 < x < 70 , 0 < t < 75 .Finally, the solution V 2 8 (x, t) in Fig. 10 indicates the singular bell-shaped wave solution within the specified interval 0 < x, t < 1 .The graphs with similar shapes for the other two solutions V 2 9 (x, t) and V 2 10 (x, t) are omitted to maintain simplicity.

Results comparison
In this part, we compare our achieved results with other existing results in the literature with the other methods.Those suggested equations presented above were examined by earlier researchers using a variety of methods.The improved Bernoulli sub-equation function approach with beta-derivative has been used to solve the abovementioned equations and we attain several remarkable outcomes that are both more wide-ranging and relevant than those of the previous researchers.A comparison between existing solutions and our obtained solutions is presented in Tables 1 and 2.
We observe that the derived solution V 1 1 (x, t) and V 1 5 (x, t) of space-time fractional NLKG equation are similar to Ege and Misirli 28 solutions.The others solution V 1 2 (x, t) − V 1 4 (x, t) and V 1 6 (x, t) attain in this study are completely new and innovative and also have not reported in the prior study.Table 1.A comparison between the solutions obtained by Ege and Misirli 28 and the solutions obtained by our study for the space-time fractional nonlinear Klein-Gordon equation using the IBSEF method.

Table 2.
A comparison between the solutions obtained by Uddin et al. 34 and the solutions obtained by our study for the space-time fractional modified regularized long-wave equation using IBSEF method.
If we insert q 0 = q 1 = q −1 = 1, in solution (4.17) then solution u 24 (x, t) becomes: Vol.:(0123456789)  34 solutions.The other solutions V 2 2 (x, t) − V 2 8 (x, t) and V 2 10 (x, t) found in this study are completely fresh and cutting-edge and have not previously been reported.
It is important to note that the achieved solutions to the space-time fractional NLKG and mRLW equations are explain the relativistic electrons in atom theory, to understand the long-wave occupancy in the ocean, including tsunamis and tidal waves, the propagation of water waves in shallow water, accounting for dispersion and dissipation effects.

Conclusion
In this research, the space-time fractional NLKG equation and mRLW equation have been studied to the potential IBSEF method to develop some fresh and more universal closed-form traveling wave solutions of those equations in the sense of beta-derivative.We achieved exponential function solutions and hyperbolic trigonometric function solutions for numerous values of free parameters including bell-shaped, kink-shaped, single soliton, periodic soliton, compaction, and anti-kink-shaped.These solutions are demonstrated by three types of diagrams namely three-dimensional plotline, plot of contour, and 2D plot by using Mathematica.The acquired solutions can describe different types of phenomena like plasma waves in complex media, the propagation of water waves in shallow water, accounting for dispersion and dissipation effects, long-wave occupancy dynamics in the ocean, including tsunamis and tidal waves, which are essential for coastal hazard assessment and maritime safety, and understanding the behavior of spinless ion in relativistic particles.The accuracy of all solutions obtained was confirmed by substituting the original equation in each case using the computational software Maple and found them correct.This approach not only yields identical solutions, but it also has the potential to identify novel solutions that have not been reported by other researchers.This study demonstrates that the proposed technique is a useful, effective, and compatible mathematical tool for solving a broad range of NLFPDEs that arise in various domains of mathematical physics and engineering.The implications are far-reaching, as the solutions acquired through this methodology not only open the doors to resolving immediate problems but also lay foundation for comprehensive research exploration.We have studied the NLFPDEs having balance number one but there are lots of NLFPDEs whose balance number is two or more.The future research might get a gorgeous way in looking for solitary wave solutions to the other NLFPDEs by the suggested technique when the balance number is two or more. https://doi.org/10.1038/s41598-023-45423-6
It is noteworthy to see that the consequent solution V 2 1 (x, t) and V 2 9 (x, t) of space-time fractional mRLW equation are analogous to Uddin et al.